利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 L+L9)8FJ
oMbCljUC
1、光疏射向光密 YSif`W!
B"v=Fr[
clear p3' +"sFU
p2K9R4
close all +n dyR
%54![-@
n1=1,n2=1.45; < r~Tj
U_No/$ b
theta=0:0.1:90; ?da 3Azp
}d(6N&;"zN
a=theta*pi/180; VUb*,/hxa
%ZK}y{u\
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); *gn*S3Is[j
.$iIr:Tc>
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); n6f|,D!?
tDo0Q/`
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
f/ZE_MN2
KX`,7-
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); w6_}]
&F
guvQISQlY
figure(1) E_oe1C:
51vK>
subplot(1,2,1); W#!\.m`5
:-)[B^0
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) !MC Wt
q}jf&xUWzH
legend('r_p','r_s','|r_p|','|r_s|') c
z|IBsa*
"^H+A-R[
xlabel('\theta_i') D }\`5L<
hd\#Vh(H
ylabel('Amplitude') ?6dtvz;K+?
,l6W|p?ZO^
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 5kF5`5+Vj
dQ[lXV[}v
axis([0 90 -1 1]) w9%gaK;
<_![~n$H
grid on B.?yHaMI[
5IMH G%W7
subplot(1,2,2); d,77L
P*SXfb"HC
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) :.cX3dP@
DQ)SMqOotw
legend('t_p','t_s','|t_p|','|t_s|') o;.-I[9h]
IQDWH/c
xlabel('\theta_i') #lO;G
k{
g:Dg?_o
ylabel('Amplitude') $sJn:
8z
&a`-NRU#
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ?,j:Y0l.L
dZW:Cf 9K
axis([0 90 0 1]) 89paR[
{x8`gP\H
grid on M Pt7 /
(GQy"IuFh
Rp=abs(rp).^2; b%kh:NV{S
+3C
S3fTq
Rs=abs(rs).^2; L6a8%%`
Y%faf.$/9
Rn=(Rp+Rs)/2; 1pV"<,t
C$bK!]a
Tp=1-Rp; V/Q6v
YX
/Yh8r1^2tZ
Ts=1-Rs; ^tah4QmUA
u7mj
Tn=(Tp+Ts)/2; *M"wH_cd
g@QpqrT
figure(2) =4zsAa
MiC&av
subplot(1,2,1); 6"DvdJ0MB
#'T|,xIr-Q
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) G>,rf
]N
3EyN"Lvp{o
legend('R_p','R_s','R_n') E8xXr>j>#
"CaVT7L
xlabel('\theta_i') |0&S>%=
4Mprc~ 7vr
ylabel('Amplitude') bfJDF(=h
vmoqsdZ/
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) K)qmJ-Gub
0h7\zoZ5
axis([0 90 0 1])
={fi&j
Pwq}
;+
grid on MwuH.# Ez
ESs)|t h
subplot(1,2,2); {e~d^^N5
&]euN~y
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) ;'h7
j*6
(p. 5J
legend('T_p','T_s','T_n') ~7ArH9k.
HZ|6&9we
xlabel('\theta_i') d1,azM
G67BQG\av
ylabel('Amplitude') BAxZR
*) wp
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) mz#(\p=T
wy^>i$TC
axis([0 90 0 1]) _])1P?.
H1l'\
grid on `S!`=26Z!
^WeT3b q
q%JV"9,
:g,r l\S7
\F>
*d!^C