利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
'T*h0xX 2hE+Om^n 1、光疏射向光密
sH\ h{^ #4ii!ev clear
){xMM Q5 H?)?(t7@ close all
Mx=L lC ) }BLT2]y0 n1=1,n2=1.45;
Oy~X@A lBfG#\rdW~ theta=0:0.1:90;
=pn(56 N t_7Z a=theta*pi/180;
J_>nn o Np4> 7Lk rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Fn*)!,) znu [i&\= rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
q%c"`u/v / za [;d4<}k tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
D8wZC'7 BxHfL8$1[$ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Wup%.yT~Ds aXyg`CDv figure(1)
:qO)^~x I=o/1:[- subplot(1,2,1);
hwdZP=X j! hdi-aTU plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
c`4i#R R*z:+p}oHy legend('r_p','r_s','|r_p|','|r_s|')
jGKas I` ]6TX ) 1
xlabel('\theta_i')
6sl2vHzA \_PD@A9 ylabel('Amplitude')
e?eX9yA7F _U/C G<n title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Y lEV @ :gacP? axis([0 90 -1 1])
7P7d[KP< g'{hp: grid on
{B+{2;Zk zHW}A
`Rz subplot(1,2,2);
:} 9Lb)Yp YztW1GvI plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
PR*qyELu SoJ'y 6 legend('t_p','t_s','|t_p|','|t_s|')
LwS>jNJx !L\P.FP7b xlabel('\theta_i')
&?j]L4% 5W~-|8 m ylabel('Amplitude')
coFQu ; i =}Xw}X+[WY title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ejI nJ o5 |P5h axis([0 90 0 1])
3QO*1P@q i>n)T grid on
r-k,4Yz g$$j:U*- Rp=abs(rp).^2;
!xo; $4 2%zJI"Ic Rs=abs(rs).^2;
Ve\=By-a| Z=\wI:TY1 Rn=(Rp+Rs)/2;
sFWH*k dP? v^QUYsar Tp=1-Rp;
Zf ub+A 1>e%(k2w% Ts=1-Rs;
%44Z7 #Mmr{4m Tn=(Tp+Ts)/2;
NA9 N#; )+xHv figure(2)
q'~F6$kv5 <764|q subplot(1,2,1);
h|S6LgB F R9*WI
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
'}eA2Q>BV Q( \2(x\ legend('R_p','R_s','R_n')
.B|a.-oA4 a}#Jcy!e xlabel('\theta_i')
"T,^>xD |37y =" ylabel('Amplitude')
I3$/ # %~A$cc title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
q"sD>Yh & eLc@w<yB axis([0 90 0 1])
K$(LiP -3 2?]LN}
grid on
z3X:.% _onEXrM subplot(1,2,2);
mfpL?N (fJ. o-LQ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
F;@A2WD :VEy\ R>W legend('T_p','T_s','T_n')
= $t L=9w
3VXS xlabel('\theta_i')
2%F!aeX wX!>&Gc. ylabel('Amplitude')
FaUc"J HD& Ag title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
21\t2<" AECaX4h+_ axis([0 90 0 1])
auKGm: * f4BD|| grid on
g9$ P J: f/Q7WXl0
Jm(sx'qPx c3NUJ~>=y 
4$=ATa;x-
