计算光波由光密入射到光疏介质时相位变化 ;n1<1M>!
@)n xX))a
clear D&0y0lxI@
0*7N=
close all {6}H}_(]
EMK>7 aks
n1=1.45,n2=1; bn$a7\X-
3jH-!M5
theta=0:0.1:90; {-?^j{O0.
JAEn
72
a=theta*pi/180; b7;`A~{9v
;;|S
QX
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); OAx5 LTd
#Xdj:T<*
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); 4I&e_b< 30
+Wgp~$o4
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); i
b$2qy
e0rh~@E
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); _4~'K?
=Rv!c+?
arp=angle(rp); 9&]g2iT P
Z]VmTB
ars=angle(rs); O>FE-0rW}e
G_m$W3 zS
atp=angle(tp); W_JFe(=3,
_4+'@u
#
ats=angle(ts); Lmp_8q-Ej
*SP@`)\D
figure(1) Q:rQ;/b0/
RQg7vv]%
subplot(1,2,1); G%Dhj)2}
]!%
p21e
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) J*I G]2'H
n*yVfI
legend('arg(r_p)','arg(r_s)') #H8% BZyV
yv),>4_6
xlabel('\theta_i') .A6i?iROe
hW~XE{<
ylabel('\phi') ,\".|m1o.
J]m[0g7O_
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Uxll<z,
^ "
axis([0 90 -3.5 3.5]) a%FM)/oI|T
!t.*xT4W
grid on e"I+5r",
cJ^:b4j
subplot(1,2,2);
)c;zNs
)c; YR}tC
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) Ja`xG{~Y7i
*PSUB{i(
legend('arg(t_p)','arg(t_s)') !=/wpsH
sfo+B$4|
xlabel('\theta_i') q
eW{Cl~
[v%j?
ylabel('\phi') kuX{2h*`
655OL)|cD6
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) )0\"8}!
\jHHj\LLr.
axis([0 90 -3.5 3.5]) 16cc9%
78u9> H
grid on V]{^}AKc
"kdmqvTHK0