基于rqpso的颗粒粒径分布反演算法-张彪.pdf

《基于rqpso的颗粒粒径分布反演算法-张彪.pdf》由会员分享，可在线阅读，更多相关《基于rqpso的颗粒粒径分布反演算法-张彪.pdf（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第 47 卷第 11 期 中南大学学报(自然科学版) Vol.47 No.112016 年 11月 Journal of Central South University (Science and Technology) Nov.2016DOI: 10.11817/j.issn.1672-7207.2016.11.040基于 RQPSO 的颗粒粒径分布反演算法张彪 1，李舒 2，许传龙 1，王式民 1(1.东南大学能源热转换及其过程测控教育部重点实验室，江苏南京，210096；2.南京市计量监督检测院科技发展部，江苏南京，210037)摘要： 针对标准量子微粒群算法 (BQPSO)，提出改进量

2、子微粒群算法 (IQPSO)和含正则化项的改进量子微粒群算法(RQPSO)，并将它们引入到粒径分布的反演中，利用光全散射法在独立模式下，通过测量可见光波段内不同波长下的光谱消光值反演几种粒径分布，其中正问题利用反常衍射近似 (ADA)计算得到估计值，测量值则通过 Mie理论计算得到。研究结果表明：与 BQPSO 相比，IQPSO 在计算效率和稳定性上得到很大提升；在粒径分布的反演中，RQPSO 提高 IQPSO的维数极限，并具有更高的反演精度、稳定性和抗噪性，为粒径分布的反演提供一种新的方法。关键词：粒径分布；量子微利群算法；正则化；光全散射法；独立模式中图分类号：TP212 文献标志码：A

3、文章编号：16727207(2016)11392207Retrieval of particle size distribution based onRQPSO algorithmZHANGBiao1,LIShu2,XUChuanlong1,WANGShimin1(1.KeyLaboratoryofEnergyThermalConversionandControlofMinistryofEducation,SoutheastUniversity,Nanjing210096,China;2.TechnologyDevelopmentDepartment,NanjingInstituteofMea

4、surementandTestingTechnology,Nanjing210037,China)Abstract: An improved quantum behavior particle swarm optimization (IQPSO) and regularized quantum behaviorparticle swarm optimization (RQPSO) were developed based on basic quantum behavior particle swarm optimization(BQPSO). Furthermore, these three

5、algorithms were introduced in retrieval of particle size distribution (PSD). Severaltypes of PSDs were retrieved by measuring the spectral extinction values in the visible spectrum, which used total lightscattering method under independent mode. In the direct problem, the anomalous diffraction appro

6、ximation was used tocalculate the estimation values, and Mie theory was used for measurement values. The results show that the efficiencyand stability of the IQPSO algorithm was proved to be more greatly improved than the BQPSO algorithm. For theretrieval of these PSDs, the RQPSO algorithm has bette

7、r performances on dimension limit, accuracy, stability and noiseimmunitythantheIQPSOalgorithm.Thus,thisalgorithmprovidesanewmethodforretrievingofPSDs.Key words: particlesizedistribution;quantumbehaviorparticleswarmoptimization;regularization;totallightscattering;independentmode收稿日期： 20160102； 修回日期：

8、20160320基金项目 (Foundation item)： 国家自然科学基金资助项目 (51506030,51376049)；江苏省自然科学基金资助项目 (BK20150622)；国家质量监督检验检 疫 总 局 科 技 计 划 项 目 (2012QK176) (Projects(51506030, 51376049) supported by the National Natural Science Foundation of China;Project(BK20150622) supported by the Natural Science Foundation of Jiangsu P

9、rovince; Project(2012QK176) supported by the StateAdministrationofQualitySupervision,InspectionandQuarantineScienceandTechnology)通信作者： 许传龙，博士，教授，博士生导师，从事多相流测试方面的研究； E-mail:万方数据第 11 期 张彪，等：基于 RQPSO 的颗粒粒径分布反演算法 3923颗粒粒径及其尺寸分布是颗粒技术最重要的参数和技术指标之一，在许多情况下，颗粒粒径不仅直接影响到产品的性能和质量，而且与能源的高效利用、环境污染的防治、生产工艺的优化、人身健康

