利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
) N"gW* *Cgd?*\ 7 1、光疏射向光密
#[#dc]D eL!G, W clear
:O}<Q L:F:ZOM6` close all
p^``hP:J NeyGIEP n1=1,n2=1.45;
WH$
Ls(' B1Iq:5nmoS theta=0:0.1:90;
wa y-Q7 [}_ar a=theta*pi/180;
uNY]%[AnJ 8)k.lPoo. rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<!|2Ru
'[#uf/~W rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
rbP3 & L _n{6/ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
/qL&)24 (Qq$ql 27 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
EH$wW l^ Dx[t?- figure(1)
ZdH WSfO)O Ck.GN<#-^P subplot(1,2,1);
6CU8BDN R GtUKr' plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
^J~}KOH O!]w J legend('r_p','r_s','|r_p|','|r_s|')
W^+b gg<. vFCp= 8h xlabel('\theta_i')
+-~;?wA } }( ~' ylabel('Amplitude')
Wu Gm~<NS $C)@GGY title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
KQ9:lJKr Cm ;N5i axis([0 90 -1 1])
{2:H`|x {Ue6DK % grid on
l d@f:Zali dadMwe_l0 subplot(1,2,2);
ci^-0l_O *j*jA/ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
&1': s|c 9wb$_j]F`# legend('t_p','t_s','|t_p|','|t_s|')
Ixec]UOS #uXOyiE xlabel('\theta_i')
D@JHi'F !7 ^He3 ylabel('Amplitude')
#/jHnRrQ wFd*6% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
d,)F #;^5 g*tLqV axis([0 90 0 1])
<zDe ;& ca- n:1 grid on
b{dzbmak Sl_zO?/PF Rp=abs(rp).^2;
z
vY DE] $c y:G Rs=abs(rs).^2;
R{rV1j#@!a lb 6s3b Rn=(Rp+Rs)/2;
G[ ^G~U\+! m X1oRhf Tp=1-Rp;
2b"5/$|6 !KlSw,&=.6 Ts=1-Rs;
=J'?>-B yR$_$N+E Tn=(Tp+Ts)/2;
&sNID4FR f: Ja figure(2)
1*jL2P]D 7cy~qg subplot(1,2,1);
z0UO<Y?9 .5yc O plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
;VH]TKkk k:qou})#4 legend('R_p','R_s','R_n')
mT>p:G 2~U+PyeNz xlabel('\theta_i')
1G;Ns] u 8{l=`y"nB ylabel('Amplitude')
qR<DQTO< yk)j;i4@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
H44&u](8{ M0e|G.S&_ axis([0 90 0 1])
w Ebs E<</ Iz8 ^? >X grid on
16Y ~5JAc C4uR5U subplot(1,2,2);
DmrfD28j~F -]A#G`' plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
S}f<@-16P y
m{/0&7 legend('T_p','T_s','T_n')
A(p D ` X6'PP xlabel('\theta_i')
Rx"Qwi, \U ]."c4S_)| ylabel('Amplitude')
9Gws Q \ Uh Tr<(@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
@/(7kh + j1C.# -P[ axis([0 90 0 1])
9tEKA|8 Ve1] ECk grid on
|Z^g\l.j{ G0> Wk#or 
mO>
M =2A CEq]B:[IC 
tuU XW5!/
