计算光波由光密入射到光疏介质时相位变化
U]tbV<m % k~Qb"6n2 clear
{WQH _PD RUJ close all
4A6D>ChB'E gKy@$at& n1=1.45,n2=1;
)nmL gsg q"xIW0 Pc theta=0:0.1:90;
h fT HP 35 I y\ a=theta*pi/180;
}1E'a>^| Y [Jt+p ] rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
/cHd&i,> t)P5bQ+$u9 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Tap=K|b ]
Q;{[U!\ : tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
],CJSA!5F ;AMbo`YK[ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
eAKK uML ,CqJ (( arp=angle(rp);
CA%p^ 4Q e FDh J ars=angle(rs);
.J :04t1 y+xw`gR: atp=angle(tp);
V^fSrW] >^}nk04 ats=angle(ts);
Ncr38~;w a)J3=Z- figure(1)
6(Za}H B;xw @:H subplot(1,2,1);
?w'a^+H 4/ YEkD plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
W:D' k^u @V{s'V legend('arg(r_p)','arg(r_s)')
AZ'
"M{wiI cpz'upVOZ xlabel('\theta_i')
`L p3snS T \%{zz_( ylabel('\phi')
%qA@)u53 JTbg8b title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
&"GHD{ix {TT@Mkz_QC axis([0 90 -3.5 3.5])
)6mx\t k^3 ?Z2a grid on
?^]29p_ t}m6]; subplot(1,2,2);
0tbximmDb me] O plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
iC-WQkQY K..L8# SC legend('arg(t_p)','arg(t_s)')
DVCO(
fz Md a~@)7$ xlabel('\theta_i')
bg.f'; C tI50z khaB ylabel('\phi')
1Xy{&Ut\ :NB|r title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
hKq <e%oVH Z}bUvr XP axis([0 90 -3.5 3.5])
+zodkB~) vd 7%#sHH& grid on
V[D[MZ Q*wub9 
