计算光波由光密入射到光疏介质时相位变化
Ac<V!v71 $ 6atr-Pb clear
2H fP$ . ],R rk]1 close all
Jis{k$4 Rj9M E,u n1=1.45,n2=1;
=UA- &x@ {'G@- +K theta=0:0.1:90;
Hk8:7"4Q bmna*!l ^M a=theta*pi/180;
i>r4R z! 9 a2Ga rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
YFOK%7K )zI<C=])" rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
b")&"o)G2W p #'BV'0bl tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
t*y4)I !gR 1 9a"@WB@ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
~}w(YQy=y uF9p:FvN 8 arp=angle(rp);
q rO] t\ c3&F\3 ars=angle(rs);
ja70w:j a d|]F^DDuI atp=angle(tp);
r Y|'<$wvg ?mS798=f ats=angle(ts);
A{\7 HV 5 PTIC2 figure(1)
h=YY>
x Hkg^ subplot(1,2,1);
z,}1K! *aS|4M- plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
p=\DZU~1 -w=rNlj legend('arg(r_p)','arg(r_s)')
|uV1S^ !A C\dQ6(3}\ xlabel('\theta_i')
k!t5>kPSQ ~el#pf~ ylabel('\phi')
@ )m[: n D4G*K*z,w4 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2{!'L'km n99:2 r_ axis([0 90 -3.5 3.5])
lYv : ^i&/ k grid on
_A kc7" zA9N<0[]o subplot(1,2,2);
AUcq\Ys 7xB#) o53 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
JM -Tp!C> 7!hL(k [ legend('arg(t_p)','arg(t_s)')
c)#b*k,lw< M:6H%6eT xlabel('\theta_i')
X< K[`
=I +che Lc ylabel('\phi')
_a8^AG IE: x&q`3 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
*58<.L| o DPs xw axis([0 90 -3.5 3.5])
%;^[WT`,
zN #$eyt grid on
N'Ywn}!js C"k8 M\RW? 
