融合黄金正弦与sigmoid连续化的海鸥优化算法

文章目录

融合黄金正弦与sigmoid连续化的海鸥优化算法
- 1.海鸥优化算法
- 2. 改进海鸥优化算法
- - 2.1 Sigmoid 函数收敛因子
  - 2.2 全局最优更新连续化
  - 2.3 黄金正弦指引机制
- 3.实验结果
- 4.参考文献
- 5.Matlab代码
- 6.Python代码

摘要：针对海鸥算法（SOA）在迭代寻优过程中容易陷入局部最优、收敛速度慢以及寻优精度低等缺陷，提出一种黄金正弦引导与Sigmoid连续化的海鸥优化算法（GSCSOA）。在海鸥迁徙阶段，使用Sigmoid函数作为非线性收敛因子引导海鸥搜寻过程，使得算法前期保持更强的全局寻优能力，后期更快收敛。在海鸥扑食阶段，引入禁忌搜索的思想，使得海鸥始终向着置信度更高的区域移动，并且在一次迭代中最优位置持续变化，从而提高寻优精度。之后使用黄金正弦机制指引种群位置跟新，缩小搜索范围，提高局部寻优能力。

1.海鸥优化算法

基础海鸥优化算法的具体原理参考，我的博客：https://blog.csdn.net/u011835903/article/details/107535864

2. 改进海鸥优化算法

2.1 Sigmoid 函数收敛因子

