计算光波由光密入射到光疏介质时相位变化 W m\HZ9PN
hvCX, ^LoJ
clear }@%A@A{R
c+$a lw L~
close all ?CIMez(h
h}r64<Y2{
n1=1.45,n2=1; s)X'PJ0&Bs
a<-NB9o~v
theta=0:0.1:90; l?GN& u
cnJ(Fv_F$
a=theta*pi/180; }tue` ">h
j x B
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); +Qy0K5Ee
L5$r<t<
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); k+
[V%[U
uj#bK
7
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); :k9n
9
Ve\^(9n
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); VBV y3fnj
.: gZ*ks~
arp=angle(rp); 6$]@}O^V
OSACH0h
ars=angle(rs); uFPJ}m[>5
_"Yi>.{]
atp=angle(tp); 6i AHus-
&, E^ y,r
ats=angle(ts); ;s{k32e
z Ic%>?w
figure(1) USB U?WDt
<NRW^#g<x
subplot(1,2,1); g=W1y
}!-BZIOlO
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) PNVYW?l
S5B12P
legend('arg(r_p)','arg(r_s)') W
&wqN
cb|cY Co5
xlabel('\theta_i') L[v-5 u)
M.B0)
ylabel('\phi') l0AVyA4RFV
sXe=4`O
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) e<F>u#d
@~WSWlQW
axis([0 90 -3.5 3.5]) ~Q<h,P
7 nr+X Os
grid on )Pr*\<C ld
!<` }m E!:
subplot(1,2,2); up`.#GWm
:
&! >.Y
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) < zUU`
)0F\[Jl}
legend('arg(t_p)','arg(t_s)') M PSoRA: h
t<sy7e='
xlabel('\theta_i')
d3% 1 P)
lJZ-*"9 V
ylabel('\phi') }~/u%vI@M5
} <G"w 5.<
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) F"2rX& W
o Efy{54
axis([0 90 -3.5 3.5]) `2}H$D
v(Vm:oK,
grid on !a
%6nBo
`{1` >5