10、的保障等等都密切相关 13，因此，准确、快速、稳定的颗粒粒径测量技术对国防、气象、煤矿、动力、化工等方面有重要意义。光散射法近 30 年来随着激光和计算机等技术的迅速发展而得到了普遍应用，是一种重要的燃烧诊断手段。它通过测量颗粒的散射光强来反演颗粒粒径分布，能够实现非接触式在线测量，具有适用性广、测量范围宽、响应快、易于实现自动化等特点 46。该方法的最大难点在于如何通过衰减光谱的反演获得准确的粒径分布，在理论上归结于第一类Fredholm 积分方程的求解。这是一个典型的不适定问题，直接求解具有很大的困难，因此，粒径分布反演算法的研究受到了众多研究者的重视 79。光全散射粒径测量方法分为独立模

11、式算法和非独立模式算法，独立模式算法事先无须假定粒径分布，通过求解离散线性方程组得到粒径分布；非独立模式算法需要假设被测颗粒系的颗粒尺寸分布满足某个已知的分布函数，通过一定的优化方法确定分布函数中的待定参数。事实上，在绝大多数实际应用中，往往不知道被测颗粒系的粒径分布规律，或者颗粒系的尺寸分布无法简单地用某个分布来描述，这就使非独立模式求得的结果变得不可靠 10。目前已经发展了多种独立模式下的粒径分布反演算法，每种算法都有各自的应用背景和局限性。传统算法如： Phillips-Twomey 算法、 Chahine迭代算法、共轭梯度法等，这些算法存在目标函数及导数计算复杂、抗噪声能力差及多峰粒径

12、重建困难等问题。近年来，遗传算法、模拟退火法等智能算法也已用于解决粒径测量问题，与传统方法相比，智能算法具有很好的全局搜索和抗噪声能力，但进化速度慢， 反演结果不稳定等缺点突出1112。本文作者利用光全散射法在独立模式下反演了几种粒径分布，其中正问题利用反常衍射近似 (ADA)计算得到估计值，测量值则通过 Mie 理论计算得到；反问题方法采用量子微粒群算法，并在标准量子微粒群算法的基础上提出了改进提出改进量子微粒群算法，在目标函数中利用Markov 随机场建立正则化项，提升了算法的反演能力，并进行了数值反演验证。1 正问题模型1.1 光全散射法测量原理光全散射法以光的散射理论为基础，当一束光强

13、为 I0，波长为 的平行单色光照射到厚度为 L的悬浮待测颗粒系时，由于颗粒对入射光的吸收和散射作用， 穿过颗粒系透射光的光强将减弱。根据 Lambert-Beer定律，如果假设颗粒系为服从一定粒径分布范围的多分散球形粒子系，并且颗粒间满足不相关单散射的条件 (忽略多次散射效应 )，则多分散球形颗粒系在波长为 时的消光值可以表示为 13maxmin ext0 ( , , )( ) 3ln ( )d( ) 2 DD Q m DI L N f D DI D (1)式中： I0为入射激光强度； I为透射光强； 0( )/ ( )I I 为波长为 的消光值，它可以通过人工实验测量得到；N为待测颗粒系的颗

14、粒总数； ( )f D 为颗粒系的体积频度分布； Dmax和 Dmin分别为颗粒粒径分布的上、下限；ext( , , )Q m D 为消光系数，它是波长 、 介质复折射率 m以及颗粒粒径 D的函数。将粒径分布区间划分为 M个子区间后，积分项可以变成离散项求和： ext0 1 ( , , )( ) 3ln ( )( ) 2 M jj jjj Q m DI L N c f DI D (2)式中： cj为数值积分系数； M 为整个待测粒径范围Dmin,Dmax内划分的子区间个数； Dj为各子区间的等效粒径。 由于消光系数 ext( , , )jQ m D 可以通过 Mie理论或者一些近似方法计算得到

