计算光波由光密入射到光疏介质时相位变化
+R\vgE68 I0'[!kBF| clear
wBInq~ K_ [Cb` { close all
90#
;?# {w8 NN-n n1=1.45,n2=1;
'Vr$MaO d7&eLLx theta=0:0.1:90;
GCttXAto $%%os6y2v a=theta*pi/180;
SR8qt z/V YqV8D&I rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Y(+^;Y3U x%< rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
.p ls! ?hR7<0 2 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
.,i(2^ 4_D@ST% ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
`s]4AKBO y?)}8T^ arp=angle(rp);
$v`afd y ]u-0 2g ars=angle(rs);
pCu!l#J Ln.ZVMZ; atp=angle(tp);
m$ LVCB K T.?X p:z ats=angle(ts);
NJ MJ @O}7XRJ_8 figure(1)
/?6gdN 8*SP~q subplot(1,2,1);
JTqq0OD} EQe5JFR plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
m))<!3 vNW jH!' legend('arg(r_p)','arg(r_s)')
|3 {&@7 fRvAKz|rL xlabel('\theta_i')
P0i V<T4^ ZCVl5R(mZ ylabel('\phi')
a1.|X i'/z mF09U(c i title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3p*-tBOO -z: &*= axis([0 90 -3.5 3.5])
[,~TaP}m 7dX/bzUVz8 grid on
#g5^SR|qE m9md| yS subplot(1,2,2);
+<.\5+ q+a.G2S plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
kL S(w??T e# /kNHl legend('arg(t_p)','arg(t_s)')
_}mK!_` 4.K'\S xlabel('\theta_i')
ue 8qIZH jko"MfJ ylabel('\phi')
?`zgq>R}w[ 3|rn] yZ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Fvxu >BK >?YNW axis([0 90 -3.5 3.5])
3 ,);0@I Ze !92 g grid on
BwJuYH7QJ$ m= b~i^@ 
