计算光波由光密入射到光疏介质时相位变化 ?~BC#B\>o
dk[!V1x4\
clear k:t ]s_`<
lA|
5E?
close all V,lOt4b
Z]>O +
n1=1.45,n2=1; x gVeN["
&#{Z( h.de
theta=0:0.1:90; ]#)()6)2v
_<n~ n]%
a=theta*pi/180; XRM_x:+]
;w{tv($$
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); lk/n}bx
#^u$
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); s Ep"D+f
(9''MlG d%
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); 2Q/x@aT,h
M'?,] an
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); 2V-
16Q'%
> c1qpk/
arp=angle(rp); GFj{K
| k&Ck
ars=angle(rs); CM/H9Kz.
>N^Jj:~l
atp=angle(tp); 1OG x>J6
cvn@/qBq*t
ats=angle(ts); bn|I> e
]R+mKUZ9
figure(1) N]>=p.#j
Ci@o|Y }tP
subplot(1,2,1); hQ,ch[j'
\j@OZ
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) dN$0OS`s[
C)um9 }
legend('arg(r_p)','arg(r_s)') Nyku4r0
Z;h t
xlabel('\theta_i') 0IsPIi"7
%f&/E"M
ylabel('\phi') QyEn pZ8?a
/ C:Y94B-z
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) v,FU^f-'
:(/~:^!
axis([0 90 -3.5 3.5]) ISYXH9V
29;?I3<
*
grid on VACQ+
lp.ldajN
subplot(1,2,2); nD_ GL
F%@A6'c
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) g#"zQv ON
KoOz#,()
legend('arg(t_p)','arg(t_s)') @:oXN]+
_
$njUX SQ;
xlabel('\theta_i') u0JB\)(-/h
: y'Ah#
ylabel('\phi') W>${zVu
iFd+2S%
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) /K+r?
]kf
sQYkQ81
axis([0 90 -3.5 3.5]) DwM)r7<Ex
!{ y@od@T
grid on 2Z+Wu3#
ew>XrT=Zm
