利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
32z2c:G Ev ]oPCeA 1、光疏射向光密
JJ4w]Dd4 awU&{<,=g clear
E>i sl" /K :H2?J close all
',m!L@7M5 r< OqI*7 n1=1,n2=1.45;
M~l\rg8 4L<;z' theta=0:0.1:90;
7Sl"q=> }DFZ9,gQ a=theta*pi/180;
KmpKyc[ J*D3=5& rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
7-("pp YX= Ti>2N rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
->r udRQ .1F41UyL tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
w~ O)DhC 1k!$#1d< ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
n'9&q]GN| ;T3}#Q*qC figure(1)
rYO ~/N PwC^
]e subplot(1,2,1);
oD 3Q{ e jhB+ ] plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
(V<pz2\ Yv" -_ legend('r_p','r_s','|r_p|','|r_s|')
>uR;^ B5m u85?f xlabel('\theta_i')
:RDQP i Jb-F*_y ylabel('Amplitude')
%9b TfX" C *]XQ1F4 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`teaE7^Wm oH1]-Nl$ axis([0 90 -1 1])
JlE b @<z#a9 grid on
1O+$" 5H j$Vtd & subplot(1,2,2);
^w*&7.Z N 4w&g- plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
G5 *_ cv9-ZOxJ legend('t_p','t_s','|t_p|','|t_s|')
CO{AC~ 1?{w~cF} xlabel('\theta_i')
] 69 z-; no9=K4h` ylabel('Amplitude')
pykRi#[UrX MrhJk title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
q ],/%W xbs X -F axis([0 90 0 1])
K-n]m#U4o <5npVm grid on
`3L?x8g ^3ysY24 Q Rp=abs(rp).^2;
LsaR w-4.c A?H.EZ Rs=abs(rs).^2;
ni-4 ~k [cT7Iqip
Rn=(Rp+Rs)/2;
$o^N_`l uZ+vY F^ Tp=1-Rp;
)w0K2&)A N[wyi&m4 Ts=1-Rs;
Atod& qH -9yW f8 ; Tn=(Tp+Ts)/2;
9` G}GU]@} ,S-zY\XB figure(2)
Vm%0436wOY crU]P $a subplot(1,2,1);
DHh30b$c X -_0wR plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
X_#,5 t=7 )P9&I.a8 legend('R_p','R_s','R_n')
J>^KQ ^i6`w_ / xlabel('\theta_i')
7F8>w 7Y] ,e+S7 YX ylabel('Amplitude')
Z'_EX7r wu19Pg?F title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
jo ~p#l.' 6Lz:J:Q) axis([0 90 0 1])
gkld}t*U U_Am Riy grid on
#RP7?yGM, !\|L(Paf subplot(1,2,2);
B8 R&Q8Q T4x[
\v5d plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
O],]\M{GL 9FmX^t$T legend('T_p','T_s','T_n')
9P# <T7 @mu=7_$ U xlabel('\theta_i')
,{sCI/ tkf^sGgNO ylabel('Amplitude')
RhI>Ak;- )}4xmf@g l title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
S3&lkN5 Fes /8*- axis([0 90 0 1])
RyZy2^0< Q# sLIZ8= grid on
<9aa@c57 |H4f&& Wd 
H05U{vR Rx .dM_S 
0uS6F8x@
