计算光波由光密入射到光疏介质时相位变化
B L^?1x Z@sDxYt9 clear
*BdKQ/Dk naOCa close all
zH)M ,+P *m$lAWB5D n1=1.45,n2=1;
)( CZK&< T- ~l2u|s theta=0:0.1:90;
Yf&P|Iiw gdr"34%vbM a=theta*pi/180;
$f>h_8cla 4{z z-4= rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
cJ^{iOQ+ I j w{g% rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
NQ&\t[R[ *h([ai"1- tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
E@}
NV|90 S Lsw '< ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
"/ 9EUbca Fi!BXngbd arp=angle(rp);
|^:qJ;dO P -/w#f&Y+]8 ars=angle(rs);
Sj(5xa[ .Tm m atp=angle(tp);
^ osXM ` WW B m*?U ats=angle(ts);
Acix `-< Vf*Z }' figure(1)
F9}j iCom `AcUxnO subplot(1,2,1);
W gTD
O3 MJ?fMR@ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
Z~S%|{&Br a]S0|\BkN legend('arg(r_p)','arg(r_s)')
:z-UnC||j * a ?qV xlabel('\theta_i')
xCTPsw]s [C-4*qOaa2 ylabel('\phi')
,%=SO 82W )hy(0 D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
_GbwyfA
n# gfr+`4H >v axis([0 90 -3.5 3.5])
1-=ZIHW Y W9+.Dc` grid on
jL6ZHEi#d7 iVSN>APe subplot(1,2,2);
:5W 8S6[o t@vVE{` plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
G(;hJ'LT T1*%]6&V | legend('arg(t_p)','arg(t_s)')
iwVsq_[] L ABaK60.O[O xlabel('\theta_i')
? h |& kRq :3N&&] ylabel('\phi')
Abc%VRsT @ ,^c?v title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
1Qk]?R/DN CqVeR';2 axis([0 90 -3.5 3.5])
%
|^V ) ;R[w}#Sm grid on
tv 7"4$T k}&7!G@T 
