计算光波由光密入射到光疏介质时相位变化 `Oxo@G*@}W
4c"x&x|
clear <[9{Lg*D
\-A=??@H
close all k)+2+hX&>
ZMs$C3
n1=1.45,n2=1; ,dhSc<:LT
R7xKVS_MP
theta=0:0.1:90; iC]=S}
{<f_,Nlc
a=theta*pi/180; # fvt:iE
QH:PClW![
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); -*;-T9
Rlvb@aXgy
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); |)4Fe/!cJ
WW5AD$P*
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); >l|ao&z>bm
kjWY{7b!
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); j. 1@{H
bB01aiUw@l
arp=angle(rp); 7GN>o@ t
.L;M-`^
ars=angle(rs); i"eUacBz/-
MXy~kb&
atp=angle(tp); y7[D9ZvZ
:by EXe;3
ats=angle(ts); X%iiz
j6zZ! k
figure(1) M<'AM4
{FV_APL9_
subplot(1,2,1); {.vU;
r`? bYoz
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) ^Fmp"[q
d@o1<Q
legend('arg(r_p)','arg(r_s)') rGPFPsMQ]
1A>>#M=A
xlabel('\theta_i') ,{{e'S9cy
P.G`ED|K!Y
ylabel('\phi') DI/yHs
mJjd2a"vi
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) %`r Z]^H
FT0HU<." 1
axis([0 90 -3.5 3.5]) F(jvdq
e;QPn(
grid on +k@$C,A
nP9zTa
subplot(1,2,2); 8t{-
E038p]M!
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) wWaO"N]
kB41{Y -
legend('arg(t_p)','arg(t_s)') sn.Xvk%75
8T3j/D<r
xlabel('\theta_i') fU@{!;|Pz
gQXB=ywF
ylabel('\phi') lO?dI=}]
r!DUsE
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 2(5HPRQ
;xp^FKP
axis([0 90 -3.5 3.5]) xp+Z%0D
y'gIx*6B@
grid on V\C$/8v
Qw&It