文章编号:1674-7070(2011)04-0360-05 应用粒子群算法优化地位指数曲线模型参数研究 摘要 以优势高和地位指数的估测误差最 小为目标函数,采用粒子群优化算法求 解地位指数曲线模型的参数.结合实例 与免疫算法比较,结果表明:粒子群优化 算法求解的参数使模型的总体误差更 小,精度更高,拟合效果更理想,更加科 学合理,同时也提高了幼林的估算精度. 研究的结果为森林经营中生长模型参数 的求解以及相关研究提供了新的应用思 路,也拓宽了粒子群优化算法在林业科 学中的应用. 关键词 粒子群算法;地位指数;参数优化 中图分类号TP39l;¥711 文献标志码A 收稿日期2011-01—18 资助项目福建省自然科学基金项目(2009. J05043;2011JO5044);福建农林大学青年教师 基金项目(2010019) 作者简介 胡欣欣,女,博士,从事计算机在林业科学 中的应用研究.xinxinhu@fjau.edu.cn 1福建农林大学计算机与信息学院,福州 350002 胡欣欣 王李进 0 引言 在森林经营管理中,准确地评价森林的立地质量,估计森林的生 产力,是一项十分重要的基础工作.地位指数是常用的评价方法,其 模型是评价森林立地质量和预估各年龄优势高的有用工具.Richards 方程…、Sloboda方程 J、Amateis方程 以及Schumacher方程 等常 用于拟合地位指数曲线模型.文献[5]在Schumacher方程的基础上建 立了能显示表达优势高和地位指数的多形地位指数曲线模型.上述 模型采用优势高的理论值和实际值的残差平方和最小,或者兼顾地 位指数的理论值和实际值的残差平方和最小为目标函数,求解其参 数.在优化拟合地位指数曲线模型的参数时,常用单纯形算法、遗传 算法和免疫算法进行求解,并取得较好的效果,为立地质量评价提供 了依据.单纯形算法是一种利用目标函数值的简单和良好局部搜索 优化算法,所求的解对线性化后的模型是最优的,但对原模型不一定 是最优的 J.遗传算法在一定条件下具有全局收敛性,但算法的交 叉、变异和选择算子都是在概率意义下随机进行的,虽然保证了种群 的群体进化,但在一定程度上不可避免出现退化现象,同时参数设置 也较为复杂_6j.免疫算法与遗传算法相似,也是从随机生成的初始解 群出发,采用复制、交叉、变异等算子进行操作,产生比父代优越的子 代,循环执行,逐渐逼近最优解,但算法为了建立最优可行解集合,获 得解特征样本,需要进行大量计算 . 粒子群优化算法(Particle Swarm Optimization,PSO)与遗传算法类 似,是一种基于进化的优化工具.同遗传算法比较,PSO的优势在于简 单容易实现并且没有许多参数需要调整.目前,PSO在林业中得到初 步应用 J,但用于求解地位指数曲线模型参数未见相关报道. 因此,本文采用粒子群优化算法求解地位指数曲线模型的参数, 旨在为建立简便实用且满足精度要求的多形地位指数曲线模型提供 一条新的途径,同时也拓宽粒子群算法在林业科学中的应用. 1粒子群优化算法 粒子群优化算法最早是由Kennedy等 提出的,它的基本原理源 于对鸟群捕食行为的仿真,是一种基于群智的计算智能方法.算法简 单易于实现,同时又有深刻的智能背景,既适合科学研究,又特别适 合工程应用.因此,PSO算法一提出,立刻引起了演化计算等领域学者 曲未信垂 乎学报:自然科学版,2011,3(4):360—364 361 Journal of Nanjing University of Information Science and Technology:Natural Science Edition,2011,3(4):360_364 们的广泛关注_6].目前,已经广泛应用于函数优化、 程等,拟合的关键是参数的求解.以往在求解上述模 神经网络训练、模糊系统控制以及其他遗传算法的 型的参数时,采用的目标函数是各年龄时的优势高 应用领域. 理论值和实际值的残差平方和最小,见式(3). PSO算法求解优化问题时,问题的解对应于搜 索空间中一只鸟,被抽象为没有质量和体积的“粒 Q=∑(日一百) . (3) 其中,日表示优势高实际值,Ⅳ表示优势高理论值. 采用这种目标函数求解参数并不一定满足在给定优 子”,并将其延伸到Ⅳ维空间.每个粒子都有自己的 位置和速度(决定飞行的方向和距离),还有一个由 被优化问题决定的适应值(fitness).每个粒子知道 势高和年龄时所确定的地位指数的误差最小.为了 自己到目前为止发现的最好位置(P )和现在的位 置,除此之外,每个粒子还知道到目前为止整个群体 中所有粒子发现的最好位置(g ),可看做是粒子 同伴的经验.各个粒子记忆并追随当前的最优粒子 在解空间中搜索,这样,每次迭代的过程不是完全随 机的,如果找到较好解,将会以此为依据来寻找下一 个解. PSO算法通常采用随机化的方式为粒子产生初 始位置和速度(随机初始解).假设d维搜索空间的 各个粒子的位置和速度分别为X=( , ,…, ) 和V=( , ,…,V ),在此之后的每一次迭代中,粒 子通过跟踪2个最优解来更新自己,第1个就是粒子 本身所找到的最好解,即个体极值点(P ),另一个 是整个粒子群目前找到的最好解,称为全局极值点 (g ).在找到这2个最优值时,粒子根据式(1)和 式(2)来更新自己的速度和新的位置. t“=O)V:+C1r1(Pb :一 :)+c2r2(gb。 :一 :),(1) tii = =X,‘+i Vi+ t ・. (2) 式中,tO是惯性系数,其主要作用是产生扰动,以防 止算法的早熟收敛,c。和c 是加速系数(或称学习 因子),分别调节向个体最好粒子和全局最好粒子方 向飞行的最大步长,若太小,则粒子可能远离目标区 域,若太大,则会导致突然向目标区域飞去,或飞过 目标区域.合适的C,和c 可在加快收敛速度的同时 还能不易陷入局部最优,通常令c =c:=2.r 和r。 为0和1之间均匀分布的随机数 』. 粒子群算法发展到现在有很多种变形及其改进 算法.为了有效地控制粒子的飞行速度使算法达到 全局探测与局部开采两者的有效平衡,本文采用带 压缩因子的粒子群算法¨。。求解目标函数. 2地位指数曲线模型 2.1求解参数的目标函数描述 拟合地位指数曲线模型常用的方程有Richards 方程、Sloboda方程、Amateis方程以及Schumacher方 预估各年龄时的优势高有意义,文献[5]对式(3)进 一步改进,使得兼顾地位指数的理论值和实际值的 残差平方和最小,见式(4). Q=∑((日一 ) +(s,一S,) ). (4) 其中,s,表示地位指数实际值,S,表示地位指数理 论值, 和 由Schumacher方程推导获得,分别见 式(5)和(6). 日=6(1一(TO/ ))×SI( 2 【卜 2)T/On, (5) SI=(H×b T/O ’) 2+(1-b2’vo/ ”. (6) 其中,b。和b 为待求解的参数, 为年龄, 为基准 年龄. 2.2模型评价 模型评价主要是对模型的误差分析,本文采用 均方根误差,即剩余标准差(简称标准差)和平均相 对误差绝对值(简称平均误差)作为评价模型误差的 指标,见式(7)和(8). E一。= √(厂— ( —— y———— —) )—一 (7) =e( l y I/ye) (8) 其中,Y 是实测值,Y 是理论值,n是观察样本集的 大小. 3参数优化 应用带压缩因子的粒子群算法求解地位指数曲 线模型的参数,其实质是将式(4)作为待求解的优化 函数.当优化函数的适应度为最小时,函数的解即为 所求模型的最理想的参数.模型参数优化的具体步 骤如下. 步骤1.初始化粒子群种群规模,随机初始化种 群中每个粒子的位置和速度. 步骤2.按式(4)评价每个粒子的适应度,初始 化个体极点(P )和全局极点(g 。 ). 步骤3.按式(1)和式(2)更新粒子的速度和 位置. 胡欣欣,等.应用粒子群算法优化地位指数曲线模型参数研究 362 HU Xinxin,et a1.Applying particle swarm optimization algorithm to optimize the site index Qu[Ne mtMel 步骤4.更新Pb ̄st和g 步骤5.若满足停止条件,迭代停止,输出最优 解,否则转步骤3. 数的2个误差指标都小于文献[5]中的优化算法. 表3给出了各龄阶对优势高和地位指数的估测 误差的检验.基于免疫算法的优势高和地位指数的 4 应用实例 4.1实验数据与处理 标准差以及平均误差在高龄阶上微弱优于粒子群算 法,但低龄阶上粒子群算法明显优于免疫算法.基于 粒子群算法的优势高和地位指数的标准差分别低于 本文的实验数据引用文献[5]中的基础数据,其 0.5和1.0,而免疫算法有6、8龄阶超过0.5和1.0. 中基准年龄为20 a,具体见表1. 粒子群优化模型的平均误差都不超过1 0%. 基于表1的数据,按模型参数优化的具体步骤 表4给出了各地位指数对优势高和地位指数的 求解模型参数b。和b .其中,种群规模取50,最大迭 估测误差的检验.基于免疫算法的优势高和地位指 代次数取100,c =c =1.495, :0.729.算法满足 数的标准差以及平均误差在较高地位指数上微弱优 停止条件,其最优解为b,=2.441 2,b =0.862 4,目 于粒子群算法,但低地位指数上粒子群算法明显优 标函数最小值为30.698 0. 于免疫算法.基于粒子群算法的优势高和地位指数 4.2结果分析 的标准差分别低于0.5和0.9,而免疫算法在8.05 表2给出了模型的整体误差分析,以及与文献 地位指数上,优势高和地位指数分别达到0.8和 [5]研究结果的比较.表2中的标准差和平均误差是 1.64.粒子群优化模型的平均误差都不超过7%,免 将表1的所有数据综合在一起计算得出的总标准差 疫算法有2个地位指数的优势高的平均误差超过 有1个地位指数的预测地位指数平均误差超 和总平均误差.结果表明:粒子群算法的目标函数值 7%, 明皿优于免疫算法的目标函数值;优势高和地位指 讨7%.表1 各地位指数和各龄阶下的优势高平均值 Table 1 Average value of dominant height under each site index and each age gradation 5结束语 翥嘉 霎 票 翥薯 粒子群优化算法的优势在于简单容易实现并且 效果更理想,同时也较好地提高了幼林时的估算精 曲未馐垂 乎学报:自然科学版,2011,3(4):360—364 363 Journal of N ing University of Information Science and Technology:Natural Science Edition,201l,3(4):360-364 表3在各龄阶下不同算法估测优势高和地位指数的误差 Table 3 Errors about dominant height and site index under each age gradation estimated by PSO and immune algorithm 表4在各地位指数下不同算法估测优势高和地位指数的误差 Table 4 Errors about dominant height and site index under each site index estimated by PSO and immune algorithm 度,使得模型更加科学合理和更具推广应用价值.研 boda multi—shape site index model of masson pine I J』. Journal of Biomathematics,2002,17(4):489-493 究的结果为森林经营中生长模型参数的求解以及相 『j J 关研究提供了新的应用思路,也拓宽了粒子群算法 在林业科学中的应用. 陈信旺,庄晨辉,江希钿.福建柏人工林地位指数曲线 模型的研究[J].华东森林经理,2005,19(1):7一l0 CHEN Xinwang.ZHUANG Chenhui.JIANG Xidian.Re— search on status index curve model of FOKIENIA H0D. 参考文献 Refefences GINSII plantation l J 1.East China Forest Management. 2005.19(1):7.10 骆期邦,吴志德,蒋菊生,等.Richards函数拟合多形 地位指数模型的研究[J].林业科学研究,1989,2 (6):534—539 [4] 洪伟,吴承祯,闫淑君.广义Schumacher模型的改进 及其应用[J].应用生态学报,2004,15(2):241-244 HONG Wei.WU Chengzhen.YAN Shujun.Modification and its application of generalized Schumacher model f J i. LUO Qibang,WU Zhide,JIANG Jusheng,et a1.A study on the establishment of polym。rphic site index model by r 1 Chinese Journal f oApplied Ecology,2004,15(2):241—244 adopting Richards function[J].Forest Research,1 989,2 (6):534-539 江希钿,庄晨辉,陈信旺,等.免疫进化算法在建立地 位指数曲线模型中的应用[J].生物数学学报,2007, 22(3):515_519 JIANG Xidian,ZHUANG Cheuhui.CHEN Xinwang.et at.The application of immune evo[utionar' ̄algorithm in [2] 吴承祯,洪伟,林成来.马尾松人工林Sloboda多形地 位指数模型研究[J].生物数学学报,2002,17(4): 489-493 setting up the site index curve model f J f.Journal of Bio— mathematics,2007,22(3):515_519 WU Chengzhen,HONG Wei,LIN Chenglai.Study on Slo一 [6] 钟一文,宁正元,蔡荣英,等.一种改进的离散粒子群 胡欣欣,等.应用粒子群算法优化地位指数曲线模型参数研究 364 HU Xirtxin,et a1.Applying particle swaFnl optimization algorithm to optimize the site index cHive model 优化算法[J].小型微型计算机系统,2006,27(10): 1893 1896 ZHONG Yiwen,NING Zhengyuan,CAI Rongying,et a1. 科学版,2010,2(4):320—323 WANG Lijin,HU Xinxin,NING Zhengyuan.Projection pursuit cluster model based on particle swaFnl optimiza— An improved discrete particle swarm optimization algo— tion and its application[J].Journal of Nanjing Universit)r of Information Science&Technology:Natural Science E— rithm j J 1.Journal of Chinese Computer Systems,2006, 27(10):1893—1896 焦李成,尚荣华,马文萍,等.多目标优化免疫算法理 论和应用[M].北京:科学出版社,2010 JIAO Licheng,SHANG Ronghua,MA Wenping,et a1. dition,2010,2(4):320 323 Kennedy J,Eberhart R.Particle swarm optimization l C If Proceedings of IEEE Internationa Conference on Neu. ral Networks.1995,4:1942—1948 Theory and application of multi—objective optimization Clerc M,Kennedy J.The particle swarm—explosion,stabil— ity,and convergence in a multidimensional complex space immune algorithm[M j.Beijing:Science Press,2010 王李进,胡欣欣,宁正元.基于粒子群优化的投影寻踪 聚类模型及其应用[J].南京信息工程大学学报:自然 [c]//IEEE Transactions on Evolutionary Computation, 2002,6(1):58-73 Applying particle swarm optimization algorithm to optimize the site index curve model HU Xinxin WANG Lijin 1 College of Computer&Infomatrion Science,Fujian Agriculture and Forestyr University,Fuzhou 350002 Abstract The parameters of the site index curve model were solved by using particle swarm optimization(PSO) algorithm,with the minimum error in estimation of dominant height and site index as the target function.An applica— tion example is introduced to test the proposed method,and the estimation result is compared between PSO and im— mune algorithm.Comparison result shows that the parameters solved by PSO can decrease the overall error,increase the precision,improve the fitting effect of the site index curve model,thus increase the estimation precision of young forest.This research is hoped to provide new idea for parameter solving of grow model in forest management and re— lated research,and expand the application of PSO in forestry science as wel1. Key words particle swarm optimization;site index;parameter optimization