利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
^i{B8]2, @o+T<}kW X 1、光疏射向光密
S,5 >g07-` lLHHuQpuj clear
DytOS}/^9 3zb;q@JV close all
%PG::b ZkYc9!an Y 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));
N DB ]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*G 4 { 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)
@Q x|!% (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 y y.3-( axis([0 90 -1 1])
~ u8}s4 n<CJx+U grid on
5 d(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 HcHfw Lin0 ylabel('Amplitude')
2O9dU 5b X@9_ukdpu title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
yixW>W} -5v c0"?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 l bg6n:@ Rn=(Rp+Rs)/2;
B[4y ( Im bq[Q Tp=1-Rp;
{, APZ`q| D{'Na 5( Ts=1-Rs;
dhK$ XG s ^V8FH Tn=(Tp+Ts)/2;
K!I]/0L ^# 3$C?d figure(2)
4N$s vA tx5_e [ subplot(1,2,1);
).oqlA ! A$XjzT R 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 R X2=
iO" axis([0 90 0 1])
'hya#r C&( {f^30F w grid on
[PX'Jer g+k 6pi* 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)u 4_ 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 
zf o.S[R@ \clWrK 
OT *W]f
