计算光波由光密入射到光疏介质时相位变化
bb@3%r|_< t>Ot)d clear
R9-Ps qmF 40M/Gu: close all
gO9\pI 2 s|q]11r+H n1=1.45,n2=1;
~T'!.^/ D.ajO^[ theta=0:0.1:90;
@pn<x"F5' ws_/F a=theta*pi/180;
gf]k@-) J+]W*?m rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
$+%e Lx* +{7/+Zz rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
@D3|Ak 1 K'DRX85F tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Y/Dah* Kh Wt9=9 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
?T ( @< T U~nW>WJ+. arp=angle(rp);
.! ><qV g k2@|fe ars=angle(rs);
{~=Z%Cj2 Q o4l=o Y:' atp=angle(tp);
0Xe?{!@a ]NRQM8\ ats=angle(ts);
wpW3%r;9 =eB^( !M figure(1)
xBevf &t P E$?:^ausu subplot(1,2,1);
zmEg 4 v'I :OFL@byS plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
kH9fK80 4!,`|W1 legend('arg(r_p)','arg(r_s)')
!I UH 5 {YT@$K]w, xlabel('\theta_i')
Lfj]Y~*z ,!{/Y7PmJ ylabel('\phi')
TQjM3Ri=V 8h=Rfa9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Obl,Qa:5 0.4Q-?J axis([0 90 -3.5 3.5])
ta)'z@V @g YYzl"<)c grid on
Pd^v - }[ 9[7Gxmf subplot(1,2,2);
hOTqbd} {rE]y C^ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
E5EAk6 017(I:V?(: legend('arg(t_p)','arg(t_s)')
< lxE^M Dg_ AoC xlabel('\theta_i')
b+DBz }L4
}@'Zt6+tS ylabel('\phi')
I" Gr <?r /"LcW"2;N title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
k5X& |L/ Kz`g Q |S axis([0 90 -3.5 3.5])
{;0+N -U 0**.:K <i grid on
E qc,/ {WYHT6Z 
