利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
BQ gK< _ |t58n{V.O 1、光疏射向光密
)h%tEY$AJ ~f 8:sDJ clear
d P/1E6*m $C{-gx+: close all
ge&!G O hus9Zv4 n1=1,n2=1.45;
\-Q6z 8 +=L^h9F theta=0:0.1:90;
5.U|CL j;&su=p" a=theta*pi/180;
?ieC>cr ~5~Cpu2v7 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
eC$ Jdf `Y< FR rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
c38EN f gCI{g. [I! tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
vSJ#
}& =?1B|hdo ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
CGZ3-OW@E 1n>(CwLG" figure(1)
7MwS[N%# eh6= - subplot(1,2,1);
/V E|F Ts /,^AG2]( f plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
a`; nB E bb`8YF+?' legend('r_p','r_s','|r_p|','|r_s|')
Au{J/G<W@ \2y [Hy? xlabel('\theta_i')
SH1)@K- ;0:[X+"( ylabel('Amplitude')
9 EV. ![ %v<BE
tq title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)s> R~7 m Bg$eiGTB axis([0 90 -1 1])
}6 Mo C0 Z3z"c
B grid on
lWk/vj<5 $3^Cp_p6 subplot(1,2,2);
)PjU=@$lI ^_G@a, plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
# 9 p|aS\ T3h 1eU legend('t_p','t_s','|t_p|','|t_s|')
Z6IWQo,)Rh K5}0!_)G xlabel('\theta_i')
ug S L=V.@? ylabel('Amplitude')
P$*9Z@ jAy 0k
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
-nDY3$U/ EM,= R axis([0 90 0 1])
$Ahe Vps@@ gvoYyO#cm grid on
|Ju d*z c]R27r E Rp=abs(rp).^2;
U}DLzn|w %`} nP3 Rs=abs(rs).^2;
zQaD&2 q W}@IUCRs Rn=(Rp+Rs)/2;
j?1wP6/NP . :~E.b Tp=1-Rp;
h`f $]_c }mpFo 2 Ts=1-Rs;
lESv zJo?,c Tn=(Tp+Ts)/2;
6o4Y]C2W{1 63/a 0Yn figure(2)
_fS4a134R O~m Q\GlW subplot(1,2,1);
z! :0%qu &cwN& XBY plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
ehTv@2b UjS+ Ddp legend('R_p','R_s','R_n')
c2-oFLNP= Yz#E0aTTA xlabel('\theta_i')
D7=Irz!O\7 D>0(*O ylabel('Amplitude')
.q1y)l-^Z -F_c Bu81V title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
nDyA][ Iw?*y.z| axis([0 90 0 1])
*Qyu
QF nw(R=C grid on
YQb503W"d~ Xk\IO0GF subplot(1,2,2);
.5);W;`X r6GXmr plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
uIBN
!\j >@Pw{Zh$ legend('T_p','T_s','T_n')
c-_1tSh}
i[I&m]N xlabel('\theta_i')
Z_fwvcZ?05 T?RN} @D ylabel('Amplitude')
A@9\Qd | P`b"x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
$B7c\MR
j Xu
E' %;: axis([0 90 0 1])
l fF RqZ ?XP4kjJ grid on
w]nt_xj J5n6K$ .d 
<$K=3&:s8q U k*(C( 
'J3yJ{
