利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
]saf<?fzr L{$ZL & 1、光疏射向光密
<a *X&P />fy@nPl| clear
B\&;eZY'G +r)'?zU close all
N7 j -A}*Aa'\ n1=1,n2=1.45;
~o`I[-g) }fw;{&s{z theta=0:0.1:90;
nf!RB-orF <Gzy*1 Q& a=theta*pi/180;
+~J?/ .)^@[yrkz rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
jh.W$.Oq iK%<0m rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
:u#Ls,OZz Uh}n'Xd#{} tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
KU[eY} oo7&.H Wf ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
!]Z> T5$ a#W:SgE?Y figure(1)
DsY-JBDvoz D"fE )@Q@Y subplot(1,2,1);
ann!" s_ ) F 6#n&2 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
v=?U{{xQ j.4oYxK!s/ legend('r_p','r_s','|r_p|','|r_s|')
SV .\B DG1
>T xlabel('\theta_i')
j[YzBXd
V pm Xx2T#= ylabel('Amplitude')
ay#cW., $DC*&hqpt title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Byey Uw F.?`<7 axis([0 90 -1 1])
{e[S?1t=l >L/Rf8j & grid on
u4B, |_MK d
-6[\S# subplot(1,2,2);
phSF. WC !_s | h@ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
&|P@$O > |E!()j= legend('t_p','t_s','|t_p|','|t_s|')
O
z%K* 8%JxXtWW` xlabel('\theta_i')
G5Y5_r6Gu 8WV1OIL ylabel('Amplitude')
Q^/5hA hu\HK81m title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(r`+q[ 49}yw3- axis([0 90 0 1])
OlhfBu)~ A|tee@H*0 grid on
` =dD6r 0CI?[R\ Rp=abs(rp).^2;
VE^NSk Oa& ?*E Y~'I Rs=abs(rs).^2;
exV6&bdu i0,{*LD%^ Rn=(Rp+Rs)/2;
?UQVmE& W^0F(9~!( Tp=1-Rp;
8r Xq-V_u &'Qz Ts=1-Rs;
Y+< C[Fiq uOc>~ITPS Tn=(Tp+Ts)/2;
]4\^> ( gO ?-0 figure(2)
;6t>!2I>C oT&JQ,i[2Q subplot(1,2,1);
];6c/#2x k);z}`7 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
8+
eZU<\B( _)2.#L legend('R_p','R_s','R_n')
7Upm m\7-/e2 a xlabel('\theta_i')
M$-4.+G u#%Ig 3 ylabel('Amplitude')
xPb;_~ m{!BSl title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
|K'{R'A UA{sUj+? axis([0 90 0 1])
M ! gX4 nFW^^v< grid on
u_@%}zo?5* w>; :mf subplot(1,2,2);
^j '; 4' \l{*1lQ` plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
NmST1pMk n )} J< legend('T_p','T_s','T_n')
x `PIJE dJ?VN!B 0 xlabel('\theta_i')
_@;2h`q ? I"Ji_4QV ylabel('Amplitude')
Q'JE DH\ v Q+}rHf`[ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
s6Bt)8 A ,.]e~O4R axis([0 90 0 1])
sxgR;gf6 seHwn'Jn grid on
KU-z;}9s cAyR)Y!I 
@<CJbFgJp jCW>=1:JGY 
H24g+<Tv
