胡春安,熊昱然

(江西理工大學(xué)信息工程學(xué)院,江西 贛州 341000)

1 引言

以模擬自然界中動物種群行為為基礎(chǔ)設(shè)計的算法被稱為元啟發(fā)式算法[1,2]。這些行為包括繁殖、捕獵和遷徙等。通過對這些行為進行數(shù)學(xué)建模,來解決一些分布式問題[3]。近年來,群優(yōu)化算法[4]逐漸被應(yīng)用于諸如圖像處理、路徑規(guī)劃和數(shù)據(jù)挖掘等方面,并取得了相應(yīng)的研究成果。目前,群智能算法發(fā)展現(xiàn)狀已經(jīng)從簡單的粒子群演變到一些具有一定感知、預(yù)警和自主搜尋能力的種群。比較常見的群智能算法有粒子群優(yōu)化算法[4]、蟻群算法[5]、麻雀搜索算法[6]和布谷鳥算法[7]等。群算法的基礎(chǔ)理論提出后主要的發(fā)展方向有2個:一個是在算法的基礎(chǔ)上進行改進,優(yōu)化其中的種群行為策略;另一個是將算法應(yīng)用在實際工程案例中,獲得更好的優(yōu)化效果。

哈里斯鷹優(yōu)化HHO(Harris Hawk Optimization)算法是由Heidari等人[8]在2019年提出的一種元啟發(fā)算法,源于觀察哈里斯鷹和它們的獵物(兔子)之間的追逐逃跑行為。哈里斯鷹捕捉獵物的主要策略是突襲。在這個智能策略里,鷹群從幾個不同的方向協(xié)同攻擊檢測到的正在逃跑的兔子。由于獵物存在一定的逃生能力,哈里斯鷹可能在幾分鐘內(nèi)在獵物附近進行多次短距離的快速攻擊。鷹群能根據(jù)環(huán)境和獵物的逃跑狀態(tài)表現(xiàn)出幾種不同的追逐風(fēng)格。當最強的鷹(領(lǐng)導(dǎo)者)俯身沖向獵物并且迷路時,就會發(fā)生策略的轉(zhuǎn)換,由鷹群中另一名成員繼續(xù)追擊。

與其他的元啟發(fā)式算法類似,傳統(tǒng)的哈里斯鷹優(yōu)化算法本身也存在在開發(fā)階段搜索容易陷入局部最優(yōu)、收斂精度低等問題[9]。為了使算法有所改善,郭雨鑫等人[10]提出了多策略優(yōu)化的哈里斯鷹優(yōu)化算法,在探索階段引入了柯西變異,局部搜索階段引入了自適應(yīng)權(quán)重粒子來更新位置,解決原始哈里斯鷹優(yōu)化算法易陷入局部最優(yōu)和收斂精度低的問題。馬一鳴等人[11]提出了一種改進的哈里斯鷹優(yōu)化定位算法,在提升原算法性能的基礎(chǔ)上保留其尋優(yōu)機制,對基于最大似然估計的適應(yīng)度函數(shù)進行改進,在優(yōu)化過程中達到更優(yōu)的適應(yīng)度值,從而提高算法的尋優(yōu)精度。王歸新等人[12,13]通過結(jié)合2種較新的群智能算法(改進灰狼優(yōu)化算法與哈里斯鷹優(yōu)化算法),從而獲得了二者易實現(xiàn)、收斂速度快和全局搜索能力強的優(yōu)點。

目前的HHO改進策略雖然在一定程度上改進了算法的尋優(yōu)能力,但只是針對其中的某一種更新策略進行更改或者改變逃逸能量的變化方式,并沒有有效解決尋優(yōu)的盲目性。針對以上問題,本文提出了一種融合分布式估計策略、混沌局部搜索策略和精英池策略的多策略改進哈里斯鷹優(yōu)化算法。

2 哈里斯鷹優(yōu)化算法

哈里斯鷹優(yōu)化算法通過構(gòu)建數(shù)學(xué)模型來模擬現(xiàn)實中哈里斯鷹在不同情況下的捕獵策略,整個算法分為探索階段和開發(fā)階段。

2.1 探索階段

在HHO算法中,哈里斯鷹個體是候選解,每個迭代步驟中的目標獵物為所求問題的最優(yōu)解,并根據(jù)2種策略來發(fā)現(xiàn)獵物,如式(1)所示:

(1)

