计算光波由光密入射到光疏介质时相位变化 Jg%sl&65
:I&iDS>u1
clear v4Nb/Y
*|` ' L
close all huN(Q{fj
1BWuFYB
n1=1.45,n2=1; bNL E=#ro
W]y$6P
theta=0:0.1:90; {fX4
[NKWudq
a=theta*pi/180; (/hF~A
/>n0&~k[h
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); DO%YOv
&1Az`[zKGW
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); ?8~$du$
N Z.aI{
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); \64(`6>
134wK]d^
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); }B]FHpi
\SMH",u
arp=angle(rp); %Z}dY~:
ZW\h,8%
ars=angle(rs); %/>_o{"hw
m+vwp\0
atp=angle(tp); +osY
iP5
5-&P4
ats=angle(ts); :;|x'[JoE?
,p 'M@[
figure(1) V];RQWs
_VM J q9.
subplot(1,2,1); zR5D)`Ph
vhUuf+P*
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) B!-W765Y
g.9MPN
legend('arg(r_p)','arg(r_s)') LLU>c]a
LpF6e9V\Wp
xlabel('\theta_i') 1yE~#KpH
''WX
ylabel('\phi') q$HBPR4h
kW(8i}bg
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) hA~}6Qn
DSnsi@Mi
axis([0 90 -3.5 3.5]) JHMj4Zkp
G-:DMjvN
grid on OZv&{_b_
uZ'5&k96T
subplot(1,2,2); k5PzY!N
hpw;w}m
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) dkVVvK
xbmOch}j6
legend('arg(t_p)','arg(t_s)') R'80 {
y\%4Dir
xlabel('\theta_i') 5N[Y2
1-b,X]i
ylabel('\phi') )c!f J7o:
a<HM|dcst
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) B$lbp03z
*/S,CV
axis([0 90 -3.5 3.5]) dA_s7),
/evh .S
grid on =lffr?#&B
%Ajf|Go0/G