计算光波由光密入射到光疏介质时相位变化
h.\9a3B:r S"Q$ Ol" clear
0CYm%p8 ! D/jS4'$vA close all
p^LUy LG` Jk.Ec )w n1=1.45,n2=1;
L;},1
\ w:}RS.AK theta=0:0.1:90;
}b#KV?xgW qYMTud[Vf a=theta*pi/180;
olC@nQ1c* <*Kj7o{Qn rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Ef;OrE"" |7jUf$Q\p rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
!2('Cq_^ +^c;4-X
0 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Y dgaZJs t._W643~ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
mn=G6h
T}W /CtR|~ w L arp=angle(rp);
D/CSR= b r+BPz%wM=O ars=angle(rs);
OG_2k3v @x>J-Owd]J atp=angle(tp);
'w+T vOB _fH.#C ats=angle(ts);
+G3&{#D
? ,qgph^C figure(1)
dpO ZqhRs. 29?{QJb subplot(1,2,1);
;[-dth m CFScT plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
nQc ]f* Xi'y-cV
^ legend('arg(r_p)','arg(r_s)')
f$x\~y<[ 9HlRf6S xlabel('\theta_i')
F z_SID Q0f7gY1-% ylabel('\phi')
sskwJu1 ,X&lVv# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/++CwRz@Gm ?hh 4M axis([0 90 -3.5 3.5])
t)n!]; ] 7C=.'Y grid on
-.| V S|y Tp.0@aC subplot(1,2,2);
[N}:Di,S 3,8>\yf` plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
R 2uo ZA, 'aQ"&GX@ legend('arg(t_p)','arg(t_s)')
Si#b"ls' 1&~u:RUXe xlabel('\theta_i')
:, $:@ LYPjdp2>"o ylabel('\phi')
(yv&&J c uI7 d?s title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
XV+s 5 C u/e-m / axis([0 90 -3.5 3.5])
vu[+UF\G 'W 5r(M4U grid on
lzz rzx^ `MAluu+b 