其中,t表示為當前迭代次數(shù),X(t+1)表示第t+1次迭代時鷹的位置;Xrabbit(t)表示第t次迭代時獵物的位置;X(t)表示第t次迭代時鷹的位置;r1,r2,r3,r4,q是(0,1)的隨機數(shù);LB和UB分別表示搜索空間的下界和上界;Xrand(t)表示第t次迭代時種群中隨機選擇的鷹的位置;Xm(t)表示第t次迭代時鷹群位置的平均值,如式(2)所示:

(2)

其中,Xi(t)表示第t次迭代時第i只鷹的位置,N表示鷹的總數(shù)。

2.2 過渡階段

根據(jù)獵物的逃逸能量,HHO算法可以從探索階段轉(zhuǎn)移到開發(fā)階段。在獵物逃跑的過程中,獵物的能量是遞減的,本文把逃逸能量抽象成式(3):

(3)

其中,E表示當前獵物的逃逸能量;T表示最大迭代次數(shù);E0表示初始能量,是(-1,1)的隨機數(shù)。

在每次迭代過程中,當逃逸能量|E|≥1時,進入探索階段;逃逸能量|E|<1時,進入局部開發(fā)階段。

2.3 開發(fā)階段

HHO提出了4種可能的策略來模擬哈里斯鷹的攻擊手段。用r來表示獵物的逃跑概率。

Case 1:軟圍攻。當r≥0.5且|E|≥0.5時,在這種攻擊方式中,哈里斯鷹群將會對它進行軟包圍,使獵物更加疲憊,然后進行突襲。該行為的抽象數(shù)學(xué)模型如式(4)所示:

X(t+1)=ΔX(t)-E|JXrabbit(t)-X(t)|,

ΔX(t)=Xrabbit(t)-X(t)

(4)

其中,ΔX(t)表示第t次迭代時獵物位置與哈里斯鷹距離之差;J表示獵物在逃跑過程中的隨機跳躍強度,是一個(0,2)的隨機數(shù),每次迭代中都會隨機變化。

Case 2:硬圍攻。當r≥0.5且|E|<0.5時,獵物此時逃逸能量低,哈里斯鷹直接進行最終突襲。在這種情況下用式(5)更新當前位置:

X(t+1)=Xrabbit(t)-E|ΔX(t)|

(5)

Case 3:累速俯沖式軟圍攻。當r<0.5且|E|≥0.5時,獵物仍然有足夠的能量成功逃脫,在突襲之前,哈里斯鷹建立了一個軟圍攻策略,將Levy飛行集成進了HHO中,這個階段的位置更新策略如式(6)～式(8)所示:

Z=Y+S·LF(D)

(6)

Y=Xrabbit(t)-E|JXrabbit(t)-X(t)|

(7)

(8)

其中,D表示當前問題的維度,S表示D的隨機行向量,Y表示累速俯沖式軟圍攻公式計算出的位置向量,Z表示經(jīng)過Levy飛行擾動過的位置,F(·)表示求當前適應(yīng)度值的函數(shù),LF(·)表示Levy飛行生成的擾動值。

Case 4:累速俯沖式硬圍攻。當r<0.5且|E|<0.5時,獵物沒有足夠的能量逃跑,在突襲捕捉獵物之前哈里斯鷹會進行一次硬圍攻。在這次圍攻中,哈里斯鷹試圖縮小它與獵物的距離,因此按照式(9)～式(11)所示的策略更新:

Z=Y+S·LF(D)

(9)

Y=Xrabbit(t)-|EXrabbit(t)-Xm(t)|

(10)

(11)

3 改進的哈里斯鷹優(yōu)化算法

通過分析HHO算法可知,HHO也存在開發(fā)能力不足、種群多樣性下降、容易陷入局部最優(yōu)等缺點,這是由于算法對搜索空間開發(fā)和探索的不平衡造成的。為了提高算法的性能,平衡HHO算法的開發(fā)和探索能力,本文提出了一種具有3種改進策略的HHO算法變體,叫做多策略改進的哈里斯鷹優(yōu)化算法MHHO(Harris Hawk Optimization algorithm using Multi-strategy)。首先,在HHO中引入混沌局部搜索策略,提高算法的開發(fā)能力。利用混沌映射的優(yōu)點,圍繞當前個體進行局部搜索,從而找到更好的個體。其次,為了增強種群多樣性,提出了精英備選池策略。此外,為了修改進化方向,本文采用了分布估計策略,通過對優(yōu)勢種群信息進行采樣,更好地引導(dǎo)種群方向,從而提高算法收斂速度。最后,采用貪婪策略充分保留優(yōu)勢個體,加快收斂速度。

