计算光波由光密入射到光疏介质时相位变化
ND*5pRzvp wuCODz@~ clear
\Uiw:
, $1YnQgpT close all
S3w? X n7vi@^lf( n1=1.45,n2=1;
uP:'e8 gueCP+a_ theta=0:0.1:90;
pB`<4+"9 Flsf5 Tr0 a=theta*pi/180;
ZC"p^~U_e[ ~9=g" v rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
0:nyOx(; Nsb13mlY rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
'tekne /OEj]DNY tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
S:wmm}XQ p)`JVq,H/B ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
G9;WO* :7gIm|2"] arp=angle(rp);
idHBz*3~ps SyI#Q[f'_ ars=angle(rs);
o 0B`~7( zrcSPh atp=angle(tp);
6v-h!1p{u UWz<~Vy ats=angle(ts);
HjrCX>v umY4tNe]$ figure(1)
5^qI6
U GKX#-zsh79 subplot(1,2,1);
$4m{g"xL #-{4F?DA]y plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
+]G;_/[2 c8h
9 legend('arg(r_p)','arg(r_s)')
V<b"jCXI 72aj4k]^ xlabel('\theta_i')
xGjEEBL uTQ/_$
ylabel('\phi')
8D-g%Aj- 8d-_'MXk3 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]O^!P,l)"
W<t,Ivg axis([0 90 -3.5 3.5])
(fC [Y :2La, grid on
f -bVcWI 6u7>S? subplot(1,2,2);
F[LBQI`zq fU,sn5zZ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
s{!F@^a dI3U*:$X
legend('arg(t_p)','arg(t_s)')
#8cY,%<S] ydw')Em xlabel('\theta_i')
IL\#!|> p tMysYT' ylabel('\phi')
.- {B o@}Jd0D4 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:*1bhk8~ }{s<!b axis([0 90 -3.5 3.5])
7^=O^!sa uGOvZO^v grid on
YoJN.],gf 0]&~ddL