计算光波由光密入射到光疏介质时相位变化 D||0c"E
!2 t7s96
clear bQ*yXJ^8
y
4i3m(S
close all e-dkvPr
ebA:Sq:w
n1=1.45,n2=1; }geb959
yY VR]H H
theta=0:0.1:90; #Y5k/NPg
c7mKE`
a=theta*pi/180; sH6;__e
$N?8[
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); 12`u[O}\}-
^Iw$ (
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); .IW`?9O$E
KVZB`c$<t
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); [ CQR
ysnW3q!@
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); MCN>3/81
snaAn?I4
arp=angle(rp); D!`;v Z\>
}nr jA0WN
ars=angle(rs); sVyV|!K
0F[ f%2j
atp=angle(tp); LS~at.3zX
]\*^G@HA2
ats=angle(ts); Z6I|Y5#H
[c{\el9H
figure(1) H07\z1?.K
bq/Aopfr
subplot(1,2,1); K P]ar .
R4p Pt
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) Tpl]\L1v-
. `Rt
legend('arg(r_p)','arg(r_s)') @&\Y:aRO%i
oC5 h-4~
xlabel('\theta_i') Hz TmNm)
ACyK#5E
ylabel('\phi') Y4k2=w:D
?.e,NHf
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) G)am ng/
3;>ls~4
axis([0 90 -3.5 3.5]) 8}|et~7!
C8m 9H8Qm
grid on ` aVp #
*(9 Tl]w
subplot(1,2,2); >jt2vU@t.
[
U:C62oK,
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) +d3|Up8=
</ [.1&S+\
legend('arg(t_p)','arg(t_s)') C hF~
k%6CkC w
xlabel('\theta_i') ABCm2$<
(%6fMVp
ylabel('\phi') /p
!A:8
PYCN3s#Gi
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) =F+v+zP7P
NWf=mrS8@$
axis([0 90 -3.5 3.5]) c_>Gl8J
ed6@o4D/kf
grid on ?Bf>G]zx
s|\)Y*B`
