利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
H,W8JNPs T+x
/ J]A 1、光疏射向光密
x Q';$& CDF;cM"td clear
eIy:5/ s $vK,Gugcx close all
+1I 7 K|M }Ml BmD n1=1,n2=1.45;
B uso
`G V FM!K$_ theta=0:0.1:90;
DE7y\oO] (f^/KB= a=theta*pi/180;
}z#M!~ !Pz#czo rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
}
xA@3RT n&o"RE 0~0 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
C~ A`h=A< N
p*T[J tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
H)h^|A/vO BW6Ox=sr< ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Pr>05lg St(jrZb figure(1)
p^}`^>OL i#^Y QCy subplot(1,2,1);
6&5D4
V {]U
\HE1w plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
~ES%=if~Y %I9f_5BlT8 legend('r_p','r_s','|r_p|','|r_s|')
~(^pGL3< q.<)0nk xlabel('\theta_i')
u,nn\>Y
qo u\4YZ ylabel('Amplitude')
r73W. & Qd\ ='*:! title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
$=,pQ q +} 1zw< axis([0 90 -1 1])
/F9Dg<#a 1=5HQ~|[TO grid on
Qe7=6< oemN$g &7 subplot(1,2,2);
*!TQC6b$ y$bY
8L plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Q"U%]2@= fVgN8b|&' legend('t_p','t_s','|t_p|','|t_s|')
]cv|d c= F-b]>3r xlabel('\theta_i')
nS h~ mP 9_d# F'#F ylabel('Amplitude')
f8SO:ihXL ]" e'z title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:!D m,PP% L C##em=Y axis([0 90 0 1])
T
iL.py, ZA. S X|m grid on
[ `1` E1X ab2Cn|F Rp=abs(rp).^2;
! [1a P, *k; bkd4x Rs=abs(rs).^2;
P7z Uf [<{r~YFjWW Rn=(Rp+Rs)/2;
@[?ZwzY:9 vf@j d}? Tp=1-Rp;
!W8=\:D[ kr~n5WiAZ Ts=1-Rs;
;,6C&|n]w E9>z. vV
Tn=(Tp+Ts)/2;
u-y?i` ~E((n figure(2)
n"aF#HR?0d X<. l(9$ subplot(1,2,1);
3u[8;1}7Q ny qX\m - plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
$#+D :W)az J^CAQfcx legend('R_p','R_s','R_n')
RCYv 2=m>Q 82ixv<B xlabel('\theta_i')
9 Xl#$d5 QICxSk ylabel('Amplitude')
j;E$7QH[ #9r}Kr =P title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Yb`b / BMR z9OpMA axis([0 90 0 1])
F (."nUrf z8'zH> grid on
4G' E< ab =H23eOS_# subplot(1,2,2);
_eq$C=3Ta w0Nm.=I- plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
u6?9#L( `:~Wu/Ogr legend('T_p','T_s','T_n')
PKntz7 rjWtioZEa xlabel('\theta_i')
4v _Hh<% XD8 I.q ylabel('Amplitude')
)\:IRr" 2jC:uk title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
@ :PMb Ub DpA)Vdj axis([0 90 0 1])
]i6* $qgma !*~QB4\2b grid on
elhP!"G M_UmnqN1C 
:Ch XzZ &}Wi@;G]2 
+x`p WH]2
