mie散射系数的新算法

《mie散射系数的新算法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、Mie散射系数的新算法王少清任中京　张希明　何芳江海鹰(山东建筑材料工业学院，济南　250022)提要：介绍了一种计算Mie散射系数的新方法，给出了计算实例。关键词：米氏散射，Mie系数，Mie计算NewalgorithmofMiescatteringparametersWangShaoqingRenZhongjing　ZhangXiming　HeFangJiangHaiying(ShandongInstituteofBuildingMaterials，Jinan　250022)Abstract：AnewalgorithmofcalculatingMieParameterswasintro

2、duced.Somecalculationresultsdonewiththisalgorithmweregiven.Keywords:Miescattering，Mieparameters，Miecalculation1　引言Mie理论是球形颗粒对单色光的散射场分布的严格解析解[1]，目前在环保、动力、气象、天文、两相流及粉体颗粒尺寸分布测量等工程技术领域中有广泛的应用。利用单一颗粒或颗粒群光散射场的测量数据。可以反推得散射颗粒或颗粒群的很多物理性质，如颗粒的尺寸、颗粒的折射率等[2]。但反推必须事先计算出各种尺寸的颗粒在各种复折射率下的散射场分布数据。1968年Dave[3]最先发表了

3、完整的Mie散射计算方法，以后Lentz[4]和Wiscomb[5]又针对部分计算提出了新算法。国内也有人发表了他们自己的算法[6、7、8]。但总的看来，这些算法均有各自的局限性。尤其是当颗粒尺寸或折射率的虚部值较大时，往往计算速度过慢或产生溢出和不收敛等现象。本文介绍作者发展的Mie散射新算法。该算法的特点是不受颗粒尺寸及折射率的限制，不会产生溢出和不收敛的现象，且具有较快的计算速度。2　Mie散射系数的计算公式Mie散射计算的中心问题是计算Mie散射系数an和bn，其表达式为[9]其中为颗粒的尺寸参数，定义为α=πd/λ，d为颗粒直径，λ为入射光在颗粒周围介质中的波长，而m为颗粒在周围

4、介质中的相对复折射率，即m=m1+im2(m2<0),式中i为虚数单位。而Ψn(Z)和ξn(Z)(Z既表示α，又表示mα)的表达式为Ψn(Z)=（πZ/2）Jn+1(Z)ξn(Z)=Ψn(Z)+iΧn(Z)Χn(Z)=̶(πZ/2)(̶1)n-1J-(n-1)(Z)=̶(πZ/2)Nn+1(Z)3计算中产生溢出的原因计算Mie散射系数须先计算Ψn和Χn。一般采用递推的方法。递推又分为向前递推(即从n=0开始)与向后递推(即从n=N开始至n=0,N为预先设定值)。实验表明，向前递推总是快于向后递推。Ψn和Χn的初值为分析以上两式可知，当m2≠0时,若颗粒尺寸d很大,或复折射率的虚部值m2很大

5、，将使得乘积m2d很大，可使两式中的项exp(-m2α)=exp(-πm2d/λ)的值超过计算机的数据限，从而产生溢出，这是产生溢出的重要原因。另外，在递推过程中，不恰当的算法也可能造成溢出。4　an和bn的新算法为解决上述问题，作者提出了新的算法。将an和bn的公式变形如下:令其中Lnr、Lnj分别表示Ln(m)的实部与虚部。将（2）式代入（1）式，并用anr、anj和bnr、bnj分别表示的实部与虚部。如此可推得在上述四个公式中采用比值的形式是非常重要的，这样可避免递推过程中当ai与bi较大时乘法运算可能产生的溢出。这是本算法的一个重要特点。在以上四式中由于均为实变量函数，计算不产生溢

6、出。关键是和的算法如何处理，才能保证计算中不产生溢出。在Lentz的算法中是采用连分式计算Ln的值，其精度的保证是由在大量计算基础上得出一个截断项数N与参数a与m的经验公式而实现的。这样的经验公式，一则有实用上的局限性，再则也会带来截断误差。文献（6）对此经验公式做了改进，但仍陷于a=1~100，m1=1~2，m2=0~1的范围。下面介绍本文作者发展的关于Ln的新算法。该算法的特点是不受a及m值的限制，不会产生溢出或不收敛等病态现象，且具有较快的计算速度。令上面导出的(3)-(20)诸式，构成了Mie系数an和bn的完整算法。由于an和bn是从n=1开始计算，利用初值公式(16)-(20)

7、即可算得任意级数的an和bn的值，故没有舍入误差的问题。从(16)式可见，因为y=m2ɑ≤0,故无论m2和ɑ取何值均不会产生溢出，再加上(3)的各式中采用了比的形式，又避免了计算过程中的溢出，这就从根本上解决了溢出的问题。5　计算实例利用以上算法编制了计算颗粒散射场强度和消光系数的计算机程序。当波长为λ的单位振幅的平面自然光入射颗粒时，颗粒的散射光强为[9]Ξ为计算机在双精度下的最小数据限。图1示出散射光强的一组计算实例