2019智能优化算法：哈里斯鹰算法-附代码
@[toc]
摘要：2019 年 Heidari 等人提出哈里斯鹰优化算法(Harris Hawk Optimization, HHO)，该算法有较强的全局搜索能力，并且需要调节的参数较少的优点。
1.算法原理
哈里斯鹰优化算法是一种模拟哈里斯鹰捕食行为的智能优化算法，主要由 3 部分组成：搜索阶段、搜索与开发的转换和开发阶段。
1.1 搜索阶段
哈里斯鹰随机栖息在某个地方，通过 2 种策略找到猎物：
其中，分别为当前和下一次迭代式时个体的位置，
为迭代次数，
为随机选出的个体位置，
为猎物位置，即拥有最优适应度的个体位置，
都是[0,1]之间的随机数。
用来随机选择要采用的策略，
为个体平均位置，表达式为：
其中， 为种群中第
个个体的位置，
为种群规模。
1.2 搜索与开发的转换
HHO 算法根据猎物的逃逸能量在搜索和不同的开发行为之间转换，逃逸能量定义为：
其中，是猎物的初始能量，为 [-1,1] 之间的随机数，每次迭代时自动更新，t为迭代次数，T 为最大迭代次数。当
时进入搜索阶段，
当时进入开发阶段。
1.3 开发阶段
定义r为[0,1] 之间的随机数，用于选择不同的开发策略。当且
时，采取软围攻策略进行位置更新：
其中， 表示猎物位置与个体当前位置的差值，
为 [0, 2] 之间的随机数。
当且
时采取硬围攻策略进行位置更新：
当且
时，采取渐近式快速俯冲的软包围策略进行位置更新：
其中， 为适应度函数，
为 2 维随机向量，元素为[0,1] 之间的随机数，
是莱维飞行的数学表达式。
当且
时，采取渐近式快速俯冲的硬包围策略进行位置更新：
算法步骤：
步骤 1：种群初始化。根据搜索空间每一维的上界和下界，初始化每个个体。
步骤 2：计算初始适应度。将适应度最优的个体位置设为当前猎物位置。
步骤 3：位置更新。先通过更新猎物逃逸能量，然后根据逃逸能量和生成的随机数执行搜索或开发行为中对应的位置更新策略。
步骤 4：计算适应度。计算位置更新后的个体适应度，并与猎物适应度值进行比较，若位置更新后的个体适应度值优于猎物，则以适应度
值更优的个体位置作为新的猎物位置。
重复步骤 3 和步骤 4，当算法迭代次数达到最大迭代次数时。输出当前猎物位置作为目标的估计位置。
2.算法结果
3.参考文献
[1] HEIDARI A A, MIRJALILI S, FARIS H, et al. Harris hawks optimization: algorithm and applications[J]. Future Generation Computer Systems, 2019, 97: 849-872.
4.Matlab代码
https://mianbaoduo.com/o/bread/aJiak5o=
文献复现：
[1]汤安迪,韩统,徐登武,谢磊.混沌精英哈里斯鹰优化算法[J/OL].计算机应用:1-10[2021-01-29].http://kns.cnki.net/kcms/detail/51.1307.TP.20210114.0947.032.html.
https://mianbaoduo.com/o/bread/YZaakp5v