15、，则可以通过测量多个波长下的消光值 0( )/ ( )I I 来反演颗粒系的体积频率分布( )jf D 。1.2 Mie 理论计算公式Mie 理论是一种球形颗粒散射特性的严格计算方法， Mie 散射利用光的电磁波性质，应用 Maxwell 方程对散射颗粒形成的边界条件进行求解，以得到光散射的物理量。可用于计算在单色平行光照射下任意尺寸和任意成分的球形颗粒的散射场 1。利用 Mie 理论推导得到的颗粒消光系数 ext( , , )jQ m D 的计算公式为ext 2 12 (2 1)( )n nnQ Re n a b (3)式中： /D ，为无因次粒径参数； an和 bn为 Mie系数，是与 B

16、essel 函数和 Hankel 函数有关的函数，( ) ( ) ( ) ( )( ) ( ) ( ) ( )n n n nn n n n nmx x m mx xa mx x m mx x (4)( ) ( ) ( ) ( )( ) ( ) ( ) ( )n n n nn n n n nm mx x mx xb m mx x mx x (5)万方数据中南大学学报 (自然科学版 ) 第 47 卷3924 1/2( ) /2 ( )n nx x J x (6)(2)1/2( ) /2 ( )n nx x H x (7)式中： r im m mi 为颗粒相对于周围介质的复折射率； 1/2( )nJ

17、 x 和 (2)1/2( )nH x 分别为半奇阶的第 1 类Bessel 函数和第 2 类 Hankel 函数。1.3 ADA 计算公式由于 Mie 理论公式计算繁琐，特别是计算中涉及到无穷级数求和问题，占用大量内存资源和计算时间， 不适合作为在线检测的正问题求解模型。反常衍射近似 (ADA)假设颗粒的消光主要是吸收以及透射光与衍射光之间的干涉引起的，它的计算速度要远快于 Mie理论的计算速度。最初的 ADA 只能适用于满足 1和 1m 1 条件的粒子， ZHAO 等 14引入边界效应，提出了修正 ADA 算法，扩展了的适用范围。对于球形粒子，标准的反常衍射近似衰减因子adQ 和边界效应 e

18、dgeQ 的表达式如下：ad 21 exp( ) 1 exp( )4 2 i iQ Re i (8)2/3edge 02 /Q c (9)式中： 2 ( 1)m ， 0 0.996 130c 。当粒径为 0.110m 时，修正消光系数为ext ad edge ad11 2/ 1/ 1 ( 1)Q Q Q m Q (10)2 反问题模型2.1 标准量子微粒群算法 (BQPSO)量子微粒群算法是由 SUN 等15提出。量子微粒群算法将 PSO 系统看成是一个量子系统，每个粒子具有量子行为，量子的状态由波函数 决定， 2 为粒子的位置 1 2( , , , , , )i i i ij iNx x x

19、 xX 的概率密度。在第 t 次迭代中，粒子 i 在 N 维搜索空间内以粒子的局部吸引因子 1 2( , , , , , )i i i ij iNq q q qQ 为中心在领域内搜索。 用 蒙 特 卡 洛 法 模 拟 后 ， 位 置 i X1 2( , , , , , )i i ij iNx x x x 可由下式表示： 1 ln(1/ )2ij ij ijx q L R (11)1 1 2 2 g,1 1 2 2ij jij CR p C R pq CR C R (12)m,2ij j ijL p x (13)m, 1Mj ijip p M (14)式中： R， R1和 R2为 0,1区间内服

