计算光波由光密入射到光疏介质时相位变化
noz1W ] * Wp?0CP clear
.)bNi*& Qa:[iF close all
='+I dn# 5 :r
"G Z n1=1.45,n2=1;
:/$_eg0A W{Z 7= theta=0:0.1:90;
=w`uZ;l$Q ?VRsgV '$ a=theta*pi/180;
:HrFbq \k"C t zoX rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
uF}B:5 3A bd% <
Jg+ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
a -F I`Dv aqQ+A:g tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
|dqESl,2 T2 rBH]5 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
(@ !K tW ;34p
[RT arp=angle(rp);
/|H9Gm ]s)Y">6 ars=angle(rs);
?GhMGpd Mq 8L_OH atp=angle(tp);
*pnaj\ W4k$m 2 ats=angle(ts);
84e8z { 3< 6h~ek ) figure(1)
*ej< 0I { f 9kd&#O& subplot(1,2,1);
'PFjZGaKR W| zPV` plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
o^"OKHU,S0 +Q);t, legend('arg(r_p)','arg(r_s)')
kF ,ME5% H7Z`a QC xlabel('\theta_i')
qbS6#7D u=]*,,5< ylabel('\phi')
q I~ *G3 -Hw3rv3o title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
5|pF*8* :wg =H axis([0 90 -3.5 3.5])
8 qw{e`c ,~1k:>njY~ grid on
_Ds,91<muQ /2/aMF(J subplot(1,2,2);
bE2O[B oUN\tOiS+ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
a.?U $F lP]Y ^Gz legend('arg(t_p)','arg(t_s)')
V'Kie d+ h. ftl2> xlabel('\theta_i')
| V{ Q MVvB d 3 ylabel('\phi')
I"lzOD; eI F$'u` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
$2i@@#g8 +0Q +0: axis([0 90 -3.5 3.5])
l ic-68T e`7>QS ;. grid on
,5} w]6bCr X;)/<:mX 
