计算光波由光密入射到光疏介质时相位变化
{Yv
|C)O c-z=(Z clear
G sE ?<3 / >.& close all
m <aM b D Eqk9Exk` n1=1.45,n2=1;
&> }MoB =@w};e#D theta=0:0.1:90;
S 3btx9y{ YVD%GJ a=theta*pi/180;
rS)7 D - stSl* rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
6L'cD1p u ~8}"X] 4 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
~O |j*T At +on9&= tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
%Xkynso~ y)Ip\.KV\ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
i|.!*/qF \#G`$JD arp=angle(rp);
$5%tGF h Ya304Pjd ars=angle(rs);
T-f+<Cxf AUzJ:([V atp=angle(tp);
'00DUUa .Uha %~% ats=angle(ts);
&{ntx~ Eq yd# 4b`8U` figure(1)
?8 YHz x\I9J4Q subplot(1,2,1);
q\d'}:kfu oV,>u5:B plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
pd>EUdbrp& n3~xiQ' legend('arg(r_p)','arg(r_s)')
)k0bP1oGS V u;tU. xlabel('\theta_i')
(O /hu3 |Z# ) 1K ylabel('\phi')
*kZJ |F9/7 z\5+ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+% '0; T~`m'4"+c axis([0 90 -3.5 3.5])
AP/tBC eM 6i=m1Yk grid on
gLd3,$ Ei X(g<rz1J] subplot(1,2,2);
5x4(5c5^ ;h0?o*i_ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
,f} s!>j .CwMxuW legend('arg(t_p)','arg(t_s)')
7}TjOWC
yZdM4` xlabel('\theta_i')
1=W>zC i#(T?=VPcy ylabel('\phi')
C O 5?UgA
%|l *=v title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
rf^1%Zo: ]>*Z 1g; axis([0 90 -3.5 3.5])
|w:7).P `Z/"Dd;F^3 grid on
A3M)yW q 6ZCt xs! 
