计算光波由光密入射到光疏介质时相位变化
C3Q[L}X\ m'k.R
j clear
WY#A9i5Ge HDZB)'I close all
o CdWf63D ~YKe:K+&z n1=1.45,n2=1;
BpZE :zZK%} G< theta=0:0.1:90;
TaeN?jc5 L O<R<zz a=theta*pi/180;
+uKh]RP RUUV" y rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
c1M/:*?% fO CLN$x^ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
lN#W I~?D^ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
(:Rj:8{ wgxr8;8`q ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
T;qP"KWZ SYRr| Lg arp=angle(rp);
R!b<Sg VV O C- : ars=angle(rs);
_Q,`Qn@|BD NrT!&>M atp=angle(tp);
dO}6zQ\ Y~#m-y ats=angle(ts);
ec+&K ?T P2pdXNV figure(1)
9`Bmop .6aC2A]es subplot(1,2,1);
@igr~hJ <dl:';@a- plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
reJw&t }Q ^.)oQo SE legend('arg(r_p)','arg(r_s)')
Bee`Pp2 Jt5 \ xlabel('\theta_i')
@dei} !e m/uBM6SXx ylabel('\phi')
;"46H'>! 2\80S[f title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7{>mm$^|V @ -CZa^g axis([0 90 -3.5 3.5])
_ K["qm{X_ VWf&F`^B( grid on
<y[LdB/a =vJ:R[Ilw subplot(1,2,2);
S?ELFq(g V
.+ mK|) plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
f Ne9as *P2_l
Q= legend('arg(t_p)','arg(t_s)')
Z6C!-a XE8>& & X xlabel('\theta_i')
KSkT6 _< j=b?WNK ylabel('\phi')
ScOiOz:Ha -P#PyZEH&I title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
z6
T3vw L kYcAY$w axis([0 90 -3.5 3.5])
}2c)UQD8
+l' l*< grid on
Kv(2x3(" [Z~h!} 
