利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 2@ S}x@^
GDLw_usV
1、光疏射向光密 8 lQ}- 8
%`]+sg[i
clear x/,;:S
Y j oe|
close all oc1BOW z
dN2JOyS
n1=1,n2=1.45; :^7w
sVyV|!K
theta=0:0.1:90; 0F[ f%2j
Z],j|r Wy6
a=theta*pi/180; LS~at.3zX
]\*^G@HA2
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); .n<vhLDQn
[c{\el9H
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); H07\z1?.K
bq/Aopfr
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); >oW]3)$4S
3aW<FSgP
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); YTTyMn
G9:XEEN
figure(1) . `Rt
@&\Y:aRO%i
subplot(1,2,1); oC5 h-4~
Hz TmNm)
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) 7
$y;-[E[
l<6G Z
legend('r_p','r_s','|r_p|','|r_s|') ceUe*}\cr
XPd mz !,b
xlabel('\theta_i') 01r%K@ xX\
x9 YQd69
ylabel('Amplitude') 5%} e j)@
$d*9]M4
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 8w5}9}xF
]oP1c-GEk
axis([0 90 -1 1]) ? i _ACKpw
</ [.1&S+\
grid on C hF~
k%6CkC w
subplot(1,2,2); 2+9VDf2
zsp%Cz7T
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) l*l*5hA
E@G Yl85fI
legend('t_p','t_s','|t_p|','|t_s|') >pF* unC;
UF+Qx/4h0
xlabel('\theta_i') g~DuK|+
3^/w`(-{@
ylabel('Amplitude') M@ U >@x;
? c#s}IH
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) vhBW1/w&F
y C xYFi
axis([0 90 0 1]) E0ED[d,
gGrVpOzBj
grid on 0he3[m}Nr
X.b8qbnq[
Rp=abs(rp).^2; gH//@`6
iVFOOsJ@
Rs=abs(rs).^2; >ai,6!
{;{U@Z
Rn=(Rp+Rs)/2; VM$n|[C~
t'U=K>7
Tp=1-Rp; DP{nvsF
9L>?N:%5
Ts=1-Rs; O=jLZ2os
#
55>?
Tn=(Tp+Ts)/2; BA h'H&;V
YYQvt
figure(2) +(*HDa|
=+iY<~8
subplot(1,2,1); t 'eaR-
cQEUHhRg!
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) !+SL=x y!{
RRja{*R
legend('R_p','R_s','R_n') L"b&O<N o
0x<ASfka
xlabel('\theta_i') {T8;-H0H
lsax.uG5x
ylabel('Amplitude') $9G&
wH>{
KZ 5%q.
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 'C5id7O&
':n`0+Eh
axis([0 90 0 1]) |S!R Q-CF
o898 pg
grid on j:%,l cF
&GLDoLk6[
subplot(1,2,2); ,[;O'g?,g
w-Ph-L/
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) %r|sb=(yT
Q+Bl1xl
legend('T_p','T_s','T_n') $9YQ aN%
9Jwd *gevV
xlabel('\theta_i') 3H,x 4L5j
wa[L[ mw
ylabel('Amplitude') TdPd8ig8{
wa@Rlzij>
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) {8.Zb NEJ
EGl<oxL*R2
axis([0 90 0 1]) "lt <$.
{dF@Vg_n
grid on qxI $F
"w? 0f["
%V=%ARP|
:of([e|u6
<{W{
Y\_A>
