利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
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 gC(@]% 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- PEg]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 yZN~A: legend('t_p','t_s','|t_p|','|t_s|')
e)N<r 4j8$&~/ xlabel('\theta_i')
ANgt\8 PT
}J.Dwx ylabel('Amplitude')
MkhD*\D
/ xkU8(= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
M6bM`wHH> ~KV{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 xcn~KF8 Rn=(Rp+Rs)/2;
>rJ**y EeT69o Tp=1-Rp;
"LJV}L IsShAi Ts=1-Rs;
%H;}+U]Z 5423Ky< Tn=(Tp+Ts)/2;
l~w^I|M^C JWYe~ figure(2)
:U *8S\$ 6yO5{._M subplot(1,2,1);
#p7gg61 R\XS5HOE( plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
0kNKt(_ 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(OeH7 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
'S9o!hb'@ E?czolNl 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]OZ+^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