利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 6P@K]jy& n
p,z>:3M
1、光疏射向光密 02EX_tt),
Zq33R`
clear U~BR8]=G
tOVTHx3E]
close all {=?[:5
92Gfxld\
n1=1,n2=1.45; /Zs_G=\>
d1.@v;
theta=0:0.1:90; 56YqYu.
j9c:SP5
a=theta*pi/180; el$@^Wy&$
GTbV5{Ss
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); U= GJuixy
5I[:.o0
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); Il\{m?Y
Y/UvNb<lK
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); V1Ft3Msq
93Gj#Mk
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); Z/ml,4e
B7qi|Fw
figure(1) hB<(~L?A]
bO$KV"*!
subplot(1,2,1); l-r$czY
xC
+>R1)
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) IX
6 jb"
?2o+x D2
legend('r_p','r_s','|r_p|','|r_s|') ^3vI
NF
eCJtNPd
xlabel('\theta_i') ;xh.95BP`
9qi|)!!L
ylabel('Amplitude') xv>8rW(Np5
i0\]^F
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) '_DB0_Dp
1M)88&
axis([0 90 -1 1]) 6E^m*la%
Zd)LVc[
grid on WUYU\J&q3
AWFq5YMSI
subplot(1,2,2); !+Cc^{
+5O^{Ce6
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) vfT
@;`
}|/<!l+;$
legend('t_p','t_s','|t_p|','|t_s|') D?xR>Oo)
m_ wvi
xlabel('\theta_i') mV}8s]29
`\z )EoI
ylabel('Amplitude') &!:mL],
`v(!IBP|
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) w=nS*Qy2
4Af7x6a;
axis([0 90 0 1]) ,3GB9
ZBK)rmhMx
grid on b:O_PS5h
0"\js:-$
Rp=abs(rp).^2; &-IkM%_A9
/i!/)]*-
Rs=abs(rs).^2; =mwAbh)[7n
u&`rK7J
Rn=(Rp+Rs)/2; w?fq%-6f*
FD~uUZTM
Tp=1-Rp; =yJc pj
bjwl21;{
Ts=1-Rs; Y!c7P,cZ+3
n|.>41bJ
Tn=(Tp+Ts)/2; hG<W*g
LChwHkRHJI
figure(2) yJ; ;&
DM)Re~*
subplot(1,2,1); HN;f~EQT
_w <6o<@
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) Y,4?>:39J
/sB,)>X
legend('R_p','R_s','R_n') !_>/ r
j=Q ?d]
xlabel('\theta_i') -q-BP}r3
S[/D._5QD%
ylabel('Amplitude') 44B9JA7u
LmR OG-9
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) maQDD*
U6nC
<3f
F
axis([0 90 0 1]) }&Xf<6
V22Br#+
grid on ;,KT+!H$
7bM
H
subplot(1,2,2); u]3VK
~^~+p
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) .'1]2/ad
bHs},i6
legend('T_p','T_s','T_n') 'A/{7*,
m0Uk*~Gz
xlabel('\theta_i') XPi5E"
?,Wm|xY
ylabel('Amplitude') riglEA[^
I'R|B\
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) srU*1jD)
:7qJ[k{g
axis([0 90 0 1]) hZO=$Mm4p
OBM&N
grid on C !Lu`y
?+{qmqN
3=lQZi<]%
pN;T t+}
>ufN[ab