利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 `W S
IKpx~
1、光疏射向光密 ts$UC $
/YR$#&N2
clear d0,F'?.0|
dPX>A4wp
close all `Pbn
>p:fWQ6
n1=1,n2=1.45; 2p[3Ap
|mA*[?ye@
theta=0:0.1:90; yln.E vJjD
|{"7/~*[
a=theta*pi/180; B 1.@K }
Lj,%pz J
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); T ua
@w+
l!'iLq"K(
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); K.zs;^
}H.vH
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ((q(Q9(F
:]jtV~E\
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); {`
t9{EO#o'k
figure(1) ajr8tp'
x{5*%}lX8
subplot(1,2,1); wUndNE
rP_)*)
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) z<*]h^!3
/Qy0vAvJ
legend('r_p','r_s','|r_p|','|r_s|') HPB1d!^
#1k,t
xlabel('\theta_i') h^,a 1'
#YdU,y=B
ylabel('Amplitude') | w -W=v
OwUbm0)h^V
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) =G6@:h=
nX'.'3
axis([0 90 -1 1]) !y.7"G*
_aYQ(FO
grid on y
6<tV.
;<H2N0qJ(
subplot(1,2,2); KG@hjO
(""&$BJQ|
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) eH6cBX#P.
RqR X
legend('t_p','t_s','|t_p|','|t_s|') +>j_[O5Y
vD t?N9
xlabel('\theta_i') g^FH[(P[G
?=&*6H_v
ylabel('Amplitude') )&K%Me
g>f394j
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) r9[J3t*({~
{rQ`#?J}^?
axis([0 90 0 1]) x`&W[AA4
5B6:pH6e
grid on `N+ P,
Vv]mME@
Rp=abs(rp).^2; |n;7fqK
A:(uK>5{Kk
Rs=abs(rs).^2; d|3[MnU[a
#q%/~-Uk
Rn=(Rp+Rs)/2; NT 'Y h
l]gfT&
Tp=1-Rp; Tec6]
:
X@rAe37h+
Ts=1-Rs; NG ~sE&,7
KC'{>rt7
Tn=(Tp+Ts)/2; va\cE*,@ns
-[z;y73]t
figure(2) dL>ZL1.$
A7se#"w
subplot(1,2,1); %|Vq"MW,I
XQ>m8K?\d
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) n7vi@^lf(
uP:'e8
legend('R_p','R_s','R_n') gueCP+a_
pB`<4+"9
xlabel('\theta_i') Flsf5 Tr0
ZC"p^~U_e[
ylabel('Amplitude') H`sV\'`!}
qmhHHFjQ
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) \TjsXy=:)
"Z
<1Msz
axis([0 90 0 1]) 3~ylBJJ
hz!.|U@,{<
grid on 0t8-oui
[||$1u\%
subplot(1,2,2); *=rl<?tX
{>#Ya;E
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) -4.+&'
+m_quQ/ys
legend('T_p','T_s','T_n') K\#+;\V
cOSUe_S0w[
xlabel('\theta_i') "b qB@)
F{v+z8nW
ylabel('Amplitude') lq74Fz&(
o}BaZ|iZ2
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) WE\V<MGS/
+NQw^!0qy
axis([0 90 0 1]) t2>Vj>U
'Q# KjY
grid on wml`3$"cf
5=eGiF;0\
r!+)U#8
|2Q;SaI^\
MOXDR