计算光波由光密入射到光疏介质时相位变化
N[Sb#w`[/ X=~V6m clear
h qmSE '8 =Q(vni83< close all
@Bs0Avj . u3ZCT" ! n1=1.45,n2=1;
$/_ qE "tu*YNP\Q theta=0:0.1:90;
&~-~5B|3" ^# e~g/ a=theta*pi/180;
op7FZHs } /e`v6 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
]gEhE H4<Nnd\ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*P2[qhP2 Qw)9r{f tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
K: r\ {# 9 ~kKrDLW+ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
{6
.o=EyM{ kzO& 24 arp=angle(rp);
ULBg {e?l8 A6UO0lyu ars=angle(rs);
RuAlB* .ve *Vp atp=angle(tp);
RT*5d;l0 !.Zt[ g} ats=angle(ts);
-5>NE35Cto l1 +[ figure(1)
`/Y+1 aD Gs V4ZZ subplot(1,2,1);
<@, $hso7: E 9 Y\X plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
tRfm+hqRZ 4]XI"-M^D legend('arg(r_p)','arg(r_s)')
6)*xU|fU >nn Y:7m xlabel('\theta_i')
3071:W X
K>&$<5{ ylabel('\phi')
G '#41>q+ 3/goCg title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
k#)Ad*t =-jD~rN4;P axis([0 90 -3.5 3.5])
8 cVzFFQP V@ :20 m grid on
8+'C_t/0i z, f=}t[.Y subplot(1,2,2);
cT' w= P-Su5F plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
E{Vo'!LY SUdm 0y legend('arg(t_p)','arg(t_s)')
RKkGITDk K|^wc$ xlabel('\theta_i')
Ruaur] U;^ [$Aq ylabel('\phi')
f7<pEGb "{BqtU*. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ax<\jW< pGWA\}' axis([0 90 -3.5 3.5])
R p.W,)i f_6`tq m % grid on
]]uHM}l [ygF0-3ND 
