计算光波由光密入射到光疏介质时相位变化
VqS1n x5k6"S"1, clear
Q
822 # 8COGe=+o close all
5s<.qDc " t,Z O n1=1.45,n2=1;
X]s="^ <W/YC 2b theta=0:0.1:90;
A bB+<0 _+<AxE 9\ a=theta*pi/180;
0+k= gO +<3e@s& rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
:JK+V2B$H Dk}txw}# rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
)H{OqZZYD nX<yB9bXDg tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
\cQ+ 9e) `m\ ?gsw7 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
dZ Ab' : (]$&.gE.F arp=angle(rp);
Zig3WiD& /Kh Y,G'Z ars=angle(rs);
v>5TTL~? !X1
KOG atp=angle(tp);
Lt {&v ^y MpJ] 1 ats=angle(ts);
JQSczE3 Hqb-)8 ~ figure(1)
tY60~@YO& I9YMxf>nI subplot(1,2,1);
d V3R) o:@A% *jg plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
G4"n`89LK Hm_&``=' legend('arg(r_p)','arg(r_s)')
p e$WSS J %Z yt ;p2 xlabel('\theta_i')
.19_E Q >+ UM. Se(kS ylabel('\phi')
Pk:b:(4 D\ P-|} title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
-_f-j 2K 2_- axis([0 90 -3.5 3.5])
>n5Kz]]% I} \ `l+ grid on
1+Bj` ACP l(~NpT{=V subplot(1,2,2);
:D:J_ {HJ p 9eTrFDy? plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Wbi12{C XpdjWLO]C< legend('arg(t_p)','arg(t_s)')
Jg@eGs\* dbfI!4 xlabel('\theta_i')
kj`h{Wc[) F ZfhiIf ylabel('\phi')
vcSb:(' {QI deB[ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
LP} j 0)n '^BTa6W}m axis([0 90 -3.5 3.5])
<KZ J ;z9U_ grid on
4x" je 5- "aK~@+ 
