利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 EzzzH(!j
B:+}^=
1、光疏射向光密 98o;_tU'
}wj*^>*
clear ukq9Cjs
uyjZmT/-
close all
CcAsJX~_
kDO6:sjR7
n1=1,n2=1.45; 8q_3*++D
rYN`u
theta=0:0.1:90; ^b>E_u
8* 6J\FE<p
a=theta*pi/180; A (;J
iWQBo>x
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); |(XV '-~
Wu.od |t0
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); !!Mp;h'}-
?_n.B=H`8
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ;hd> v&u#
Zb);08X
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); G.B^C)g uu
]~kqPw< R
figure(1) Mt@K01MI%
,,BNUj/:
subplot(1,2,1);
s .GTY@t
w[4SuD
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) VF A1p)n
D8[&}D4
legend('r_p','r_s','|r_p|','|r_s|') ,EwJg69
K=`*cSU>
xlabel('\theta_i') Qyd3 e O_
l*% voKZG
ylabel('Amplitude') 1`2lTkg
hw~c S7
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) S$muV9z2=
y7,~7f!N2
axis([0 90 -1 1]) 5\QNGRu"
:MJBbrV
,
grid on e$F7wto
{"{J*QH
subplot(1,2,2); WxO *{`T!
Z.<1,EKi=
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) 9h*$P:S;1v
6LvUi|~"<
legend('t_p','t_s','|t_p|','|t_s|') ] p0m6}B
G*J( 4~Yw}
xlabel('\theta_i') X8; 03EW;
GIZNHG
ylabel('Amplitude') dF1Bo
fCEd
:Kr
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 8
LaZ5
-P'>~W,~
axis([0 90 0 1]) rM>&! ?y+
f<kL}B+,Og
grid on JmN;v|wF:c
XTZWbhNF
Rp=abs(rp).^2; C:
e}}8i
bdcu O)3
Rs=abs(rs).^2; VT'$ lB%IK
xYwbbFGrG
Rn=(Rp+Rs)/2; c1 :op@t
ZU "y <
Tp=1-Rp; }T(|\
X
e,1Jxz4QH
Ts=1-Rs; 9 lA YCsX
Yq5}r?N
Tn=(Tp+Ts)/2; /9GqEQsfM
F1%'
zsv
figure(2) w_|WberU
$L`7(0U-
subplot(1,2,1); %Yd}},X_E
K9 0Zf
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) < W&~tVv
|:
nuT$(
legend('R_p','R_s','R_n') AvV.faa
88G[XkL$2
xlabel('\theta_i') UdLC]
Ho"FB|e
ylabel('Amplitude') IOmIkx&`GP
\t!+]v8f8
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) da[l[b;
%LVk%kz
axis([0 90 0 1]) 4Nq n47|>e
]L_HnmD6
grid on gNZ^TeT
4 6e;UUf!d
subplot(1,2,2); S2fw"1h*x
w"l8M0$m
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) q5OW1%
c~C :"g.y
legend('T_p','T_s','T_n') T>?sPq
/6S/a*`<X
xlabel('\theta_i') ,;D74h2F
$6'xRUx X
ylabel('Amplitude') .R-:vU880
`I'=d 4
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) z1V 0WDVm
PDpIU.=!0
axis([0 90 0 1]) w=^*)jZ8
UT^t7MY#O
grid on %@BQv 4oJ
}xdI{E1 q)
H%%#^rb^
M#|T Qa N
C ]r$
