利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
uRKCvsi sX 86~HkHliv 1、光疏射向光密
'z2} qJJ) _tL*sA>[~) clear
7ju7QyR !J[3U
close all
BSG_),AH hZ.Sj~> 7` n1=1,n2=1.45;
R)t"`'6| rSB"0 W7 theta=0:0.1:90;
{SW104nb&# $g
_h9L a=theta*pi/180;
3~BL!e, %$I\\q q>{ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
J}TS-j 0 Y*3qH ] rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
/H@k;o tsU.c"^n tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
s'ntf $ # @G! ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
g||{Qmr=1 '@wYr|s4 figure(1)
=+97VO(w]G e6k}-<W*q subplot(1,2,1);
'+Dn~8Y+9 xzy7 I6X plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
] ;^A8? 0kpRvdEr- legend('r_p','r_s','|r_p|','|r_s|')
`Qv7aY abWmPi xlabel('\theta_i')
on(F8%]zE 9C$b^wHd ylabel('Amplitude')
(} "r 5 WO)rJr!C title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
T^aEx.`O}` "4H&wHhT! axis([0 90 -1 1])
,$mnD@) 2<yi8O \ grid on
Dz]&|5'N +e%9P%[+ subplot(1,2,2);
z*l3O~mZ U{RW=sYB~9 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
;) 5d
wq 7h<Q{X<A legend('t_p','t_s','|t_p|','|t_s|')
$/5Jc[Ow /Bid :@R xlabel('\theta_i')
2K!3+D" K/+5$S jF ylabel('Amplitude')
bec n$R gf2l19aP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
B1JdkL 3h ,4jkTQ*@2 axis([0 90 0 1])
CwTx7
^qa Q{UR 3U'Q grid on
fr;>`u[; za1MSR Rp=abs(rp).^2;
4VSlgoz ?7@B$OlU Rs=abs(rs).^2;
,<zGvksk :qc@S&v@] Rn=(Rp+Rs)/2;
*O#%hTYq h$ D Fp Tp=1-Rp;
EJ.oq*W!*J 7qA0bU ee5 Ts=1-Rs;
X~sl5? MmI4J$F Tn=(Tp+Ts)/2;
(8qMF{ KIC5U50J figure(2)
JaX T
B"e efy65+~GG subplot(1,2,1);
JWr:/? #g MMh B= plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
wxkCmrV ,icgne1j legend('R_p','R_s','R_n')
?!U=S=8 A)kx,,[ xlabel('\theta_i')
2xhwi.u ^rfY9qMJr8 ylabel('Amplitude')
=\B{)z7@6D Y|_ #yb title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}tF/ca:XPQ ` ky<
* axis([0 90 0 1])
&dh%s Fy OKNA36cU' grid on
c$rkbbf~V X]U,`oE)9 subplot(1,2,2);
8V >j-C );_g2=:# plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
F^ 7qLvh o%i^t4J$e legend('T_p','T_s','T_n')
!jN}n)FSq `utv@9 _z xlabel('\theta_i')
#:s*)(Qn q&Y'zyHLP ylabel('Amplitude')
klxVsx%I{G mTz %;+|L title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
nBI?~hkP3 rmu5K$pl axis([0 90 0 1])
NZ djS9 J 05@SG': grid on
% o bR2% X^ckTIdR 
t.$3?"60~ "~(qp_AI 
OSQZ5 :g|
