智能优化算法：鹈鹕优化算法-附代码

智能优化算法：鹈鹕优化算法

1.鹈鹕优化算法简介

鹈鹕优化算法(Pelican Optimization Algorithm,POA)是2022年由Pavel Trojovský和Mohammad Dehghani 提出的，该算法模拟了鹈鹕在狩猎过程中的自然行为。

2.鹈鹕优化算法基本原理

2.1灵感来源和鹈鹕在狩猎时的行为

鹈鹕很大，喙很长，喉咙里有一个大袋子，用来捕捉和吞咽猎物。这种鸟喜欢群体和社交生活，并以数百只鹈鹕为一群生活。鹈鹕的外观如下：它们的体重约为2.75至15公斤，高度约为1.06至1.83米，翼展约为0.5至3米。鹈鹕的食物主要是以鱼类为主，不常吃的是青蛙、海龟和甲壳类动物；如果它很饿，它甚至会吃海鲜。鹈鹕经常组团一起捕猎，鹈鹕在确定猎物的位置后，从10-20米的高度俯冲到猎物身上。当然，有些物种也会在低海拔地区找到它们的猎物。然后，它们在水面上展开翅膀，迫使鱼进入浅水区，以便它们能很容易地抓住鱼。当捕获鱼后，大量的水进入鹈鹕的喙，在吞下鱼之前将头部向前移动以除去多余的水。

鹈鹕在捕猎时的行为和策略是一个智能的行为过程，使得这种鸟类成为熟练的猎手。该鹈鹕优化算法设计的主要灵感来自于上述策略的建模。

2.2算法的数学模型

在鹈鹕优化算法中，模拟了鹈鹕在攻击和狩猎时的行为和策略，以此来更新候选解。该狩猎过程分为两个阶段：逼近猎物（勘探阶段），水面飞行（开发阶段）。

2.2.1 初始化

鹈鹕种群初始化数学描述如下：
$$x_{i,j}=l_j+rand\cdot \left(u_j-l_j\right),i=1,2,\dots ,N,j=1,2,\dots ,m\text{\hspace{1em}(1)}$$
式中：$x_{i,j}$为第i个鹈鹕的第j维的位置；N为鹈鹕的种群数量；m为求解问题的维度；$rand$是[0,1]范围内的随机数；和$u_j$和$l_j$分别是求解问题的第j维的上下边界。

在鹈鹕优化算法中，鹈鹕种群可以用以下种群矩阵表示：
$$X={\left[\begin{array}{c}{X}_1\\ ⋮\\ X_i\\ ⋮\\ X_N\end{array}\right]}_{N\times m}={\left[\begin{array}{ccccc}{x}_{1,1}& \cdots & x_{1,j}& \cdots & x_{1,m}\\ ⋮& \ddots & ⋮& & ⋮\\ x_{i,1}& \cdots & x_{i,j}& \cdots & x_{i,m}\\ ⋮& & ⋮& \ddots & ⋮\\ x_{N,1}& \cdots & x_{N,j}& \cdots & x_{N,m}\end{array}\right]}_{N\times m}\text{\hspace{1em}(2)}$$
式中：$X$为鹈鹕的种群矩阵；$X_i$为第i个鹈鹕的位置；

在鹈鹕优化算法中，求解问题的目标函数可以用来计算鹈鹕的目标函数值；鹈鹕种群的目标函数值可以用目标函数值向量表示：
$$F={\left[\begin{array}{c}{F}_1\\ ⋮\\ F_i\\ ⋮\\ F_N\end{array}\right]}_{N\times 1}={\left[\begin{array}{c}F\left(X_1\right)\\ ⋮\\ F\left(X_i\right)\\ ⋮\\ F\left(X_N\right)\end{array}\right]}_{N\times 1}\text{\hspace{1em}(3)}$$
式中：$F$为鹈鹕种群的目标函数向量；$F_i$为第i个鹈鹕的目标函数值。

2.2.2 第一阶段：逼近猎物（勘探阶段）

在第一阶段，鹈鹕确定猎物的位置，然后向这个确定的区域移动。对鹈鹕逼近猎物策略进行建模，使得POA算法可以对搜索空间进行扫描，进而发挥POA算法在搜索空间中的不同区域的勘探能力。在POA算法中，需要重视的一点是：猎物的位置在搜索空间中是随机生成的，这样增加了POA算法在解决精确搜索问题上的勘探能力。对上述概念和逼近猎物策略进行数学建模，如下：
$$x_{i,j}^{P_1}=\{\begin{array}{ll}{x}_{i,j}+\text{ rand }\cdot \left(p_j-I\cdot x_{i,j}\right),& F_p<F_i\\ x_{i,j}+\text{ rand }\cdot \left(x_{i,j}-p_j\right),& \text{ else }\end{array}\text{\hspace{1em}(4)}$$
式中：$x_{i,j}^{P_1}$为基于第1阶段更新后第i个鹈鹕的第j维的位置；rand是[0,1]范围内的随机数；$I$为1或2的随机整数；$p_j$为猎物的第j维的位置；$F_p$为猎物的目标函数值。

在POA算法中，如果目标函数值在该位置得到改善，则接受鹈鹕的新位置。在这种类型的更新中，也被称为有效更新，该算法不能移动到非最优区域。这个过程可以如下公式描述：
$$X_i=\{\begin{array}{cc}{X}_i^{P_1},& F_i^{P_1}<F_i\\ X_i,& \text{ else }\end{array}\text{\hspace{1em}(5)}$$
式中：$X_i^{P_1}$为第i个鹈鹕的新位置；$F_i^{P_1}$为基于第一阶段更新后的第i个鹈鹕的新位置的目$X_i^{P_1}$标函数值。

2.2.3 第二阶段：水面飞行（开发阶段）

在第二阶段，当鹈鹕到达水面后，它们在水面上展开翅膀，将鱼向上移动，然后把猎物放在它们的喉咙袋里。鹈鹕水面飞行的这种策略可以使得它们在被攻击区域内捕获更多的鱼。对鹈鹕的这种行为过程进行建模，可以使得POA算法收敛到狩猎区域更好的位置，这样增加了POA算法的局部搜索能力和开发能力。从数学的视角来看，该算法必须检查鹈鹕位置附近的位置，以便算法可以收敛到更好的位置。鹈鹕在狩猎过程中的这种行为在数学建模为：
$$x_{i,j}^{P_2}=x_{i,j}+R\cdot \left(1-\frac{t}{T}\right)\cdot (2\cdot rand-1)\cdot x_{i,j}\text{\hspace{1em}(6)}$$
式中：$x_{i,j}^{P_2}$为基于第2阶段更新后第i个鹈鹕的第j维的位置；rand是[0,1]范围内的随机数；$R$为0或2的随机整数；t为当前迭代次数；T为最大迭代次数。

算法流程图如下：

3.实验结果

4.参考文献

[1] Trojovský P, Dehghani M. Pelican Optimization Algorithm: A Novel Nature-Inspired Algorithm for Engineering Applications. Sensors (Basel). 2022 Jan 23;22(3):855. doi: 10.3390/s22030855. PMID: 35161600; PMCID: PMC8838090.

5.Matlab代码