计算光波由光密入射到光疏介质时相位变化
p7QZn.,=u .Y6v#VI clear
sU>IE To ^BA
I/WP close all
Bz/ba * +,,(8=5 g n1=1.45,n2=1;
lp(2"$nQ rt^~
I \V theta=0:0.1:90;
_eGYwBm SZH`-xb!+5 a=theta*pi/180;
~u&gU1}
{&+M.Xn rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
NF0_D1Goi 'I>USl3 hI rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
T +vo)9w x+cL(R tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
d.? }>jl NK q I x ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*QQeK# $s t*hy"e{*a arp=angle(rp);
=mXC,<] Y\9}LgIvr ars=angle(rs);
h^g0|p5 h/ n( atp=angle(tp);
loIb}8 N# }A9t ats=angle(ts);
+ydd"` 1 !_p
figure(1)
H$Kc~# = 'i% r subplot(1,2,1);
WkXgz6 P x|m9?[
!_ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
HQ@ g6 joI) 6c legend('arg(r_p)','arg(r_s)')
>Lo\?X ~ VgVDTWs7 xlabel('\theta_i')
~ vA{I%z5~ "}v.>L<P ylabel('\phi')
0Fb ];:a OTF/Pu$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'^ TeV= /K|(O^nw axis([0 90 -3.5 3.5])
p{5m5x sQ`G'<! grid on
@dv8 F
"v 2 {?]W/&fS subplot(1,2,2);
f|,Kh1{e @ZU$W9g plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
7_\F$bp` O2 >c|=# legend('arg(t_p)','arg(t_s)')
u{DE OhtI4 s $Vv xlabel('\theta_i')
+51heuu[o cTGd< ylabel('\phi')
36{GZDGQ Wu
0:X*>}p title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+vy fhw4 b7>'ARdbzX axis([0 90 -3.5 3.5])
['o ueOg Mw0>p5+ cy grid on
*,JE[M -8 ^qtB 
