计算光波由光密入射到光疏介质时相位变化
,dMi+c`ax /! G0 g%k clear
8?pZZtad [6N39G$ close all
xF+x I6 bvTkS EN n1=1.45,n2=1;
. Dxr c 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_B 5@9e# vLS 6Gb't rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
}e pN<DL 4`]1W,t tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
m q9&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);
4ebGA g ?_ ;Yj&7k1 figure(1)
d a9 *>+[ i[J', subplot(1,2,1);
+[+ Jd)Z RL9 P:]
^ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
|J Q05nb 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"`4 6 0kfw8L on subplot(1,2,2);
u :$x ,Q B %/Pn
2 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
$'Z!Y;U e i`;I"o Y4 legend('arg(t_p)','arg(t_s)')
lvlH5Fc nFSa~M xlabel('\theta_i')
lL v0lf DS fK Ux& ylabel('\phi')
arJ[ .f9s 8:bNFgJD title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Mtr~d 19 _F\3 2 axis([0 90 -3.5 3.5])
qMkP/BjV l{4\Wn Va grid on
&~8} y+z y:L|]p}huE 
