计算光波由光密入射到光疏介质时相位变化 ~'N+O K
BOv ^L?)*Z
clear p tlag&Z
p ^Y2A
close all S2EV[K8 #
y 7z)lB y\
n1=1.45,n2=1; [Xww`OUsh
zEDN^K '
theta=0:0.1:90; 9[h8Dy
3|D .r-Q
a=theta*pi/180; 1z_1Hl
U5[r&Y
D
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); @>+`1C
d3NER} f4V
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); AJ
z 1
>J{e_C2ZS
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); _ o==
l=C| 4@
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); =KHX_ib
#JR$RH
arp=angle(rp); E5Lq-
60l!3o"p!
ars=angle(rs); B?ipo,2~{
XpPcQIM*
atp=angle(tp); (O/W`qo
le6eorK8
ats=angle(ts); YD#L@:&gv
!ZY1AhGZ
figure(1) f;!L\$yKy
/sSM<r]5j
subplot(1,2,1); t-Ble
6n
H'NNS:J
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) %!R\-Vej
V!s#xXD }
legend('arg(r_p)','arg(r_s)') Y6 a9S`o
k=X)ax t1
xlabel('\theta_i') +No Ve#
&D&U!3~(
ylabel('\phi') pL: r\Y:R
TJ|do`fw >
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) $ |VdGRZ1
gp+@+i>b+[
axis([0 90 -3.5 3.5]) 10_ >EY`
R tEx
WTc
grid on <&HHo>rl
'tt4"z2
subplot(1,2,2); f DPLB[
~|ha9 1
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) k< RaC=
9TEAM<b;
legend('arg(t_p)','arg(t_s)') f n8|@)J
pxCQ=0k
xlabel('\theta_i') 4}0Ry\
6
/ 1eeNbd
ylabel('\phi') + ij xv
me90|GOx+
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) $5DlCN
:bkmm,%O
axis([0 90 -3.5 3.5]) b]*OGp4]5
sIR fC<
/P
grid on {uxTgX
.]N` ]3$=
