计算光波由光密入射到光疏介质时相位变化
)bPwB.} kq W >;AMun clear
2|w(d T|ZT&x$z close all
T JLz^%t *E+) mB"~ n1=1.45,n2=1;
5j,)}AYO H*; J9{ theta=0:0.1:90;
G@ed2T r 3p fG a=theta*pi/180;
{%b >/r ,&z_ 2m rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
si%f.A # 2z ArAch rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
%+xh nPvR tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
~sMn/T*fv izmL8U
?t ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
0z.` C @hnT <e arp=angle(rp);
$P 4hN b bZERh:%o ars=angle(rs);
LTBqXh [C
P V5\2 atp=angle(tp);
Kze\|yJ - uliND ats=angle(ts);
89+m?H] K mY dU`j figure(1)
''v_8sv ~ EE*/vX subplot(1,2,1);
n3~xiQ' :m@(S6T m plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
dRas9g Woe sE:NiR legend('arg(r_p)','arg(r_s)')
v\$XhO K |F9/7 z\5+ xlabel('\theta_i')
m' z <d mbJ#-^}V ylabel('\phi')
z }u u+XZdV title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~`8`kk8 ?%*Zgk!l7 axis([0 90 -3.5 3.5])
J=zh+oLCV _U#u e grid on
@qg=lt|(F &[23DrI8 subplot(1,2,2);
L{<E'#@F vV8 y_ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
EQu M|4$ix n8R{LjJ2@ legend('arg(t_p)','arg(t_s)')
W# -M| zx0{cNPK 5 xlabel('\theta_i')
$;$_N43 B>|@XfPM ylabel('\phi')
|w:7).P `Z/"Dd;F^3 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'rq
[P ", [A[vR7&S axis([0 90 -3.5 3.5])
A GPZd9 9b()ck-\F# grid on
R &T(S LxsB.jb- 
