计算光波由光密入射到光疏介质时相位变化 b|DU
*GBV[D[G,
clear L/-SWid)
FZ+2{wIV^
close all 8Nyz{T[
+pYwc0~
n1=1.45,n2=1; F9K`N8wlu
LTo!DUi`
theta=0:0.1:90; eaDZ^Z
Er
"N=$=Dy>
a=theta*pi/180; `^hA &/1
y] D\i5Xv
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); X=)L$Kd7
a6./;OC
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); bO/r1W
m[2[9bQ0
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ||pOiR5
qp6'n&^&
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); e.DN,rhqI
wZ\93W-}
arp=angle(rp); XIbxi
\EoE/2"<
ars=angle(rs); \:1$E[3v
bF_0',W
atp=angle(tp); IO"P /Q
R$
+RTG:E
ats=angle(ts); !]g[u3O
/]Y#*r8jRi
figure(1) oYf+I
E Zf|>^N
subplot(1,2,1); 8\'tfHL
=UK:83R(
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) G`0{31us
<|Lz#iV37
legend('arg(r_p)','arg(r_s)') Tb?X KO,
[9Rh" H;h
xlabel('\theta_i') )z74,n7-
t!/~_}eD J
ylabel('\phi') k`N*_/(|n
4To$!=
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) "'Fvt-<^S7
dazML|1ow
axis([0 90 -3.5 3.5]) JB'tc!!*
h'h8Mm
grid on `V V>AA5
J9NuqV3
subplot(1,2,2); v+Y^mV`|
>m1b/J3#
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) R+uZi~
L,#YP#O,j
legend('arg(t_p)','arg(t_s)') lN5PKsGl
Ce<z[?u
xlabel('\theta_i') RPeH [M^
p7`9
d1n
ylabel('\phi') N-F&=u}
D60quEe3%
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) {{P 3Z[
s,[I_IiPf
axis([0 90 -3.5 3.5]) MDIPoS3BRa
j1{`}\e
grid on CckfoJ 9
O
]
!tK