计算光波由光密入射到光疏介质时相位变化 b7Jk{x #u
M%pxv6?""{
clear g0GC g
VeeQmR?u-
close all O`5PX(J1&
`V Y -3
n1=1.45,n2=1; :pj #t$:!
K.4t*-<`[
theta=0:0.1:90; vggyQf%
n ,}\;Bp
a=theta*pi/180; ~c~N _b
C-' n4AY^
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); QxG:NN;jW
p9rnhqH6
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); h'YC!hjp
|1C=Ow*"
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); $ ,Y\
vp7J';
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); B'"(qzE-kM
hi4#8W
arp=angle(rp); !PJD+SrG
>utm\!Gac
ars=angle(rs); k44s V.G4L
C1_':-4
atp=angle(tp); [F{q.mZj
m [7@l
ats=angle(ts); q66!xhp;?
dlk xA^
figure(1) TOmq2*,/
6&/n/g
subplot(1,2,1); de{@u<Y Zb
5/4N Y
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) 0.C[/ u[
@)=\q`vV
legend('arg(r_p)','arg(r_s)') #AHX {<
}tue` ">h
xlabel('\theta_i') j x B
+Qy0K5Ee
ylabel('\phi') wh8 h1I
Z9TmX
A@
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) L :UJur%
[yj-4v%u`
axis([0 90 -3.5 3.5]) VBV y3fnj
$j=c;+W
grid on kU ^*hd ]
OSACH0h
subplot(1,2,2); fui;F"+1
Za,rht
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) 1t7 vP;
2;X{ZLo
legend('arg(t_p)','arg(t_s)') g2Pa-}{
b #\i]2b:
xlabel('\theta_i') AmYqrmJ
2MS-e}mi
ylabel('\phi') PNVYW?l
b5Wt L+Z
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 1=Npq=d
,m ^q >
axis([0 90 -3.5 3.5]) "Z xM,kI
'Pk1 4`/
grid on vb^/ DMhz
z*,P^K 0T
