计算光波由光密入射到光疏介质时相位变化 %PdYv _5
8U:dgXz
clear rHBjR_L.2
JR<-'
close all Bpo68%dx89
TIh zMW\/K
n1=1.45,n2=1; z slEUTj)
h5J wB<8
theta=0:0.1:90; EM
w(%}8w
A^@ <+?
a=theta*pi/180; u\geD
@d ^h/w
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); Hpg;?xAT
Y <k,E
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); :QL p`s
RZh)0S >J
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); N_Ld,J%g
[=F
|^KL
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); "s<l Lgi
,wo" (E!4e
arp=angle(rp); +* {5ORq=
fO(S +}
ars=angle(rs); ( #iM0{
>D4Ez
atp=angle(tp); SfL`J Ni)
.
\0=1P:
ats=angle(ts); I8]NY !'cW
.%Q Ea_\
figure(1) %ys}Q!gR
",V5*1w
subplot(1,2,1); HYmUxheN2
32P ]0&_O
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) ^tcBxDC"]
1+}Ud.v3VW
legend('arg(r_p)','arg(r_s)') 2I 7 `
Bic {
H
xlabel('\theta_i') RAbq_^Q
Q b{5*>
ylabel('\phi') )_K@ ?rWS
W(4?#lA2W
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) %{ abRBny
4F6o
axis([0 90 -3.5 3.5]) =T}uQ$X
QXj(U&#rp
grid on i 4rF~'h@
7zZ|=W?&{
subplot(1,2,2); Rc P5].^T
{=7i}xY]T
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) |UO;St F
y visoZX
legend('arg(t_p)','arg(t_s)') T=d vc }
):ZumG#o
xlabel('\theta_i') " <a|Q ,!
7!nAWlQ&-E
ylabel('\phi') r/L]uSN
53 05N!
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) eJp-s" %
y<d#sv(s
axis([0 90 -3.5 3.5]) gj
X1b2
'FFc"lqj
grid on <U pjAuG8
Fsj[J E
