利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
v7-'H/d. B@4#y9`5 1、光疏射向光密
AYi$LsLhO >@X=E3 clear
}MX`WW0\]Z o.Q9kk? L close all
1ysA ~2 -'nx7wnj2 n1=1,n2=1.45;
}LRAe3N%8 *,. {Xf theta=0:0.1:90;
+SJ aE] $ C" s a.#} a=theta*pi/180;
)sQbDA|p l20q(lb rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
M]W4S4&Y= m4iR
'~L} rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
P$D r6; F,D & tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
u,C-U!A G$Dg*< ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*c d9[ ~ `8Om*{xg figure(1)
Ua:EI!` cr!I"kTgD subplot(1,2,1);
}X$>84s>[P D&'".N,} plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
^b%AwzHH} +,%x&L&I legend('r_p','r_s','|r_p|','|r_s|')
LP87X-qkjW >?kt3.IQ!X xlabel('\theta_i')
l`1ZS8 [. h m"B kOA ylabel('Amplitude')
u]}Xq{ZN kTzZj|l^\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
l"- D@]" Ro1l:P)C` axis([0 90 -1 1])
CN >q`[! T}zi P grid on
t512]eqhb( ( R0 subplot(1,2,2);
7(oxmv}#Q eBP
N[V plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
$% jV%k ubgq8@; legend('t_p','t_s','|t_p|','|t_s|')
Ji7A9Hk )~Q$ tM` xlabel('\theta_i')
ZpwB"%e$ k'{lo _ ylabel('Amplitude')
~"LOw_BRh Wl*\kQ}U title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
IX3r$}4 yXkgGY5 axis([0 90 0 1])
t%/5$<!b m.2 grid on
Pukq{/27 ]2(
%^#qBG Rp=abs(rp).^2;
pI h@!C W}EO]A%f.\ Rs=abs(rs).^2;
"=Ziy4V ]-0
&[@I4@ Rn=(Rp+Rs)/2;
q'W`t>2T ]CTu | Tp=1-Rp;
_)p% &"d4J?io` Ts=1-Rs;
k1]?d7g$w t4,6`d?C Tn=(Tp+Ts)/2;
`T }e3l aZmSCi:&' figure(2)
c|I{U[(U T,;6q!s= subplot(1,2,1);
ptUnV3h (0Y6tcV]R plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
(=7e~'DC Ha;^U/0| legend('R_p','R_s','R_n')
l_$~~z ~ .lqo>Ta
y xlabel('\theta_i')
`m3C\\ 9; ~uV.jh ylabel('Amplitude')
6n <:ph,h; i 1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
YsLEbue [<bfwTFsl axis([0 90 0 1])
tV%M2 D xS (3 B;
V grid on
%1xo|6hm- e PlEd'Z subplot(1,2,2);
`y
m^0x8 sexnO^s plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
#j'O rD p bAQf3 legend('T_p','T_s','T_n')
w0VJt <e* f |aO9w xlabel('\theta_i')
+ d]} fok OjTE ylabel('Amplitude')
v?5Xx{ym <(;"L<?D<C title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
y~''r%] miEfxim axis([0 90 0 1])
:,6dW?mun6 Yw4c`M yL grid on
U*1rA/"n @UE0. R< 
0 }-&v+ Hk(w\
9kWyO:a_(
