利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
"'p+qbT8 ~ dI&> CL 1、光疏射向光密
odT7Gq \kJt@ [w% clear
Z?j='/u>@ D 5wR?O close all
m^u&g&^ {p_vR/yN n1=1,n2=1.45;
QIz N#;g Gy Xs{* theta=0:0.1:90;
JU,ROoz( ';T5[l, a=theta*pi/180;
e.Q K% qYiK bzy rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`h{mj|~ u|*|RuY rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
$j2)_(<A%Q H5MAN,` tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
=
(h;L$ 6)2M/( ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
nX aX= bFxJ| figure(1)
{gHscj;SM zZ0V6T} subplot(1,2,1);
)0V]G{QN 4bjp*1 *] plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
yzzre>F Uht:wEr legend('r_p','r_s','|r_p|','|r_s|')
$P9$ ,w4 ~g6[ [ xlabel('\theta_i')
V>j` :q[n1
O[Ch ylabel('Amplitude')
8VGXw;(Y,d NGZ>: title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(G!J== bjX$idL axis([0 90 -1 1])
_rh.z_a7w t|"d#5' grid on
=C2C~Xd T>o# *{qn subplot(1,2,2);
t[ocp;Q 9f[[%80 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
9+|,aG s Q*{
2 legend('t_p','t_s','|t_p|','|t_s|')
`g1~ya(MC f^)nZ:~ xlabel('\theta_i')
Ci:QIsu* <xNM@!'\h ylabel('Amplitude')
NistW+{< pD/S\E0@t title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
->{WO+6( t!C-G+It axis([0 90 0 1])
bcs(# yln.E vJjD grid on
8R D)yRJ DZZt%n8J Rp=abs(rp).^2;
:UciFIa ri`|qy6! | Rs=abs(rs).^2;
>f/g:[ (k45k/PAP Rn=(Rp+Rs)/2;
-lEh}r $sJfxh
r Tp=1-Rp;
/Qy0vAvJ + k:?;ZG Ts=1-Rs;
X2Mj|_#u .ji%%f Tn=(Tp+Ts)/2;
`6su_8Hno =G6@:h= figure(2)
/+YWp>6LU IuD<lMeJJ subplot(1,2,1);
*ni|I@8 39#>C~BOl plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
eH6cBX#P. {wySH[V legend('R_p','R_s','R_n')
, 88}5)b[ ipG+qj/= xlabel('\theta_i')
g>f394j /FY_LM ylabel('Amplitude')
wC~Uy% M*%Z5,Tc title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
pm+[,u!i Y!zlte|P axis([0 90 0 1])
Zlhr0itf }J1#UH_E grid on
>b0}X)Z+U &CgD smJo# subplot(1,2,2);
cqDnZ`|6 #%~wuCn<K plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
K9c:K/H }tN"C 3)@ legend('T_p','T_s','T_n')
&LE/hA TD'1L:mv xlabel('\theta_i')
C'Q} Z_ }/=_ ylabel('Amplitude')
HzV3O-Qz] U<$ |ET' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
gwyz)CUkL pnWDsC~) axis([0 90 0 1])
"b qB@) NeYj[Q~xy grid on
)9s[-W,e kjH0u$n
+]G;_/[2 :EZTJu
5bk5EE`