Mie散射Matlab源代码

Mie散射的Matlab源代码，适合研究Mie散射的研究员参考使用。

function result = Mie（m x）

% Computation of Mie Efficiencies for given
% complex refractive-index ratio m=m‘+im“
% and size parameter x=k0*a where k0= wave number in ambient
% medium a=sphere radius using complex Mie Coefficients
% an and bn for n=1 to nmax
% s. Bohren and Huffman （1983） BEWI:TDD122 p. 103119-122477.
% Result: m‘ m“ x efficiencies for extinction （qext）
% scattering （qsca） absorption （qabs） backscattering （qb）
% asymmetry parameter （asy=） and （qratio=qb/qsca）.
% Uses the function “Mie_abcd“ for an and bn for n=1 to nmax.
% C. M妕zler May 2002.

if x==0 % To avoid a singularity at x=0
result=[real（m） imag（m） 0 0 0 0 0 0 1.5];
elseif x>0 % This is the normal situation
nmax=round（2+x+4*x^（1/3））;
n1=nmax-1;
n=（1:nmax）;cn=2*n+1; c1n=n.*（n+2）./（n+1）; c2n=cn./n./（n+1）;
x2=x*x;
f=mie_abcd（mx）;
anp=（real（f（1:）））; anpp=（imag（f（1:）））;
bnp=（real（f（2:）））; bnpp=（imag（f（2:）））;
g1（1:4nmax）=[0; 0; 0; 0]; % displaced numbers used for
g1（11:n1）=anp（2:nmax）; % asymmetry parameter p. 120
g1（21:n1）=anpp（2:nmax）;
g1（31:n1）=bnp（2:nmax）;
g1（41:n1）=bnpp（2:nmax）;
dn=cn.*（anp+bnp）;
q=sum（dn）;
qext=2*q/x2;
en=cn.*（anp.*anp+anpp.*anpp+bnp.*bnp+bnpp.*bnpp）;
q=sum（en）;
qsca=2*q/x2;
qabs=qext-qsca;
fn=（f（1:）-f（2:））.*cn;
gn=（-1）.^n;
f（3:）=fn.*gn;
q=sum（f（3:））;
qb=q*q‘/x2;
asy1=c1n.*（anp.*g1（1:）+anpp.*g1（2:）+bnp.*g1（3:）+bnpp.*g1（4:））;
asy2=c2n.*（anp.*bnp+anpp.*bnpp）;
asy=4/x2*sum（asy1+asy2）/qsca;
qratio=qb/qsca;
result=[real（m） imag（m） x qext qsca qabs qb asy qratio];
end;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1570  2007-02-25 16:01  Mie\Mie.m

     文件       1398  2007-02-25 16:03  Mie\MieExample.m

     文件        925  2007-02-25 16:02  Mie\Mie_abcd.m

     目录          0  2007-02-25 16:02  Mie

----------- ---------  ---------- -----  ----

                 4111                    5