利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
|> mx*G 4N&4TUIM 1、光疏射向光密
Dk$[b9b NbPv>/r clear
tz{W69k+ tE"aNA#= close all
@"[xX}xK; )@"iWQ 3K n1=1,n2=1.45;
(<RZZ{m ,1-n=eTQ theta=0:0.1:90;
s?h=%; T[ zqqu7.` a=theta*pi/180;
o' U:: [gK (x% rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
c#l W ? +k=BD s rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
h *J=F0KM zNe>fZ
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
%au2kG, 6|q \ M ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
u(W%snl Oy>u/g~ figure(1)
g8<Ja (J 0OJBC~?{\ subplot(1,2,1);
q}vz]L&o ` Mj>t( plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
? OrRTRW :Osw4u]JXd legend('r_p','r_s','|r_p|','|r_s|')
`?Wy;5- nt$ V H xlabel('\theta_i')
7G N>o@ t .L;M-`^ ylabel('Amplitude')
i"eUacBz/- MXy~kb& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
y7[D9ZvZ <=m@Sg{ o axis([0 90 -1 1])
mj\]oWS7 d j6zZ! k grid on
M<' AM 4 {FV_APL 9_ subplot(1,2,2);
{.vU; r`? bYoz plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
kNW}0CDgs 3yT7;~vPj legend('t_p','t_s','|t_p|','|t_s|')
\[m{ &%^G ,{{e'S9cy xlabel('\theta_i')
Yvky =RM u} y)'eH ylabel('Amplitude')
eBw6k09C+ R $vo title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
TL%2?'G ?!PpooYK axis([0 90 0 1])
<B,z)c #
tN#_<W grid on
`/WX!4eR, $w,&h:.p Rp=abs(rp).^2;
d9' gH#f? u;{,,ct Rs=abs(rs).^2;
AQCU\E v; !f Rn=(Rp+Rs)/2;
^zdZ"\x u 19 d!#g Tp=1-Rp;
Q& p'\6~ 7WuhY Jbf
Ts=1-Rs;
Aqi9@BH i+`N0!8lY Tn=(Tp+Ts)/2;
y _M te xW`,@a } figure(2)
-Xm/sq(i)% &<Gs@UX~w subplot(1,2,1);
8Ja 't8 6rBXC <Z plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
|&O7F;/_ 3`V #ImV> legend('R_p','R_s','R_n')
:i;iSrKy q-3,p . xlabel('\theta_i')
^Q)&lxlxpx ^
+e5 M1U= ylabel('Amplitude')
e;Iz K]kP 2m$\]\kCUv title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
zUw=e}?: Vn4y^_H axis([0 90 0 1])
]{mz %\ s>i`=[qFc grid on
Uc j
eB D_n(T ') subplot(1,2,2);
]`p*ZTr)\ /hx|KC&:e plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
]yvHb)X ,!m][ legend('T_p','T_s','T_n')
*3,Kn}ik p3sR>ToJ xlabel('\theta_i')
_]g?3Gw7! *G]zN "Y ylabel('Amplitude')
;AL keUR[ $-tgd<2h title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
STfcx] L d nZA+Pa axis([0 90 0 1])
_A[k&nO!&J q4'Vb grid on
Sq t '} rKK{*%n 
4P1<Zi+< ;rD
M%S@ 
W7k0!Grrl
