计算光波由光密入射到光疏介质时相位变化
;cIs$ $d.UF!s clear
-%"Kxe dC;@ Fn close all
?fcQd6-} n?vrsqmZ n1=1.45,n2=1;
RajzH2j+> 0 }jB/Z_T theta=0:0.1:90;
eZs34${fN gS ~QlW V a=theta*pi/180;
Ak[} s|,) 09pnM|8A rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
SN O'* ? k_.%(ZE rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
bCHA!zO mu0ER 3o tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
;_)&#X,?( <)+9PV<w ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
n8 #i L `~QS3zq arp=angle(rp);
+s.r!?49+ `qp[x%7^ ars=angle(rs);
dX-Xz g }7E2,A9_" atp=angle(tp);
9/TF # <7cm[ ats=angle(ts);
hl?G _%a Zo}wzY~x>I figure(1)
GZQy~Uk~ >8+:{NW subplot(1,2,1);
yZ=wT,Y hb'S !N5 m plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
x;2tmof=L uFNVV;~RFI legend('arg(r_p)','arg(r_s)')
&wr0HrE\ $+qJ#0OE$ xlabel('\theta_i')
8C3oj I5PI;t + ylabel('\phi')
'v iF8?_ 8V@ /h6-e, title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:hB/|H*= <5
G+(vP axis([0 90 -3.5 3.5])
)J#7:s]eo :786Z,') grid on
Kf#9-.}? +m kub}<a subplot(1,2,2);
@DG $ fvg jqiT plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
C5GO?X2 G$&SlJZEk legend('arg(t_p)','arg(t_s)')
zI77#AUM d[ (KgX9 xlabel('\theta_i')
_Vr>/f &y"e|aE ylabel('\phi')
R@ Ch3l@ -E_lwK title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<hazr KUn ~Qj f-| axis([0 90 -3.5 3.5])
jy#'o adS? s
8O"U% grid on
@ *'$QD, [;yEG$)K 
