利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
kr?| >6 ? *}J_STM 1、光疏射向光密
;1%a:#5 Z K vh] clear
q
mB@kbt hD*?\bBs0 close all
vEsSqzc [=Z{y8#:J n1=1,n2=1.45;
N9Ml&*%oX{ RBXoU'. theta=0:0.1:90;
]Oy<zU 5f5ZfK3<i a=theta*pi/180;
eBB
D9 SI 7d)aDc*TjW rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
z -D pLV |#wz)=mD rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
f"Iyo:Wt ?<;<# JN tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`9-Zg??8r wOOPWwk ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
b ~gF,^w `Nn?G figure(1)
kQ8 WO|bA 5IK@<#wE subplot(1,2,1);
3fPv71NVtt zMKL: Um" plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
~'3% Qr `2PT 8UM legend('r_p','r_s','|r_p|','|r_s|')
Q)>'fZ) Vz*'^=(o& xlabel('\theta_i')
5*$Zfuf >y[S?M ylabel('Amplitude')
GN0'-z6Uy v knFtpx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,whNh V*}zwm s6 axis([0 90 -1 1])
7 %"7Rb^@ BP$#a
# grid on
R{_ IrYk K }BX6dA subplot(1,2,2);
([~`{,sv t^hkGYj!2 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
&u-Bu;G.e "=.|QKC1` legend('t_p','t_s','|t_p|','|t_s|')
Oe`t! &v G.8b\E~ xlabel('\theta_i')
$P3nP=mf [2V/v ylabel('Amplitude')
]||=<!^kn Hea<!z PH title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"[yiNJ"kt T*yveo &j axis([0 90 0 1])
Kf#!IY][ gN~y6c:N grid on
dL(|Y{4 k qw? X{ Rp=abs(rp).^2;
]?{lQ0vw'w TzGm562o% Rs=abs(rs).^2;
HU9y{H 6l'y Rn=(Rp+Rs)/2;
U I C? S 8
-A7 Tp=1-Rp;
$:!T/ *p* "6Uj:9 Ts=1-Rs;
-QDgr`%5 B_glyC Tn=(Tp+Ts)/2;
(B<AK4G Ir #V2]$ figure(2)
#Ca's'j&f !`o:+Gg@ subplot(1,2,1);
{LHe 6# X>7]g670@ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
[t{ #@X _9|@nUD legend('R_p','R_s','R_n')
bK9~C" k MXk. 2 xlabel('\theta_i')
J5_Y\@ /
L M ylabel('Amplitude')
f*{;\n (.t 9t.yP;j\Y title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
=@=R)C4f* 3Lki7QW` axis([0 90 0 1])
jDXGm[U rq["O/2 grid on
2Q|*xd4B^ ^jjJM| a subplot(1,2,2);
D*'M ^k|1 x9A
ZS#e)[ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
O>M*mT M }, &, Dt legend('T_p','T_s','T_n')
'8kjTf#g<l \R qh|T<D xlabel('\theta_i')
|#:d C # .y9rM{h}b ylabel('Amplitude')
umZ
g}|C_ /3 Ix,7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Cmx2 /N -u9yR"n\} axis([0 90 0 1])
tO>OD# !riMIl1 grid on
' oZn<c` K6*UFO4}i 
y{Y+2}Dv/ 2vdQ&H4 
Jo aDX ,