20、从均匀分布的随机数； C1和 C2为加速系数； pij为粒子 i的个体历史最优位置在第 j维上的坐标； g,jp 为群体历史最优位置在第 j 维上的坐标； 为吸引扩散系数，在 1.781时可以保证量子微粒群的全局收敛。2.2 改进量子微粒群算法 (IQPSO)对于一般的优化过程来说，初始粒子位置相对分散有利于群体的寻优，在标准量子微粒群算法中，初始化中是利用随机数产生位置的，这样有可能在初始化的过程中，产生位置相近的粒子，本文将典型的logistic 映射应用于产生混沌的信号，这样初始化得到的粒子位置会相对分散。 (0) ()( 1) ( ) 1 ( ) ; 0,1, , 1r randr k

21、 r k r k k M (15)式中： 为混沌系数，当 9.3 时， logistic 映射将处于一个完全的混沌状态。 在标准量子微粒群算法中，局部吸引因子 iQ 是 1个非常重要的参数， 1 个不利的局部吸引因子容易让算法陷入局部收敛，因此，本文提出了一种服从正态分布的局部吸引因子，它以群体历史最优位置 gP 为平均值，以 m gP P 为标准偏差的正态分布随机数。可以通过下式来定义： g m g ,i N Q P P P (16)如同其他群体智能优化算法一样，标准量子微粒群算法在进化后期容易丧失种群的多样性，这样也不利于优化，本文引进了基因变异机制，它可以帮助粒子逃出局部最优位置而增加种

22、群的多样性。位置更新可由下式来确定： msgn( ) (1 ) ln(1/ )i i iR r R X Q P X (17)33,0,R Rr R (18)式中： sgn()为符号函数； R 为在 1,1 区间内服从均匀分布的随机数； R 为服从标准正态分布的随机数；R3为 0,1区间内服从均匀分布的随机数； 为变异因子，为 1 个小于 0.5 的正数。2.3 含正则化项的量子微粒群算法 (RQPSO)当反演参数过多时，标准量子微粒群算法和改进万方数据第 11 期 张彪，等：基于 RQPSO 的颗粒粒径分布反演算法 3925量子微粒群算法容易达到维数极限，导致反演结果不稳定而失效，这时通常需要

23、对目标函数进行光滑化处理。本文在改进量子微粒群算法的基础上，利用Markov随机场理论认为每个待测粒径范围 Dmin,Dmax的子区间上的体积频率密度仅仅与它邻近子区间的体积频率密度有关，与其他子区间内的体积频率密度无关，整个待测粒径范围内的体积频率密度为一个Markov 随机场，因而建立起目标函数的正则化项，其中目标函数的表达式如下： obj 2 2F M E AE (19)3 2 1 0 0 0 0 0 02 5 2 1 0 0 0 0 01 2 6 2 1 0 0 0 00 0 0 00 0 0 00 0 0 00 0 0 0 1 2 6 2 10 0 0 0 0 1 2 5 20 0

24、0 0 0 0 1 2 3 A (20)式中： M 为测量值； E 为估计值； A 为光顺矩阵； 2为向量的 2-范数。3 结果与讨论3.1 正问题验证为了验证 ADA 的可靠性，本文通过 Mie 理论的计算结果来进行检验。在可见光范围内，光全散射法的最佳粒径测量范围为 0.110.0m1，本文粒子的复折 射 率 参 照 典 型 飞 灰 粒 子 ， 它 的 范 围 分 别 为r 1.18,1.92m 和 i 0.01,1.13m 6。因此，分别利用Mie理论和 ADA计算复折射率为 1.44 1.00m i 的粒子在波长为 0.45 m 下消光因子和复折射率为1.235 0.01m i 的粒子

25、在波长为 0.73 m 下消光因子，计算结果如图 1 所示。从图 1 可以看出： ADA 的计算结果与 Mie 理论的计算结果十分吻合，但 ADA 的计算效率是 Mie 理论的 1758 倍。因此， ADA 适合作为粒径反演中正问题的求解方法，而 Mie 理论的计算结果可以作为粒径反演中的测量值。3.2 函数反演为了比较标准量子微粒群算法 (BQPSO)和改进量子微粒群算法 (IQPSO)的性能，文中通过选用 Sphere函数、 Rastrigin 函数和 Schaffer 函数 3 种标准测试函(a)m=1.44+1.00i， =0.45m;(b)m=1.235+0.01i， =0.73m图

