利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
[8oX[oP aMJ2bu 1、光疏射向光密
q]SH'Wd ltNY8xrdGN clear
:()K2<E |)*! &\Ch close all
kV!1k<f C#3&,G W n1=1,n2=1.45;
#MiO4zXgd [
<k&]Kv theta=0:0.1:90;
IH5^M74b B+pLW/4l a=theta*pi/180;
kqq1;Kd 1|:;~9n<t rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
u^WZsW ,6"[vb#* 3 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
If8
^ 6KPj ZC< tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
+OtD@lD`! }h8U.k?v ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
T
eu.i L6ypn)l figure(1)
>enP~uW[# Kq+vAp). subplot(1,2,1);
\nL@P6X 4\uq$.f- plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
QXEZ?gx AgUjC legend('r_p','r_s','|r_p|','|r_s|')
nB5^ =2vZqGO30 xlabel('\theta_i')
niYD[Ra\xP !FB2\hiM ylabel('Amplitude')
!d[]Qt%mA / sPa$D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{>TAnb?n _) 4zm axis([0 90 -1 1])
%>y!N!.F Hd
gABIuX grid on
f 3:dn7 Xd1+?2 subplot(1,2,2);
fWDTP|DV ft iAty0n plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
?iHcY, sAZL,w legend('t_p','t_s','|t_p|','|t_s|')
5Z#(C# vB5mOXGN q xlabel('\theta_i')
rm|,+ { AU9:Gu@M/ ylabel('Amplitude')
%L-{4Z!"sI uOZSX.o^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
* )HVK& ' ma"M? aM axis([0 90 0 1])
FR]uCH v^[!NygShs grid on
&[yYgfsp <\d2)I v Rp=abs(rp).^2;
>km$zfM2- <YCR^?hJSi Rs=abs(rs).^2;
V warU(* G,(Xz"`, Rn=(Rp+Rs)/2;
<N=ow "rD eq9qE^[Z& Tp=1-Rp;
'}F9f? M[_Ptqjb Ts=1-Rs;
xq%BR[1 D? H|O[ Tn=(Tp+Ts)/2;
{WeRFiQ?- = 2} bQW figure(2)
9R&.$5[W(s M,Lq4 bz subplot(1,2,1);
Q
H>g-@ FE1En plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
'p%w_VbI 7CX5pRNL legend('R_p','R_s','R_n')
-}8r1jQH; -YipPo "a xlabel('\theta_i')
f
uH3C~u7< c5[ ~2e ylabel('Amplitude')
?9*[\m?- \z.p [;'ir title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
5yroi@KT 5iGz*_
m axis([0 90 0 1])
oiz]Bd ItAC=/(d grid on
V9`jq $ /5\{(= 0 subplot(1,2,2);
Q"o JhxS h{?f
uoZj% plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
S7R*R} }F3}"Ik'L legend('T_p','T_s','T_n')
F-Ku0z]){? H LE%f; xlabel('\theta_i')
owO &[D/ iX>)6)uJ ylabel('Amplitude')
obgO-d9l LM!@LQAMY title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
j?! /#' a]I~.$G
axis([0 90 0 1])
/j\.~=,_ ?.nD!S @ grid on
n+Ng7 tZr_ {F @ 
eHyIFoaC/ 1L3 $h0i 
',6d0>4 *
