计算光波由光密入射到光疏介质时相位变化
B'yN &3 Vy}:Q[ clear
h JF Q/( jq.@<<j|$ close all
YI%7#L7C YLPiK n1=1.45,n2=1;
$23="Jcl c0Q`S"o+ theta=0:0.1:90;
ucoBeNsHx ik &lo M_ a=theta*pi/180;
C,tl p D3XQ>T [*q rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
XHN?pVZ7 , wX/cUyZ
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
m?[F)<~a c?q#?K
aF tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
1-w1k ^e !m_'<=)B4~ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
3D<P
[.bS Yn J=&21 arp=angle(rp);
!vImmhI!I W!IK>IW" ars=angle(rs);
L x9`y t6 Lg Bs<2 atp=angle(tp);
3kKXzIh oWXvkDN
ats=angle(ts);
L0+@{GP? {_k 6 t figure(1)
\BJnJk!% vtL) subplot(1,2,1);
WU=Os8gR 3*8#cSQ/6o plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
^GRd;v =-@ +?mZ_sf8w legend('arg(r_p)','arg(r_s)')
r43dnwX -$e\m]
}Z xlabel('\theta_i')
&m Y<e4 X_%78$N-a` ylabel('\phi')
Ui? t@. !_x-aro3< title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,{E'k+ YdeSJ(: axis([0 90 -3.5 3.5])
a 2).Az ss,6;wfX grid on
A`OU} 'v?L >TL0hBaaR subplot(1,2,2);
B'~.>, fg N|7._AR2 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
N b B`6@ r R \`,Q'3 legend('arg(t_p)','arg(t_s)')
c-s ~q/ 0 'L+9T5 xlabel('\theta_i')
*B3` #t E; RI.6y ylabel('\phi')
{a(YV\^y|H Xq J@NgsY title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^-=,q.[7 w l9E axis([0 90 -3.5 3.5])
a<vCAFQ T}4 RlIZF grid on
: [AW y YF80mnJz 
