计算光波由光密入射到光疏介质时相位变化
f` J|>Vk G/(*foT8SE clear
)E~ _rDTl ppFYc\&= close all
.[,6JU% @'gl~J7 n1=1.45,n2=1;
j{"[Ec ^ $wJi9D6 theta=0:0.1:90;
0[9I0YBJ /{R3@,D[] a=theta*pi/180;
oZ-FF' A$<>JVv rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
lR}%)3_k @G(xaU'u rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
\k4pK &b @'HT;Q!\Vd tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Q"'V9m7
i VF= Z` ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
HtUG#sc&`{ "& Mou arp=angle(rp);
G
* @@K }9=2g`2Q ars=angle(rs);
Yw[{b eo ^9'$Oa,* atp=angle(tp);
r w\D>} \ ngLJ@TP- ats=angle(ts);
4Cd#sQ |!xpYT: figure(1)
Rn;VP:H M Vrg3{@$ subplot(1,2,1);
1 KB7yG-#6 \n;g2/VjO plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
'z-D%sCA Wvb Eh|y legend('arg(r_p)','arg(r_s)')
Wxs>o sq GmAj< /~ xlabel('\theta_i')
6e/7'TYwT E\iJ P^n ylabel('\phi')
6y9C@5p}B 4(8<w cL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
9fMSAB+c% {NV:|M ! axis([0 90 -3.5 3.5])
ssT@<Tk^4 @`Wt4< grid on
J[6VBM.Y c"qPT jY subplot(1,2,2);
6J"(xT )^";BVY plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
wn1,
EhHt Mlwdha0 legend('arg(t_p)','arg(t_s)')
ke^d8Z. ,S0UY):( A xlabel('\theta_i')
?DRR+n _ =+4 _j ylabel('\phi')
/:KQAM0 jOv~!7T title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:~i+tD Tm\OYYyk axis([0 90 -3.5 3.5])
=R2l3-HA= F: ,#? grid on
*73AAA5LKa kJ__:rS(T_ 