Ps(x)=A×Z+Pbs(x)Z=∣Pbs(x)+2×A×rd×(Pbs(x)−Pb(x)∣×x×y×z(11)\begin{aligned} &P_s(x)=A \times Z+P_{\mathrm{bs}}(x)\\ &Z=\mid P_{\mathrm{bs}}(x)+2 \times A \times r d \times\left(P_{b s}(x)-P_b(x) \mid \times x \times y \times z\right. \end{aligned}\tag{11} Ps(x)=A×Z+Pbs(x)Z=∣Pbs(x)+2×A×rd×(Pbs(x)−Pb(x)∣×x×y×z(11)

式(11)由式(10)转换而来, 由式(11)可知：当收敛因子 AAA 的值较大时, 全局搜索能力强, 当 AAA 的值小时, 局部寻优能力强。SOA算法使用线性惯性权重, AAA 的值随着迭代次数增加从 2 线性下降至 0 。这种线性惯性权重虽然能在一定程度上平衡全局搜索和局部搜索, 但是实际的搜索过程是一个极为复杂的非线性过程, 线性权重会降低算法寻优的性能。
本文提出的 GSCSOA 算法使用 Sigmoid 函数作为非线性收敛因子 AfA_fAf, 替换原有的收敛因子 A∘AfA \circ A_fA∘Af 的表达式如式(14) 所示, 首先将 xxx 缩放到 (−5,5),Af(-5,5), A_f(−5,5),Af 从接近 2 的值非线性下降至接近 0 的值, 改进后的位置更新由式(12)计算。
Ps(x)=Af×Z+Pbs(x)(12)S(x)=21+e−x(13)Af(x)=fc×S(−x−12Maxitea Max⁡itea ×10)=fc1+exp⁡(x−12Maxitea Mitea ×10)(14)\begin{gathered} P_s(x)=A_f \times Z+P_{\mathrm{bs}}(x)&(12)\\ S(x)=\frac{2}{1+e^{-x}}&(13)\\ A_f(x)=f_c \times S\left(-\frac{x-\frac{1}{2} M a x_{\text {itea }}}{\operatorname{Max}_{\text {itea }}} \times 10\right)=\frac{f_c}{1+\exp \left(\frac{x-\frac{1}{2} M a x_{\text {itea }}}{M_{\text {itea }}} \times 10\right)}&(14) \end{gathered} Ps(x)=Af×Z+Pbs(x)S(x)=1+e−x2Af(x)=fc×S(−Maxitea x−21Maxitea ×10)=1+exp(Mitea x−21Maxitea ×10)fc(12)(13)(14)
本文使用 Sigmoid 函数作为收敛因子，在前期保持种群较大的个体自由度，增强全局寻优能力。在后期的个体自由度迅速下降，局部寻优能力加强。相比于线性惯性权重，Sigmoid 函数作为非线性收敛因子来平滑过度迁徙和扑食过程，能更优地平衡全局搜索和局部搜索，使寻优过程变成非线性。因此这种改在进理论上可以提高种群寻优精度，同时加速寻优速度。

2.2 全局最优更新连续化

SOA 算法每次迭代中，对于种群更新后的位置全部接受，迭代完成后种群才进行全局最优查找。这样会存在两个缺点：一是种群多样性过早下降。由式(11)可知，SOA 算法更新后的种群分布以上次搜索的全局最优为中心，呈螺旋状形态。随着迭代次数的增加，种群越接近中心，个体多样性锐减。二是全局最优更新非连续。在一次迭代中，全局最优在每个搜索代理位置更新后都有可能产生，SOA 算法是在种群位置更新完成后才更新全局最优位置，这种更新会导致 P s (x)不连续。使的本次迭代与下一次迭代之间不连续，非连续化的更新使得最新产生的全局最优不能及时的用来指导后面的搜索，降低了算法的效率和收敛速度。

GSCSOA 算法引入贪心策略(15)对每个搜索代理的位置进行一次篮选, 引导种群向最优位置靠近。这一策略一方面保留了较大适应度的搜索代理, 使得平均适应度较高, 由式 (1)知搜索代理自身的经验也在寻优中发挥作用, 因此有利于提高算法整体的寻优性能。另一方面, 抛弃了本次迭代中的一部分个体, 保留了历代搜索中的一部分较高适应度的个体, 使得每次迭代计算的搜索代理来自不同的搜索历史, 从而种群的多样性得到丰富, 个体分布较为分散, 避免陷入局部收敛, 解决了 SOA 遇到的第一个问题。
为解决第二个问题, GSCSOA 在每个搜索代理位置更新后使用式(16)更新全局最优，在每次迭代中全局最优可以多次变化。及时地使用最新全局最优来指导后面的更新, 避免了无效的计算, 提升算法效率。这种机制使 Pbs(x)P_{b s}(x)Pbs(x) 的更新相对连续化, 每次迭代之间更加平滑地过度。
Pbs(x+1)={Ps(x+1),if fit (Ps(x+1))≥fit(Pbs(x))Pbs(x),else (16)P_{\mathrm{bs}}(x+1)=\left\{\begin{array}{l} P_{\mathrm{s}}(x+1), \text { if fit }\left(P_{\mathrm{s}}(x+1)\right) \geq f i t\left(P_{b \mathrm{~s}}(x)\right) \\ P_{\mathrm{bs}}(x), \text { else } \end{array}\right. \tag{16} Pbs(x+1)={Ps(x+1), if fit (Ps(x+1))≥fit(Pb s(x))Pbs(x), else (16)

2.3 黄金正弦指引机制

在 GSCSOA 算法中, 通过融合搜索代理自身的历史经验和种群的全局最优经验对搜索代理的位置进行二次更新, 如式(17)所示。
Ps(x)=xit∗∣sin⁡(r1)∣+r2∗sin⁡(r1)∗∣θ1∗Pbs(x)−θ2∗Ps(x−1)∣(17)\begin{gathered} P_s(x)=x_i^t *\left|\sin \left(r_1\right)\right|+r_2 * \sin \left(r_1\right)^* \\ \left|\theta_1 * P_{b s}(x)-\theta_2 * P_s(x-1)\right| \end{gathered}\tag{17} Ps(x)=xit∗∣sin(r1)∣+r2∗sin(r1)∗∣θ1∗Pbs(x)−θ2∗Ps(x−1)∣(17)
其中, r1r_1r1 是 (0,2Π)(0,2 \Pi)(0,2Π) 中的一个随机数, r2r_2r2 是 (0,Π)(0, \Pi)(0,Π) 中的一个随机数。 θ1\theta_1θ1 和 θ2\theta_2θ2 是黄金正弦系数, 由式(18)获得, τ\tauτ 是黄金分割比例, a、ba 、 ba、b 是参数, 初始值分别为 Π\PiΠ 和 Π\PiΠ, 通过黄金分割算法 [15]更新。
θ1=a∗(1−τ)+b∗τ,θ2=a∗τ+b∗(1−τ)(18)τ=1−52=0.618033..(19)\begin{gathered} \theta_1=a^*(1-\tau)+b^* \tau, \theta_2=a^* \tau+b^*(1-\tau) &(18)\\ \tau=\frac{1-\sqrt{5}}{2}=0.618033 . .&(19) \end{gathered} θ1=a∗(1−τ)+b∗τ,θ2=a∗τ+b∗(1−τ)τ=21−5=0.618033..(18)(19)
黄金正弦指引机制通划分种群范围得到最佳搜索区域, 将其用于指导搜索代理的寻优过程, 使得 GSCSOA 算法在两个方面有较大改善。一是海鸥的迁徙和和扑食操作都在最优范围内进行, 搜索代理更容易寻找到全局最优位置, 从而提高种群搜索精度。二是搜索范围以黄金比例缩小, 种群迅速最佳区域靠近, 从而提高收玫速度。

结合以上改进方法, 本文改进的海鸥算法伪代码如下:
算法 1GSCSOA1 \mathrm{GSCSOA}1GSCSOA
输入：海鸥种群 Ps(s=1,2…N)P_s(s=1,2 \ldots N)Ps(s=1,2…N) 。
输出：全局最优 PbsP_{b s}Pbs 。
初始化算法超参数 fcf cfc, 根据种群规模 NNN, 最大迭代次数 Tmax T_{\text {max }}Tmax , 搜索空
间维度 ddd 初始化种群;
计算种群每个搜索代理的适应度, 寻找种群最优位置 Pbs(x)P_{b s}(x)Pbs(x);
While (t<Tmax⁡)//Tmax⁡\left(t<T_{\max }\right) \quad / / T_{\max }(t<Tmax)//Tmax 为最大迭代次数
For in range (N)//N(N) \quad / / N(N)//N 为搜索代理数量
使用式 (14) 计算 Sigmoid 收敛因子 AfA_fAf, 并使用式 (12)更新搜索代理 sss 的位置 Ps(x)P_s(x)Ps(x);
使用式 (15)(16)更新搜索代理 sss 的位置 Ps(x)P_s(x)Ps(x) 和全局最优位置 Pbs(x)P_{b s}(x)Pbs(x);
End for
For in range (N)(N)(N)
按式 (17)对所有搜索代理进行位置二次更新; ///黄金正弦引导
End for
End while
输出全局最优 Pbs(x)P_{b s}(x)Pbs(x);

3.实验结果

4.参考文献

[1]王宁,何庆.融合黄金正弦与sigmoid连续化的海鸥优化算法[J/OL].计算机应用研究:1-8[2021-08-06].https://doi.org/10.19734/j.issn.1001-3695.2021.05.0176.

5.Matlab代码

6.Python代码