3.1 混沌局部搜索策略

混沌局部搜索策略通過搜索每一個解決方案的附近區(qū)域,從而找到更好的解決方案,因此該策略能夠有效提高算法的開發(fā)能力。另一方面,混沌映射具有隨機性和遍歷性的特點,所以利用混沌映射能夠進一步提高局部搜索策略的有效性。在MHHO中,混沌局部搜索策略只應(yīng)用于種群的優(yōu)勢種群。混沌局部搜索策略的公式如式(12)所示:

(12)

混沌映射的方法有很多種[14],最常見的有Logistic映射、Tent映射和Sin映射。但是,Logistic映射的遍歷性較差,對初始參數(shù)的設(shè)置較敏感,映射點在邊緣位置的密度大,在中間位置的密度相對較小,與本文想要達到的種群均衡性的預(yù)期不符合。而Tent映射值的分布比較平坦、均勻,映射的均勻性相比于Logistic映射的均勻性更好,能夠提高算法的尋優(yōu)速度?？紤]結(jié)合完全隨機化的特點,MHHO將Tent混沌映射生成Ct,同時在[0,1]產(chǎn)生分布較均勻的初始值。Tent混沌映射分布值如圖1所示,Tent映射表達式如式(13)所示:

Figure 1 Tent chaos mapping

(13)

3.2 精英備選池策略

在早期迭代階段,HHO通過跟隨最優(yōu)個體更新其余個體位置信息,雖然這有利于算法加快收斂,但會減弱種群的多樣性,容易導(dǎo)致算法陷入局部最優(yōu)[15]。此外,在迭代后期,使用Levy捕食策略時,同樣僅跟隨最優(yōu)個體。為了保證算法開發(fā)和探索的平衡,本文提出了一種精英備選池策略,將當前最好的3個個體放入一個集合,如式(14)所示:

Xesp={Xesp1,Xesp2,Xesp3}

(14)

其中,Xesp1、Xesp2和Xesp3是種群中最好的3個個體。

食物源每次都從這3個個體中隨機選擇。通過使用精英備選池策略,食物源的位置從最優(yōu)個體變?yōu)樽詈玫?個個體中的一個,這在一定程度上避免了因最優(yōu)個體陷入局部最優(yōu)而導(dǎo)致算法早熟收斂。為了平衡算法的開發(fā)和探索,本文算法還將全局最優(yōu)個體放入精英備選池,以確保每一個個體有機會向最優(yōu)個體靠近,確保算法的收斂效率。因此,精英備選池策略的最終數(shù)學(xué)模型描述如式(15)所示:

Xesp={Xesp1,Xesp2,Xesp3,Xbest}

(15)

避免局部最優(yōu)一直是群智能算法的一個重大難題。本文通過精英池策略能夠有效地避免個體陷入局部最優(yōu)。最終的精英池由幾個最好的個體隨機組成,這保留了選擇優(yōu)勢個體的可能性,也給當前種群提供了更大的選擇范圍。由于精英池內(nèi)的個體是隨著每次迭代進行更新的,所以不必擔(dān)心精英池內(nèi)的個體陷入局部最優(yōu)導(dǎo)致整個算法陷入局部最優(yōu)的情況發(fā)生。

為了更直觀地比較原始HHO算法與使用精英池策略改進后HHO算法的區(qū)別,本文選取CEC2017中F3、F4、F9和F10作為測試函數(shù)。

Table 1 Test functions

從圖2中可以很明顯地看出,在4個測試函數(shù)上,精英池HHO算法比原始HHO算法有更好的收斂性和收斂速度。

3.3 分布估計策略

在原始的HHO中,迭代過程中種群的更新主要由迄今為止得到的最佳解引導(dǎo),因此種群的多樣性在算法過程的后期迅速減少,導(dǎo)致HHO很容易陷入局部最優(yōu)。分布估計策略的核心是加權(quán)協(xié)方差矩陣,具有較強的局部最優(yōu)規(guī)避能力。因此,本文采用加權(quán)極大似然估計法對分布模型進行估計,選取表現(xiàn)較好的前一半作為優(yōu)勢總體,以充分利用有希望的種群信息來估計更好的進化方向。該策略的數(shù)學(xué)模型描述如下:

(16)

(17)

Cov=

(18)

(19)

(20)

