利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
Z)(C7,Xu [XQoag;! 1、光疏射向光密
\9BIRY ` lBTgI"n=eK clear
@ B3@M ~c3CyOab close all
o0Hh&:6!M P>kS$U) n1=1,n2=1.45;
#,qikKjt2 O ox5${#^ theta=0:0.1:90;
d=wzN3 ;- *pv hkJ g( a=theta*pi/180;
\Jv6Igu +B '9!t4 2 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
.x1EdfHed/ $0K9OF9$ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
VxfFk4 c{#lKD<7 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Xf9VW}`*8 $X-,6* ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
G#CWl ),= W?/7PVGv5h figure(1)
]]%CO$`T [ nS'0i&<{1 subplot(1,2,1);
;mo}$^49* #,
vN plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
88)0Xi|]KP T cSj `- legend('r_p','r_s','|r_p|','|r_s|')
Jb {m "mT~_BsD xlabel('\theta_i')
y)#Ib*? 9^s
sT>&/ ylabel('Amplitude')
v?h#Ym3e< fwxyZBr title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
oChcEx% hlkf|H axis([0 90 -1 1])
it>FG9hVo Jp5~iC2d grid on
]Wv\$JXI FQ(=Fnqn subplot(1,2,2);
Nc]oA Y 7,jh44(\= plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
(g ~&$&pa tp*AA@~ legend('t_p','t_s','|t_p|','|t_s|')
h+ELtf ;Cyt2]F xlabel('\theta_i')
jB?Tua$,s bEy j8=P; ylabel('Amplitude')
}[l`R{d5q> tM ]qR+ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
RpS 'Tz} <(@m913| axis([0 90 0 1])
M*+_E8Lh W/{HZ< :. grid on
cD2}EqZ 9 Y@^M U->+ Rp=abs(rp).^2;
S4`u NB#Ht Lf rS:g Rs=abs(rs).^2;
\\x``* Alp9]
0( Rn=(Rp+Rs)/2;
,uE WnZ"4 }3^t,>I=,6 Tp=1-Rp;
aLuxCobV a{69JY5 Ts=1-Rs;
CQzJ_aSJ ( hMeE@Q0 Tn=(Tp+Ts)/2;
H`aqpa"C )UZ
's>O figure(2)
ShHm7+fV
E@P8 -x'i subplot(1,2,1);
hq$:62NYg [ZOo%"M _Y plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
FrB}2 >K;p+( <6 legend('R_p','R_s','R_n')
!1)lGjMW y$#mk3(e~t xlabel('\theta_i')
oC.:mI SG ;]Vr ylabel('Amplitude')
{&" rv<p lJu2}XR iU title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~ %k<N/B {Qg"1 +hhM axis([0 90 0 1])
xhj
A!\DS wp-*S}TT grid on
z';p27 5 xv9SQ,n< subplot(1,2,2);
TR5"K{WDx BRFA%FZ, plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
G;qC& 7T rM< c;iQ legend('T_p','T_s','T_n')
92bvmP*o4 P[?~KNS:/ xlabel('\theta_i')
s==gjA e: ;:hyW,J ylabel('Amplitude')
#ZRQVC; b; X'IW &^kI title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0^5SL/2 DL_2%&k/ axis([0 90 0 1])
|u<qbl c*9Rz D#Zj grid on
;Q:^|Fw!F `<8~tS/. w 
<hkg~4EKc IFH%R>={ 
,/C<GFae
