旅鼠优化算法(ALA)-2025年SCI顶刊新算法—附MATLAB完整版代码

智能信息时代的到来见证了跨各个学科的复杂优化问题的激增。尽管现有的元启发式算法在许多场景中都表现出了功效，但它们仍在与某些挑战作斗争，例如过早收敛、探索不足以及在高维、非凸搜索空间中缺乏鲁棒性。这些限制凸显了对能够在保持计算效率的同时更好地平衡探索和利用的新颖优化技术的需求。针对这一需求，我们提出了人工旅鼠算法（ALA），这是一种受生物启发的元启发式算法，它在数学上模拟了自然界中旅鼠的四种不同行为：长距离迁徙、挖洞、觅食和躲避捕食者。
具体而言，长距离迁移和挖洞行为致力于高度探索搜索域，而觅食和躲避捕食者行为则在优化过程中提供了开发。于2025年1月发表在中科院2区（IF>10）.

2.2 数学模型

2.2.1 初始化

ALA 是一种基于种群的算法，要求在进入迭代过程之前初始化所有搜索代理的位置。候选解集 $\vec{Z}$ 是一个由 $N$（种群大小）行和 $Dim$（给定问题的维度数）列组成的矩阵，位于上下界之间，如公式 (3) 所示。每次迭代中的最佳位置被认为是到目前为止获得的最优解或接近最优解。通过公式 (4) 计算每个维度的决策变量 $z_{i,j}$：

$$\vec{Z}=\left[\begin{array}{ccccc}{z}_{1,1}& z_{1,2}& \cdots & z_{1,Dim-1}& z_{1,Dim}\\ z_{2,1}& z_{2,2}& \cdots & z_{2,Dim-1}& z_{2,Dim}\\ ⋮& ⋮& ⋮& ⋮& ⋮\\ z_{N-1,1}& z_{N-1,2}& \cdots & z_{N-1,Dim-1}& z_{N-1,Dim}\\ z_{N,1}& z_{N,2}& \cdots & z_{N,Dim-1}& z_{N,Dim}\end{array}\right]$$
然后得到初始化如下
$$z_{i,j}=L{B}_j+rand\times (U{B}_j-L{B}_j),i=1,2,\cdots ,N,j=1,2,\cdots ,Dim$$

其中 rand 是范围在 0-1 之间的随机值，$L{B}_j$ 表示第 $j$ 维的下限，$U{B}_j$ 是第 $j$ 维的上限。

2.2.2 长距离迁移（探索）

在第一种行为中，当食物稀缺时，鼹鼠会随机进行长距离迁移。在这里，鼹鼠将根据当前位置和种群中随机个体的位置探索搜索空间，以找到资源丰富的栖息地，从而获得更好的生存条件和资源。同时，值得注意的是，鼹鼠迁移的方向和距离由于各种因素（如生态环境）而不稳定。为了模拟这种行为，提出了以下公式：

$${\vec{Z}}_i(t+1)={\vec{Z}}_{best}(t)+F\times BM\times \left(\vec{R}\times \left({\vec{Z}}_{best}(t)-{\vec{Z}}_i(t)\right)+\left(1-\vec{R}\right)\times \left({\vec{Z}}_b(t)-{\vec{Z}}_i(t)\right)\right)$$

其中 ${\vec{Z}}_i(t+1)$ 表示第 $i$ 个搜索代理在 $(t+1)$ 次迭代中的位置，${\vec{Z}}_{best}(t)$ 表示当前最优解。$F$ 用于改变搜索方向，通过公式 (7) 计算，这有助于避免局部最优并为搜索代理提供更高的机会来严格扫描问题域。$BM$ 表示随机数向量，用于描述布朗运动，其具有动态和均匀的步长，以探索搜索空间中的某些潜在区域。标准布朗运动的步长通过正态分布的概率密度函数获得，方差为 1，均值为 0，如公式 (6) 所示。$\vec{R}$ 是一个 $1\times Dim$ 的向量，其元素是均匀分布在区间 [-1, 1] 中的随机数，通过公式 (8) 生成。该向量用于表示种群中当前最优个体和随机个体之间的相互作用。${\vec{Z}}_i(t)$ 表示第 $i$ 个搜索代理的当前位置。${\vec{Z}}_b(t)$ 表示从种群中随机选择的搜索个体的位置，$a$ 是介于 1 和 $N$ 之间的整数索引。图 4 展示了二维中长距离迁移行为的数学建模原理。

