【论文复现|智能算法改进】融合黑寡妇思想的蜣螂优化算法
目录
- 1.算法原理
- 2.改进点
- 3.结果展示
- 4.参考文献
- 5.代码获取
1.算法原理
【智能算法】蜣螂优化算法(DBO)原理及实现
2.改进点
ICMIC混沌映射
z n + 1 = sin ( α z n ) , α ∈ ( 0 , + ∞ ) (1) z_{n+1}=\sin(\frac{\alpha}{z_n}),\alpha\in(0,+\infty)\tag{1} zn+1=sin(znα),α∈(0,+∞)(1)
α 是映射系数, 为了进一步探索映射参数的影响, 采用李雅普诺(lyapunov)指数(λ)来衡量映射效果。李雅普诺夫指数表征了一个系统随着时间运动发生改变的程度,当指数小于 0,说明运动系统体积缩小,系统趋于稳定;而当指数大于 0,则说明运动系统体积增大,系统趋于混沌。
λ = lim n → ∞ 1 n ∑ n = 0 n − 1 ln ∣ d f ( z n , α ) d z ∣ (2) \lambda=\lim_{n\to\infty}\frac{1}{n}\sum_{n=0}^{n-1}\ln|\frac{df(z_n,\alpha)}{dz}|\tag{2} λ=n→∞limn1n=0∑n−1ln∣dzdf(zn,α)∣(2)
论文中给出李雅普诺夫指数图
复现一下这张图
黑寡妇算法
论文提出使用黑寡妇算法替换滚球行为算法。 黑寡妇算法将黑寡妇蜘蛛在蛛网中的运动模型转化为线性和螺旋两种形式,位置更新公式:
x i ( t + 1 ) = { x ∗ ( t ) − m x r 1 ( t ) if α ≤ 0.3 x ∗ ( t ) − cos ( 2 π β ) x i ( t ) α > 0.3 (3) x_i(t+1)=\begin{cases}x^*(t)-mx_{r_1}(t)&\text{if}~\alpha\leq0.3\\x^*(t)-\cos(2\pi\beta)x_i(t)&\alpha>0.3\end{cases}\tag{3} xi(t+1)={x∗(t)−mxr1(t)x∗(t)−cos(2πβ)xi(t)if α≤0.3α>0.3(3)
信息素(pheromone)是黑寡妇种群的另一个特质,在交配的过程中起到了重要的作用。黑寡妇算法将信息素定义为:
p h e r o m o n e ( i ) = f i t n e s s max − f i t n e s s ( i ) f i t n e s s max − f i t n e s s min (4) pheromone(i)=\frac{fitness_{\max}-fitness(i)}{fitness_{\max}-fitness_{\min}}\tag{4} pheromone(i)=fitnessmax−fitnessminfitnessmax−fitness(i)(4)
当信息素过小时,个体将会被新个体替代,位置更新:
x i ( t ) = x ∗ ( t ) + 1 2 [ x r 1 ( t ) − ( − 1 ) σ ∗ x r 2 ( t ) ] (5) x_i(t)=x^*(t)+\frac{1}{2}[x_{r_1}(t)-(-1)^\sigma*x_{r_2}(t)]\tag{5} xi(t)=x∗(t)+21[xr1(t)−(−1)σ∗xr2(t)](5)
自适应 T 分布改进因子
采用以迭代次数自由度参数的t 分布变异扰动因子,对小蜣螂觅食行为进行扰动。t分布即是学生分布,决策参数为自由度参数 m,其概率密度函数为:
p ( x ) = Γ ( m + 1 2 ) m π ∗ Γ ( m 2 ) ∗ ( 1 + x 2 m ) − m + 1 2 , − ∞ < x , ∞ (6) p(x)=\frac{\Gamma(\frac{m+1}{2})}{\sqrt{m\pi} ^*\Gamma(\frac{m}{2})}*(1+\frac{x^2}{m})^{^{-\frac{m+1}{2}}},-\infty<x,\infty \tag{6} p(x)=mπ∗Γ(2m)Γ(2m+1)∗(1+mx2)−2m+1,−∞<x,∞(6)
当 m=1 时, t 分布即是柯西分布, m 趋于无穷, t 分布即是高斯分布。因此自适应 t 分布变异扰动因子具有柯西、高斯分布两种特性,使得蜣螂算法在迭代前期具有较好的全局开发能力,在迭代后期具有良好的局部探索能力,并提高算法的收敛速度。位置更新:
x i ( t ) = X b + t { exp [ 4 ∗ ( t M ) 2 ] } ∗ X b (7) x_i(t)=X^b+t\{\exp[4^*(\frac{t}{M})^2]\}^*X^b\tag{7} xi(t)=Xb+t{exp[4∗(Mt)2]}∗Xb(7)
流程图
3.结果展示
采用CEC2005测试集测试对GWO,SSA,DBO,BWDBO算法进行测试对比:
CEC2005-F1
CEC2005-F1 箱型图
CEC2005-F1 Friedamn排名
**CEC2005-F2 **
**CEC2005-F3 **
**CEC2005-F3 **
从适应度曲线,箱型图,Friedamn排名来看,改进蜣螂优化算法DWBDO效果不错😊
4.参考文献
[1] 孔令崧,石颉,孙浩,等.融合黑寡妇思想的蜣螂优化算法[J/OL].微电子学与计算机,1-15[2024-06-01].
5.代码获取
智能算法,论文复现,算法应用(机器学习、二维&三维路径规划、UAV路径规划、布局优化、调度优化、VRP问题等),定制算法可以联系我~
资源清单:https://docs.qq.com/sheet/DU1V0QWtSeHJnY0JU?u=989d90f9b14449ec9419aa7b51473c1b&tab=BB08J2