智能优化算法：纵横交叉算法-附代码

智能优化算法：纵横交叉算法-附代码

摘要：纵横交叉算法(Crisscross optimization algorithm，CSO) 是一种全新的基于种群的随机搜索算法。有两种不同的交叉方式，能较好地解决一般智能算法存在的局部最优问题。

1.算法原理

纵横交叉算法在迭代过程中每一代都会进行横向交叉和纵向交叉两种交叉方式，从而使得种群中某些陷入维局部最优的维有机会跳出迭代。纵横交叉算法每次交叉后得出的解称为中庸解 ( $M{S}_{hc},M{S}_{vc}$ ) ，通过引入竞争算子，使得这两种交叉方式有机的结合起来：每次交叉操作之后都会进入竞争算子，与父代进行竞争，只有比父代更优秀的粒子会被保留下来进入下次迭代，得出的解称为占优解 ( $D{S}_{hc},D{S}_{vc}$) 。

1.1 横向交叉操作

横向交叉类似于 GA 的交叉操作，是在种群中两个不同个体粒子相同维之间进行的一种算数交叉。假设父代个体粒子 $X(i)$ 和 $X(j)$ 的第 $d$ 维进行横向交叉，则它们产生子代的公式如下:
$$M{S}_{hc}(i,d)=r_1\ast X(i,d)+(1-r_1)\ast X(j,d)+c_1\ast (X(i,d)-X(j,d))\text{\hspace{1em}(1)}$$

$$M{S}_{hc}(j,d)=r_1\ast X(j,d)+(1-r_1)\ast X(i,d)+c_1\ast (X(j,d)-X(i,d))\text{\hspace{1em}(2)}$$

式中： $r_1,r_2$ 是 [0, 1] 之间的随机数； $c_1,c_2$ 是 [ - 1, 1] 之间的随机数； $X(i,d),X(j,d)$分别是父代种群中个体粒子 $X(i)$ 和 $X(j)$ 的第 $d$维； $M{S}_{hc}(i,d)$和 $M{S}_{hc}(j,d)$ 分别是$X(i,d)$ 和 $X(j,d)$ 通过横向交叉产生的第 $d$ 维子代。

不同个体粒子的相同维度的上下限范围是相同的，从社会学的角度来看，式 (1) 中的 $r_1\ast X(i,d)$是粒子 $X(i)$ 的记忆项，是粒子本身的当前最优值。$(1-r_1)\ast X(j,d)$ 是粒子 $X(i)$ 和 $X(j)$ 的群体认知项，表示了不同粒子间的相互影响。这两项通过惯性权重因子$r_1$ 较好的结合在一起。 $c_1$ 是学习因子，第三项$c_1\ast (X(i,d)-X(j,d))$ 可以增大搜索区间，在边缘寻优。横向交叉操作完成后，得到的中庸解 $M{S}_{hc}(i,d),M{S}_{hc}(j,d)$ 必须分别与父代粒子 $X(i),X(j)$ 的适应度比较，只有适应度更好的中庸解才可以保留下来，成为占优解 $D{S}_{hc}$ ，参与下一次迭代。

1.2 纵向交叉操作

纵向交叉是种群中一个粒子的两个不同的维之间进行的一种算数交叉。由于不同维元素的取值范围不同，因而交叉前须将两维进行归一化处理，并且为了使已陷入维局部最优而停滞的那一维跳出局部最优而又不破坏另外一维的信息，每次纵向交叉操作只产生一个子代粒子，只对其中一维进行更新。

假定粒子 $X(i)$ 的第 $d_1$ 维和第 $d_2$ 维是参与纵向交叉，根据公式 (3) 产生中庸解 $MSvc(i,d_1)$ 。
$$M{S}_{vc}(i,d_1)=r\ast X(i,d_1)+(1-r)\ast X(i,d_2),i\in N(1,M),d_1,d_2\in N(1,D)\text{\hspace{1em}(3)}$$
式中： $r\in [0,1]$ ； $M{S}_{vc}(i,d_1)$ 是个体粒子 $X(i)$ 的第 $d_1$ 维和第 $d_2$ 维通过纵向交叉产生的第 $d_1$ 维后代。

式 (3) 中的第一项是粒子 $X(i)$ 的第 $d_1$ 维的记忆项，第二项是粒子 $X(i)$ 的第 $d_1$ 维和第 $d_2$ 维相互影响，通过惯性权重因子 $r$ 结合在一起。这样得到的中庸解$M{S}_{vc}(i,d_1)$不但包含父代粒子 $X(i)$ 的第 $d_1$ 维的信息，还依一定概率含有 $X(i)$ 的第 $d_2$ 维信息，并不会破坏 $X(i)$的第 $d_2$ 维信息。中庸解 $M{S}_{v}c(i,d_1)$ 与父代粒子$X(i)$ 比较适应度，较好的保留下来作为占优解 $D{S}_{vc}$ ，进行下一次迭代。

CSO 的操作流程可总结如下：
1）根据具体的问题初始化种群
2）执行横向交叉，然后执行竞争算子，得到 DSh 作为纵向交叉的父代
3）执行纵向交叉，然后执行竞争算子，得到 DSv 作为下一次横向交叉的父代
4）判断是否达到最大迭代次数或者指定的结束条件，如果是，则终止算法，否，
则转到步骤 2）。

2.算法结果

3.参考文献

[1] 孟安波, 胡函武, 刘向东. 基于纵横交叉算法优化神经网络的负荷预测模型[J]. 电力系统保护与控制, 2016, 44(07):102-106.

4.Matlab代码

5.Python代码

个人资料介绍