利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
`EW_pwZPA "oGM> @q=B 1、光疏射向光密
N X #/1= 9Z_OLai
clear
m">
=QP R@{/$ p: close all
C%7 ,#}[U/ z4%F2Czai& n1=1,n2=1.45;
hF%~iqd g 5YsV p theta=0:0.1:90;
1VM5W!} 3'NL1d u a=theta*pi/180;
+!Gr`&w*) 5 WNRo[`7 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
D0_x|a vrEaNT$J- rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
C36.UZoc K*i1! "w tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
r H_:7#.E |~vQ0D
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
'C8=d(mR=m g"AfI figure(1)
>Ti2E+}[M 9^h% }> subplot(1,2,1);
~;4k UJD wk 7_(gT`0 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Xv(9 Yh S wuC tg= legend('r_p','r_s','|r_p|','|r_s|')
m6ws #%|[ :[rKSA]@ xlabel('\theta_i')
uTloj . 8q LgB
ylabel('Amplitude')
u! FSXX< .7 ^-*HT} title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
s
!vROJ Yxq Q g axis([0 90 -1 1])
lTd2~_ <UL|%9 =~ grid on
4E(5Ccb -"tgEC\tD subplot(1,2,2);
NB#*`|qt hd
BC ^n plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
>43yty\
~F6gF7] z legend('t_p','t_s','|t_p|','|t_s|')
?B!ZqJ # ^W05Z!} xlabel('\theta_i')
JX<W[P>M @{a-IW 3 ylabel('Amplitude')
YGmdiY:;1 j7 3@Yi% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
P&^7wud-sb E.bbIV6mQ axis([0 90 0 1])
9>>}-;$ ;i?!qB>baX grid on
1}M.}G2u/ )t=u(:u] Rp=abs(rp).^2;
L=FvLii. cb,sb^- Rs=abs(rs).^2;
j}* +-.YF # Kr.!uD Rn=(Rp+Rs)/2;
xA E@cwg kw}1 CXD Tp=1-Rp;
<\EfG:e (:x"p{ Ts=1-Rs;
}]?G"f
t K jP'b! 4 Tn=(Tp+Ts)/2;
*?"{T;4u~O e;[8 GE.
figure(2)
3) 0~: AAY UXY! subplot(1,2,1);
wKbymmG @&B!P3{ f plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
9r#{s Y #L$ I %L" legend('R_p','R_s','R_n')
AwrK82 .p}Kl$K ] xlabel('\theta_i')
Q0_W<+` Cw5K* ylabel('Amplitude')
+N@F,3yNa &/?j MyD@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uy:=V } p 5^/[] * axis([0 90 0 1])
} O wCmv/m grid on
& 2igX?60 .x8$PXjPG subplot(1,2,2);
[ 1GEe D:EF@il plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
H\bIO!vb - S-1<xR legend('T_p','T_s','T_n')
TMsoQ82 %MNV 5UA[w xlabel('\theta_i')
;# j 82 i`'^ zR(`i ylabel('Amplitude')
Ti'kn{
Zv Fo~v.+^? title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^W'[l al. ]L^M7SKE6 axis([0 90 0 1])
]S geZ07 AoeW<}MO grid on
efR$s{n! /)TEx}wk 
R ?s\0 Bdepvc}[# 
,[ L$
