计算光波由光密入射到光疏介质时相位变化 d+| !6
@ %B!$\]
clear x@RA1&c
E;Ftop
close all ,f""|X5
B0SmE_u_N
n1=1.45,n2=1; |~vQ0D
'C8=d(mR=m
theta=0:0.1:90; g"AfI
>Ti2E+}[M
a=theta*pi/180; 9^h%}>
vpw&"?T
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); Pw0 KQUs
QZq9$;>dW
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); v\tbf
uNGxz*e
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); CoN/L`.SN
F!cAaL1
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); KO;6 1y:
x;]{ 8#-z
arp=angle(rp); =
y,avR
;Z~.54Pf{d
ars=angle(rs); 0mi[|~x=
40sLZa)e
atp=angle(tp); PvBbtC-9b
w+(wvNmNEK
ats=angle(ts); s7.*o@G
5K-)X9z?
figure(1) (dt_ D
=}KbE4D+8
subplot(1,2,1); %{_
YJXpO
xa*gQ%+F
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) ti2_kYq
"8
mulE,
legend('arg(r_p)','arg(r_s)') 491I
e%Xf*64
xlabel('\theta_i') CXFAb1m
;I@L
ylabel('\phi') ll09j Ef
F|K4zhK
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) io'Ovhf:
V-r<v1}M
axis([0 90 -3.5 3.5]) vaZZzv{H
{4q:4i
grid on 0>MI*fnY"
Bb"4^EOZ,
subplot(1,2,2); 3pH`]m2
<~*Ol+/
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) j]D = \
!QspmCo+
legend('arg(t_p)','arg(t_s)') jch8d(`?d
<%7
V`,*g/
xlabel('\theta_i') sB/s17ar
\8aF(Y^H
ylabel('\phi') >7q,[:(gs
iRtDZoiD'
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) zL}hFmh
Z!eq /
axis([0 90 -3.5 3.5]) @&B!P3{f
9r#{s Y
grid on #L$ I%L"
AwrK82