利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
8TnByKZz W:16qbK 1、光疏射向光密
jVpk) ;vC 3 _tO clear
]gQ 4qu5 7) Rx- close all
TBYRY)~f U8(Rye$ n1=1,n2=1.45;
,&IBj6%Y #ab=]}2W_g theta=0:0.1:90;
5? s$(Lt~ ctL@&~*nY a=theta*pi/180;
]l'W=_XDg w(9.{zF|vQ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
u F T&r| ]S~Z8T-[ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
YFCP'J"Z &V&0kp@+ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
#` mo5 V4OhdcW{ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
dSb|hA}@ ZAPT5 figure(1)
v&B*InR?+ *_ z5Pa`A subplot(1,2,1);
y$At$i>u mq?5|` plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
yjVPaEu]aU D/Y .'P:j legend('r_p','r_s','|r_p|','|r_s|')
p_jDnb# g(Jzu' xlabel('\theta_i')
u9FXZK7 qlSc[nEk ylabel('Amplitude')
|Y!#` !0~$u3[b title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
XX,iT~+- cRX0i;zag axis([0 90 -1 1])
z1ltc{~Z pCNihZ~ grid on
)dJaF#6j Qc z7IA subplot(1,2,2);
$@t-O or; I#kK! m1Q plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
+!V*{<K 92GO.xAD? legend('t_p','t_s','|t_p|','|t_s|')
dm& F1NkT 8Z!+1b xlabel('\theta_i')
OZ1+` 4 v K^k1]!W= ylabel('Amplitude')
02} &h oQ 5g0(J~ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{b>tX)Tep a2*WZc` axis([0 90 0 1])
Xz0jjO, U%ce0z grid on
-L6 rXQV@j N
G1]!Vz5 Rp=abs(rp).^2;
)kl| 5 i H>e?FDs0*R Rs=abs(rs).^2;
d+%1q 8h~v%aZ1 Rn=(Rp+Rs)/2;
4em7PmT =8j;!7 p Tp=1-Rp;
TdP_L/>|J s@f4f__(] Ts=1-Rs;
#H(|+WEu (
=->rP Tn=(Tp+Ts)/2;
,*r"cmz C8y 3T/G figure(2)
~ -Rr[O=E VY@uQ#&A subplot(1,2,1);
*"0Yr`)S L._I"g5 H9 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
cn=~}T@~Z \w^i SK- legend('R_p','R_s','R_n')
` &|Rs \i "I1xU xlabel('\theta_i')
=R!=uml( 52tIe|K wL ylabel('Amplitude')
6!/e_a ";9cYoK RY title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}+QgRGQ ,>2ijk# axis([0 90 0 1])
J & +s B @UaaWh grid on
FgNO# % R* E/E subplot(1,2,2);
4>{q("r, ' G -]> plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
NyJU?^f&v Et%s,zeA{2 legend('T_p','T_s','T_n')
oKz|hks[6 =XJ
SE+ 7 xlabel('\theta_i')
?.IT!M}DR pQ*9)C ylabel('Amplitude')
$uawQf+S r `i<XGPJ% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
u=jF\W9 F^IYx~: axis([0 90 0 1])
c&3
]%urL vd SV6p.d grid on
9]VUQl9gh e1S |&W8 
IQoz8!guh: Jj|HeZ1C f 
$/5Jc[Ow
