计算光波由光密入射到光疏介质时相位变化
y'`7zJ rfVQX<95=/ clear
7+f6?
)nf%S+KV close all
]/Nt >)C7IQ/ n1=1.45,n2=1;
'.z7)n %XN;S29d5W theta=0:0.1:90;
R98YGW_
dT [;IE Z/ZX a=theta*pi/180;
#tA9`! @!oN]0`F; rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`XE>Td>Bs D+;4|7s+ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
\?t8[N\_[( G{6@]72 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
wxcJ2T d H mK$E&,OkA ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
i/-IjgM"- Uv[a
~' arp=angle(rp);
y
_ap T<P @C@9Tw2Y ars=angle(rs);
p@~Y[a = g7?[}?]3"p atp=angle(tp);
N KgEs (*r2bm2FPO ats=angle(ts);
2~/`L=L Z*w({k7] figure(1)
]vB^% Dti-*LB1 subplot(1,2,1);
0
N7I:vJ 'fK=;mM plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
p_BG#dRM Ht+ro Y legend('arg(r_p)','arg(r_s)')
?;](;n#lU uq9mq" xlabel('\theta_i')
!NNPg?Y 0^Ldw)C" ylabel('\phi')
i#YDdz ;DX{+Z[ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
X~m57bj -SD:G]un
axis([0 90 -3.5 3.5])
Ay6T*Nu` z^gz kXx7 grid on
z5({A2q b/*QV0( subplot(1,2,2);
An(gHi;1$ FEhBhv|m plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
o7+<sL 1f^oW[w& legend('arg(t_p)','arg(t_s)')
zx"EAF{ hU( xlabel('\theta_i')
0e"KdsA:<U =3hJti9[ ylabel('\phi')
?~.9:93 1c"s+k]9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Bz,D4E$ Us "G X_ axis([0 90 -3.5 3.5])
@.cord` A<H]uQ> grid on
W!8$:Ih_Z *]q`:~u2