利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
6 lN?) <uQ )\um "l*\c 1、光疏射向光密
=t@: F '&RZ3@}+ clear
F<h &3 sVE>=0TVP close all
I h5/=_n iSOD&J_ n1=1,n2=1.45;
_u#/u2< Tf9&,!> V theta=0:0.1:90;
W H{cJ7wCL 'wCS6_K a=theta*pi/180;
D9}d]9]$ !}iL O0 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
%T3j8fC{s "g
x5XW& rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
6>@( /mh* lf?Z{^ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
:jZ*,d%1={ # rh 0r` ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
zd?bHcW/h cFRSd
}p= figure(1)
{;]uL`abi? j09mI$2y67 subplot(1,2,1);
1j oc<EI p5lR-G plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
O{;M6U8C\ JA }S{ legend('r_p','r_s','|r_p|','|r_s|')
F@> w&A ~K VFe-#"0ZO xlabel('\theta_i')
X ulh.: N} 1.hOE>A% ylabel('Amplitude')
G/x3 wR W)6U6 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
hXV4$Dai > xxXPvM<` axis([0 90 -1 1])
NTv#{7q d_4n0Kh0 grid on
n`F QgC uKLOh<oio subplot(1,2,2);
rnzsfr-|(2 5pNvzw plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
8.Pcr< +~R.7NE% legend('t_p','t_s','|t_p|','|t_s|')
sRkz
WMl kcg)_]~6 xlabel('\theta_i')
EG Q1l i'B "wc`fg"3 ylabel('Amplitude')
cH qvkN` DI $ mD { title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~G|{q VO7A LM'*OtpDG axis([0 90 0 1])
kMwIuy ^L*VW
gi9 grid on
j8D$/ 73!
x@Duh Rp=abs(rp).^2;
Nl f&]^4(0 0C9QAJa Rs=abs(rs).^2;
D{&+7C:8. 0ER6cTo-t Rn=(Rp+Rs)/2;
uK"$=v6| (HTk;vbZm Tp=1-Rp;
d'**wh, .@x"JI> ; Ts=1-Rs;
2vW,.]95M @=aq&gb Tn=(Tp+Ts)/2;
;6\Ski0=l D6pEQdX` figure(2)
{ra Esb-X h!JjN$ subplot(1,2,1);
CwCo"%E8} z9uEOX&2\ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
% (O^as \8C*O{w legend('R_p','R_s','R_n')
-Z\UYt <O`q3u'l xlabel('\theta_i')
1W6n[Xg V)@nRJ g ylabel('Amplitude')
aaqjE
J&]
XLr.j title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
HzO0K=Z=R0 ]i_) :@ axis([0 90 0 1])
R!M|k%( j:vD9sdQ grid on
Do{*cSd 8Vg`;_ - subplot(1,2,2);
'!|E+P- N;Gf,pE plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
\gPNHL* { &JurZ legend('T_p','T_s','T_n')
,y{fqa4 (WU~e !} xlabel('\theta_i')
{(zL"g46 S)AE ylabel('Amplitude')
N?u2,h- zDQ\PZ~ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
tIp\MXkTQ& h 19.b:JT axis([0 90 0 1])
jWb\"0) 9#=IrlV 4 grid on
xBx ?>nN a!6OE"?QQ 
y O9pEO|W -<VF6k< 
I^C
]6D{
