利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
B>GE 9y5 jBU4F~ 1y 1、光疏射向光密
Z\HX~*,6 T}"6wywM clear
6(G?MW. :{ oZ ~< close all
'UMXq~RMe Ti'}MC+0 n1=1,n2=1.45;
9nH?l{As OtGb<v<_H theta=0:0.1:90;
|\U5) ,m z}{afEb a=theta*pi/180;
(iR ide !OT-b>*w rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
BM(]QUxRd D5:{fWVsV/ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
NcP.;u;` (4\d]*u5-c tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
KyzFnVH3) om]4BRe ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
h;-yU.(w j S]><rm figure(1)
ak:ibV _z54Ycr4H subplot(1,2,1);
1c;6xc,ub XXh6 ^@H= plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
;s~ X v$=QA:!U legend('r_p','r_s','|r_p|','|r_s|')
PI*82,f3dE MNg^]tpf xlabel('\theta_i')
#Jna6 J7v|vj I ylabel('Amplitude')
hRq3C1 mR 1 ]}#)- title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
z.CywME<)t t w(J ZDc axis([0 90 -1 1])
B}eA\O4}I fd #j Y} grid on
3mSXWl^? 8kbY+W%n subplot(1,2,2);
Iu{kP yx &bhq`> plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
E( h<$w8s 0Hs|*:Y1D legend('t_p','t_s','|t_p|','|t_s|')
Lgz$]Jbl8 ZbT/$\0(6 xlabel('\theta_i')
Y $u9%0q|? .@"q$\ ylabel('Amplitude')
=E&OuX-R 0A:n0[V:] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
1.D-FPK 2Sq_Tw3^ axis([0 90 0 1])
KIeT!kmDl xBFJ} v grid on
m'-|{c !m(4F(!"h Rp=abs(rp).^2;
Jd,)a#<j +m}Pmi$ Rs=abs(rs).^2;
eh* 6cQ.0 L'i-fM[# Rn=(Rp+Rs)/2;
81%8{yn!$" &t(0E:^TRU Tp=1-Rp;
hSGb-$~F Sb".] >^ Ts=1-Rs;
8wH.et25k ,&HR(jTo Tn=(Tp+Ts)/2;
'<iK*[NW $ XBAZ<"hd figure(2)
X~ Rl 6/, ]=73-ywn] subplot(1,2,1);
#(LfYw.P1V z5Qs @dG plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
'rb'7=z5 = \ ,
qP legend('R_p','R_s','R_n')
8 iC:xcN3 I /x iT xlabel('\theta_i')
V+wH?H= }iBC@`mg( ylabel('Amplitude')
?lCd{14Mkh 9O,,m~B title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
d6 _C"r ^b?2N/m@ axis([0 90 0 1])
|f 2A89 ~H@+D}J? grid on
$<?X7n^ H ;@!?I subplot(1,2,2);
^ D Aa%u @fn 6<3 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
\X}8 q `QP
~ legend('T_p','T_s','T_n')
Vj^dD9: $. e) xlabel('\theta_i')
d }0qJoH4 jPu5nwvUV> ylabel('Amplitude')
q7]>i! A 89eq[ |G_ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/r$&]C:Fi 6517Km 4- axis([0 90 0 1])
M?:\9DDd whkJ pK(
grid on
s^C;> iK}v`xq 
2*[U n( ]if;A ) ' 
@&R1wr1>I5
