计算光波由光密入射到光疏介质时相位变化
hYMIe]kJ o"~ODN" L clear
z"* X/T 9PjL
4A close all
ez=$ ]cln })!d4EcZf n1=1.45,n2=1;
+]uW|owxo hO(8v&ns3 theta=0:0.1:90;
Hy5_iYP5 {[G2{ijRz a=theta*pi/180;
JIvVbI Kdh(vNB> rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
bhe~ekb SQ| pH" rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
x? Z)q4 jT"P$0sJAd tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ZOc1 v j ATCFdtNc ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
@%&;V( ";}Lf1M9 arp=angle(rp);
dl;~-'0 F +PIZ% ars=angle(rs);
8v@6 &ras@ EW* !_| atp=angle(tp);
pg~vteq5 I Gv_s+O-* ats=angle(ts);
tagkklJ~ +qq,;npi figure(1)
fpf1^ TZ 8^M5k%P subplot(1,2,1);
$'e;ScH } Uki )3( plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
r9z_8#cR txQyHQ) @ legend('arg(r_p)','arg(r_s)')
_ _cJ+%e N
?Jr8 xlabel('\theta_i')
Yao>F--? +Q_(wR"FS ylabel('\phi')
6l&m+!i :kh l}| title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
@Tb
T 9 \i;zpN\ axis([0 90 -3.5 3.5])
6cz%>@ ;KJJK#j grid on
c nvx TI< bc L>S$B subplot(1,2,2);
&;L4Cj$ q n4\6\0jq6 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
4NN-'Z>a j[NA3Vj1P legend('arg(t_p)','arg(t_s)')
xal,j* ,OWdp<z xlabel('\theta_i')
Hn)K;?H4 d,[.=Jqv[ ylabel('\phi')
sj a;NL *}R5=r0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
;e;lPM{+ S!.aBAW axis([0 90 -3.5 3.5])
v:1l2Y)g > "?HbR9 grid on
& tQHxiDX ;5^ grr@,4 
