利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
Uc4L|: l]klV+9t 1、光疏射向光密
TLL[F;uZ 9snyX7/!L clear
J%O4IcE LN3dp?;_{ close all
NV:XPw/ o YI=p3l n1=1,n2=1.45;
s*~jvL 9bcyPN theta=0:0.1:90;
Hb AMoow! M_.,c Vk a=theta*pi/180;
DneSzqO"o b=QGbFf rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
2}#wdJ` KutgW#+40 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
3_eml\CY []N$;~R7 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`@.s!L(V *ZSp9g"Z ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
C=_-p"O# sr($Bw figure(1)
w KMk|y> vyx\N{ subplot(1,2,1);
53+rpU_ ]E8<;t)# plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
qpwh #^2 acju!,G legend('r_p','r_s','|r_p|','|r_s|')
IweK!,:>dN Q3O .<9S xlabel('\theta_i')
SnF[mN' 674oL, ylabel('Amplitude')
ch8w' 0|> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Dx# @D# 3;l>x/amk axis([0 90 -1 1])
_}9R} 6ewOZ,"j"4 grid on
ZM^;%( ^FMa8;'o subplot(1,2,2);
n
ua8y(W 56
kgL;$h plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
e%c5OZ3~ ~$ qJw?r
legend('t_p','t_s','|t_p|','|t_s|')
N[bf.5T -r'seb5 xlabel('\theta_i')
f]A6Mx6 XM@i|AK
M0 ylabel('Amplitude')
2GKU9cV*` #bZ=R title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
&b6@_C9 xv% USm axis([0 90 0 1])
u`Qcw|R+ wfTv<WG,.E grid on
'=5_u tP(bRQ> Rp=abs(rp).^2;
Fx:en|g p}zk&` Rs=abs(rs).^2;
-Fc# "XsY~ Rn=(Rp+Rs)/2;
%+B-Z/1} 6')SJ*|yS Tp=1-Rp;
O*/-I
pM z==}~|5 Ts=1-Rs;
cB F%])! Wk6&TrWlY Tn=(Tp+Ts)/2;
x&/Syb +Y]*>afG figure(2)
V;]VwsZ" e27CbA{_w subplot(1,2,1);
uvv-lAbjw C^=gZ
6m plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
skx=w<YO6] -XIvj'u legend('R_p','R_s','R_n')
~)IJE+e>} rkc%S5we xlabel('\theta_i')
~I|R}hS uWtS83i ylabel('Amplitude')
V+~{a:8[pq e.ym7L]$O title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
KP-z Bo\v-97 axis([0 90 0 1])
U105u.#7 [Q_|6Di grid on
EjE`S_i= 5f@YrTO[@ subplot(1,2,2);
4m!3P"$ H08YMP>dc plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
;p!hd}C 1gej$G@ legend('T_p','T_s','T_n')
>t2)Z|1 mf[79:90^ xlabel('\theta_i')
/_\W*@ E uOqDJM'RM ylabel('Amplitude')
j=% -b] C\@YH] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}M@Jrq+7 aMg f6veM axis([0 90 0 1])
Mf13@XEo !rGI), grid on
-mo
'
$1 rB|:r\Z(jG TFNU+ 8vo7~6yy d6
EJn/