计算光波由光密入射到光疏介质时相位变化
/-4%ug tD$ #vS>^OyP clear
jliKMd<? (iS94}-) close all
e 94csTh= Y+G4 : n1=1.45,n2=1;
2+?M(=4 xSLN theta=0:0.1:90;
\{~x< <qFd i.by Hz?/ a=theta*pi/180;
WnIh (
0 ].1R~7b rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
;CU3CLn dB1bf2'b# rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
L:%h]- q& P" tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
zfxxPL' vwT?Bp ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
s=KA(4p F!Nx^M 1 arp=angle(rp);
&Wn!W U :IQWl C ars=angle(rs);
+i
K.+B HF0J>Clq atp=angle(tp);
UDxfS4yI KWd]?e) ats=angle(ts);
Xf Y]qQP ]i{-@Ven figure(1)
$osDw1C t4 aa5@r subplot(1,2,1);
, {]>U'- ="XxS|Mq3 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
==Y^~ab;K rVZk G,Q legend('arg(r_p)','arg(r_s)')
G `TO[p]q [Si`pPvl xlabel('\theta_i')
)C <sj E pPKo ylabel('\phi')
[dUW3}APV ,@GI3bl title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/y
NU0/ %"{SG p axis([0 90 -3.5 3.5])
! 5 ]/2 6LUB3;g7 grid on
M< Eg<* l$z-' subplot(1,2,2);
6 R"& !.ZF HbV[L)zYG plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
%/~Sq?f-9@ RD,` D! legend('arg(t_p)','arg(t_s)')
wl}Q|4rZ +AXui|m n xlabel('\theta_i')
6$`8y,TMSt * y u|]T ylabel('\phi')
X(N!y"z OBu$T& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
i]YH"t8GY \iO
,y: axis([0 90 -3.5 3.5])
VYik#n>|Gp d TS 7l02 grid on
$FS
j^v] I+ydVj(Op 