26、 1 ADA与Mie理论的计算结果对比Fig. 1 ComparisonofcalculationresultsbyADAandMietheory数对这 2 个算法进行研究，其具体参数如表 1 所示。在这 2 种算法中粒子的总数都为 Np=50，最大迭代次数都为 Nt=50 000，吸引扩散系数 随着迭代从0.9 线性下降到 0.4，在改进量子微粒群算法中混沌系数为 3.9，变异因子为 0.05，本文中若未作特别说明，算法的参数设置保持不变。为了对比这 3 种算法的计算效率和鲁棒性，将每种算法都独立执行100 次，计算结果如表 2 所示。从表 2 可以看出：随着测试函数维数的增加，这 2 种算

27、法的迭代次数均增加，标准偏差也在增大；相比于标准量子微粒群算法而言，改进量子微粒群算法在计算效率和稳定性上得到了很大的提升，说明了改进量子微粒群算法在性能上具有很大的优势。在后续的算例中，本文都采用改进量子微粒群算法作为反问题的求解方法。万方数据中南大学学报 (自然科学版 ) 第 47 卷3926 表 1 3个测试函数的具体参数Table 1 Parametersofthreetestfunctions函数 表达式 维数 N 搜索空间 精度Sphere 21N ii x 10 100,100N 10100Rastrigin 21 10cos(2 ) 10N i ii x x 10 5.12,5

28、.12N 10100Schaffer 22 2 2 2 21 2 1 2sin 0.5 1 0.001( ) 0.5x x x x 2 100,100N 1016表 2在固定精度下3种算法执行100次的迭代次数Table 2 Iterationsofthreealgorithmsin100independentrunsunderfixedaccuracy函数 Dim BQPSO IQPSOSphere 1 609.232.7 603.833.12 934.832.7 787.633.33 1220.839.8 936.933.34 1504.244.2 1089.033.15 1800.642

29、.2 1284.930.46 2086.348.4 1516.037.47 2398.261.2 1811.643.58 2729.752.5 2180.654.49 3120.860.0 2636.466.410 3493.455.1 3150.284.5Rastrigin 1 635.840.3 598.029.42 1086.397.1 791.032.83 1579.0298.4 961.347.24 2399.3978.4 1155.865.65 3374.22006.4 1398.1125.26 5268.94280.5 1672.1184.47 6935.66902.1 2080

30、.2202.28 8443.06495.9 2560.2286.89 15281.912710.1 3171.7406.510 16031.311891.7 3959.6736.0Schaffer 2 983.1336.2 341.6539.63.3 单峰分布反演为了考察量子微粒群算法在粒径反演中的效果，本文利用改进量子微粒群算法 (IQPSO)和含正则化项的改进量子微粒群算法 (RQPSO)反演 1 个服从单峰R-R 分布的颗粒系，表达式为 2.55 1 2.552.55( ) exp0.45 0.45 0.45D Df D (21)已知颗粒系的粒径范围为 0.011.00m，颗粒的复折射率

31、为 1.235 0.01m i 。首先将粒径范围均匀的分成 10 个子区间，通过在可见光波段内选择 10 个波段来测量粒子系的消光值，反演出 10 个子区间内的体积频率密度，反演结果如图 2(a)所示，其中算法的搜索范围为 max min0,( )/ D D M 。为了衡量反演结果的优劣，定义误差为 1/22ret ori1 1/22ori1( ) ( )( )M i ii M iif D f Df D (22)式中： fret和 fori分别为颗粒体积频率密度的反演值和真值。 它们的反演误差分别为 IQPSO,10 0.211 835 2 和RQPSO,10 0.137 687 8 ，从图

