计算光波由光密入射到光疏介质时相位变化
,dMi+c`ax /! G0 g%k clear
8?pZZtad [6N39G$ close all
xF+x I6 bvTkSEN n1=1.45,n2=1;
.Dxrc Vp4] theta=0:0.1:90;
Bn:sN_N C}E
ea~ a=theta*pi/180;
nEr, jd~f ?C\9lLX rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
d_B5@9e# vLS6Gb't rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
}epN<DL 4`]1W,t tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
mq9&To! .?:~s8kB ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
xc3Q7u!| mn` Ae= arp=angle(rp);
?!Gt.
fb ![jP)WgF ars=angle(rs);
nJC}wh2d# aQ\O ]gCE atp=angle(tp);
Vd~k4 4{E=wg^p ats=angle(ts);
4ebGAg ?_ ;Yj&7k1 figure(1)
d a9 *>+[ i[J', subplot(1,2,1);
+[+Jd)Z RL9P:]
^ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
|JQ05nb a}\JA`5;)Z legend('arg(r_p)','arg(r_s)')
=h{2!Ah7
X >O=V1 xlabel('\theta_i')
K2Zy6lGOZ |{Q,,<C ylabel('\phi')
^;bkU|(`6 24fWj?A| ^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+a;j>hh 9
|Y?#oZ1 axis([0 90 -3.5 3.5])
>sq9c/}X K.~U%v} grid on
mH"`46 0kfw8Lon subplot(1,2,2);
u:$x,Q B%/Pn
2 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
$'Z!Y;Ue i`;I"oY4 legend('arg(t_p)','arg(t_s)')
lvlH5Fc nFSa~M xlabel('\theta_i')
lLv0lf DS fKUx& ylabel('\phi')
arJ[.f9s 8:bNFgJD title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Mtr~d 19 _F\32 axis([0 90 -3.5 3.5])
qMkP/BjV l{4\Wn Va grid on
&~8}y+z y:L|]p}huE