计算光波由光密入射到光疏介质时相位变化 Rl&nR$#
_Q 'f^Kj
clear g:Q:cSg<
PeIKx$$Kl{
close all 85e*um^
eA *We
n1=1.45,n2=1; fr'DV/T
b;F aTm@
theta=0:0.1:90; n-g#nE c:
+p [O|[z
a=theta*pi/180; W[ R`],x`
wrc1N?[bn
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); Fi/`3A@68
&@FufpPw/
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); z%BX^b$Hj
Jd~M q9(
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); &dPUd ~&EL
Cn0s?3Fm
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); 1]9l
SE!E7
:|8!w
arp=angle(rp); ;66{S'*[
*pDS%,$xe
ars=angle(rs); $&!|G-0'
#l h'
!
atp=angle(tp); 3,EtyJ3[Bh
-BSO$'{7
ats=angle(ts); Khl0 ~
2I| lY>Z
figure(1) M)U)Sc zHO
g6+5uvpd
subplot(1,2,1); N f)SR#;
[aG
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) z K I1
R #i`H(N
legend('arg(r_p)','arg(r_s)') 0d Q\Y]b
QyY<Z i;6
xlabel('\theta_i') YI;MS:Q j
qA[cF$CIl)
ylabel('\phi') )c?nh3D
8)2M%R\THn
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) z` eMb
+z0s)HU>j
axis([0 90 -3.5 3.5]) xB]^^ NYE=
OI8} v
grid on 0x<G\ l4
GHo
mk##0E
subplot(1,2,2); ;U$Rd,T4S
y2mSPLw
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) 52 NI{"
LM 1Vsh<
legend('arg(t_p)','arg(t_s)') x8x-b>|$&<
Jl6lZd(Np
xlabel('\theta_i') L4ct2|w}ul
\j-:5M#m
ylabel('\phi') Bj"fUI!dK
<:&{ c-f/
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) lauq(aD_C
ZD7qw*3+
axis([0 90 -3.5 3.5]) ,b5vnW\
N7KG_o%
grid on ^.
P RNq8nmxC