32、2(a)可以看出： 2 种算法反演出的分布都基本复合分布规律；从误差水平上看， RQPSO 算法比 IQPSO 算法的误差要小，说明加上正则化后，反演结果更接近真实分布。 为了考察 2 种算法的维数极限，将粒径分布范围的子区间划分成 20 份，反演的误差水平分别为IQPSO,20 0.598 462 1 和 RQPSO,20 0.105 985 5 ，反演结果如图 2(b)所示。从图 2(b)可以看出： IQPSO 算法的反演结果已经严重偏离真实分布，而 RQPSO 算法的反演结果能很好地符合真实分布。 为了考察 RQPSO 算法能力的提升程度，将粒径分布范围的子区间分别划分成 30 份和 5

33、0 份，它们的反 演 误 差 水 平 分 别 为 RQPSO,30 0.118 601 8 和RQPSO,50 0.149 625 6 ，反演结果分别如图 2(c)(d)所示。从图 2 可以看出：反演的分布都基本符合真实分布。说明 RQPSO 的反演能力得到了很大提升，不仅提高了反演的维数极限，而且提高了反演的精度。3.4 双峰分布反演为了进一步的考察含正则化项的改进量子微粒群算法的性能，本文利用其反演了 1 个服从双峰 R-R 分布的颗粒系，表达式为万方数据第 11 期 张彪，等：基于 RQPSO 的颗粒粒径分布反演算法 3927子区间数量/ 份：(a)10 ；(b)20 ；(c)30 ；(

34、d)50图 2单峰分布的反演结果Fig. 2 Retrievalresultsofmonomodaldistribution5.9 1 5.95.9( ) 0.3 exp0.25 0.25 0.25D Df D 10.0 1 10.010.0(1 0.3) exp0.7 0.7 0.7D D (23)颗粒系的粒径范围和复折射率与上述颗粒系相同，分别将粒径范围划分成 20 份和 30 份，它们的反演 误 差 水 平 分 别 为 RQPSO,20 0.153 405 0 和RQPSO,30 0.124 550 8，反演结果如图 3 所示。从图3 可以看出：反演的分布都基本符合真实分布，进一步证明了

35、 RQPSO 算法的可靠性。3.5 误差分析为了考察测量误差对反演结果的影响，本文以单峰分布为例，将粒径范围划分成 30 份，分别加入 0，2%， 5%和 10%的测量误差后，反演误差水平分别为0 0.063 856 4 ， 2 0.227 241 5 ， 5 0.244 149 3 和 10 0.346 607 7 ，反演结果如图 4 所示。从图 4可以看出：随着测量误差的增大，反演的误差水平有所增加，但反演的分布基本符合真实分布，说明了含正则化的改进量子微粒群算法具有很好的抗噪性能。 子区间数量 /份： (a)20； (b)30图 3双峰分布的反演结果Fig. 3 Retrievalres

36、ultsofbimodaldistribution万方数据中南大学学报 (自然科学版 ) 第 47 卷3928测量误差 /%： 1 Original； 2 0； 3 2； 4 5； 5 10。图 4含测量误差时的反演结果Fig. 4 Retrievalresultswithmeasurementerror4 结论1) 针对标准量子微粒群算法，提出了改进量子微粒群算法和含正则化项的改进量子微粒群算法，通过函数反演证明了改进量子微粒群算法相对于标准量子微粒群算法具有更高的效率和稳定性。 2) 将改进的量子微粒群算法引入到粒径分布的反演当中，利用光全散射法在独立模式下反演了服从单峰分布和双峰分布的颗

37、粒系粒径分布，得出含正则化项的改进量子微粒群算法提高了改进量子微粒群算法的维数极限，并具有更高的精度和可靠性。 3) 含正则化项的改进量子微粒群算法具有很好的抗噪性能，为独立模式下的粒径分布反演提供了一个新的方法。参考文献：1 唐红 . 光全散射法颗粒粒径分布反演算法的研究 D. 哈尔滨 :哈尔滨工业大学电气工程及自动化学院 ,2008:18.TANG Hong. Study of inversion algorithm of particle sizedistribution using total light scattering methodD. Harbin:Harbin Instit

