利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
1(VskFtZF XOZ@ek)LY 1、光疏射向光密
ta SYR$V J w 3L+7V, ! clear
/jU4mPb;\D f*[Uq0? close all
ehX4[j6 ( >om.FM n1=1,n2=1.45;
%( 1y {RH)&k&% theta=0:0.1:90;
PiX( Ase M[Jy?b) a=theta*pi/180;
`]2y =f<{X
({t6Cbw rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`b5pa `\4 q[.,i{2R} rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
e5 sQl1 o PA m* tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
VC NQ}h[D 74~ %4 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
,Ct1)%
`[4{]jX+< figure(1)
T!H }^v s9?H#^Y5u subplot(1,2,1);
sCAWrbOe> ?CuwA-j plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
z`y^o*qc] R ?kyJ4S legend('r_p','r_s','|r_p|','|r_s|')
]*AQT 7PH v}"DW? xlabel('\theta_i')
TP)}1 @ /. @"wAw: ylabel('Amplitude')
{X&H wb-_CQ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0a's[ >-'A nA#dXckoc axis([0 90 -1 1])
.u 9,w ncij)7c)u grid on
)L7h:%h# ~@VyJT% subplot(1,2,2);
$)M 5@KT yUFT9bD plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
D3;#: kCU (Hi`Q legend('t_p','t_s','|t_p|','|t_s|')
$+[
v17lF s\
YHT.O? xlabel('\theta_i')
iXuSFman vHx[:vuq: ylabel('Amplitude')
b(: U]>J kt hy9<!$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
-Y/c]g V3> JZH` axis([0 90 0 1])
Wr\A ->+ .d%CD`8! grid on
i~EFRI@ X?OH//co Rp=abs(rp).^2;
Z {*<G x r/mKuGa] Rs=abs(rs).^2;
|]x>|Z?/u xU;;@9X Rn=(Rp+Rs)/2;
IkJ-*vI6 {3 *Zx"e![ Tp=1-Rp;
D1 f}g B)u*c]<qU Ts=1-Rs;
nK9?|@S*' 5X,|Pn Tn=(Tp+Ts)/2;
Ex$i8fO( aqN6.t figure(2)
WBY_%RTx EHX/XM subplot(1,2,1);
YS+|n%? ,ft KRq plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
5? 1:RE(1 tsN,yI]-VA legend('R_p','R_s','R_n')
zP|^) h5 <K zEn+ xlabel('\theta_i')
|'N)HH>; q jml wVw ylabel('Amplitude')
)V?:qCuY> :[.**,0R title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
$ir F p5\B0G<m axis([0 90 0 1])
M)j.Uu `0Bk@B[> grid on
yJ?S7+b \*5${[ subplot(1,2,2);
6h:2,h
pE nv/'C=+L plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
fKOC-% w }GL@?kAGR5 legend('T_p','T_s','T_n')
M.?[Xpa VQwF9Iq]` xlabel('\theta_i')
VH7nyqEM 6>&h9@ ylabel('Amplitude')
fm1yZX?` 6g &E v' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+ Un(VTD 3
G_0DS axis([0 90 0 1])
aGq1 YOD[$ r6gfxW5 grid on
GY wU3`{ ?Es(pwJB 
xY>@GSO1 u4vyj#V 
FE,BvNBZ
