利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 y3pr(w9A
Avi8&@ya
1、光疏射向光密 NRZ>03w
(f?&zQ!+
clear Dv[ 35[Yh
i*Ee(m]I
close all |csR"DOqz
.+Fh,bNYK
n1=1,n2=1.45; x@480r
pTk1iGfB
theta=0:0.1:90; "+:~#&r
,MdK "Qa>
a=theta*pi/180; ^PI8Bvs>j
7cMSJM(]G
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); _T[m YY
ay,E!G&H
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); ({#M*=&"
J,j!
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); I:G4i}mA
h?ZxS
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); iLD:}yK
b{wj4
figure(1) p@+r&Mg%W"
ds" q1
subplot(1,2,1); BV01&.<|
&o$Pwk\p/
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) x-HN]quhe
9WhZ=
Xk
legend('r_p','r_s','|r_p|','|r_s|') p!xCNZ(m
OFe?T\dQn
xlabel('\theta_i') +>wBGVvS
dr iw\
ylabel('Amplitude') yxz"9PE/P
8RJ^e[?o(
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) !*@sX7H
a@&P\"k
axis([0 90 -1 1]) /"%(i#<)xs
k2cC:5Xf3
grid on $D)Ajd;
vMB`TpZ
subplot(1,2,2); .V Cfh+*J#
c;dMXv
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) T,D(Xh
F6\Hqv
legend('t_p','t_s','|t_p|','|t_s|') A\v(!yg
O W|5IEC
xlabel('\theta_i') F+3}Gkn
efu'PfZ`&
ylabel('Amplitude') M'D l_dx-
z[`OYwsW
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) bYKe5y=
M#7w54~b?M
axis([0 90 0 1]) ',Q|g^rF]
#{BHH;J+
grid on awic9uMH
Ob#d;F
Rp=abs(rp).^2; M)JKe!0ad1
q&RezHK l
Rs=abs(rs).^2; \6)]!$F6:
Z/f%$~Ch
Rn=(Rp+Rs)/2; RU_wr<
88l\8k4r
Tp=1-Rp; t.`&Q|a
L"
GQQ
Ts=1-Rs; ?8GggJC
v{[:7]b_=
Tn=(Tp+Ts)/2; 4Lb!Au|Y
Jb (CH4|7
figure(2) >3MzsAH\
%qYiE!%&
subplot(1,2,1); &x YO6_.
KW1b #g%Z
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) %A^V@0K3
O;dtz\
legend('R_p','R_s','R_n') [ <,i}z
FP_q?=~rFs
xlabel('\theta_i') (/a#1Pd&
^.HvuG},O
ylabel('Amplitude') 6B=: P3Y
!5}u \
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ,|RN?1 ?U
H6t'V%Ys
axis([0 90 0 1]) iX|K4.Pz{
.;$Ub[
grid on TF1,7Qd
S<Os\/*
subplot(1,2,2); js..k*j
=G,wR'M
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) R ~ZcTY[8
?-Zl(uX
legend('T_p','T_s','T_n') PptVneujI
a-bj! Rs
xlabel('\theta_i') tg%#W`
sW^e D;
ylabel('Amplitude') m0Geq.
Q_r}cL/A
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) W2Ik!wEe&
8= kwc
axis([0 90 0 1]) }Be;YIhG
!*eDT4a
grid on yt@7l]I
\5iMr[s
1)o6jGQ
QJ,[K_
\]8F_K