计算光波由光密入射到光疏介质时相位变化
'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=80 s0 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<R7 J 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);
3Vc T7y*{P rWoe
?g ars=angle(rs);
OgEUq' ' S 9G8aea/ atp=angle(tp);
rDIhpT)a :%/\1$3P ats=angle(ts);
_kn]#^ucCe R=\v3 m 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 -&tr k xlabel('\theta_i')
Dp!;7e s| S\B5&W ylabel('\phi')
"3@KRb4f d O+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 
