智能优化算法：猎食者优化算法-附代码

智能算法研学社（Jack旭）

已于 2022-05-25 10:49:39 修改

阅读量5.9k

点赞数 4

分类专栏：智能优化算法文章标签：算法线性代数几何学

于 2022-03-12 15:05:38 首次发布

本文链接：https://blog.csdn.net/u011835903/article/details/123444104

版权

智能优化算法专栏收录该内容

149 篇文章 1193 订阅

订阅专栏

智能优化算法：猎食者优化算法

文章目录

智能优化算法：猎食者优化算法

摘要：猎食者优化算法（Hunter–prey optimization，HPO），是于2022年提出的一种新型智能优化算法，该算法通过模拟动物猎食的过程，来对问题进行寻优，具有收敛速度快，寻优能力强的特点。

1.算法原理

1.初始化

与其他优化算法一样，猎食者种群在搜索范围内随机初始化位置。
$x_{i} = \operatorname{rand}(1, d) . *(u b-l b)+l b \tag{1}$
其中， $x_i$ 是猎食的位置， $l b$ 是问题变量的最小值(下界)， $u b$ 是问题变量的最大值(上界)， $d$ 是问题变量的数量(维度)。

2.猎食者搜索

对于猎食者的搜索机制，式(2)给出了其数学模型：
$x_{i, j}(t+1)=x_{i, j}(t)+0.5\left[\left(2 C Z P_{p o s(j)}-x_{i, j}(t)\right)+\left(2(1-C) Z \mu(j)-x_{i, j}(t)\right)\right] \tag{2}$
其中， $x (t)$ 是当前猎食者位置， $x (t + 1)$ 是猎食者的下一次迭代位置，是 $P_{pos}$ 猎物的位置， $\mu$ 是所有位置的平均值， $Z$ 是由式(3)计算的自适应参数：
$\begin{array}{l} P=\vec{R}_{1}<C ; I D X=(P==0) \\ Z=R_{2} \otimes I D X+\vec{R}_{3} \otimes(\sim I D X) \end{array} \tag{3}$
其中， $\overrightarrow R1$ 和 $\overrightarrow R3$ 是[0,1]内的随机向量， $P$ 是 $\overrightarrow R1<C$ 的索引值， $\overrightarrow R2$ 是[ 0 , 1 ] 内的随机数， $I D X$ 是满足条件( $P = = 0$ )的向量 $\overrightarrow R1$ 的索引值， $C$ 是探索和开发之间的平衡参数，其值在迭代过程中从1减小到0.02，计算如下：
$t\left(\frac{0.98}{M a x I t}\right) \tag{4}$
其中， $i t$ 是当前迭代次数， $M a x I t$ 是最大迭代次数。计算猎物的位置( $P_{pos}$ )，以便首先根据式(5)计算所有位置的平均值( $\mu$ )，然后计算每个搜索个体与该平均位置的距离。
$\mu=\frac{1}{n} \sum_{i=1}^{n} \vec{x}_{i} \tag{5}$
根据式(6)计算欧几里得距离：
$D_{e u c(i)}=\left(\sum_{j=1}^{d}\left(x_{i, j}-\mu_{j}\right)^{2}\right)^{\frac{1}{2}}\tag{6}$
根据式(7)，距离位置平均值最大的搜索代理被视为猎物( $P_{pos}$ ):
$\vec{P}_{\text {pos }}=\vec{x}_{i} \mid i \text { is index of } \operatorname{Max}(\mathrm{end}) \operatorname{sort}\left(D_{\text {euc }}\right) \tag{7}$
如果每次迭代都考虑到搜索代理与平均位置( $\mu$ )之间的最大距离，则该算法将具有延迟收敛性。根据狩猎场景，当猎食者捕获猎物时，猎物会死亡，而下一次，猎食者会移动到新的猎物位置。为了解决这个问题，考虑一种递减机制，如式(8)所示：
$\text { kbest }=\operatorname{round}(C \times N) \tag{8}$
其中 $N$ 是搜索代理的数量。改变式(7)，将猎物的位置计算为式(9)：
$\vec{P}_{\text {pos }}=\vec{x}_{i} \mid i \text { is sorted } D_{\text {euc }}(\text { kbest }) \tag{9}$
在算法开始时， $k b e s t$ 的值等于 $N$ 。最后一个距离搜索个体的平均位置( $\mu$ )最远的搜索个体被选择为猎物，并被猎食者捕获。假设最佳安全位置是最佳全局位置，因为这将使猎物有更好的生存机会，猎食者可能会选择另一个猎物。式(10)用于更新猎物位置：
$x_{i, j}(t+1)=T_{p o s(j)}+C Z \cos \left(2 \pi R_{4}\right) \times\left(T_{p o s(j)}-x_{i, j}(t)\right)\tag{10}$
其中， $x (t)$ 是猎物的当前位置； $x (t + 1)$ 是猎物的下一次迭代位置； $T_{pos}$ 是全局最优位置； $Z$ 是由式(3)计算的自适应参数； $R_4$ 是范围[ − 1 , 1 ]内的随机数； $C$ 是探索和开发之间的平衡参数，其值在算法的迭代过程中减小，并由式(4)计算； $c o s$ 函数及其输入参数允许下一个猎物位置在不同半径和角度的全局最优位置，并提高开发阶段的性能。为了选择猎食者和猎物，结合式(2)和(10)提出了式(11)：
$x_{i}(t+1)=\left\{\begin{array}{ll} x_{i}(t)+0.5\left[\left(2 C Z P_{\text {pos }(j)}-x_{i}(t)\right)+\left(2(1-C) Z \mu(j)-x_{i}(t)\right)\right] & \text { if } R_{5}<\beta \\ T_{\text {pos }}+C Z \cos \left(2 \pi R_{4}\right) \times\left(T_{\text {pos }}-x_{i}(t)\right) & \text { else } \end{array}\right .\tag{11}$
其中， $R_5$ 是[ 0 , 1 ] 范围内的随机数， $\beta$ 是一个调节参数。