在每一次迭代完成之后,使用貪婪策略保留父代群體和子代群體中的最優(yōu)NP個個體形成新的群體,這樣進一步提高了MHHO算法的全局收斂能力。

3.4 改進算法流程圖

MHHO的流程圖如圖3所示。

Figure 3 Flow chart of MHHO algorithm

3.5 時間復(fù)雜度分析

本節(jié)對HHO算法和MHHO算法的時間復(fù)雜度進行分析。對于原始的HHO算法而言,種群規(guī)模為NP,搜索空間維度為D,最大迭代次數(shù)為T,需要優(yōu)化的目標函數(shù)為f(x),由此可知HHO算法的時間復(fù)雜度為:O(NP)+O(T·NP·D)+O(T·NP)。

MHHO算法每次只對優(yōu)秀種群使用分布估計策略產(chǎn)生出M個變異新解。MHHO會計算M+NP/2個個體的適應(yīng)度值,從得到的適應(yīng)度值中挑選出前NP/2個個體進入下一輪更新。計算M+NP/2個個體的適應(yīng)度值的時間復(fù)雜度為O(T·(M+(NP/2)))。將這些得到適應(yīng)度值的個體進行排序的時間復(fù)雜度為:O(T·(M+(NP/2))log(M+(NP/2))),所以MHHO算法的總時間復(fù)雜度為:O(NP)+O(T·(M+(NP/2)))+O(T·(M+(NP/2))log(M+(NP/2)))+O(T·NP·D)。

4 仿真實驗和結(jié)果分析

4.1 測試函數(shù)和初始化參數(shù)

為了充分驗證MHHO算法的優(yōu)越性,本文使用IEEE CEC2017單目標測試函數(shù)對算法進行測試。CEC2017測試集由16個測試函數(shù)組成。F1是一個只有一個全局最優(yōu)解的單峰函數(shù),用來驗證算法的局部搜索能力。F2～F4是多峰函數(shù),主要用于測試算法跳出局部最優(yōu)的能力。F9～F12、F17和F18是混合函數(shù),F19、F20、F25～F28是復(fù)合函數(shù),可以用來測試算法在現(xiàn)實世界中解決復(fù)雜優(yōu)化問題的潛力。表2給出了函數(shù)及其全局最優(yōu)點。

Table 2 CEC2017 test functions

為了保證公平性,對于經(jīng)典測試函數(shù),所有算法都采用相同的維度,最大迭代次數(shù)設(shè)置為300,種群數(shù)量設(shè)置為50,分別對每個函數(shù)獨立求解30次。對于CEC 2017測試集,維度(Dim),最大迭代次數(shù)tmax和種群數(shù)(NP)分別被統(tǒng)一設(shè)置為30,600和500。所有測試函數(shù)獨立運行51次。本文實驗在一臺裝有Intel i5-8500U處理器和8 GB內(nèi)存的計算機上進行,并使用Matlab R2020a進行編程。

4.2 實驗結(jié)果與分析

本節(jié)將最近提出的5種算法與MHHO進行比較。這些最先進的算法包括蝴蝶優(yōu)化算法BOA(Butterfly Optimization Algorithm)、樽海鞘群算法SSA(Salp Swarm Algorithm)、算術(shù)優(yōu)化算法AOA(Arithmetic Optimization Algorithm)、尋路者算法PFA(PathFinder Algorithm)和被囊群算法TSA(Tunicate Swarm Algorithm)[16-20]。所有算法的參數(shù)設(shè)置與原文獻的一致,如表3所示。

Table 3 Algorithm parameters setting

本文通過數(shù)值分析、收斂分析、穩(wěn)定性分析、Wilcoxon符號秩檢驗、Friedman檢驗對MHHO的性能進行全面評估,結(jié)果如表4所示。

Table 4 Statistical results of seven algorithms on CEC2017 test functions

從表4可知,對于單峰測試函數(shù)F1,MHHO性能優(yōu)于所有對比算法的,雖然MHHO無法取得最優(yōu)解,但在所有對比算法中表現(xiàn)最好,表明了MHHO具有最強的開發(fā)能力;對于多峰測試函數(shù),MHHO在2個測試函數(shù)(F2和F3)上表現(xiàn)最好,在F4上PFA取得了最好的結(jié)果,MHHO均排在第2位,MHHO在多峰函數(shù)上的表現(xiàn)說明改進算法能夠保持足夠好的種群多樣性,避免了陷入局部最優(yōu);對于混合函數(shù)和復(fù)合函數(shù),每個算法各有優(yōu)劣:MHHO在F9～F12、F17、F20、F26～F28上獲得了最優(yōu)解,PFA則在F18和F25上求得了更好的解,BOA在F19上表現(xiàn)優(yōu)于其他對比算法的?？偟膩碚f,MHHO在復(fù)合函數(shù)和混合函數(shù)上取得的結(jié)果均排在前兩位,很好地表明了MHHO解決現(xiàn)實世界中復(fù)雜優(yōu)化問題的潛力。值得注意的是,MHHO在所有測試函數(shù)上的表現(xiàn)都優(yōu)于HHO,表明本文提出的改進算法性能更好。