$$f_{BM}(x;0,1)=\frac{1}{\sqrt{2\pi }}\times \exp \left(-\frac{x^2}{2}\right)$$
其次：
$$F=\{\begin{array}{ll}1& \text{if }[2\times rand+1]=1\\ -1& \text{if }[2\times rand+1]=2\end{array}$$
其中
$$\vec{R}=2\times rand(1,Dim)-1$$

图 4 展示了二维中长距离迁移模型的示意图

2.2.3 挖洞（探索）

鼹鼠的第二种行为涉及在栖息地中挖洞以形成复杂的隧道，提供安全的庇护所和食物储存。鼹鼠将根据当前位置和种群中随机个体的位置随机挖新洞。这种设计有助于它们快速逃离捕食者并更有效地寻找食物。图 5 展示了挖洞的简单示例。该行为通过公式 (9) 模拟：

$${\vec{Z}}_i(t+1)={\vec{Z}}_i(t)+F\times L\times \left({\vec{Z}}_{best}(t)-{\vec{Z}}_i(t)\right)$$

其中 $L$ 是与当前迭代次数相关的随机数。${\vec{Z}}_b$ 表示从种群中随机选择的搜索个体，$b$ 是介于 1 和 $N$ 之间的随机整数索引。$L$ 和 ${\vec{Z}}_b$ 用于描述鼹鼠在挖新洞时的相互作用。$L$ 的值计算如下：

$$L=rand\times \left(1+\sin \left(\frac{t}{T_{max}}\right)\right)$$

图 5 展示了二维中挖洞模型的示意图

2.2.4 觅食（开发）

在第三种行为中，鼹鼠在其栖息地的洞穴中广泛而任意地移动，依靠其敏锐的嗅觉和听觉来定位食物来源。鼹鼠通常在其栖息地内建立一个相对较小的觅食区域，该区域取决于食物的丰富度和可用性。为了尽可能多地摄入食物，鼹鼠将在其觅食区域内随机游荡。图 6 展示了鼹鼠的随机觅食行为。螺旋缠绕机制被认为可以模拟这一阶段，如下所示：

$${\vec{Z}}_i(t+1)={\vec{Z}}_{best}(t)+F\times spiral\times rand\times {\vec{Z}}_i(t)$$

其中 spiral 表示觅食过程中随机搜索的螺旋形状，通过公式 (12) 和 (13) 计算：

$$spiral=radius\times \left(\sin \left(2\times \pi \times rand\right)+\cos \left(2\times \pi \times rand\right)\right)$$
其中
$$radius=\sqrt{\frac{Dim}{{\sum }_{j=1}^{Dim}{\left(z_{best,j}(t)-z_{i,j}(t)\right)}^2}}$$

其中 radius 表示觅食范围的半径，即当前位置与最优解之间的欧几里得距离。

图 6 展示了二维中觅食模型的示意图

2.2.5 逃避天敌（开发）

在最后一个阶段，建模重点放在鼹鼠在遇到危险时的逃避和保护行为上。洞穴为鼹鼠提供了庇护所。一旦检测到敌人，鼹鼠利用其出色的奔跑能力逃回洞穴。同时，鼹鼠还进行欺骗性机动以逃避捕食者的追击。图 7 展示了鼹鼠逃避捕食者的行为。相应的数学表达式在公式 (14) 中描述：

$${\vec{Z}}_i(t+1)={\vec{Z}}_{best}(t)+F\times G\times Levy(Dim)\times \left({\vec{Z}}_{rest}(t)-{\vec{Z}}_i(t)\right)$$
其中
$$G=2\times \left(1-\frac{t}{T_{max}}\right)$$

其中 $G$ 是鼹鼠的逃避系数，表示其逃避能力，随着迭代次数的增加而降低，参见公式 (15)。$T_{max}$ 表示最大迭代次数。$Levy(\cdot )$ 是莱维飞行函数，用于模拟鼹鼠在逃避过程中的欺骗性机动。莱维飞行函数表示如下：

$$Levy(x)=0.01\times \frac{u\times \sigma }{|x{|}^{\frac{1}{\beta }}},\sigma ={\left(\frac{\Gamma \left(1+\beta \right)\times \sin \left(\frac{\pi \beta }{2}\right)}{\Gamma \left(\frac{1+\beta }{2}\right)\times \beta \times 2^{\frac{\beta -1}{2}}}\right)}^{\frac{1}{\beta }}$$

其中 $u$ 和 $v$ 是区间 [0-1] 中的随机值，$\beta$ 是常数，取值为 1.5。

图 7 展示了二维中捕食者规避模型的示意图

2.2.6 从探索到开发的过渡

在 ALA 中，四种搜索策略与鼹鼠自身的能量水平密切相关。在初始阶段，鼹鼠主要参与探索过程以定位有前景的区域，而在搜索的后期阶段，它们倾向于实施局部开发以获得更优越的全局最优解。为了在探索和开发之间保持最佳平衡，设计了一个能量因子，随着迭代次数的增加而减少。当鼹鼠有足够的能量时，它们选择性地迁移或挖洞；否则，它们觅食和逃避捕食者。能量因子的计算公式如下：

$$E(t)=4\times \arctan \left[1-\frac{t}{T_{max}}\right]\times \ln \left(\frac{1}{rand}\right)$$

图 8 展示了 500 次迭代中 $E$ 的变化曲线。可以看出，随着迭代次数的增加，$E$ 逐渐减少到 0 并出现波动。我们将阈值设为 1，因此鼹鼠在迭代过程中几乎有相同的机会进行探索和开发。设：

$$\delta =\left(1-\frac{t}{T_{max}}\right)$$

然后

$$E(t)=4\times \arctan \delta \times \ln \left(\frac{1}{rand}\right)$$

$E(t)>1$ 的概率通过公式 (20) 计算：

$$P(E(t)>1)={\int }_0^1{\int }_0^{\frac{1}{1-\delta }}\frac{1}{4{\pi }^2}\frac{1}{{\cos }^2\left(\frac{1}{4\rho }\right)}d\rho d\delta \Rightarrow y=\frac{1}{4\arctan \delta },{\int }_{-\infty }^{-\frac{1}{\delta }}{e}^u\frac{1}{4{\rho }^2{\cos }^2\left(\frac{1}{4\rho }\right)}dy\approx 0.5060$$

因此，$E>1$ 的概率约为 0.5，这有助于优化器从探索平稳过渡到开发。

参考文献

[1] Xiao Y, Cui H, Khurma R A, et al. Artificial lemming algorithm: a novel bionic meta-heuristic technique for solving real-world engineering optimization problems[J]. Artificial Intelligence Review, 2025, 58(3): 84. https://doi.org/10.1007/s10462-024-11023-7

完整代码

https://github.com/suthels/-/edit/main/README.md