计算光波由光密入射到光疏介质时相位变化 } bv+^#
XnZ$ %?$
clear 7K\v=
UHZ&7jfl
close all TQT3]h6
5G}4z>-]F)
n1=1.45,n2=1; s#s">hMrI
E+lR&~mK=
theta=0:0.1:90; ks0Q+YW
k~ 8-E u1
a=theta*pi/180; ?>h
~"D#
o'#& =h$_
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); * @4 @eQF
;#85 _/
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); zWC| Qe
~#sD2b` 0
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); KEOk%'c,
rkjnw@x\
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); i!%bz
bFG~08Z ,d
arp=angle(rp); 3 J#LxYK
2t\0vV2)/O
ars=angle(rs); ='h2z"}\Bn
!q2zuxq!R
atp=angle(tp); y0 d=
5skN'*oG
ats=angle(ts); *D?_,s
*]]Zpa6
figure(1) #dZs[R7h
|{ N{VK
subplot(1,2,1); ZM.'W}J{ *
S ~Gse+*
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) K W .4 9
(Mk7"FC7
legend('arg(r_p)','arg(r_s)') gK rUv0&F
p / ITg
xlabel('\theta_i') Mh)? A/e
ab6I*DbF
ylabel('\phi') }^&S^N 7
\l59/ZFan
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) :*Wq%Y=
Wlc&QOfF
axis([0 90 -3.5 3.5]) MK BDWLCB
1 $&@wG
grid on K0.a U
`==l 2AX
subplot(1,2,2); ME)Tx3d
dbby.%
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) WE) *~5
zy,SL
|6:
legend('arg(t_p)','arg(t_s)') Jy(G
A
P{o)Ir8Tt
xlabel('\theta_i') q94;x |63
?AyxRbk
ylabel('\phi') L8`v
yT7$6x
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 8US#SI'x
]h S:0QE
axis([0 90 -3.5 3.5]) ,6>3aD1w~q
OijuOLt
grid on eAlOMSL\
>Z!!` 0{
