计算光波由光密入射到光疏介质时相位变化 g/$R uT2U
U+B{\38
clear F7C+uG Ts
6rE8P#
close all
rwSR
l =P)$O|=w
n1=1.45,n2=1; KI\bV0$p<
H'$H@Kn]-
theta=0:0.1:90; GBS+ 4xL|
M*T!nwb
a=theta*pi/180; Qo3Enwap=
t
4M-;y
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); (g[WZB3x
<6.`(isph
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); wU'+4N".
N?rE:0SJ
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); h]i vXF*
PQ_A^ 9 5
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); W +GBSl
%b_0l<+
arp=angle(rp); 2H8\ P+
4yTg H0(T
ars=angle(rs); dhP")@3K;p
g*_ n|7pB
atp=angle(tp); PS!or!m
k&[6Ld0~56
ats=angle(ts); T%"wz3~
Jc]k\U
figure(1) 4,YL15.
q 3v5gz^t
subplot(1,2,1); (u@X 5O(a
Y0Tw:1a
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) _ A{F2M
A ^hafBa
legend('arg(r_p)','arg(r_s)') )iC@n8f7o
I%Z=O=
xlabel('\theta_i') V#Y"0l+~
#_3ZF"[zq
ylabel('\phi') eq"Xwq *
l*CCnqE
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Y HS/|-
' qT\I8%
axis([0 90 -3.5 3.5]) ][//G|9
iM1E**WCtv
grid on 3H>\ hZ
L"c.15\
subplot(1,2,2); :T@} CJ
fpO2bD% $8
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) ]sP
%ib7)8Ki0
legend('arg(t_p)','arg(t_s)') vq3:N'
FH%: NO
xlabel('\theta_i') h'vBWtMa
3)e{{]6
ylabel('\phi') 211T}a
[T [] U
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) fX\y/C
Dequ'
axis([0 90 -3.5 3.5]) i3L2N~:V
2zv:j7
grid on JXt_
^e\$g2).
