计算光波由光密入射到光疏介质时相位变化
2]y Hxo/6 U2)?[C1q{ clear
]VR79l 7z_ZD0PxPc close all
Fh/ psd I<Cm$8O? n1=1.45,n2=1;
T/0cPn0> :%gM
Xsb theta=0:0.1:90;
PWeWz(]0Z4 t3 3\f<e a=theta*pi/180;
O2":)zU. *~~J1.ja> rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
I s|_ Ey.%:
O-Dv rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Scug
wSB X(O:y^sX} tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
a ]:xsJ~ _%3p&1ld ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
c'XSs i%GiWanG arp=angle(rp);
2%v6h gu VuO ars=angle(rs);
fRxn,HyV n2dOCntN> atp=angle(tp);
<00nu'Ex1v g:.LCF ats=angle(ts);
qC=9m[MI 9@?|rj e9 figure(1)
nXk9
IG( $> ;| subplot(1,2,1);
nI4oQE :L6%57 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
?^J%S, :fDzMD legend('arg(r_p)','arg(r_s)')
M;W&#Fz% M1]w0~G xlabel('\theta_i')
i03= Af3 GDs/U1[* ylabel('\phi')
O 3?^P"C lKf kRyO_S title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7L!}F;yT mhM;` d l axis([0 90 -3.5 3.5])
wz@[rMf >Hmho' grid on
j+>[~c;0) t\]kVo) subplot(1,2,2);
W4qnXD1n fLeHn,*," plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
c}!`tBTm 2"k|IHs1 legend('arg(t_p)','arg(t_s)')
RameaFX8 3DxgfP%n xlabel('\theta_i')
|9F-ZH~6 aO}p"-' ylabel('\phi')
e\O625 9':Hh' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`9k\~D=D~ unRFcjEa axis([0 90 -3.5 3.5])
\II^&xSF F&#I[]# grid on
a^^OI|? dQFUQ 
