计算光波由光密入射到光疏介质时相位变化
;9z|rWsF >(a35 b$ clear
2QaE&8vW `r1}:`.m, close all
gP
QOv vz4(
k/ n1=1.45,n2=1;
>M2~p& Si BQ Pmo1B theta=0:0.1:90;
uluAqDz` ]4H)GWHKg a=theta*pi/180;
` qhT $4g {4-) rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
DK?aFSf\ Y
]()v rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
x [{\Aw>$. LH4#p%Pb% tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
?]D &D:Z?I -b^dK)wR~ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
ly`
A,dh iG54 +] arp=angle(rp);
Qg/FFn^Kg* M6p\ QKi ars=angle(rs);
s_y8+BJaV _\\Al v. atp=angle(tp);
MBt\"b# t As46:< !2 ats=angle(ts);
eX#.Zt] Ext C\(X; figure(1)
rn)Gx2 5 &/U fXK r subplot(1,2,1);
^o5;><S] 4:rwzRDY plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
oH=?1~ e L-`V^{R] legend('arg(r_p)','arg(r_s)')
sA? 8i:]O: ZJ%NZAxy xlabel('\theta_i')
2|bt"y-5r < ?B3^z$ ylabel('\phi')
l[i4\ CT '=$`NG8 l title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:-jbIpj' n8Qv8 axis([0 90 -3.5 3.5])
3 zh:~w_ B$rhsK % grid on
E& cC2(w f>6{tI 5X subplot(1,2,2);
eXKEx4rU *xxk70Cb plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
@LOfqQ$FE m"~ddqSMT legend('arg(t_p)','arg(t_s)')
}=EJM7sM|k } iKjef#J xlabel('\theta_i')
:7 LA /j sf2%WPK
ylabel('\phi')
By@65KmR" gf>H-718F title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2!-Q!c`y Mbxl{M
> axis([0 90 -3.5 3.5])
mQ`atFz:Z S'HA] grid on
R*bx&..< !Np7mv\7 
