利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
o8Gygi5 15,JD 1、光疏射向光密
:aK?Dt Z 8! rdqI clear
!
5NuFLOf ZZ7qSyB s? close all
d&[RfZ` <^'{=A>
n1=1,n2=1.45;
m1U:&{:^ dH!z<~ theta=0:0.1:90;
* <\K-NSL fEX=csZ86 a=theta*pi/180;
l6y}>] qh:Bc$S rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Aeb(b+= -cM1]soT rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
p,goY F?? MDU#V tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
&CQO+Yr$l 0Gc@ AG{ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
-}9^$}PR N,c!1: b figure(1)
DK\XC%~m ?9kC[4G subplot(1,2,1);
3o %vV* {d'-1 z"q plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
N+=|WeZ ,|{`(y/v
legend('r_p','r_s','|r_p|','|r_s|')
M'PZ{6; R7kkth xlabel('\theta_i')
&ASR2J motK}G ylabel('Amplitude')
[mA-sl] eF,F <IJT{ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
f9W:-00QD #GE]]7:Na axis([0 90 -1 1])
<qwf"Ey - 2M~KlYl grid on
=~GP;=6 z# ^fS
| subplot(1,2,2);
.9 r YBy +e-F`k plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
c Mgd ]>t~Bcn m legend('t_p','t_s','|t_p|','|t_s|')
u]P| 9{*{Ba xlabel('\theta_i')
#;]#NqFX U!aM63F3 ylabel('Amplitude')
D%p* G5Bg3 r D].=.?1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uiVN z8H ) *~A|[ axis([0 90 0 1])
h Ma; \ k 9 {&g.+ grid on
)l7XZ_gw' qq/ > E*~ Rp=abs(rp).^2;
QB*,+ u4 !6 KX^j- Rs=abs(rs).^2;
6ZqU:^3 {^WK#$] Rn=(Rp+Rs)/2;
c ZY y+ &-3 e 3) Tp=1-Rp;
Xp:A;i9 b=,B Le\ Ts=1-Rs;
#i bwD:{ BNfj0e 5b Tn=(Tp+Ts)/2;
m,k 0 h% T/_u;My; figure(2)
ppyy0E^M 42NfD/"g+s subplot(1,2,1);
}QFL [';o -c"! plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
'J,UKK\5 L4>14D\ legend('R_p','R_s','R_n')
o,*m,Qc q Gk.7wf% xlabel('\theta_i')
FDMQ Lx f V<QpC5 ylabel('Amplitude')
: _8K8Sa &C9IR,& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
B\J[O5}, Kh]es,$D axis([0 90 0 1])
f%]@e9dD ZISIW! grid on
\ND]x]5d UPhO =G subplot(1,2,2);
~H}Z;n]H 9@*pC@ I) plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
f 2WVg;Z SZWNN#w60? legend('T_p','T_s','T_n')
!> sA.L&= M hMiSsZ xlabel('\theta_i')
N[_T3( Y 25`vE( ylabel('Amplitude')
2{ F-@}= E`)e
;^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,} t%7I REh"/d axis([0 90 0 1])
* ~P B /TMV Pnvz. grid on
xA3_W +~5Lo' ^ 
V IzIl\<aM T<uX[BO-a 
+6WjOcu
