蚁狮优化算法( Ant Lion Optimizer,ALO)-Matlab源码

获取更多资讯，赶快关注上面的公众号吧！

蚁狮优化器Ant Lion Optimizer (ALO)

2015年，Seyedali Mirjalili提出了一种新型的受自然启发的算法——蚁狮优化算法(Ant Lion Optimizer ,ALO)，该算法模模仿了自然界中蚁狮在捕猎蚂蚁时的智能行为。关注公众号后发送“蚁狮”或“ALO”获取Matlab源码下载链接！

启发

蚁狮属于蚁蛉科，为脉翅目类昆虫。蚁狮的生命周期包括两个主要阶段:幼虫和成虫，自然总寿命可达3年，主要发生在幼虫(成年期为3-5周)。

它们的名字源于它们独特的捕猎行为和它们最喜爱的猎物。一只蚁狮幼虫沿着一条圆形的路径在沙子里挖出一个圆锥形的坑，并用它巨大的下颚把沙子扔出去。图1（a）显示了几个不同大小的锥状坑。挖好陷阱后，幼虫藏在圆锥体的底部(作为一个坐等捕食者)，等待昆虫(最好是蚂蚁)被困在坑里，如拖1（b）所示。锥的边缘很急剧，昆虫可以很容易地掉到陷阱的底部。一旦蚁狮意识到猎物在陷阱里，它就会试图抓住它。然而，昆虫通常不会立即被捕捉，它们会试图逃离陷阱。在这种情况下，蚁群会聪明地将沙子扔向洞穴边缘，将猎物滑入洞穴底部。当猎物被抓住时，它会被拉到土壤下面被吃掉。吃掉猎物后，蚁狮会把剩下的猎物扔到坑外，并修复坑为下次捕猎做准备。

图1 左（a），右（b）。锥状陷阱和蚁狮的捕猎行为

另一个观察到的关于蚁狮生活方式的有趣行为是陷阱的大小和两件事的相关性：饥饿程度和月亮形状。当蚁狮感到饥饿或者月圆时，往往会挖出更大的陷阱。它们已经进化并适应了这种方式来提高它们的生存机会。人们还发现，安狮子不会直接通过观察月亮的形状来决定陷阱的大小，但它有一个内部的月球时钟来做出这样的决定。

ALO算子

ALO算法模拟了蚁狮和陷阱中的蚂蚁之间的交互，为了模拟这种交互作用，蚂蚁需要在搜索空间中移动，而蚁狮则可以捕猎它们。由于蚂蚁在自然界中寻找食物时是随机移动的，因此通过选择一个随机行走来模拟蚂蚁的运动：
$$X(t)=\left[0,\text{ cumsum }\left(2r\left(t_1\right)-1\right),\text{ cumsum }\left(2r\left(t_2\right)-1\right),\dots ,\text{ cumsum }\left(2r\left(t_n\right)-1\right)\right]\text{\hspace{1em}(1)}$$

其中$\text{cumsum}$计算了累积和，$n$为最大迭代次数，$t$表示随机游走的步数（即迭代），$r(t)$通过以下定义的随机函数：

