计算光波由光密入射到光疏介质时相位变化
B/G d(S`@q W>M~Sk$v clear
N~pIC2Woo }X;U|]d close all
+%N
KQ'49I Pv<FLo%u< n1=1.45,n2=1;
o{ *ay$vA] )'
xETA theta=0:0.1:90;
=2y8 CgLj pium$4l2# a=theta*pi/180;
x'Pi5NR E kC UT ^ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
19(Dj&x
u@ p? rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
9kKnAf4Z 7%OKH<i\2< tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
qQ7w&9r.M z"DkFvA ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
)a=/ 8ofe i g 0u^BC arp=angle(rp);
>EXb|vw
]WJfgN4 ars=angle(rs);
'f/Lv@]a PiLJZBUv atp=angle(tp);
%HtuR2#ca |m,VTViv;i ats=angle(ts);
^TXf sQs R*1kR|*_) figure(1)
j1Yq5`ia ,]Z p+>{
subplot(1,2,1);
#>@z
2K7 C9Bh@v%90^ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
^ _p%Yv =A[5 =
k> legend('arg(r_p)','arg(r_s)')
Qq 0l* )mX 2aNT#J"_ xlabel('\theta_i')
AB.(CS=i v7trr W} ylabel('\phi')
<c,~aq#W' 77gysd\( title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
vxN0,l j<tq1?? [b axis([0 90 -3.5 3.5])
z&nZ <ih
ARh6V&Hi- grid on
7z~Ghz Z& !!]"I subplot(1,2,2);
=G-N`
39 FE5Q?*Ea plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Zu<]bv w,.qCp T$_ legend('arg(t_p)','arg(t_s)')
O>DNC-m)i{ E\)eu1Hw4B xlabel('\theta_i')
wpPn}[a SKXBrD=- ylabel('\phi')
=N.!k Vkl qFR dg V>8 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
?~ULIO ' 2%rLoL$Y2+ axis([0 90 -3.5 3.5])
M<8ML!N0;t 8IY19>4'5J grid on
BQgoVnQo_c L hJ a)jFQ 
