计算光波由光密入射到光疏介质时相位变化 2V#6q,2
O@U?IF$
clear Wz=OSH7"f
mE $dO3
close all ryD%i"g<
1/j$I~B
n1=1.45,n2=1; \4~AI=aw,T
eO5ktEo J
theta=0:0.1:90; :=g.o;(/N
&/F_ *=VE
a=theta*pi/180; vP]9; mQ
T nxKR$Hoh
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); U ^#?&u
_\AUQ{
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); X B I;Lg
e9@(/+
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); _x\m|SF_g
q(9%^cV6
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); a| x8=H
8GB]95JWwp
arp=angle(rp); =:K@zlO:
N=fz/CD)I
ars=angle(rs); uCUu!Vfeg
.j 'wQ+_
atp=angle(tp); X\A]"su
JieU9lA^&B
ats=angle(ts); b~w KF0vq
i.@*t IK
figure(1) h%b hrkD
Cg6;I.K
subplot(1,2,1); qpgU8f
>ZCo 8aK
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) I)[B9rbe
&c^7O# j
legend('arg(r_p)','arg(r_s)') u~Lu<3v
S:97B\ u`
xlabel('\theta_i') $%}>zqD1
mTd<2 Hy
ylabel('\phi') Q;gQfr"c7
' o=E!?
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 2 ]Fu
1
gE= Wcb!
axis([0 90 -3.5 3.5]) Q x.jCy@
j;0ih_Z@4W
grid on HWbBChDF
pGk"3.ce
subplot(1,2,2); mVrK z
( -2R{! A
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) 24|<<Xn
Gc}0]!nrW9
legend('arg(t_p)','arg(t_s)') r%[1$mTOR
^Ezc y?
xlabel('\theta_i') {{zua- F
c,AZ/t
ylabel('\phi') e[p^p!a
T g\ hx>
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) VD4S_qx
3:O +GQ*
axis([0 90 -3.5 3.5]) 7 s-`QdWX
{.#j1r4J`
grid on uXP-
J ]>
{jUvKB_x
