计算光波由光密入射到光疏介质时相位变化 hcM 0?=
Px)VDs=k
clear E7y<iaA{~
E-&=I> B5
close all %#"uK :(N
6nxf <1
n1=1.45,n2=1; Lt=#tu&d
dB< \X.
theta=0:0.1:90; !+CRS9\D
wgI$'tI
a=theta*pi/180; E]"ePdZZ/
=q>eoXp
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); <L:v2 8c
RYmk6w!w
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); S1zV.]
1\_4# @')
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); "jb?P$
m<r. sq&;
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); sL[,J[AN;
<A+Yo3|7
arp=angle(rp); -s4qm)\
7?B]X%
ars=angle(rs); Ks9"U^bPs
b\H~Ot[i
atp=angle(tp); Mx[tE?!2
/q(+r5k \
ats=angle(ts); Rl<~:,D
~h0SD (
figure(1) ~M,nCG^4
Qz[~{-<
subplot(1,2,1); JF!!)6!2#
N',]WZ }
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) PK!=3fK4\F
e6#^4Y/+`
legend('arg(r_p)','arg(r_s)') "l2_7ZXsPT
4*d_2:|u
xlabel('\theta_i') bV`Z o(z
Q!|. ,?V
ylabel('\phi') k45xtKS>d
O*!+D-
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) a +,)rY9
Vh4z+JOC
axis([0 90 -3.5 3.5])
X1y1
$UO7AHk
grid on ym8\q:N(R
q<.^DO~$L
subplot(1,2,2); q'd6\G0 }
f4]nz:2
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) a!xKS8-S==
aW$7:<A{
legend('arg(t_p)','arg(t_s)') v!K %\h2A
Vz51=?75
xlabel('\theta_i') Kx$?IxZ
z`m-Ca>6
ylabel('\phi') My
Af~&Y+
4s.wQ2m
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Xy=|qu
w ,9$*=k
axis([0 90 -3.5 3.5]) p*n$iroy_{
4|7L26,]5
grid on 2 u/(Q >#
3-~_F*%ST