$$r(t)=\{\begin{array}{ll}1& \text{ if rand }>0.5\\ 0& \text{ if rand }⩽0.5\end{array}\text{\hspace{1em}(2)}$$

其中$rand$为[0,1]内服从均匀分布的随机数。

在优化过程中将蚂蚁的位置保存至以下矩阵中：
$$M_{Ant}=\left[\begin{array}{ccccc}{A}_{1,1}& A_{1,2}& \dots & \dots & A_{1,d}\\ A_{2,1}& A_{2,2}& \dots & \dots & A_{2,d}\\ :& :& :& :& ⋮\\ :& ⋮& :& :& ⋮\\ A_{n,1}& A_{n,2}& \dots & \dots & A_{n,d}\end{array}\right]\text{\hspace{1em}(3)}$$

其中$M_{Ant}$为用于存储每只蚂蚁位置的矩阵，$A_{i,j}$为第$i$只蚂蚁的第$j$个变量（维度）的值，$n$为蚂蚁数量，$d$为变量个数。这里的蚂蚁有点类似于PSO中的粒子，矩阵$M_{Ant}$中保存了所有蚂蚁的位置（所有解的变量）。

采用如下的适应度函数评估每只蚂蚁的目标值：
$$M_{OA}=\left[\begin{array}{c}f\left(\left[A_{1,1},A_{1,2},\dots ,A_{1,d}\right]\right)\\ f\left(\left[A_{2,1},A_{2,2},\dots ,A_{2,d}\right]\right)\\ ⋮\\ ⋮\\ f\left(\left[A_{n,1},A_{n,2},\dots ,A_{n,d}\right]\right)\end{array}\right]\text{\hspace{1em}(4)}$$
其中$M_{OA}$中保存了每只蚂蚁的适应度值，$f$为目标函数。

除了蚂蚁，我们认为蚁狮也隐藏在搜索空间的某个地方，同样它也有对应的位置矩阵和适应度值矩阵：

$$M_{\text{Antlion}}=\left[\begin{array}{ccccc}A{L}_{1,1}& A{L}_{1,2}& \dots & \dots & A{L}_{1,d}\\ A{L}_{2,1}& A{L}_{2,2}& \dots & \dots & A{L}_{2,d}\\ :& :& ⋮& :& :\\ ⋮& ⋮& ⋮& :& ⋮\\ A{L}_{n,1}& A{L}_{n,2}& \dots & \dots & A{L}_{n,d}\end{array}\right]\text{\hspace{1em}(5)}$$
其中$M_{\text{Antlion}}$为每只蚁狮的位置矩阵，$A{L}_{i,j}$为第$i$只蚁狮的第$j$个维度的值，$n$为蚁狮数量，$d$为变量个数（维度）。

$$M_{0AL}=\left[\begin{array}{c}f\left(\left[A{L}_{1,1},A{L}_{1,2},\dots ,A{L}_{1,d}\right]\right)\\ f\left(\left[A{L}_{2,1},A{L}_{2,2},\dots ,A{L}_{2,d}\right]\right)\\ ⋮\\ ⋮\\ f\left(\left[A{L}_{n,1},A{L}_{n,2},\dots ,A{L}_{n,d}\right]\right)\end{array}\right]\text{\hspace{1em}(6)}$$

$M_{0AL}$为保存每只蚁狮适应度值的矩阵，$f$为目标函数。

优化时需要满足以下条件：

• 蚂蚁使用不同的随机游走在搜索空间中移动。
• 随机游走应用于蚂蚁的所有维数。
• 随机游走受到蚁狮陷阱的影响。
• 蚁狮可以建立与适合度成比例的坑(适合度越高，坑越大)。
• 坑大的蚁狮更容易捉到蚂蚁。
• 每只蚂蚁都可以在每次迭代中被一只蚁狮和精英(最优的蚁狮)捕获。
• 自适应地减小随机游走范围来模拟蚂蚁滑向蚁狮。
• 如果蚂蚁优于蚁狮，这意味着它会被蚁狮抓住并拖到沙子下面。
• 蚁狮会根据最近捕获的猎物重新定位，并筑起一个坑来改善每次捕猎后捕获另一个猎物的变化。

随机游走

随机游走均基于式(1)。蚂蚁在每一步的优化中都通过随机游走来更新它们的位置。但是，由于每个搜索空间都有一个边界(变量范围)，式（1）不能直接用于蚂蚁位置的更新。为了使随机游走保持在搜索空间内，使用下式对其进行归一化（最小-最大归一化）：
$$X_i^t=\frac{\left(X_i^t-a_i\right)\times \left(d_i-c_i^t\right)}{\left(d_i^t-a_i\right)}+c_i\text{\hspace{1em}(7)}$$

其中$a_i$为第$i$个变量随机游走的最小值，$d_i$是最大值。$c_i^t$是第$t$次迭代中第$i$个变量的最小值，$d_i^t$是第$t$次迭代中第$i$个变量的最大值。

在每次迭代中都需要使用式（7），以保证搜索空间内出现随机游走。

困于蚁狮坑

如上所述，蚂蚁的随机游走会受到蚁狮陷阱的影响。为了对这一假设进行数学建模，提出以下方程：
$$c_i^t={\text{Antlion}}_j^t+c^t\text{\hspace{1em}(8)}$$

$$d_i^t={\text{Antlion}}_j^t+d^t\text{\hspace{1em}(9)}$$
其中$c^t$为第$t$次迭代中所有变量的最小值，$d^t$为第$t$次迭代中所有变量的最大值，$c_i^t$为第$i$只蚂蚁所有变量的最小值，$d_i^t$为第$i$只蚂蚁所有变量的最大值，${\text{Antlion}}_j^t$表示第$t$次迭代中选择的第$j$只蚁狮的位置。

式（8）和（9）表达了由向量$c$和$d$所定义的超球面中所选择的蚁狮附近的蚂蚁随机游走，如图3所示。

图3 蚂蚁在蚁狮陷阱内的随机游走

构造陷阱

为了模拟蚁狮的狩猎能力，使用了轮盘赌。如图3所示，假设蚂蚁只被困在一个选定的蚁狮中。ALO算法需要利用轮盘赌算子在优化过程中根据蚁狮的适应度选择蚁狮。这种机制为更适应的蚁狮捕捉蚂蚁提供了很大的机会。

蚂蚁滑向蚁狮

根据目前提出的机制，蚁狮可以根据自身的适合度设置陷阱，而且蚂蚁需要随机移动。然而，一旦蚁狮意识到有一只蚂蚁被困在陷阱里，它们就会把沙子从洞中央往外喷。这种行为可以将试图逃跑的被困蚂蚁滑选。通过自适应地减小蚁群随机游走超球面半径，建立了蚁群随机游走超球面的数学模型：
$$c^t=\frac{c^t}{I}\text{\hspace{1em}(10)}$$

$$d^t=\frac{d^t}{I}\text{\hspace{1em}(11)}$$
其中$I$为比率，其中$c^t$为第$t$次迭代中所有变量的最小值，$d^t$为第$t$次迭代中所有变量的最大值。

在式（10）和（11）中，$I=10^w\frac{t}{T}$，其中$t$为当前迭代，$T$为最大迭代次数，$w$是根据当前迭代定义的常量（$t>0.1T$时$w=2$，$t>0.5T$时$w=3$，$t>0.75T$时$w=4$，$t>0.9T$时$w=5$，$t>0.95T$时$w=6$）。$w$可以调整利用的精度水平，如图4所示。该方程缩小了蚂蚁位置的更新半径，模拟了蚂蚁在坑内的滑动过程，从而保证了在搜索空间内进行利用。

图4 自适应的上下界

捕捉猎物和重建坑

捕猎的最后阶段是当一只蚂蚁到达洞穴底部时被蚁狮的下巴夹住。过了这个阶段，蚁狮就把蚂蚁拉进沙子里并吃掉它的身体。为了模拟这个过程，假设当蚂蚁比对应的蚁狮更适应(进入沙子里)时，捕获猎物就会出现。然后，蚁狮需要将自己的位置更新到被猎杀蚂蚁的最新位置，以增加它捕捉新猎物的机会。对此有如下公式：
$${\text{ Antlion }}_j^t=An{t}_i^t\text{ if }f\left(An{t}_i^t\right)>f\left({\text{Antlion}}_j^t\right)\text{\hspace{1em}(12)}$$
${\text{ Antlion }}_j^t$表示第$t$次迭代中第$j$只蚁狮的位置，${\text{ Ant }}_i^t$表示第$t$次迭代中第$i$只蚂蚁的位置

精英

精英主义是进化算法的一个重要特征，它使进化算法能够在优化过程的任何阶段保持获得的最佳解。在本研究中，保存到目前为止在每个迭代中获得的最佳蚁狮，并将其视为一个精英。因为精英是最适合的蚁狮，它应该能够在迭代过程中影响所有蚂蚁的运动。因此，假设每只蚂蚁同时围绕通过轮盘赌随机选择的蚁狮和精英进行随机游走，有如下：
$$An{t}_i^t=\frac{R_A^t+R_E^t}{2}\text{\hspace{1em}(13)}$$

其中$R_A^t$为在第$t$次迭代中围绕根据轮盘赌选择的蚁狮的随机游走，$R_E^t$为在第$t$次迭代中围绕精英的随机游走，$An{t}_i^t$表示第$t$次迭代中第$i$只蚂蚁的位置。

ALO算法

根据上述提出的算子，现在可以定义ALO优化算法。该算法被定义为一个三元组函数，用于逼近优化问题的全局最优，如下所示：
$$ALO(A,B,C)\text{\hspace{1em}(14)}$$

其中$A$用于生成随机初始解，$B$用于对初始解进行操作，当满足终止条件时$C$返回true。函数$A,B,C$的定义如下：

∅ ⟶ A { M A n t , M O A , M A n t l i o n , M O A L } (15) \emptyset \stackrel{A}{\longrightarrow}\left\{M_{A n t}, M_{O A}, M_{A n t l i o n}, M_{O A L}\right\}\tag{15} ∅⟶A​{MAnt​,MOA​,MAntlion​,MOAL​}(15)

{ M A n t , M A n t l i o n } → B { M A n t , M A n t l i o n } (16) \left\{M_{A n t}, M_{A n t l i o n}\right\} \stackrel{B}{\rightarrow}\left\{M_{A n t}, M_{A n t l i o n}\right\}\tag{16} {MAnt​,MAntlion​}→B{MAnt​,MAntlion​}(16)

{ M A n t , M Antlion } → C { true, false } (17) \left\{M_{A n t}, M_{\text {Antlion}}\right\} \stackrel{C}{\rightarrow}\{\text {true, false}\}\tag{17} {MAnt​,MAntlion​}→C{true, false}(17)

ALO算法的伪代码定义如下：

图5 ALO算法伪代码