利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
^i{B8]2, @o+T<}kW X 1、光疏射向光密
S,5>g07-` lLHHuQpuj clear
DytOS}/^9 3zb;q@JV close all
%PG::b ZkYc9!anY n1=1,n2=1.45;
Vgs( feGs z[Kxy1, theta=0:0.1:90;
]{OEU]I@ &7 }!U a=theta*pi/180;
%E7+W{?*1 nzDS rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
NDB ]8C Z*kGWL rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
\n850PS [;VNuF tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
=JyYU*G4 {e2 ( ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
a#1LGH7E8 CF^7 {g(y_ figure(1)
)J_!ZpMC ]Bs ? subplot(1,2,1);
%^"T z,f vjJ!d#8 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
@Qx|!% (FMYR8H*( legend('r_p','r_s','|r_p|','|r_s|')
N_L&!%s tGzYO/Zp xlabel('\theta_i')
TE*$NxQ 2 Ta!.oC[
ylabel('Amplitude')
AEi@t0By 'N5qX>Ob title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
q4 'x'8 q yy.3-( axis([0 90 -1 1])
~u8}s4 n<CJx+U grid on
5d(A( O6OP{sb subplot(1,2,2);
%>i7A?L !j4C:L3F plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
;FZ\PxN Sct-,K%i legend('t_p','t_s','|t_p|','|t_s|')
5DyN=[b ^V_acAuS^ xlabel('\theta_i')
j1YE_U HcHfwLin0 ylabel('Amplitude')
2O9dU 5b X@9_ukdpu title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
yixW>W} -5vc0"?E axis([0 90 0 1])
5vbnO]8 K;6K!6J:[ grid on
"MC&!AMv 'xStA Rp=abs(rp).^2;
u{H,i(mx? e^&YQl Rs=abs(rs).^2;
RO wbzA)]r lbg6n:@ Rn=(Rp+Rs)/2;
B[4y(Im bq[Q Tp=1-Rp;
{, APZ`q| D{'Na5( Ts=1-Rs;
dhK$XG s^V8FH Tn=(Tp+Ts)/2;
K!I]/0L ^#3$C?d figure(2)
4N$svA tx5_e[ subplot(1,2,1);
).oqlA! A$XjzTR plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
~g|e?$j U"m!f*a legend('R_p','R_s','R_n')
jcq(=7j `t!iknOQ$ xlabel('\theta_i')
Mh+'f 93 #Z$6>
Xt ylabel('Amplitude')
@z/]!n\~ qZ\zsOnp title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/=e[(5X|O RX2=
iO" axis([0 90 0 1])
'hya#rC&( {f^30Fw grid on
[PX'Jer g+k6pi* subplot(1,2,2);
XIjSwR kYJ mH{cGu? plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
(
?/0$DB huKz["]z[ legend('T_p','T_s','T_n')
Plq[Ml9
=r-Wy.a@ xlabel('\theta_i')
mu{%%b7|^ JyB>,t) ylabel('Amplitude')
lZ)u4_ 2{V| title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c(Ha"tBJ l?FNYvL axis([0 90 0 1])
--^D)n b$$XriD] grid on
\~?s= LT wj'iU&aca
zfo.S[R@ \clWrK
OT *W]f