利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
8[Qw8z5- W=b5{
6 1、光疏射向光密
urN& ."c k^L (q\D clear
k~gQn:.Cx &sllM close all
RHBEC@d[} *kIJv?%_} n1=1,n2=1.45;
&sKYO<6K } 8e-{S~@W theta=0:0.1:90;
bw[!f4~ 1TVTP2&Rd a=theta*pi/180;
OW.ckYt% IS!+J.2 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
(>jME m&x0,8 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
7}#vANm (PNvv/A tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
FxUH ?%w a2i
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
l^vq'<kI s)N1@RBR figure(1)
OO$<Wgh ^NCH)zK]v subplot(1,2,1);
AV'> tQ/ w\6{ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
wS5hXTb" dfrq8n] legend('r_p','r_s','|r_p|','|r_s|')
-py.Y Z J+;.t&5R xlabel('\theta_i')
L.GpQJ8u XI0O^[/n{ ylabel('Amplitude')
JvUKfsn u{ 87HV D Di title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"<&F=gV saV3<zgx axis([0 90 -1 1])
OVd"'|&6_ hsl8@=_ B grid on
;?y?s'>t& ?;*mSQA`J subplot(1,2,2);
55;xAsG $v^F>*I1 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
,4\vi| |%tR#!&[:g legend('t_p','t_s','|t_p|','|t_s|')
v-l):TL+= Y,8M[UIK xlabel('\theta_i')
F| PYDC FCI T+ 8K ylabel('Amplitude')
>GjaA1, 9+/<[w7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
N(
/PJJ~ fLy s$*^)^ axis([0 90 0 1])
x=H*"L= T.{I~_ grid on
%xg"e
O2x <1@_MY o Rp=abs(rp).^2;
:l6s ESr Y, 0O&' > Rs=abs(rs).^2;
{
~Cqb7 4Q]+tXes Rn=(Rp+Rs)/2;
[<%yU y O,bj_CW x Tp=1-Rp;
/yI~( 8bO Bb1d H/8 Ts=1-Rs;
GoF C!nx jR48 .W Tn=(Tp+Ts)/2;
$EtZ5 ?qS h)Yq C$A-s figure(2)
! g}9xIL 0h; -Yg subplot(1,2,1);
zX5p'8- ]v lQNd? plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
")LF;e c XY!b=9 legend('R_p','R_s','R_n')
j$Kubg(I5 ];< [Cln% xlabel('\theta_i')
#ucOjdquq u?aq'
"t ylabel('Amplitude')
}
HvVL}7 ufm#H#n)#X title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7lh%\ Z|j8:Ohz axis([0 90 0 1])
E-Mp|y /V \Q3m?)X=Gd grid on
>#V8l@IH 7tyn?t0 n subplot(1,2,2);
O 4}cv 8Mp plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
FUHjY (C. 1'<] legend('T_p','T_s','T_n')
1 ltoLd\{ ;/YSQt)rc> xlabel('\theta_i')
HFf|
>&c& UL{Xe&sT ylabel('Amplitude')
xXyzzr1[ 7g=Ze~aq title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
M"P$hb'F a1Gy I axis([0 90 0 1])
Xe%n.DW m R!,RZ?|v grid on
zL> nDnL 4 MLp5Y\8* 
xxedezNko "r|O / 
OCX?U50am
