利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
! &NrbiuN 9C0 #K\ 1、光疏射向光密
+.OdrvN4) $L?KNXHAF! clear
b6p'%;Y/ m^=El7+ close all
'4Fwh]Ee lk[Y6yE n1=1,n2=1.45;
R<(x WH h72CGA| theta=0:0.1:90;
Z*Gf`d: C, GZ a=theta*pi/180;
Vxu V`P lf P.QF9% rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
-6~ .;M 5 NzTF2ve( rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Ip:54 V; CPn tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
C/'w )*S:C ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Am_>x8z u6 Lx3 figure(1)
)%3T1
D/ 7 )r L<+ subplot(1,2,1);
4^M"V5tDx q b/}&J7+ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
H-U _ eZN"t~\rX legend('r_p','r_s','|r_p|','|r_s|')
7GWOJ^) 7(N+'8 xlabel('\theta_i')
5j6`W?|q PP >6 ylabel('Amplitude')
j49Uj}:j d7
H *F title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
R&J?X Q :dAd5v2f axis([0 90 -1 1])
xLDD;Qm, Y)+q[MZ R grid on
9 fYNSr 7t:tS7{} subplot(1,2,2);
$2?j2}M |K06 H
?6X plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
+a 1iZ bh y!.jpF'uI legend('t_p','t_s','|t_p|','|t_s|')
[ nP s [!-gb+L xlabel('\theta_i')
1 {
%y (?` P
<+0sh ylabel('Amplitude')
HOF$(86zqA eY5mwJ0K title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'.(Gg%*\. E<[
s+i X axis([0 90 0 1])
]b!n ;{5 .'gm2 grid on
HdJ g U5OX.0 Rp=abs(rp).^2;
'#Au~5 ?YLq
iAA Rs=abs(rs).^2;
(^m]
7l yJq< &g Rn=(Rp+Rs)/2;
yXJ25Axb 49$<:{ ~ Tp=1-Rp;
VNHce H 7|DG1p9C Ts=1-Rs;
Y \-W` ehr-o7]( Tn=(Tp+Ts)/2;
Wye* ~t EIg~^xK figure(2)
5mD8$% \8 L(VFzPkY% subplot(1,2,1);
^H3N1eC,`F ` 3qf}=Z` plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
alaL/p{O K)7T]z` legend('R_p','R_s','R_n')
ZH_$Q$9 ,,gMUpL7_8 xlabel('\theta_i')
G?1GkR L7-BuW}& ylabel('Amplitude')
Mw/9DrE7/ IR6 W'vA title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c\{N:S> &v{#yzM axis([0 90 0 1])
$^ee~v;m4 j Q5 F} grid on
u
]e-IYH s5nw<V9$] subplot(1,2,2);
H9/!oI1P? ^ ` y7JXI: plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
b &:v6#i [a2]_]E% legend('T_p','T_s','T_n')
pCs3-&rI3 5H0qMt P xlabel('\theta_i')
9\_AB .Z: "GO!^ZG] ylabel('Amplitude')
G%
tlV&In {EoYU\x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/iU<\+ H MA `nFkVK axis([0 90 0 1])
.\R9tt} -tj#BEC[H( grid on
|nefg0`rk i 1vz{Tc 
I zpE|8l )QE6X67i 
xE:jcA
d$}
