利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
+R:(_:7 [0("Q;Ec[j 1、光疏射向光密
P@~yx#G X1_ 5KH clear
vtJJ#8a]
%^6F_F_jS close all
AFt s( ,|/f`Pl n1=1,n2=1.45;
m+$VVn3Z}
K 3l95he theta=0:0.1:90;
W\V.r$? v *|HY>U. a=theta*pi/180;
n~Lt\K: J] r^W)O rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
7F.4Ga; ;722\y(Y rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
ZS o) 7 _[L o4_ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
f*
wx< >/6 _ ^ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^xk'Z |JsZJ9W+J figure(1)
;< 4a*;IO %B?=q@!QWn subplot(1,2,1);
RT8 ?7xFc ,<X9 Y2B plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
1k^oS$UT AQ^u legend('r_p','r_s','|r_p|','|r_s|')
;5AcFB q ZZK#,Qb xlabel('\theta_i')
cm +P]8o%{ j 5 h-dK ylabel('Amplitude')
B `sAk
% 62NsJ<#> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ue>D 7\8 :rP=t , axis([0 90 -1 1])
\GU<43J2uo UC$ppTCc? grid on
$<O D31T o{[q Zc_% subplot(1,2,2);
D)}v@je"yP ^=*;X;7 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
!p/goqT~dY -t U'yKhn legend('t_p','t_s','|t_p|','|t_s|')
9j Gu}V o o#N+Y?O xlabel('\theta_i')
dQG=G% W qxJ\ye+'* ylabel('Amplitude')
rS Ni@; >.D4co> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
G KeU%x {\\T gs axis([0 90 0 1])
- !
S_ryL ^kSqsT" grid on
K.yb
^dg5 Kf-JcBsrT Rp=abs(rp).^2;
$Z> 'Jp 3r1*m
+ Rs=abs(rs).^2;
VS|2|n1<6 ,]/X\t5]D Rn=(Rp+Rs)/2;
/Gfw8g\} :MDKC /mC Tp=1-Rp;
'O-"\J\ M'l ;: Ts=1-Rs;
nT)vNWT= o<!?7g{ Tn=(Tp+Ts)/2;
.o}v#W+st @ [v~y"tE} figure(2)
U`s{Jm >5SSQ\ 2~a subplot(1,2,1);
k|f4Cf, t ZB<on<.) plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
9MqGIOQ${j 8a"% 0d# legend('R_p','R_s','R_n')
S`]k>'
l [a<SDMR xlabel('\theta_i')
-D~%|).' ]J]h#ZHx ylabel('Amplitude')
v(%*b,^
l9H!au= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+qdEq _ m Uo ix axis([0 90 0 1])
'/%H3A#L z([</D? grid on
FXU8[j0P_G pI<f) r subplot(1,2,2);
m~ee/& T 7 Fs ay+a plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
dUdT7 ixo hK|Ul]q I legend('T_p','T_s','T_n')
6D_D' ;o @`Su0W+. xlabel('\theta_i')
+8ZF"{y w@fi{H (R ylabel('Amplitude')
Fv`,3aNB c Q_Hp
<D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Y u2Bkq+ "-V"=t' axis([0 90 0 1])
Nmh*EAJSy ]')RMg zM* grid on
~4cC/"q$X R0-j5&^jju 
+X]vl=0 K\c#ig 
C$) onk
