粒子群算法概述

粒子群优化算法（PS0）ParticleSwarmOptimization由Kennedy和Eberhart于1995年提出．群体迭代，粒子在解空间追随最优的粒子进行搜索.简单易行粒子群算法:收敛速度快设置参数少已成为现代优化方法领域研究的热点．粒子群算法发展历史简介PSO算法特点PSO算法属于进化算法的一种，和模拟退火算法相似，它也是从随机解出发，通过迭代寻找最优解，它也是通过适应度来评价解的品质，但它比遗传算法规则更为简单，它没有遗传算法的“交叉”(Crossover)和“变异”(Mutation)操作，它通过追随当前搜索到的最优值来寻找全局最优。这种算法以其实现容易、精度高、收敛快等优点引起了学术界的重视，并且在解决实际问题中展示了其优越性。粒子群算法是一种并行算法。PSO应用PSO存在的问题（1）PSO存在的问题（2）粒子群算法的基本思想粒子群算法的思想源于对鸟群捕食行为的研究．模拟鸟集群飞行觅食的行为，鸟之间通过集体的协作使群体达到最优目的，是一种基于SwarmIntelligence（群智能）的优化方法。因此，也称粒子群优化算法或鸟群觅食算法。马良教授在他的著作《蚁群优化算法》一书的前言中写到：大自然对我们的最大恩赐！“自然界的蚁群、鸟群、鱼群、羊群、牛群、蜂群等，其实时时刻刻都在给予我们以某种启示，只不过我们常常忽略了大自然对我们的最大恩赐！......”向大自然学习——生物智能粒子群算法的基本思想设想这样一个场景：一群鸟在随机搜索食物在这块区域里只有一块食物;所有的鸟都不知道食物在哪里;但它们能感受到当前的位置离食物还有多远.已知那么:找到食物的最优策略是什么呢？搜寻目前离食物最近的鸟的周围区域．根据自己飞行的经验判断食物的所在。PSO正是从这种模型中得到了启发．PSO的基础:信息的社会共享生物学家对鸟(鱼)群捕食的行为研究社会行为(Social-OnlyModel)个体认知(Cognition-OnlyModel)粒子群特性算法介绍每个寻优的问题解都被想像成一只鸟，称为“粒子”。所有粒子都在一个D维空间进行搜索。所有的粒子都由一个fitnessfunction（适应度函数）确定适应值以判断目前的位置好坏。每一个粒子必须赋予记忆功能，能记住所搜寻到的最佳位置。每一个粒子还有一个速度以决定飞行的距离和方向。这个速度根据它本身的飞行经验以及同伴的飞行经验进行动态调整。粒子群优化算法求最优解D维空间中，有N个粒子；粒子i位置：xi=(xi1,xi2,…xiD)，将xi代入适应函数f(xi)求适应值；粒子i速度：vi=(vi1,vi2,…viD)粒子i个体经历过的最好位置：pbesti=(pi1,pi2,…piD)种群所经历过的最好位置：gbest=(g1,g2,…gD)通常，在第d（1≤d≤D）维的位置变化范围限定在内，速度变化范围限定在内（即在迭代中若超出了边界值，则该维的速度或位置被限制为该维最大速度或边界位置）min,max,[X,]ddX,max,[-V,]maxddVidvid、x粒子i的第d维速度更新公式：粒子i的第d维位置更新公式：—第k次迭代粒子i飞行速度矢量的第d维分量—第k次迭代粒子i位置矢量的第d维分量c1,c2—加速度常数，调节学习最大步长r1,r2—两个随机函数，取值范围[0，1]，以增加搜索随机性w—惯性权重，非负数，调节对解空间的搜索范围pbestid—粒子i个体经历过的最好位置，gbestd—粒子群经历过的最好位置kk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx11kkkidididxxvkidvkidx粒子速度更新公式包含三部分：第一部分为粒子先前的速度第二部分为“认知”部分，表示粒子本身的思考，可理解为粒子i当前位置与自己最好位置之间的距离。第三部分为“社会”部分，表示粒子间的信息共享与合作，可理解为粒子i当前位置与群体最好位置之间的距离。kk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx區域最佳解全域最佳解運動向量慣性向量12X=X,X,...,Xiiiid12V=V,V,...,ViiiidStudyFactor12(1)()()()()()()()ididididgdidttttvwvcrandpxcrandpx(1)()()iiitttxxvpg1111122V=V+C*r*(Pbest-X)+C*r*(gbest-X)kkkkiiiii11X=X+Vkkkiii12NX=X,X,...,Xiiii12NV=V,V,...,Viiii算法流程1.Initial：初始化粒子群体（群体规模为n），包括随机位置和速度。2.Evaluation：根据fitnessfunction，评价每个粒子的适应度。3.FindthePbest：对每个粒子，将其当前适应值与其个体历史最佳位置（pbest）对应的适应值做比较，如果当前的适应值更高，则将用当前位置更新历史最佳位置pbest。4.FindtheGbest：对每个粒子，将其当前适应值与全局最佳位置（gbest）对应的适应值做比较，如果当前的适应值更高，则将用当前粒子的位置更新全局最佳位置gbest。5.UpdatetheVelocity：根据公式更新每个粒子的速度与位置。6.如未满足结束条件，则返回步骤2通常算法达到最大迭代次数或者最佳适应度值的增量小于某个给定的阈值时算法停止。maxG粒子群优化算法流程图开始初始化粒子群计算每个粒子的适应度根据适应度更新pbest、gbest，更新粒子位置速度结束noyes达到最大迭代次数或全局最优位置满足最小界限？2維簡例Note合理解目前最優解區域最佳解全域區域粒子群算法的构成要素-群体大小mm是一个整型参数．m很小:m很大:当群体数目增长至一定水平时，再增长将不再有显但收敛速度慢．著的作用．陷入局优的可能性很大．PSO的优化能力很好，粒子群算法的构成要素-权重因子权重因子：惯性因子—、学习因子—失去对粒子本身的速度的记忆社会经验部分前次迭代中自身的速度自我认知部分基本粒子群算法粒子的速度更新主要由三部分组成：惯性因子kvkk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx111()kididcrpbestx122()kdidcrgbestx粒子群算法的构成要素-权重因子权重因子：惯性因子、学习因子社会经验部分前次迭代中自身的速度自我认知部分粒子的速度更新主要由三部分组成：kvkk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx111()kididcrpbestx122()kdidcrgbestx学习因子无私型粒子群算法“只有社会，没有自我”迅速丧失群体多样性，易陷入局优而无法跳出．粒子群算法的构成要素-权重因子权重因子：惯性因子、学习因子社会经验部分前次迭代中自身的速度自我认知部分粒子的速度更新主要由三部分组成：kvkk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx111()kididcrpbestx122()kdidcrgbestx自我认知型粒子群算法“只有自我，没有社会”完全没有信息的社会共享，导致算法收敛速度缓慢学习因子粒子群算法的构成要素-权重因子权重因子：惯性因子、学习因子社会经验部分前次迭代中自身的速度自我认知部分粒子的速度更新主要由三部分组成：kvkk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx111()kididcrpbestx122()kdidcrgbestxc1,c2都不为0，称为完全型粒子群算法完全型粒子群算法更容易保持收敛速度和搜索效果的均衡，是较好的选择．粒子群算法的构成要素-最大速度但在于维护算法的探索能力与开发能力的平衡．Vm较大时，探索能力增强，作用:Vm较小时，开发能力增强，Vm一般设为每维变量变化范围的10％~20％.但粒子容易飞过最优解．容易陷入局部最优．粒子群算法的构成要素-邻域的拓扑结构全局粒子群算法和局部粒子群算法．粒子群算法的邻域拓扑结构包括两种，一种是将群体内所有个体都作为粒子的邻域，另一种是只将群体中的部分个体作为粒子的邻域．群体历史最优位置邻域拓扑结构决定由此，将粒子群算法分为粒子群算法的构成要素-邻域的拓扑结构全局粒子群算法1.粒子自己历史最优值2.粒子群体的全局最优值局部粒子群算法1.粒子自己历史最优值2.粒子邻域内粒子的最优值邻域随迭代次数的增加线性变大，最后邻域扩展到整个粒子群。经过实践证明：全局版本的粒子群算法收敛速度快，但是容易陷入局部最优。局部版本的粒子群算法收敛速度慢，但是很难陷入局部最优。现在的粒子群算法大都在收敛速度与摆脱局部最优这两个方面下功夫。其实这两个方面是矛盾的。看如何更好的折中了。粒子群算法的构成要素-停止准则停止准则一般有如下两种：最大迭代步数可接受的满意解粒子群算法的构成要素-粒子空间的初始化较好地选择粒子的初始化空间，将大大缩短收敛时间．初始化空间根据具体问题的不同而不同，也就是说，这是问题依赖的．从上面的介绍可以看到，粒子群算法与其他现代优化方法相比的一个明显特色就是所需调整的参数很少．相对来说，惯性因子和邻域定义较为重要．这些为数不多的关键参数的设置却对算法的精度和效率有着显著影响．惯性权重1998年，Shi和Eberhart引入了惯性权重w，并提出动态调整惯性权重以平衡收敛的全局性和收敛速度，该算法被称为标准PSO算法惯性权重w描述粒子上一代速度对当前代速度的影响。w值较大，全局寻优能力强，局部寻优能力弱；反之，则局部寻优能力强。当问题空间较大时，为了在搜索速度和搜索精度之间达到平衡，通常做法是使算法在前期有较高的全局搜索能力以得到合适的种子，而在后期有较高的局部搜索能力以提高收敛精度。所以w不宜为一个固定的常数。线性递减权值wmax最大惯性权重，wmin最小惯性权重，run当前迭代次数，runmax为算法迭代总次数较大的w有较好的全局收敛能力，较小的w则有较强的局部收敛能力。因此，随着迭代次数的增加，惯性权重w应不断减少，从而使得粒子群算法在初期具有较强的全局收敛能力，而晚期具有较强的局部收敛能力。maxmaxminmax()*run1999年，Clerc引入收缩因子以保证算法的收敛性。速度更新公式为其中，收缩因子K为受φ1φ2限制的w。φ1φ2是需要预先设定的模型参数收缩因子法控制系统行为最终收敛，且可以有效搜索不同区域，该法能得到较高质量的解。1122[()()]ididididdidvKvrpbestxrgbestx1222,,424K收缩因子法PSO研究热点PSO研究热点PSO研究热点