利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
E( Z8 }n4V|f- 1、光疏射向光密
[D<1 CF /\4'ddGU clear
_Nqt21sL !0vG|C ;' close all
*tK\R&4,4s '2[albx Sc n1=1,n2=1.45;
:gep:4&u
&2{ tF theta=0:0.1:90;
B Ere*J *q-['"f a=theta*pi/180;
_N{RVeO v%iof1 T'
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
w t iny,6 IX>d`O61*g rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*1fb}C_ Aj+2;]M tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
):b$xNn }x07^4$j ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
V0P>YQq9s M!REygyx figure(1)
`Ch9~* p -)RH5WG S subplot(1,2,1);
wfq7ob4^ \cUC9/
b plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
)|DM~%$QM E: $P=%b legend('r_p','r_s','|r_p|','|r_s|')
mj KS{ r} mb Xvn xlabel('\theta_i')
J
/f
.Z JRO>S ylabel('Amplitude')
}wHW7SJ t3&LO~Ye title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
&GB:|I'%7 u (em&M axis([0 90 -1 1])
'U\<IL#U hNH'XQxO grid on
Y}#J4i0b* mg$]QnbAnH subplot(1,2,2);
i03w 1pSH, VrxQc qPr` plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
<7j8 7 ! 3 f?:M legend('t_p','t_s','|t_p|','|t_s|')
q@nP}Pv&5 j]aoR xlabel('\theta_i')
c0;t4(
&8 Z(ACc9k6:' ylabel('Amplitude')
o! {w"K t\%HX.8[;% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+KvU$9Ad> [`c^ 4 E axis([0 90 0 1])
xK6n0] A N#DYJ-~* grid on
EB2!Hp uQ3 a6WI170^1 Rp=abs(rp).^2;
bdF. qO9
3
%|86:* Rs=abs(rs).^2;
t 0|!(3 17G'jiY H Rn=(Rp+Rs)/2;
[N#, K02mk 6u7?dG'4 Tp=1-Rp;
b{]z
w pf Z^IPZF Ts=1-Rs;
/eOzXCSws ]2\VweV Tn=(Tp+Ts)/2;
] 0B2#
d ),ma_{$N figure(2)
j{9D{ !VI]oRgP subplot(1,2,1);
<oX7P69 h >s!K9 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
4bn(zyP #R.- KUW: legend('R_p','R_s','R_n')
( 8c9 /7h kT'u1q$3Vo xlabel('\theta_i')
Gq }U|Z ; o 0&`b? ylabel('Amplitude')
s#t Zg !=:$lzS^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
TG+VEL |T k+8q{5>A< axis([0 90 0 1])
^]a #7/]o li 6%) grid on
7TDy.] Wu_kx2h subplot(1,2,2);
?`#/ 8PN 30.@g[~ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
iMP]W _ !Y\D?rKZ legend('T_p','T_s','T_n')
@p NNq dRaNzK)M xlabel('\theta_i')
FcYFovS 7 El[ > ylabel('Amplitude')
/(BMG/Tb Hqn#yInA7~ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/gu%:vq (;57 Vw axis([0 90 0 1])
70;Jl).\{ lAJxr8 . grid on
'_/Bp4i cdZ~2vk 
Bxfc}vC. d~`x )B( 
$[/&74#0HX
