利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
obGhO !MrQ-B ( 1、光疏射向光密
w=CzPNRHH! J0
dY%pH# clear
l[]cUE *T#^|<.XG close all
[!>9K}z,= c:52pYf+ n1=1,n2=1.45;
qco uZO 8{]nS8i theta=0:0.1:90;
o <J6KTLv (c&%1bJ a=theta*pi/180;
, a?oGi Fc8 0HK5R rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
, O/IY i6A9|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));
:<s`) xm Eom tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
eX>* }pI SSxz1y ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
o
^ \+Ua Q-!gO figure(1)
REg M g1?9ge 1 subplot(1,2,1);
*`|xa@1v` $ _,-ES I plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
P;@j Z:<an+v|5 legend('r_p','r_s','|r_p|','|r_s|')
3z, Ci$[ "OA{[)fw" xlabel('\theta_i')
=-_hq'il /M::x+/T ylabel('Amplitude')
vG.KSA n1b:Bv4"]# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
AC'_#nPL# >d=pl}-kOQ axis([0 90 -1 1])
mhX66R |`Noj+T47I grid on
fN&uat 7 #AB5}rPEI subplot(1,2,2);
j}B86oX }IZ w6KiN plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
-|^)8 b1cVAfUP legend('t_p','t_s','|t_p|','|t_s|')
Ncsh{. 4xq| xlabel('\theta_i')
N6of$p'N Y)]C.V,~ ylabel('Amplitude')
L -:@Om! 0 }qlZFB title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<K<#)mcv 09anQHa axis([0 90 0 1])
|,5|ZpgL zRD-[Z/- grid on
A7%/ sMv )qw;KG0F Rp=abs(rp).^2;
D*[J rq, a*LfT<hmU3 Rs=abs(rs).^2;
DpvHIE:W &Jb\}c } Rn=(Rp+Rs)/2;
@y~kQ5 k U|Gy 9" Tp=1-Rp;
[:#K_EI5% -y$6gC RY Ts=1-Rs;
P_NF;v5 v G%: 3.:E" Tn=(Tp+Ts)/2;
DBB&6~;? jLt3jN figure(2)
![_0 GFbT M jQju@ subplot(1,2,1);
pvUV5^B(M -,b+tC<V)0 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
!rr,(!Ip?O 0"% dPKi legend('R_p','R_s','R_n')
h*C!b?:" L;$>SLl, xlabel('\theta_i')
N'9T*&o+ 1w(3!Ps+ ylabel('Amplitude')
=WW5 H\? vg[3\!8z[ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
RU4X#gP4Vh o.A:29KoU axis([0 90 0 1])
sAqy(oy#M J](NCD grid on
86 e13MF nL/]Q'(5 subplot(1,2,2);
mc8Q2eQat} h2f8-}fsq plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
$7DW-TA ;z.niX .fx legend('T_p','T_s','T_n')
~Ay)kv; dB[4N T xlabel('\theta_i')
~ [t#$2d} #8"oqqYi ylabel('Amplitude')
"%Y= + X: g5;NT title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
m","m J'|=J axis([0 90 0 1])
dGBjV #bNT Y06^ M?} grid on
n]'
r3 b; UDgq8v 
oH%[8!# [QgP6f]= 
IUv#nB3
