利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
';,Rq9-' K-X@3&X} 1、光疏射向光密
S~GL_#a %0&c0vT clear
^(79SOZC
7j,u&%om close all
MeCHn2zwB fC|NK+Xd` n1=1,n2=1.45;
4C]>{osv >n(Ga9E theta=0:0.1:90;
&[#iM0;)W0 Z~[E ZgIg a=theta*pi/180;
R%EpF'[~[ E'KKR1t rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
OUM^ u* ^ s.necg0 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
}[FP"# DWXxB tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<4Ev3z*;Z t?l0L1; ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
P[L] S7FTr ) y#~eYn figure(1)
zLt7jxx |K( j XZ) subplot(1,2,1);
XSBh+ )0Ww #^i.[7p plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
oZS.pi e09('SON( legend('r_p','r_s','|r_p|','|r_s|')
_DD.#YB</ A(2!.Y
2?* xlabel('\theta_i')
S W; %2 ca6kqh" ylabel('Amplitude')
4). i4]%LH (U?*Z/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
V!&O5T (~ Or:a\qQ1 axis([0 90 -1 1])
h+d \u I7C*P~32{n grid on
DMXm$PU4V 8`]1Nt!*B subplot(1,2,2);
L kK# =v 2N /4. plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
n`TXm g ) 1PjI9M legend('t_p','t_s','|t_p|','|t_s|')
,GVD.whUl n 97pxD_74 xlabel('\theta_i')
%4#Q 3YlyD rEfo)jod ylabel('Amplitude')
oU[>.Igi ZIr&_x#e title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
9V/:1I0?&0 /9 wmc2 axis([0 90 0 1])
ba`V`0p- ( @b] ?Gg grid on
}<7S% ?TY dd>
qy Rp=abs(rp).^2;
BXj] ]S2 6jaol'{SuH Rs=abs(rs).^2;
lKT<aYX u |eV'-R)s Rn=(Rp+Rs)/2;
G9qN1q~ yKb+bm&5:' Tp=1-Rp;
HQ0fY ,e93I6 Ts=1-Rs;
~u^MRe|` a 9H^e<g Tn=(Tp+Ts)/2;
(fUpj^E)p =F 9!)r figure(2)
!M*$p Qi} sngM4ikhs subplot(1,2,1);
^f?>;,<& E|~)"= plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
D.;iz>_}Y nP 5T*-~ legend('R_p','R_s','R_n')
Jb1L[sT2 Ng 3r`S"_< xlabel('\theta_i')
^Ee"w7XjD duT'$}2@> ylabel('Amplitude')
tX'2 $} Mv%"aFC title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+ _"AF| j8[RDiJ axis([0 90 0 1])
+>*! 3x+sE $AyE6j_1gX grid on
*kM^l!<g /U P1*L subplot(1,2,2);
kH|cB!?x u,&[I^WK`C plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
(_.0g}2 +eb mve \+ legend('T_p','T_s','T_n')
QRFBMq}' Mq,2S xlabel('\theta_i')
JK=0juv<E 2M$^|j:[ ylabel('Amplitude')
Gw@]w;ed tmVGJ+gz title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
b2vCr F; K?$|Y-_D^M axis([0 90 0 1])
:rzq[J^ WT_4YM\bz grid on
UVz}"TRq. Q*K31Ln 
N O'37d S |pf .l 
=b_/_b$q
