计算光波由光密入射到光疏介质时相位变化 }vOg9/[{
!\d~9H%`B
clear ^W%F?#ELN2
qzA_ ~=g
close all ar,v/l>d4N
YcwDNsk
n1=1.45,n2=1; 4guR8 elM
N }NKQ]=
theta=0:0.1:90; MaD| X_g
m-tn|m!J
a=theta*pi/180; cr ]b #z
60n P'xfR
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); B6TE9IoSb8
y4|<+9<7
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); ) :Z #!O<
qrj f
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); M=ag\1S&ZF
YU+P+m2X
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); "s]
4I2:"CK06
arp=angle(rp); $8 &Y(`
P*K"0[\n
ars=angle(rs); <A| z
[*(1~PrlO,
atp=angle(tp); mS;WNlm\
^q/$a2<4
ats=angle(ts); Nl,iz_2]
a JjUy %
figure(1) 8mRZ(B>% X
+]- '{%-zK
subplot(1,2,1); \4j_K*V
hWFOed4C
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) n/*" 2
L
aA<`
legend('arg(r_p)','arg(r_s)') +_*NY~
<"`f!k#[
xlabel('\theta_i') <P0 P *>M
fJ
_MuAv
ylabel('\phi') LE 5N2k
K
re*~ "
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ,?qJAV~>
l qh:c
axis([0 90 -3.5 3.5]) }S')!3[G
xZ%3e
sp
grid on <3N\OV2
''q;yKpaz
subplot(1,2,2); e:4,rfF1
*\}$,/m['
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) ht6}v<x.eA
/g9^g(
legend('arg(t_p)','arg(t_s)') FYE(lEjxi
NYg& 8 s.
xlabel('\theta_i') nm% qm
lf
KV%
ylabel('\phi') _7;G$\^&.
3$K[(>s
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) qR cg|']R
;g+fY 6
axis([0 90 -3.5 3.5]) ?z-nY,'^uq
vR5X
grid on {&D$U'ye
lB/ ^