38、ute of Technology. School of Electrical EngineeringandAutomation,2008:18.2 QI Hong, RUAN Liming, WANG Shenggang, et al.Applicationof multi-phase particle swarm optimization technique to retrievethe particle size distributionJ. Chinese Optics Letters, 2008,6(5):346349.3 TANG Hong. Retrieval of spheri

39、cal particle size distributionwith an improved Tikhonov iteration methodJ. ThermalScience,2012,16(5):14001404.4 菅立川 , 齐宏 , 阮立明 , 等 . 微粒群算法反演粒径分布的研究J. 节能技术 ,2007,141(25):713.JIAN Lichuan, QI Hong, RUAN Liming, et al. Inverse analysisof particles size distribution by particle swarm optimizeralgorithmJ

40、. Energy Conservation Technology, 2007, 141(25):713.5 HE Zhenzong, QI Hong, WANG Yuqing, et al. Inverseestimation of the spheroidal particle size distribution using AntColony Optimization algorithms in multispectral extinctiontechniqueJ.OpticsCommunications,2014,328(19):822.6 QIHong,NIU Chunyang,GON

41、GShuai,etal.Applicationofthehybrid particle swarm optimization algorithms for simultaneousestimation of multi-parameters in a transient conduction-radiation problemJ. International Journal of Heat and MassTransfer,2015,83(4):428440.7 孙晓刚 , 唐红 , 原桂彬 . 非独立模式算法下粒径分布反演及分类的研究 J. 光谱学与光谱分析 ,2008,28(5):1111

42、1114.SUNXiaogang,TANGHong,YUANGuibin.Studyofinversionand classification of particle size distribution under dependentmodel algorithmJ. Spectroscopy and Spectral Analysis, 2008,28(5):11111114.8 QI Hong, HE Zhenzong, GONG Shuai, et al. Inversion ofparticle size distribution by spectral extinction tech

43、nique usingthe attractive and repulsive particle swarm optimizationalgorithmJ.ThermalScience,2015,19(6):21512160.9 HE Zhenzong, QI Hong, YAO Yuchen, et al. Inverse estimationof the particle size distribution using the fruit fly optimizationalgorithmJ. Applied Thermal Engineering, 2015, 88(9):306315.

44、10 徐峰 , 蔡小舒 , 苏明旭 , 等 . 独立模式算法求解颗粒粒径分布的研究 J. 中国激光 ,2004,31(2):223228.XUFeng, CAIXiaoshu,SUMingxu,et al.Studyof independentmodelalgorithmfordeterminationofparticlesizedistributionJ.ChineseJournalofLasers,2004,31(2):223228.11 王丽 , 孙晓刚 . 基于模式搜索的光谱消光粒径分布反演算法的研究 J. 光谱学与光谱分析 ,2013,33(3):618622.WANG Li, SU

45、N Xiaogang. Research on pattern search methodfor inversion of particle size distribution in spectral extinctiontechniqueJ. Spectroscopy and Spectral Analysis, 2013, 33(3):618622.12 HE Zhenzong, QI Hong, YAO Yuchen, et al. An effectiveinversion algorithm for retrieving bimodal aerosol particle sizedi

46、stribution from spectral extinction dataJ. Journal ofQuantitativeSpectrascopyandRadiativeTransfer,2014,149(18):117127.13 QIHong, ZHANG Biao,RENYatao,et al.Retrievalofsphericalparticle size distribution using ant colony optimizationalgorithmJ.ChineseOpticsLetters,2013,11(11):112901.14 ZHAO Jianqi, HU Yinqiao. Bridging technique for calculatingthe extinction efficiency of arbitrary shaped particlesJ.AppliedOptics,2003,42(24):49374945.15 SUN Jun, FENG Bin, XU Wenbo. Particle swarm optimizationwith particl