计算光波由光密入射到光疏介质时相位变化 O[ Vet/^)
s@F&N9oh
clear ]vvYPRV76
F8f}PV]b
close all =,6H2ew
&lQ%;)'
n1=1.45,n2=1; X Q#K1Z
D=0YLQ*rP
theta=0:0.1:90; srGOIK .
?=? _32O
a=theta*pi/180; k]~o=MLmj
XX8HSw!w
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); q oO`)<
&k:xr,N=
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); DB/~Z
$w0TEO!
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); EeW
,-I
X1Yw=t~a
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); DrK]U}3fh"
Z0,jg)s A4
arp=angle(rp); x-BU$bx5
o(*\MT t?
ars=angle(rs); 2S?7j[@%i`
vO}r(kNJ
atp=angle(tp); %{u@{uG0'3
r9[S%Def
ats=angle(ts); uvId],dQ5
DrJ?bG;[
figure(1) u8Ys2KLpL
[G<ga80
subplot(1,2,1); fVbjU1N
uo#1^`P
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) !oWB5x~:P
7Ak<e tHD
legend('arg(r_p)','arg(r_s)') (RddR{m X
cQ8[XNa
xlabel('\theta_i') (95|DCL
Cv**iW
ylabel('\phi') _@DOH2 lXJ
V]7/hN-Y}
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) -D`1z?zHra
:zY4p hR
axis([0 90 -3.5 3.5]) "#1 \ uoH
@<\oM]jX
grid on Pb5yz-?
)UzJ2Pa<+_
subplot(1,2,2); }y(t')= 9
w!F >fcm
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) =:H-9
X +`Dg::
legend('arg(t_p)','arg(t_s)') OX_y"]utU
%U\,IO `g
xlabel('\theta_i') $ L*gtZ
uNyU]@R<W
ylabel('\phi') w 1/QnV
\<7Bx[/D4
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) -S sgW
Lq(=0U\"P
axis([0 90 -3.5 3.5]) Y~"tL(WfJl
69c4bT:b"
grid on yE:y[k0E
.S
k+"iH5
