计算光波由光密入射到光疏介质时相位变化
>| hqt8lY ]C>h_,EZc clear
<H| ]^An!H q^b12@.
close all
+Rq]_ sDu L] =LY n1=1.45,n2=1;
/?S,u,R }L X.gm theta=0:0.1:90;
cLIeo{H .FvIT] k- a=theta*pi/180;
fKa]F`p_h ($[@'?Z1 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
dDqT#N?Y $
3 R5p rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
6g "qwWZp 2l +t- tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
?d)FY B PBAQ
KQ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
`W u.wx [UB]vPXm$ arp=angle(rp);
&IFXU2t} #wF6Wx iG ars=angle(rs);
3MoVIf1 Rc~63![O. atp=angle(tp);
V/J- zH& df9$k0Fx ats=angle(ts);
da$ErN '{ }SG b`l figure(1)
/+f3jy:d 1P/4,D@ subplot(1,2,1);
1Pd2% m[A$Sp_"-h plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
3EyVoS6D { ?]& P legend('arg(r_p)','arg(r_s)')
a7e.Z9k! ?YnB:z*eV xlabel('\theta_i')
y {QF#&lW 8o,"G}Hjk ylabel('\phi')
<2fvEW/#v 7@ mP;K0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
O y y E0 Y NG S"3F axis([0 90 -3.5 3.5])
86~q pN <is%lx(GDX grid on
^j7]> I Wq*W+ 7=. subplot(1,2,2);
d0-T\\U +qmV|$rm M plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
%~qY\> mA6Nmq%{ F legend('arg(t_p)','arg(t_s)')
c 4L++
u# MW)=l
| G xlabel('\theta_i')
"ax"k0 E=l^&[dIl ylabel('\phi')
D<xDj#Z~1 e `Tssa+ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
qr? RU .W \ > axis([0 90 -3.5 3.5])
X=pPkgW i}Cy q grid on
{_]<mw d usI$ 
