计算光波由光密入射到光疏介质时相位变化 T^(W _S
BV:Ca 34&
clear kfZ `|w@q
At?|[%< `
close all Fa [^D~$l*
*1{A'`.=\
n1=1.45,n2=1; 4\u `M R
#qY gQ<TM!
theta=0:0.1:90; /51$o\4 S
!p!Qg1O6o
a=theta*pi/180; fmj}NV&ma
`;+x\0@<
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ;!A8A4~nu
$"W[e" Q
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); O-
QT+]
qnboXGaFu
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); KR?aL:RYb
T` h%=u|D
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); pxgf%P<7
\-pwA j?
arp=angle(rp); cM' MgX9
t3Gy *B
ars=angle(rs); T&d Njx
~6n|GxR.[
atp=angle(tp); +_J@8k
U1y!R<qlp
ats=angle(ts); f?tU5EX
En9J7es_
figure(1) s(3u\#P
\ YF@r7
subplot(1,2,1); H 4 ELIF#@
U<0Wa>3zj
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) oQ
r.cKD ?
kE'p=dXx
legend('arg(r_p)','arg(r_s)') nc:/GxP
a5/, O4Q
xlabel('\theta_i') ME]4tu
[k\VUg:P
ylabel('\phi') ,\"x#Cc f
bT2c&VPCE
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) "KK}} $>
##P zc~xSn
axis([0 90 -3.5 3.5]) '+/m t_re=
d ^RcJ3w
grid on ) u Hat#
6O{QmB0KK
subplot(1,2,2); 02;'"EmP$
>$rH,Er
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) V=+|]`
PU.j (0
legend('arg(t_p)','arg(t_s)') P{K\}+9F
`;}qjm0a
xlabel('\theta_i') >Um(gbG
#@BhGB`9Qt
ylabel('\phi') F%y#)53g
?SYmsaSr5
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) C<a&]dN/
h~sTi
axis([0 90 -3.5 3.5]) Hfj.8$
ck
`td%
grid on #B?lU"f8q^
6bj77CoB
