利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 1">]w2je:
d<l-Ldle
1、光疏射向光密 blA]z!FU
\ |9KOulr
clear r.;(Kx/M
hDs.4MZC`
close all `)R@\@jt
QZ54Osdl
n1=1,n2=1.45; to*<W,I
)04lf*ti
theta=0:0.1:90; IRQ3> 4hI
=xX\z\[A
a=theta*pi/180; R!.HS0i.
dc0Ro,
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); id.W"5+
0IEFCDeCO
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); E5 0$y:
P'6(HT>F?
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); 7a:mZ[Vh
(FAd'$lhX}
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); <XxFR
;DKwv}
figure(1) $zBG19 [%
3f:I<S7
subplot(1,2,1); HSOdqjR*
eMvb*X6
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) @.X}S"yr
M_4g%uHG
legend('r_p','r_s','|r_p|','|r_s|') (&$|R\W.
ft[g1
xlabel('\theta_i') 1W5\
+B@NSEy/+
ylabel('Amplitude') #={L!"3?e
=#<hT
s
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) <z|? C
%d9UW Q
axis([0 90 -1 1]) }x1mpPND
#7U,kTj9
grid on soZw""|v
j]<T\O>t>
subplot(1,2,2); R6Cm:4m}I
_|Ml6;1aZ
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) ]Nvtiw 6
|Tz4 xTK
legend('t_p','t_s','|t_p|','|t_s|') k[Iwxl;/
v[6 BESu
xlabel('\theta_i') |pp @
&hba{!`y
ylabel('Amplitude') Y(SgfWeK@1
~]/X,Cf
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) N)h>Ie
XI\aZ\v
axis([0 90 0 1]) 7Yxy2[
G6eC.vU]j
grid on Ik1,?A
4T9hT~cT7
Rp=abs(rp).^2; ZZE
fu=}E5ScK
Rs=abs(rs).^2; RQU5T 2,
_:'m/K3Ee
Rn=(Rp+Rs)/2; 2>O2#53ls0
=,[46 ;q
Tp=1-Rp; i:kWO7aP
P+3G*M=}
Ts=1-Rs; 0 \LkJ*i
_|TE )h
Tn=(Tp+Ts)/2; G-3.-
%T6#c7U_
figure(2) cJ8F#t
?GFxJ6!%I
subplot(1,2,1); d0 V>;Q
*-|+phim
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) Hs=!.tZ,
cT0utR&
legend('R_p','R_s','R_n') R~w(]
m4Wn$Z
xlabel('\theta_i') YF>t {|
;6b#I$-J-
ylabel('Amplitude') faO8
&
HB`pK'gz
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Xe+&/J5b
+X*`}-3
axis([0 90 0 1]) D hZtiqL#_
N0vd>b
grid on Xp} vJl
{sm={q
subplot(1,2,2); Y[~6f,?^
RaU.yCYyu
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) X^|oY]D
}=;>T)QmMO
legend('T_p','T_s','T_n') u=E?N:I~F
=|U2 }U;
xlabel('\theta_i') @i{JqHU"
n~&R_"mv(
ylabel('Amplitude') gs8L/veP
<go~WpA|r
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) @6j*XF
99ha/t
axis([0 90 0 1]) 7lVIN&.=
y{<#pS.
grid on gw*d"~A
kuH;AMdv
D@^F6am%
Pm1
"
0
/M3D[aR<d