為了分析改進算法所求解的分布特性,根據(jù)各算法51次獨立求解的結(jié)果繪制了箱式圖,如圖4所示。對于每種算法,每個框的中心標記表示51次求解函數(shù)結(jié)果的中位數(shù),框的底部和頂部邊緣表示一等分點和三等分點,符號“+”表示不在箱子內(nèi)的壞值。從圖4可以得知,在求解其中5個測試函數(shù)(F3、F12、F20、F27和F28)時,MHHO沒有異常值,這表明MHHO所求解的分布非常集中;同時,對于其它在箱式圖中存在異常值的測試函數(shù)(F1、F2、F9～F11、F17、F26),MHHO有著最小的中位數(shù),表明MHHO的解的質(zhì)量相對更優(yōu)。因此,本文提出的改進算法具有很強的魯棒性。

Figure 4 Box diagrams of CEC2017 test functions

為進一步說明算法的收斂性能,將適應(yīng)度值為以10為底的對數(shù)計量的適應(yīng)值(Fitness)作為每個算法收斂精度的評判標準。圖5展示了7種算法求解CEC2017測試集的平均誤差收斂曲線圖(迭代曲線圖),橫坐標表示算法的迭代次數(shù),縱坐標表示算法的收斂精度,也就是最終得到的適應(yīng)度值。結(jié)合以上7種算法在每個測試函數(shù)所展現(xiàn)出的進化曲線和收斂精度,實驗結(jié)果分析如下:

Figure 5 Convergence diagrams of CEC2017 test functions

對于單峰函數(shù)F1,MHHO有更好的收斂精度和更快的收斂速度。對于多峰函數(shù)F2～F4,MHHO在F2和F3上有更好的收斂精度和更快的收斂速度,PFA則在F4上具有更好的收斂精度。對于復(fù)合函數(shù)和混合函數(shù),MHHO在求解其中的9個函數(shù)(F9～F12、F17、F20、F26～F28)時,收斂速度更快,收斂精度更高。值得注意的是,MHHO在所有測試函數(shù)上的尋優(yōu)精度和收斂速度都比HHO的更好。綜上所述,MHHO在收斂精度和收斂速度方面優(yōu)于對比算法。

為了測試函數(shù)的時間消耗(單位為s),本文挑選了CEC2017中幾個有代表性的單峰、多峰以及混合函數(shù)來對改進策略進行評估。

從表5可以發(fā)現(xiàn),無論是單峰函數(shù)、多峰函數(shù)還是最復(fù)雜的混合函數(shù),MHHO的時間消耗都控制在7 s左右;HHO的時間消耗極度不穩(wěn)定,且遠高于MHHO的,這說明改進算法具有更好的性能。

Table 5 Time consumption of two algorithms on test functions

僅根據(jù)平均值比較算法性能還不充分,為了避免測試中的偶然性,本文使用Wilcoxon符號秩檢驗來驗證改進算法與對比算法是否在統(tǒng)計上存在顯著差異。表6列出了每種算法和MHHO的Wilcoxon符號秩檢驗的結(jié)果。在表6 中,“+”表示MHHO的優(yōu)化結(jié)果優(yōu)于對比算法的,“-”表示結(jié)果較差,“=”表示結(jié)果相似;符號“R+”是一個正等級值,表示MHHO比對比算法更好的程度,“R-”表示相反的結(jié)果。從表6可以看出,MHHO至少在12個測試函數(shù)上優(yōu)于所有對比算法,并且在所有測試函數(shù)(15個)上均優(yōu)于原始HHO,在統(tǒng)計學(xué)上驗證了改進算法的出色性能。

Table 6 Results of Wilcoxon signed rank test

通過CEC2017測試集對MHHO的性能進行驗證,實驗結(jié)果表明:改進策略均能有效提高算法性能;和BOA、AOA、TSA、SSA和PFA相比,MHHO無論是在收斂精度和收斂速度,還是在算法穩(wěn)定性上,均有顯著優(yōu)勢,證明了本文改進算法的優(yōu)越性。

