计算光波由光密入射到光疏介质时相位变化 !%v=9muay
KDt@Xi6||
clear Y$tgz)
zxo0:dyw7
close all ^
W/,Z`
kTCWyc
n1=1.45,n2=1; C3m](%?
kaKV{;UM
theta=0:0.1:90; P:`tL)W_
G/cE2nD
a=theta*pi/180; ^;KL`
C}})dL;(
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); CBj&8#8Z
1m$< %t.>
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); $s[DT!8N
Muhq,>!U
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); SfHs,y6
PA=.)8
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); WKHEU)'!
xt{f+c@P
arp=angle(rp); d{~5tv- H
$
N7J:Q
ars=angle(rs); h[Hn*g
AdCi*="m
atp=angle(tp); %cPz>PTW@
Hcts^zm2u
ats=angle(ts); &s]wf
|l*#pN&L
figure(1) oH[4<K>
Wa/&H$d\u@
subplot(1,2,1); X={n9*Sd8
f]P&>j|
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) L `7~~
D$
>gAv
legend('arg(r_p)','arg(r_s)') 3`W=rIMli
cqr4P`Oj
xlabel('\theta_i') - %ul9} .
}w,^]fC:
ylabel('\phi') xj1FCT2
EQ;,b4k?&g
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) im}=
AbWnDqv
axis([0 90 -3.5 3.5]) (|(#W+l~
W~TT`%[
grid on 2 g)W-M
%B;e7
UJ
subplot(1,2,2); sz5&P )X
~ jR:oN
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) OZHQnvZ
G! y~Y]e
legend('arg(t_p)','arg(t_s)') K):MT[/"
W6b5elH@
xlabel('\theta_i') h4j{44MT
$m.e}`7SF!
ylabel('\phi') G=[=[o\
O<@S,/Q4
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ?5%0zMC
OOa}+^-j
axis([0 90 -3.5 3.5]) 4 Ar\`{c>
|{KZ<
grid on Sdc
yL%6!
~"<AYJlO