计算光波由光密入射到光疏介质时相位变化
'je8k7`VA dqD;y#/ clear
*f `s%&Y]s 18.Y/nZAgQ close all
x[PEn a:yB%:2 n1=1.45,n2=1;
z [qdmx^ x/%7%_+' theta=0:0.1:90;
:'d76pM- t&R!5^R a=theta*pi/180;
+MNSZLP] QJ a4R rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
n5egKAgA gb=80s0 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
8Wdkztp/S GB<R7J tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
_\,rX\ (B>)2: T1 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
k;;nE o~6 >`hSye{ arp=angle(rp);
3VcT7y*{P rWoe
?g ars=angle(rs);
OgEUq'' S9G8aea/ atp=angle(tp);
rDIhpT)a :%/\1$3P ats=angle(ts);
_kn]#^ucCe R=\v3m figure(1)
5mI?pfm >"cr-LB subplot(1,2,1);
A$7Eo`Of ElQJ\% plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
OM20-KDc5 v[R_S legend('arg(r_p)','arg(r_s)')
e>
ar Q&u>7_, Du xlabel('\theta_i')
99F>n[5 M0O>Ljo4RN ylabel('\phi')
lir&e
9I+ 'rS'B.D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)UR1E?' [zY!'cz? axis([0 90 -3.5 3.5])
6RH/V:YY X$BXT grid on
%S{o5txo 8W-]t1O%! subplot(1,2,2);
?N4A9W9 &bB6}H( plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
\4OK!6LkI n<{aPLQ legend('arg(t_p)','arg(t_s)')
myD{sE2A -&trk xlabel('\theta_i')
Dp!;7e s| S\B5&W ylabel('\phi')
"3@KRb4f dO+kPC title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"VV914*z (.PmDBW axis([0 90 -3.5 3.5])
(F_w>w.h rwoF}} grid on
r k@UsHy DWuRJ