5 工程約束問題

工程設(shè)計問題是一個幾何形狀復(fù)雜、設(shè)計變量多、實際工程約束多的非線性優(yōu)化問題。本文通過求解實際工程問題對MHHO的性能進行評估?？紤]到這些工程設(shè)計問題都是包含不等式和等式的約束優(yōu)化問題,利用罰函數(shù)將約束優(yōu)化問題轉(zhuǎn)化為無約束優(yōu)化問題。這些問題包括壓力容器設(shè)計問題和拉壓彈簧設(shè)計問題,均使用相同的迭代次數(shù)(500)和種群(200)來解決這些優(yōu)化設(shè)計問題。每個問題獨立運行30次,并將統(tǒng)計結(jié)果與文獻中的其他算法進行比較。

5.1 壓力容器設(shè)計問題

圖6所示的壓力容器設(shè)計問題是一個典型的混合優(yōu)化問題,其目標是降低包括成形成本、材料成本和焊接成本在內(nèi)的總成本。有4個不同的變量:容器厚度Ts(x1)、封頭厚度Th(x2)、內(nèi)徑R(x3)和容器圓柱截面長度L(x4)。這個問題的數(shù)學(xué)描述如式(21)～式(24)所示:

Figure 6 Schematic diagram of pressure vessel design

(21)

g2(X)=-x2+0.00954x3≤0

(22)

(23)

g4(X)=x4-240≤0

(24)

約束條件如式(25)所示:

變量范圍:1×0.0625≤x1,

x2≤99×0.0625,10≤x3,x4≤200

(25)

各算法運行結(jié)果如表7所示。分析表7中數(shù)據(jù)可知,MHHO在求解壓力容器設(shè)計問題時的尋優(yōu)精度和穩(wěn)定性是最好的。MHHO得到的最佳方案是[0.778 2,0.384 7,40.321 0,199.981 2],說明MHHO在處理該工程問題時能獲得最佳結(jié)果。

Table 7 Comparison of optimal solutions for pressure vessel design problem

5.2 拉壓彈簧設(shè)計問題

拉伸/壓縮彈簧設(shè)計問題是一個機械工程設(shè)計優(yōu)化問題,可以用來評價算法的優(yōu)越性。如圖7所示,此問題的目標是減輕彈簧的重量。它包括4個非線性不等式和3個連續(xù)變量:導(dǎo)線直徑w(x1)、線圈平均直徑d(x2)、線圈長度或數(shù)量L(x3)。

Figure 7 Schematic design of tension spring design

該問題的數(shù)學(xué)模型可以用式(26)～式(31)所示來描述:

(26)

(27)

(28)

(29)

(30)

約束條件如式(31)所示:

變量范圍:0.05≤x1≤2,

0.25≤x2≤1.3,2.0≤x3≤15.0

(31)

各算法運行結(jié)果如表8所示。分析可知,MHHO在4項統(tǒng)計數(shù)據(jù)中均表現(xiàn)最好,說明MHHO具有最好的尋優(yōu)精度和穩(wěn)定性,MHHO的最佳方案為[0.051 600,0.354 400,11.424 800]。

Table 8 Best solution for tension spring design problem

6 結(jié)束語

哈里斯鷹優(yōu)化算法是近幾年提出的一種新穎的元啟發(fā)式算法,其原理簡單易懂,算法本身參數(shù)較少,易實現(xiàn),適合與其他的元啟發(fā)式算法進行結(jié)合,也適合解決各種高維問題和多極值問題。本文為了能使算法的性能進一步提升,提出了融合多策略改進的哈里斯鷹優(yōu)化算法。引入了混沌局部搜索策略,提高了算法的開發(fā)效率,圍繞單個個體進行擾動,從而找到更好的個體。為了增加種群多樣性,提出了精英備選池策略。最后為了修改種群的進化方向,加入了分布式估計策略,能更好地對優(yōu)質(zhì)種群進行采樣獲得其權(quán)重值,把種群向更好的方向引導(dǎo)。通過CEC2017測試函數(shù)與其他5種新型群智能算法進行比較,實驗結(jié)果表明,融合了多策略的哈里斯鷹優(yōu)化算法的性能得到大幅度提高。今后準備把改進后的HHO算法運用到更多的實際工程問題,為解決各種問題提供另一種可行的思路。