蜣螂优化算法（DBO）详解：从原理到实战

蜣螂优化算法（DBO）详解：从原理到实战

1 算法概述与生物基础

蜣螂优化算法（Dung Beetle Optimizer, DBO）是东华大学沈波教授团队于2022年提出的一种新型群体智能优化算法。该算法模拟了蜣螂（俗称屎壳郎）在自然界中的滚球、跳舞、繁殖、觅食和偷窃等行为，通过智能优化机制在解空间中进行高效搜索。

1.1 生物行为基础

蜣螂作为一种奇特的昆虫，其生活习性包含多种智能行为：

• 滚球行为：蜣螂将粪便滚成球状，并利用天体线索（太阳、月亮和偏振光）进行导航，以直线方式将粪球滚到安全地点。
• 跳舞行为：当遇到障碍物时，蜣螂会爬到粪球上方"跳舞"（旋转和停顿），重新确定前进方向。
• 繁殖行为：蜣螂将粪球滚到安全地点藏起来，雌性蜣螂在粪球内产卵，为后代提供营养和保护。
• 觅食行为：成熟的小蜣螂从地下钻出来寻找食物，它们会在最佳觅食区域活动。
• 偷窃行为：部分蜣螂会偷取其他蜣螂的粪球，这种行为增加了种群的多样性。

1.2 算法基本思想

DBO算法将优化问题的搜索过程映射为蜣螂群体的智能行为：

• 每个蜣螂的位置代表一个候选解
• 粪球的质量对应解的适应度值
• 导航行为引导群体向最优解区域移动
• 五种行为的协调平衡了全局探索和局部开发能力

DBO算法与其他群体智能算法（如粒子群优化、遗传算法）相比，具有参数少、结构简单、收敛速度快等优点，在各类优化问题上表现出色。

2 数学模型与算法原理

2.1 滚球蜣螂的数学模型

滚球行为是DBO算法的核心探索机制。当蜣螂无障碍物时，其位置更新公式为：

x i ( t + 1 ) = x i ( t ) + α × k × x i ( t − 1 ) + b × Δ x Δ x = ∣ x i ( t ) − X w ∣ \begin{aligned} x_i(t+1) &= x_i(t) + \alpha \times k \times x_i(t-1) + b \times \Delta x \\ \Delta x &= |x_i(t) - X^w| \end{aligned}xi​(t+1)Δx​=xi​(t)+α×k×xi​(t−1)+b×Δx=∣xi​(t)−Xw∣​

其中：

• x i ( t ) x_i(t)xi​(t)是第i ii只蜣螂在第t tt次迭代的位置
• α \alphaα是自然系数（取值为1或-1），模拟自然因素导致的偏离
• k kk是偏转系数（取值范围为( 0 , 0.2 ] (0, 0.2](0,0.2]）
• b bb是常数（取值范围为( 0 , 1 ) (0, 1)(0,1)）
• X w X^wXw是全局最差位置，Δ x \Delta xΔx模拟光强变化

表1：滚球行为参数说明

2.2 跳舞行为数学模型

当蜣螂遇到障碍物（以概率0.1判断）时，通过跳舞重新定向：

x i ( t + 1 ) = x i ( t ) + tan ⁡ ( θ ) × ∣ x i ( t ) − x i ( t − 1 ) ∣ x_i(t+1) = x_i(t) + \tan(\theta) \times |x_i(t) - x_i(t-1)|xi​(t+1)=xi​(t)+tan(θ)×∣xi​(t)−xi​(t−1)∣

其中θ ∈ [ 0 , π ] \theta \in [0, \pi]θ∈[0,π]是偏转角。当θ \thetaθ等于0 00,π / 2 \pi/2π/2, 或π \piπ时，位置不更新。正切函数提供了方向变化的敏感性，使蜣螂能够有效避开障碍物。

2.3 繁殖行为数学模型

繁殖行为模拟蜣螂为后代选择安全产卵区域的过程：

L b ∗ = max ⁡ ( X ∗ × ( 1 − R ) , L b ) U b ∗ = min ⁡ ( X ∗ × ( 1 + R ) , U b ) B i ( t + 1 ) = X ∗ + b 1 × ( B i ( t ) − L b ∗ ) + b 2 × ( B i ( t ) − U b ∗ ) \begin{aligned} Lb^* &= \max(X^* \times (1-R), Lb) \\ Ub^* &= \min(X^* \times (1+R), Ub) \\ B_i(t+1) &= X^* + b_1 \times (B_i(t) - Lb^*) + b_2 \times (B_i(t) - Ub^*) \end{aligned}Lb∗Ub∗Bi​(t+1)​=max(X∗×(1−R),Lb)=min(X∗×(1+R),Ub)=X∗+b1​×(Bi​(t)−Lb∗)+b2​×(Bi​(t)−Ub∗)​

其中：

• X ∗ X^*X∗是当前局部最优位置
• R = 1 − t / T max ⁡ R = 1 - t/T_{\max}R=1−t/Tmax​是动态收敛因子
• L b ∗ Lb^*Lb∗和U b ∗ Ub^*Ub∗定义产卵区域的边界
• b 1 b_1b1​和