计算光波由光密入射到光疏介质时相位变化 RLw/~
w(M i?
clear tsB}'+!v#
tnX W7ej ^
close all ARJtE@s6Y
DfOig LG*
n1=1.45,n2=1; N>j*{]OY+{
}6__E;h#J
theta=0:0.1:90; : T4ap_Ycq
F Go)] U
a=theta*pi/180; g1UQ6Oa
;$r!eFY;
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); !$-QWKD4
NUi&x +
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); #\}xy PS
WK SWOSJ
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); VM3)L>x]/
o:%;AOcl
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); @nj`T{*.
+jGUp\h%9;
arp=angle(rp); (WGEX(|
eYg0 NEq{
ars=angle(rs); 0NSCeq%;6q
?zXl Lud8
atp=angle(tp); aTLr%D:Ka
4 Gh% PUV#
ats=angle(ts); )B^T7{
mndNkK5o
figure(1) (bogAi3<F
f./j%R@
subplot(1,2,1); y|# Fu
yA<\?P s
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) %G] W Oq=q
rIj B{X{Z
legend('arg(r_p)','arg(r_s)') y( r(q
LC/%AbM
xlabel('\theta_i') )]JQlm:H
ZN`I4Ak
ylabel('\phi') qS[nf>"
V#-8[G6Ra
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) }4*~*NoQ
p:4oA<V
axis([0 90 -3.5 3.5]) `[4{]jX+<
T!H }^v
grid on s9?H#^Y5u
A@ G%*\UZ
subplot(1,2,2); qP[jtRIN
MJ@PAwv"
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) dMH}%f5;1
&m@DK>
legend('arg(t_p)','arg(t_s)') (DKQHL;
#Ic-?2Gn4<
xlabel('\theta_i') %y \
4{=^J2z
ylabel('\phi') ]A:G>K
}xy[ &-dh
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Ykj+D7rA:
Ivc/g,
axis([0 90 -3.5 3.5]) !JwR[X\f
3gba~}c)
grid on i}LVBx"K(
~0 gHh
