计算光波由光密入射到光疏介质时相位变化
U(f @zGV tKt} ]KHV clear
sg,\!' "zcAYg^U close all
zdwQpB,+^ rT ` sY n1=1.45,n2=1;
FXr^ 4B} =Q 9^|& 6 theta=0:0.1:90;
G %s O{k7 K F 'fg
R a=theta*pi/180;
qefp3&ls Z@euO ~e~ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
{Z/iYHv~#c 'f{13-# X@ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
II Ap-Y~B t#.}0Te7 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
k\O<pG[U M1eh4IVE? ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
) ' xyK ?>+uO0*S arp=angle(rp);
>IS4 1T#-1n%[k( ars=angle(rs);
Ze%S<xT!O 2;sTSGD G atp=angle(tp);
U1:m=!S;x o*204BGB ats=angle(ts);
rS>.!DiYr, jP<6J( figure(1)
p^Ey6,!8]D diNSF-wi,, subplot(1,2,1);
P1OYS\ P.>fkO1\ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
oq vu8" _i ztQ78 legend('arg(r_p)','arg(r_s)')
)$p36dWl Ia%cc
L= xlabel('\theta_i')
Vb? wwx7= |\Gkhi>; ylabel('\phi')
j?C[ids< P>q~ocq< title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
A@~9r9Uf z2V ->UK) axis([0 90 -3.5 3.5])
Srw`vql{( `}t5` :#k grid on
1;m?:|6K{ kR1
12 J9P subplot(1,2,2);
c@,1?q1bv . ?#Q(eLj plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
`%|3c CHS}tCfos> legend('arg(t_p)','arg(t_s)')
~Q"qz<WO L ntRLB' xlabel('\theta_i')
Ox
,Rk R[j'< gd. ylabel('\phi')
[|$C 2Dhw= DBy%"/c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ih("`//nP sz%]rN6$ axis([0 90 -3.5 3.5])
@ [FO;4w UK'8cz9 grid on
i*l =xW;bM -c8h!.Q$ 
