计算光波由光密入射到光疏介质时相位变化 f T7Z6$
e|ChCvk
clear - v\n0Jt
%}t<,ex(yO
close all @fDWp/
[ &IJ y
n1=1.45,n2=1; -,;r %7T
YA7h! %52)
theta=0:0.1:90; <2{CR0]u
[^iQE
a=theta*pi/180; ze%kP#c6!
/FzO9'kj
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); qu}`;\9@ld
6 fP"I_c
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); PS*=MyNa
2(_+PQ6C=
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); p&Os 5zw;|
'Q R
@G
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); B vXA9YQ3
Equj[yw% @
arp=angle(rp); !/1aot^(
Z/y&;N4
ars=angle(rs); =Gka;,n
P>*B{fi^
atp=angle(tp); a4 zq`n|3U
?*2DR:o>@
ats=angle(ts); M qy5>f)
Gp0B^^H $
figure(1) 8M{- RlR
o\n9(ao
subplot(1,2,1); k!{0ku}]
s-\.j-Sa
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) p};B*[ki
<!+T#)Qi
legend('arg(r_p)','arg(r_s)') 'qhi8=*
4$j7DJ8dj
xlabel('\theta_i') 5PaOa8=2f
m)5,ut /
ylabel('\phi') {8'f>YP
`v]|x,l+C
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ~ULD{Ov'F
(\CT
"u-
axis([0 90 -3.5 3.5]) |4=Du -e
sj"z gE)
grid on z.NJu
q
Bh'_@PHP
subplot(1,2,2); G5C=p:o{/
# :^aE|s
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) f!3$xu5
3WOm`<
legend('arg(t_p)','arg(t_s)') -i,=sZXB
H$M#+EfL
xlabel('\theta_i') [ ?dsS$Y3
29W~<E8K-
ylabel('\phi') | )"`v'8>
oxs0)B
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) fyRSg B00$
({"jL*S,q
axis([0 90 -3.5 3.5]) E%A] 8y7
>I9w|z FA
grid on 2j9+ f{ l
XZ|%9#6
