利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
6]#pPk8[Z Ay ZL( 1、光疏射向光密
s u ]x b]s.h8+v; clear
9$7tB ey K=F:GO close all
(rV#EA+6[` .du FMJl n1=1,n2=1.45;
{f;] i;!#:JX theta=0:0.1:90;
1Ao YG_ p`ai2`qC` a=theta*pi/180;
#w3cImgp2 YK Nz[x$| rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<
&[=,R0 @ ng $`<~=)\ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
iXpLcHi $CXKeWS=Q. tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
-}nTwx:|5u Gmo Y~}cg~ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
p3Uus''V4 Y <;A989D figure(1)
9l9h*P gt [ix45xu7 subplot(1,2,1);
M$j ]VZ wmh[yYWc plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
3#7D
g't .\1{>A legend('r_p','r_s','|r_p|','|r_s|')
_ _jFSa`at 6@_Vg ~=S xlabel('\theta_i')
u`Kj s}F' l n}2 ylabel('Amplitude')
-z%|
Jk NWCJ| title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
vr#_pu)f4 lTOO`g axis([0 90 -1 1])
ts r cX FL -yt grid on
rdd%"u+ G&FA~c subplot(1,2,2);
h\Gl yH~ bN-ljw0& plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
W ~sP7&sp &y-(UOqbkP legend('t_p','t_s','|t_p|','|t_s|')
gKK*`
L~ NIn# xlabel('\theta_i')
gGl}~ T@c{ 5a ylabel('Amplitude')
Z 6G>j kqJ \kd title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!
Ff/RRo L' w
} axis([0 90 0 1])
Y{~[N y E 5 "1kfB3v grid on
<0!/7*;#ZT k!Y7 Rc{" Rp=abs(rp).^2;
C^>txui8 UR;F W` Rs=abs(rs).^2;
$ K>.|\ <C0~7]XO Rn=(Rp+Rs)/2;
9\F:<Bf$# Qp9)Rc5 Tp=1-Rp;
(gd+-o4 JY4sB8 Ts=1-Rs;
fz
H$`X'M 5v:c@n Tn=(Tp+Ts)/2;
g~EJja ; /Q
Xq<NG figure(2)
~Dsz9 f
B fmSM 9 subplot(1,2,1);
"/v{B?~%! |FK ##8 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
|> mx*G =?Y %w%2 legend('R_p','R_s','R_n')
=bt/2 nPV Dk$[b9b xlabel('\theta_i')
NbPv>/r W10=SM } ylabel('Amplitude')
)%D2J C 59e q"08 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c=? =u qi !Nv$e axis([0 90 0 1])
EC *r d ~/0 t<^ grid on
G,+-}~ $_ e))fbv&V subplot(1,2,2);
.GG6wL<$? l (;~9u0sa plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
US <bM@[ /: -8 ,` legend('T_p','T_s','T_n')
:CH?,x^!@ EjZ_|Q xlabel('\theta_i')
1\G S"4~P <3aiS?i.h ylabel('Amplitude')
[kfLT::mT 7g'jg7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}A@ op+0E /b 3b0VfF axis([0 90 0 1])
QIZ }7 $]eU'!2) grid on
j<+Q Gd% sC
,[CN:b 
ySyA!Z !RX7TYf 
U&DD+4+28:
