粒子群优化算法-详细易懂-很多例子-

粒子群优化算法（PS0）ParticleSwarmOptimization智能算法向大自然学习遗传算法(GA)物竞天择，设计染色体编码，根据适应值函数进行染色体选择、交叉和变异操作，优化求解人工神经网络算法(ANN)模仿生物神经元，透过神经元的信息传递、训练学习、联想，优化求解模拟退火算法(SA)模模仿金属物质退火过程�解决最优化问题的方法传统搜索方法保证能找到最优解HeuristicSearch不能保证找到最优解�由Kennedy和Eberhart于1995年提出．群体迭代，粒子在解空间追随最优的粒子进行搜索.简单易行粒子群算法:收敛速度快设置参数少已成为现代优化方法领域研究的热点．粒子群算法发展历史简介�粒子群算法的基本思想粒子群算法的思想源于对鸟群捕食行为的研究．模拟鸟集群飞行觅食的行为，鸟之间通过集体的协作使群体达到最优目的，是一种基于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—惯性权重，非负数，调节对解空间的搜索范围kk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx11kkkidididxxvkidvkidx�粒子速度更新公式包含三部分：第一部分为粒子先前的速度第二部分为“认知”部分，表示粒子本身的思考，可理解为粒子i当前位置与自己最好位置之间的距离。第三部分为“社会”部分，表示粒子间的信息共享与合作，可理解为粒子i当前位置与群体最好位置之间的距离。kk-111idid1122v=wv()()kkididdidcrpbestxcrgbestx�區域最佳解全域最佳解運動向量慣性向量12X=X,X,...,Xiiiid12V=V,V,...,ViiiidStudyFactor12(1)()()()()()()()ididididgdidttttvwvcrandpxcrandpx(1)()()iiitttxxvpg�1111122V=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的优化能力很好，�粒子群算法的构成要素-权重因子权重因子：惯性因子、学习因子1c2c失去对粒子本身的速度的记忆1社会经验部分前次迭代中自身的速度自我认知部分基本粒子群算法粒子的速度更新主要由三部分组成：惯性因子kv�粒子群算法的构成要素-权重因子权重因子：惯性因子、学习因子1c2c社会经验部分前次迭代中自身的速度自我认知部分粒子的速度更新主要由三部分组成：kv学习因子无私型粒子群算法“只有社会，没有自我”迅速丧失群体多样性，易陷入局优而无法跳出．�粒子群算法的构成要素-权重因子权重因子：惯性因子、学习因子1c2c社会经验部分前次迭代中自身的速度自我认知部分粒子的速度更新主要由三部分组成：kv自我认知型粒子群算法“只有自我，没有社会”完全没有信息的社会共享，导致算法收敛速度缓慢学习因子�粒子群算法的构成要素-权重因子权重因子：惯性因子、学习因子1c2c社会经验部分前次迭代中自身的速度自我认知部分粒子的速度更新主要由三部分组成：kvc1,c2都不为0，称为完全型粒子群算法完全型粒子群算法更容易保持收敛速度和搜索效果的均衡，是较好的选择．�粒子群算法的构成要素-最大速度但在于维护算法的探索能力与开发能力的平衡．Vm较大时，探索能力增强，mV作用:Vm较小时，开发能力增强，mVmVVm一般设为每维变量变化范围的10％~20％.但粒子容易飞过最优解．容易陷入局部最优．�粒子群算法的构成要素-邻域的拓扑结构全局粒子群算法和局部粒子群算法．gp粒子群算法的邻域拓扑结构包括两种，一种是将群体内所有个体都作为粒子的邻域，另一种是只将群体中的部分个体作为粒子的邻域．群体历史最优位置邻域拓扑结构决定由此，将粒子群算法分为�粒子群算法的构成要素-邻域的拓扑结构全局粒子群算法1.粒子自己历史最优值2.粒子群体的全局最优值局部粒子群算法1.粒子自己历史最优值2.粒子邻域内粒子的最优值邻域随迭代次数的增加线性变大，最后邻域扩展到整个粒子群。经过实践证明：全局版本的粒子群算法收敛速度快，但是容易陷入局部最优。局部版本的粒子群算法收敛速度慢，但是很难陷入局部最优。现在的粒子群算法大都在收敛速度与摆脱局部最优这两个方面下功夫。其实这两个方面是矛盾的。看如何更好的折中了。�粒子群算法的构成要素-停止准则停止准则一般有如下两种：最大迭代步数可接受的满意解�粒子群算法的构成要素-粒子空间的初始化较好地选择粒子的初始化空间，将大大缩短收敛时间．初始化空间根据具体问题的不同而不同，也就是说，这是问题依赖的．从上面的介绍可以看到，粒子群算法与其他现代优化方法相比的一个明显特色就是所需调整的参数很少．相对来说，惯性因子和邻域定义较为重要．这些为数不多的关键参数的设置却对算法的精度和效率有着显著影响．�第九讲daili粒子群算法293.粒子群算法示例例求解如下四维Rosenbrock函数的优化问题．322211min()[100()(1)]iiiifxxxx种群大小：[30,30](1,2,3,4)ixi解算法的相关设计分析如下．编码：因为问题的维数是4，所以每个粒子的位置和max60V即算法中粒子的数量，取5m速度均4维的实数向量．设定粒子的最大速度：�第九讲daili粒子群算法30初始位置：0ix设各粒子的初始位置和初始速度为：0iv对粒子群进行随机初始化包括随机初始化各粒子的位置和速度(0)1{21.721,9.13677,6.62244,3.84079}x(0)3{29.6563,0.871811,27.8912,17.7425}x(0)2{13.5001,23.6131,17.4462,29.0515}x(0)5{28.0992,22.6482,0.675616,8.43752}x(0)4{23.6218,16.4885,22.7019,25.4033}x�第九讲daili粒子群算法31初始速度：0ix设各粒子的初始位置和初始速度为：0iv对粒子群进行随机初始化包括随机初始化各粒子的位置和速度(0)1{19.9048,29.562,22.104,5.45346}v(0)3{7.83576,55.7173,40.9177,28.255}v(0)2{20.5922,28.6944,26.3216,19.0615}v(0)5{17.561,13.5365,51.2722,56.098}v(0)4{11.6373,41.0138,17.7311,14.87}v�第九讲daili粒子群算法32初始速度：(0)(0)(0)(0)(0)12345,,,,vvvvv初始位置：(0)(0)(0)(0)(0)12345,,,,xxxxx(0)71()2.3881710fx计算每个粒子的适应值(0)72()4.4530610fx322211()[100()(1)]iiiifxxxx按照计算适应值(0)1gpx(0)75()8.5067410fx(0)74()6.5688810fx(0)83()1.3537610fx历史最优解0,(1,2,3,4,5)iiipx�第九讲daili粒子群算法33更新粒子的速度和位置：12()2(),kkkkgkvvpxpx11kkkxxv122cc01c取，,得到速度和位置的更新函数为初始速度：(0)(0)(0)(0)(0)12345,,,,vvvvv初始位置：(0)(0)(0)(0)(0)12345,,,,xxxxx群体历史最优解：(0)1gpx0,(1,2,3,4,5)iiipx个体历史最优解：�第九讲daili粒子群算法34(1)1{19.9048,29.562,22.104,5.45346}v(1)2{40.0498,3.76972,44.9573,75.6939}v更新速度，得：(1)3{14.8665,59.3694,25.667,22.1122}v初始速度：(0)(0)(0)(0)(0)12345,,,,vvvvv初始位置：(0)(0)(0)(0)(0)12345,,,,xxxxx群体历史最优解：(0)1g