智能优化算法：海鸥算法原理及Matlab代码

文章导航

海鸥算法(Seagull optimization algorithm,SOA)，是由Gaurav Dhiman和Vijay kumar于2018年提出的一种新型群体智能优化算法，该算法主要模拟了自然界中海鸥迁徙以及迁徙过程中的攻击行为(觅食行为)。

1. 算法原理

海鸥算法主要模拟了海鸥的迁徙行为和攻击行为 。迁徙行为即海鸥从一个现阶段不适宜生存的地方飞往另一个适宜生存的地方，迁徙行为影响着SOA算法的全局探索能力；攻击行为即海鸥在飞行过程中对地面、水域内食物的攻击觅食，攻击行为影响着SOA算法的局部开发能力。

2. SOA算法流程

2.1 迁徙行为(exploration ability)

根据原始文献所述，海鸥在迁徙行为即算法中个体从一个位置探索至另一个位置。在这个阶段，需要注意三点：避免碰撞、向最优个体靠拢、与最优代理保持密切联系。

为了避免与周围海鸥发生碰撞，算法采用附加变量A调整海鸥的位置。

$$C_s=A\ast P_s(x)$$其中A代表海鸥在给定搜索空间中的迁徙行为，A的大小通过$f_c$进行控制。。$$A=f_c-(x\ast f_c/Ma{x}_{iteration})$$最终A的大小根据迭代次数，从2→0呈线性下降。

在保证海鸥个体之间不会发生碰撞之后，令所有海鸥向着最佳海鸥靠拢。$$M_s=B\ast (P_{best}(x)-P_s(x))$$其中，$M_s$表示个体向最佳海鸥的收敛方向，$B$作为平衡算法探索和开发能力的重要算参数，其变化遵循：$$B=2\ast A^2\ast rand$$

在计算出每个海鸥的收敛方向后，海鸥开始向着这个位置进行移动：
$$D_s=|C_s+M_s|$$其中$D_s$即海鸥的新位置。

2.2 攻击行为(exploitation ability)

海鸥在觅食时产生的攻击行为呈螺旋状下降，通过不断改变角度和半径在空中进行移动。海鸥在三维空间中的具体位置为：
$$x=r\ast cos(k)$$$$y=r\ast sin(k)$$$$z=r\ast k$$$$r=u\ast e^{kv}$$其中，$k$是一个位于$[0,2\pi ]$的随机数，算法通过$u$和$v$控制螺旋半径$r$，$u$和$v$通常取1。

结合海鸥新位置，得到整体海鸥位置更新公式如下：
$$P_s(x)=(D_s\ast x\ast y\ast z)+P_{best}(x)$$

3. SOA算法流程如下：

step1：初始化参数
step2：计算每只海鸥的适应度值，即目标函数值
step3：根据上述公式(1-5)计算$D_s$
step4：根据上述公式(6-10)计算$P_s(x)$
step5：更新最佳海鸥的位置信息和适应度值，$ineration=iteration+1$
step6：如果$interation>Ma{x}_{iteration}$，跳转至step7，否则跳转至step3
step7：输出最佳海鸥位置和适应度值

4. 运行结果

以CEC标准测试函数——单峰函数、多峰函数进行2维和30维测试，测试结果如下图所示。


测试函数为：30维下的Sphere函数，所测结果可以看出SOA算法具有非常优秀的寻优能力！

具体Matlab代码链接如下：
https://mianbaoduo.com/o/bread/aZqamJc=

       有网友私聊我，想让我分享一下海鸥算法的改进与创新，那么他来了。大家可以去看我的第三篇文章——基于量子的鸽群优化算法，将两者结合可以发表不错的文章哦。文章链接https://blog.csdn.net/qq_38643813/article/details/109730959