利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 ^s*} 0
C AA~VEUL
1、光疏射向光密 XTJ>y@
jj2UUQ|
clear ~83P09\T%
IY&a!
close all ,d+mT^jN
%ab79RS]C
n1=1,n2=1.45; _<7FR:oBZ
ihVQ,Ct h
theta=0:0.1:90; v=>Gvl3&U
v(O.GhJ@
a=theta*pi/180; p3 o?_ !Z
._Xtb,p{
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); v2' J L(=
qu%s 7+
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); ?+\,a+46P_
A@OV!DJe]
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); Mm!saKT%
|9I; `{@
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); [ )A#9L~s=
~aG-^BAS
figure(1) *@nUas 2"
?
h%+2
subplot(1,2,1); $5r,Q{;$
)QD}R36Ic
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) [B o $?
A\)X&vR[6
legend('r_p','r_s','|r_p|','|r_s|') @U JmbD{
(T01 hR&
xlabel('\theta_i') p/~k w:I
rYQ@"o0/Y
ylabel('Amplitude') b3lpNJ J
X;:xGZ- oY
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) P} `1#$
*
U4:K@ y
axis([0 90 -1 1]) }o- P
,be?GAq
grid on @ `D6F;R
8A .7=C' z
subplot(1,2,2); $4L3y
uH
N~jQ!y
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) KT71%? P
I:<R@V<~#
legend('t_p','t_s','|t_p|','|t_s|') ;7k7/f :
8 zWPb
xlabel('\theta_i') nh )R
w.p'Dpw
ylabel('Amplitude') T P ::y
P*Sip?tdE
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ^6>|!
v!A|n3B]p
axis([0 90 0 1]) V&v~kzLr+
V JL;+
grid on ot#kU 8f
cV,Dl`1r
Rp=abs(rp).^2; q)+ n2FM
8 { %9%{
Rs=abs(rs).^2; M) XQi/
mp3_ n:R?
Rn=(Rp+Rs)/2; 6 JYOe
J6D$ i+
Tp=1-Rp; }F6b ]
Zb;$ZUWQX
Ts=1-Rs; PM!t"[@&
[Od9,XBa
Tn=(Tp+Ts)/2; \3?;[xD
h\ \2r>
figure(2) JBAK* g
M|e
n>P
subplot(1,2,1); }A/&]1GWk
T JS1,3<
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) [Vp2!"
bI6wE'h
legend('R_p','R_s','R_n') e7qMt[.
I} fcFL 8
xlabel('\theta_i') =2]rA
.t[ZXrd| 0
ylabel('Amplitude')
0 P3|1=
)Q/`o,Vm
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) LE80`t>M#
l'". }6S
axis([0 90 0 1]) K |} ]<
Z)T@`B6
grid on ` +BaDns
yi-"hT`
subplot(1,2,2); rRrW
0*/ r'
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) Qwo9>ClC
8kM0
legend('T_p','T_s','T_n') UJ)M:~O
! v VjZ
xlabel('\theta_i') AnE_<sPA
8421-c6y>
ylabel('Amplitude') `=$jc4@J
jsi# l
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) s5a OAyb*w
7a]Zws
axis([0 90 0 1]) #0<y 0uJ(y
;@YF} %!+W
grid on U2>dwn
bo ! ]
C\^< v&
1L =Qg4 H
6O@ ^`T
