利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
whxE[Xnv ~`
#t?1SP 1、光疏射向光密
m#DC;(Pn /sC$ ;l clear
|67<h5Q1 }1 Q]C"hY close all
I O|">a6 @P[Tu; 4 n1=1,n2=1.45;
0xZq?9a ND CZc_ theta=0:0.1:90;
H W D
%6dFACv a=theta*pi/180;
%w_h8 [C/h{WPC- rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
|.)dOk,o :] {+ 3A rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Cc/?-0a2! @dD70T tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
_T8#36iR IP/
zFbc ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
A)s $w#C;2k]N figure(1)
vvFXdHP +'JM:};1X8 subplot(1,2,1);
>*B/Wy @kRe0:t plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
4Poi:0oOys 97NF*-)N legend('r_p','r_s','|r_p|','|r_s|')
oq(um:m &
/UcFB xlabel('\theta_i')
cAV9.VS<L _
B",? } ylabel('Amplitude')
F+v? 2|03 |:L<Ko title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
FsTE.PT +#A >[,U axis([0 90 -1 1])
MjQKcL4%7 p4bQCI grid on
srw5&s(3X 0^?(;AK subplot(1,2,2);
218ZUg -a !! %vs
6 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
JyvXNV, @;"HslU\Q legend('t_p','t_s','|t_p|','|t_s|')
7-nwfp&|$ & 2& K9R xlabel('\theta_i')
{~.h;'m oIAP dn ylabel('Amplitude')
?tQU ZO JLp.bxx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[Z` q7ddd^ W+ Z]
Y axis([0 90 0 1])
h?p&9[e` vA10'Gx ' grid on
C6Ap
4 \=?f4*4|/ Rp=abs(rp).^2;
]HgAI$aA, vClD)Ar Rs=abs(rs).^2;
uZ&,tH/ 6!
\a8q'z Rn=(Rp+Rs)/2;
Jk!}z+X'A aF1pq Tp=1-Rp;
Ji\8(7
{8 NmthvKhH Ts=1-Rs;
FZ
DC? w&8gA[y*u Tn=(Tp+Ts)/2;
, 9^wKS!7$ U< |kA(5 figure(2)
@@9#od O B9`nV.a subplot(1,2,1);
w8:~LX.n n a])bBn plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
&EYO[~D06 \/I@&$"F legend('R_p','R_s','R_n')
n qg=I ,kfUlv= xlabel('\theta_i')
M]oO1GM +]*4!4MK6 ylabel('Amplitude')
/!rH DcR |T uk9d4] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"yxIaTZu Zmw'.hL axis([0 90 0 1])
?u{y[pI6 6":=p:PT. grid on
]fg?)z-Z wyeiz7 subplot(1,2,2);
}A;J-7g6 8Ih+^Y
a plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
M^DYzJ 1nhC! jDD legend('T_p','T_s','T_n')
NitWIj[U; #Y`GWT1== xlabel('\theta_i')
QL%&b\K ,E2Tw-% ylabel('Amplitude')
4Z5ZV! #?-2f{ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{{w5F2b((% o Mk6ZzZ,> axis([0 90 0 1])
$x# 0m 0+p
<Jc! grid on
:^]rjy/|+ ]E$h7 I 
Ucr$5^ME ?_*X\En*3 
&g%9$*gm T
