利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 aGNb Cm
J,:Wv`N:9~
1、光疏射向光密 5a&BgBO1M
MY'T%_i d
clear aC94g7)`
[ J4n%
close all @% jY
@,m 7%,
n1=1,n2=1.45; XhUVDmeUMb
9[R+m3V/`
theta=0:0.1:90; rvuasr~
-"rANP-UI
a=theta*pi/180; nK }-^ Ur
7%Ou6P$^fr
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); QXW> }GdKZ
g x8i| ]
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); QmMA] Q
-S#j Or
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); w xJoWbn
m 48Ab`
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); YJ|U| [
"B>8on8O
figure(1) "U/yq
6^lix9q7
subplot(1,2,1); B=~uJUr
CB#B!;I8v
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) D MOP*;Uk
>qBJK)LHOv
legend('r_p','r_s','|r_p|','|r_s|') Xl:. `{5L
dQ_hlx!J
xlabel('\theta_i') p3> Md?e
rv`GOta*
ylabel('Amplitude') R90# T6^
4'TssRot@h
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) = ;(L$:l~
}0?XF/e(R
axis([0 90 -1 1]) Z/T( 4
I^HwXp([
grid on qeb} ~FL"o
vR&b2G7o
subplot(1,2,2); :| !5d{8S8
AiB] A}
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) oJQ
\?~
P EMBh?)g
legend('t_p','t_s','|t_p|','|t_s|') ZL:SJ,C
?L0 |$#Iw
xlabel('\theta_i') \}~7 1y}
ym+Ezb#o
ylabel('Amplitude') upZ tVdd
0Y:)$h2 ?
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 2x<!>B
z{ & z
axis([0 90 0 1]) pXvys] @
T$<yl#FY
grid on |QD#Dx1_
v^)B [e!
Rp=abs(rp).^2; [%O f
o/Q|R+yXV
Rs=abs(rs).^2; +z:>N l
r Nurzag
Rn=(Rp+Rs)/2; P)#h4|xZ
@;x*~0GZ
Tp=1-Rp; )+DDIq
97qf3^gGd
Ts=1-Rs; ?rA3<j
=XK}eQ_d
Tn=(Tp+Ts)/2; (G/(w%#7_
g5RH:]DV
figure(2) rxqSi0p
E-( $Xc
subplot(1,2,1); S^;;\0#NK
A7,%'.k
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) B"zB =A w
,i Y:#E
legend('R_p','R_s','R_n') bt(Y@3;
^B%c3U$o
xlabel('\theta_i') CyS %1 1L
w/"vf3}(9
ylabel('Amplitude') _aBy>=2c$
%-$BtR2@o
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) |a
a\t
i7Cuc+ j8
axis([0 90 0 1]) T?QW$cU!e:
,RM8D)m\
grid on ];"40 /X
.6LRg
subplot(1,2,2); 5gII|8>rQ
"}#%h&,
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) }+bo?~2E&
Jm#p!G+
legend('T_p','T_s','T_n') O!^ >YvOh
J3~%9MCJ
xlabel('\theta_i') {Z 7ixc523
u|T]Ne
ylabel('Amplitude') #oFyi @U
S,ea[$_
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) G ;iH.rCH
0[M2LF!m
axis([0 90 0 1]) .@%L8_sMR
Kh[l} ;/F
grid on _)~1'tCs}h
?QzL#iO }h
$ v[mIR
Shn=Q
+3o0GJ
