利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 Fsj[J E
xR`2+t&t
1、光疏射向光密 |h/2'zd^-
uZI a-b
clear 9(j!#`O7&
??V["o T
close all 1Di&vpn0u
Q*(o;\s
n1=1,n2=1.45; hiHp@"l<
^T,cXpx|
theta=0:0.1:90; q:.BY}X9
.p>8oOp
a=theta*pi/180; &}1)]6q$
c.v)M\:
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); K_n%`5
c
@fc7
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); Q2?qvNZ
3/FB>w gt
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ;D:T
^4
_s8_i6 Y
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); p]4
sN
GK&Dd"v
figure(1) n\Ixv
HXI}f\6x
subplot(1,2,1); 90}B*3x
)zo ;r!eP
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) !d(V7`8
`f]O
legend('r_p','r_s','|r_p|','|r_s|') ]EQ/*ct
T1=M6iJ
xlabel('\theta_i') q3`t0eLZ
>k|[U[@
ylabel('Amplitude') e.V){}{V
aK'r=NU
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) mGvP9E"&
:f;|^(]"
axis([0 90 -1 1]) aDuanGC/V
gzF&7trN
grid on za7wNe(s
K#r`^aUc
subplot(1,2,2); E"=$p$k
Di*>PE@
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) cDg27xOUi
plfB}p
legend('t_p','t_s','|t_p|','|t_s|') p~IvkW>ln)
^ a:F*<D
xlabel('\theta_i') 'Iu(lpF&
`2B+8,{%
ylabel('Amplitude') *Y Ox`z!R
~_%[j8o&l
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) u:ISwAp
^iNR(cwgX
axis([0 90 0 1]) 0P(}e[~Z
rNc>1}DDS
grid on 7[LC*nrr
. {I7sUQ
Rp=abs(rp).^2; n:P}K?lg
zYgLGwi{
Rs=abs(rs).^2; kWFR(J&R
z61
o6mb
Rn=(Rp+Rs)/2; ?t-2oLE
|4vk@0L
Tp=1-Rp; o/R-1\Dn
X}$S|1CjO
Ts=1-Rs; G1rgp>m
^gG,}GTl
Tn=(Tp+Ts)/2; :[doYizk:
2 O%UT?R
figure(2) wr`eBPu
=E}/Z
subplot(1,2,1); *RPI$0
+\4=G@P.J
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) U|=y&a2Rb
-IE;5f#e
legend('R_p','R_s','R_n') NQmDm!-4
2;/hFwm
xlabel('\theta_i') <d4^gAfs*
eIJQ|p<v
ylabel('Amplitude') 9V~yK?
Wxjpe4
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) v!2`hqO
Oaui@q
axis([0 90 0 1]) mYCGGwD
Jy9&=Qh
grid on _>s.V`N'
fOfp.`n
subplot(1,2,2); g(1'i 1
y^ohns5{
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) Y3?kj@T`i
; ?!sU
legend('T_p','T_s','T_n') 5@Sb[za
J{H475GqiT
xlabel('\theta_i') q07>FW R
IcB>Hg5
ylabel('Amplitude') p~t$ll0s
@B+];lr/-
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) -
0zo>[c/p
.fgoEB,(
axis([0 90 0 1]) Js'|N%pi
:H~r
_>E
grid on 6 `'^$wKs
bkb}M)C
rS=6d6@
kns[b [!H
c8(.bmvF