粒子群算法（带约束处理）——Python&Matlab实现

粒子群算法思想来源于实际生活中鸟捕食的过程。假设在一个n维的空间中，有一群鸟（m只）在捕食，食物位于n维空间的某个点上，对于第i只鸟某一时刻来说，有两个向量描述，一个是鸟的位置向量，第二个是鸟的速度。假设鸟能够判断一个位置的好坏，所谓“好坏”，就是离食物更近了还是更远了。鸟在捕食的过程中会根据自己的经验以及鸟群中的其他鸟的位置决定自己的速度，根据当前的位置和速度，可以得到下一刻的位置，这样每只鸟通过向自己和鸟群学习不断的更新自己的速度位置，最终找到食物，或者离食物足够近的点。

粒子群是比较经典的有自适应过程的算法，类似的有蝙蝠算法、布谷鸟算法、蜂群算法等，某个粒子的移动会参考历史最优和当前最优，可以通过设置c1、c2表示两者的重要程度，每个粒子移动的速度也是不一样的，速度也会收到w的影响，以上三者共同决定粒子下一步到达的位置。

2.1.2 更新规则

PSO初始化为一群随机粒子(随机解)。然后通过迭代找到最优解。在每一次的迭代中，粒子通过跟踪两个“极值”(pbest，gbest)来更新自己。在找到这两个最优值后，粒子通过下面的公式来更新自己的速度和位置。

其中：f为当前适应度函数的平均值；v为粒子速度；x为粒子位置；rand()为介于（0，1）之间的随机数；c1、c2是学习因子，通常c1=c2=2；best为当前最优解；bestx为历史最优解。

公式(1)的第一部分称为**【记忆项】，表示上次速度大小和方向的影响；公式(1)的第二部分称为【自身认知项】，是从当前点指向粒子自身最好点的一个矢量，表示粒子的动作来源于自己经验的部分；公式(1)的第三部分称为【群体认知项】，是一个从当前点指向种群最好点的矢量，反映了粒子间的协同合作和知识共享。粒子就是通过自己的经验和同伴中最好的经验来决定下一步的运动。以上面两个公式为基础，形成了PSO的标准形式。**

公式(2)和 公式(3)被视为标准PSO算法

2.1.3 惩罚项

惩罚项。它的思想类似线性规划内点法，都是通过增加罚函数，迫使模型在迭代计算的过程中始终在可行域内寻优。

在粒子群算法中，每一步迭代都会更新 Pbest 和 Gbest，虽然可以将有约束问题转换为无约束问题进行迭代求解，但是问题的解 xi依然存在不满足约束条件的情况，因此需要编制一些规则来比较两

个粒子的优劣，规则如下：

1.如果两个粒子 xi 和 xj 都可行，则比较其适应度函数 f(xi)和f(xj)，值小的粒子为优。

2.当两个粒子 xi 和 xj 都不可行，则比较惩罚项 P(xi)和 P(xj）， 违背约束程度小的粒子更优。

3.当粒子 xi 可行而粒子 xj 不可行，选可行解。对于粒子的上下限约束 可以体现在位置更新函数里，不必加惩 罚项。 具体思路就是遍历每一个粒子的位置，如果超除上下限，位置则更改为上下限中的任何一个位置

2.2 程序框图

---

2.2.1 伪代码

2.2.2 程序框图