计算光波由光密入射到光疏介质时相位变化 hag$GX'2k
6v!`1}
~
clear ,t744k')
2nObl'ec
close all g =hg%gRy"
s)t@ol
n1=1.45,n2=1; <}Vrl`?h
?<,l3pwqa
theta=0:0.1:90; E~T-=ocKE
?81c 4w
a=theta*pi/180; ]c*4J\s
FUzzB94a
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ~/U1xk%
P; no?
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); ;1=1:S8
2.y-48Nz
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); kazzVK5x
uMv,zO5
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); :4w ?#
?R
'r4P,
arp=angle(rp); ~P
qM]^
M0"_^?
ars=angle(rs); zI uJ-8T"
ttQGoUkj
atp=angle(tp); MJ)RvNF
n&/
`
ats=angle(ts); ) q4[zv9
>6-`}G+|
figure(1) }v{LRRi
MchA{p&Ol
subplot(1,2,1); (lqC[:
G!##X: 6'
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) |-ALklXr
E]d.z6k
legend('arg(r_p)','arg(r_s)') d1T!+I
,qwuLBW
xlabel('\theta_i') R\f+SvE
Gz0]}]A
ylabel('\phi') wJqMa9|
>R_&Ouh:
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) U/M>?G~
HJH{nz'Lw
axis([0 90 -3.5 3.5]) (WJRi:NP?
3}1u\(Mf
grid on M_f:A
NxY#NaE:?4
subplot(1,2,2); ^\% (,KNo
|r/"
|`
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) cjY-y-vO
@HC Vmg:
legend('arg(t_p)','arg(t_s)') gQuw1
]EAO+x9
xlabel('\theta_i') >U>(`r*
j (d~aqW
ylabel('\phi') r6qj7}\
X?',n
1
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ?V=ZIGj
}X6m:#6
axis([0 90 -3.5 3.5]) #{6/ (X
eByz-,{P
grid on =nS3p6>rZ
3d]S!=4H"