利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
B$cx
'_zF ij~023$DTt 1、光疏射向光密
'HDbU#vD dLI`\e<r&[ clear
eA(c{ gA gP(" close all
S[W|=(f9 5UHxB"`C n1=1,n2=1.45;
Nm]\0m0p-
;qT7BUh(% theta=0:0.1:90;
jNA ^
(|: E-q*u(IW a=theta*pi/180;
="*8ja-K ^zr]#`@G rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
7`f',ZK% sw715"L rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
>O&:[CgEF kGSB6 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
lj!f\C}d h&6v&%S/L ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
v/yt C/WH" "m:4e`_dz figure(1)
JH0L^p &% \`Lwh subplot(1,2,1);
'
Z}/3 d p
"cUCB plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
\kGi5G] T=QV =21qn legend('r_p','r_s','|r_p|','|r_s|')
cSSrMYX2 ,,>b=r_r& xlabel('\theta_i')
" '/$ZpY ^#4?v^QNh ylabel('Amplitude')
-v(.]`Wo&; l;dZJ_Ut$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
NN9` jP2 R2af >R axis([0 90 -1 1])
fs cAG\>8 @*gm\sU4 grid on
a9GLFA8Vq Z) zWfv} subplot(1,2,2);
iPxhDn<B bstc|8< plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
5Go&+|c vJ lll]FJ1 legend('t_p','t_s','|t_p|','|t_s|')
0/."R ; XKU+'Tz xlabel('\theta_i')
#D$vH ji[O ? ylabel('Amplitude')
x^| J- *kXSl73 k title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
d4#Q<!r lr[a~ca\ axis([0 90 0 1])
pCQB<6&1N cia 4!-# grid on
E Cx_
[|3{ i~dW)7 Rp=abs(rp).^2;
l0 H,TT~2 L$GhM!c Rs=abs(rs).^2;
$ GyO+xF T7AFL= Rn=(Rp+Rs)/2;
e4V4 %Qw bqXCe\# Tp=1-Rp;
V6,D ~7 nz|;6?LCLY Ts=1-Rs;
BHE((3 ,R1`/aRy Tn=(Tp+Ts)/2;
u":D{+wC | U-I LzK figure(2)
FKd5]am C^S?W=1=w subplot(1,2,1);
u,
%mVd ~a[]4\ m; plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
pCz;km 0KyujU?sF legend('R_p','R_s','R_n')
kSH3)CC P \#[W8k<Z xlabel('\theta_i')
`~hAXnQK= }ZEh^zdz8 ylabel('Amplitude')
fv !l { G q I^$5? title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"-y\F}TE ]],6Fi+
axis([0 90 0 1])
Wiqy".Y Y (XJehdB0 grid on
ei;wT 2*Uwp; 0 subplot(1,2,2);
. *>L D F_! 6C-z plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
X-! yi J p.Sow legend('T_p','T_s','T_n')
kx 'ncxN~ 4:8#&eF xlabel('\theta_i')
J.:"yK""
/I
ylabel('Amplitude')
y U Q;t TI 4@|K^nT` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h:(Jes2 lph3"a^ axis([0 90 0 1])
!%NxSJ EA2BN} grid on
.s| n}{D_i 6g!#"=ls; 
!o*BRR* 8Uvf9 ,I' 
4x=rew>Ew
