利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
I ]a [Ngj QB3er]y0% 1、光疏射向光密
-"rANP-UI Irui{%T clear
.uSVZqJ7 _Kbj?j close all
OXnTD!m>{ N}nE?|N=5 n1=1,n2=1.45;
'< $*N T: SqENV theta=0:0.1:90;
$WXO1o(O \k{d'R#~( a=theta*pi/180;
6O_l;A[=1 s|I$c;> rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
VTwQD"oB |
{Q}:_/q rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Gk
xtGe g C(@]% tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
i*r ag0Mw \-SC-c ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
]=PkgOJD 6F5g2hBz figure(1)
nk;^sq4M : ;iW>i8 subplot(1,2,1);
9N<=,!;5~s Xck`"RU<xA plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
51tZ:-1 ! NFF!g]QN legend('r_p','r_s','|r_p|','|r_s|')
^7a@?|,q8 Ww"]3 xlabel('\theta_i')
uPxJwWXO 'uF75C ylabel('Amplitude')
SLRF\mh!L C80< L5\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
P EMBh?)g {$#88Qa\- axis([0 90 -1 1])
Wt=\hixj- PE g]z grid on
j+rY }1 %% ` subplot(1,2,2);
e ^ ,IZ{ `sDLxgwI plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
=dsEt\
j yZ N~A: legend('t_p','t_s','|t_p|','|t_s|')
e )N< r 4j8$& ~/ xlabel('\theta_i')
ANgt\8 PT
}J.Dw x ylabel('Amplitude')
MkhD*\D
/ xkU8( = title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
M6bM`wHH> ~K V{m axis([0 90 0 1])
JPHM+3v Nukyvse grid on
^tjw }sE _ ]- 4UA- Rp=abs(rp).^2;
C"$ ~w3A k vCNq2l^CW Rs=abs(rs).^2;
I ~^Xw7 x cn~KF8 Rn=(Rp+Rs)/2;
>rJ**y EeT 69o Tp=1-Rp;
"LJV} L IsSh Ai Ts=1-Rs;
%H;}+U]Z 5423Ky< Tn=(Tp+Ts)/2;
l~w^I|M^C JW Ye~ figure(2)
:U *8S\$ 6yO5{. _M subplot(1,2,1);
#p7gg61 R\XS5HOE( plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
0kN Kt( _ NbTaI{r legend('R_p','R_s','R_n')
d7mn(= & KeRC8mYp xlabel('\theta_i')
j7QK8O$XL $(+xhn(O ylabel('Amplitude')
/zb/ am1# Kh[l} ;/F title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2~l7WW+lx, $ v[mIR axis([0 90 0 1])
Shn=Q +3o0GJ
grid on
_p5#`-%mM x,]x>Up subplot(1,2,2);
9<&M~(dwT4 9(O eH7 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
'S9o!hb'@ E?czol Nl legend('T_p','T_s','T_n')
bmK !02y'JS1 xlabel('\theta_i')
P+ CdqOL c*. ylabel('Amplitude')
U._fb= dNNXMQ0" title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Du65> O s]O Z+^Z axis([0 90 0 1])
FU/:'/ L s~2o<# grid on
z-b78A/8 h@\-] zN{ 
J]XLWAM ?bt;i>O\ 
H#D:'B j29
