计算光波由光密入射到光疏介质时相位变化
i]YV { /'l{E clear
i1NY9br irFMmI b close all
lm-dW'7& w[zjerH3 n1=1.45,n2=1;
v1+3}5b'uF IEsEdw]aZE theta=0:0.1:90;
59Xi3KY :t8b3 9 a=theta*pi/180;
}@'$b<!B cdl&9-} rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
k;AD`7(= dJeN bVd rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Ln') QN v &Yi tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
.w=/+TA aQ?/%\> ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
hV8[@&Sx3 =.f-w0V arp=angle(rp);
5cL83F Qh |:q=T
~x ars=angle(rs);
e6R}0w~G (C-{B[Y atp=angle(tp);
)t0$qd ] 42{Ew8 ats=angle(ts);
\o}xF@sM5 V:In>u$QJ! figure(1)
\qdHX :Dr Wq{4 subplot(1,2,1);
@5=oeOg36 I5 o)_nc plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
+7D|4 QR c{vUR& legend('arg(r_p)','arg(r_s)')
a
oD`=I*< [')m|u~FS4 xlabel('\theta_i')
j Sh5!6O J*yf2& lI5 ylabel('\phi')
u)X=Qm) H[s(e5 6z title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ck ]Do!h ~
-4{B axis([0 90 -3.5 3.5])
YSr9VpqWV 8;b( 0^ grid on
gn8R[5:!V #|Y5,a ,{ subplot(1,2,2);
[w>$QR 3 KR d plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
]
bM)t< YIn',]p: legend('arg(t_p)','arg(t_s)')
?{P"O!I{ g`{;(/M+ xlabel('\theta_i')
kw >v:F<M R kH W
ylabel('\phi')
V]|^&A _c c(3~0Yr title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
R0P
iv: ~bM4[*Q7 axis([0 90 -3.5 3.5])
CY
4gSe? s>B5l2Q4 grid on
nHnK)9\ N prC;L*~8 
