计算光波由光密入射到光疏介质时相位变化 Skux&'N:
`q\v~FT
clear EW)r/Av:,
[a^<2V!vMn
close all D[YdPg@-
~g~`,:Qc
n1=1.45,n2=1; bhZ5-wo4%
q} g0-Da
theta=0:0.1:90; vL "n oLs
lL'K1%{+
\
a=theta*pi/180; aH50 0
t3Iij0b~
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); %a-fxV[
3`"k1W
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); k.Z? BNP
Jx4~ o{Z}c
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ] HG> O g
@$!"}xDR'
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); WCI'Kh
8Tc:TaL
arp=angle(rp); .e S* F
,fm{
krE
ars=angle(rs); w;Pe_m7\EO
_xP@kN~
atp=angle(tp); MF::At[4
1< M~ #
ats=angle(ts); ;/^O7KM -
z-S8s2.Fd
figure(1) f5n AD
qMBEJ< o
subplot(1,2,1); 6<n+p'+n
& #]||T-
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) Nn5sD3z#
F@X8a/;F-
legend('arg(r_p)','arg(r_s)') {JcMJZ3
kr]_?B(r
xlabel('\theta_i') [rhK2fr:i
+$G P(Uu,
ylabel('\phi') a'\fS7aE0l
.lppT)P
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) _3_o/I
!z?:Y#P3
axis([0 90 -3.5 3.5]) MXY!N /
S3 /Z]?o
grid on *>Zq79TG
Uw>g^[V;
subplot(1,2,2); #!M;4~Sfx
ezk:XDi4
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) 4*+ )D8
9KZLlEk5O
legend('arg(t_p)','arg(t_s)') X57\sggK
J,wpY$93
xlabel('\theta_i') sX= _|<[
Y3f2RdGl
ylabel('\phi') ^G(+sb[t
{UEZ:a
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) N%QVkuCbM
qznd '^[
axis([0 90 -3.5 3.5])
-B#>Jn#F
+P
9h%/Yk
grid on *ps")?tlC
v!? >90a
