求解全局優(yōu)化問(wèn)題的改進(jìn)灰狼算法

張 陽(yáng), 周溪召

（上海理工大學(xué) 管理學(xué)院，上海 200093）

元啟發(fā)式算法是一類在特定的資源（時(shí)間等）約束條件下求解問(wèn)題近似最優(yōu)的優(yōu)化技術(shù)，具有簡(jiǎn)單、高效和計(jì)算成本小等優(yōu)點(diǎn)。自1983 年Kirkpatrick等提出模擬退火算法（simulated annealing algorithm,SA）[1]和1992 年Holland 提出遺傳算法（genetic algorithm, GA）[2]以來(lái)，眾多學(xué)者投入到這一研究領(lǐng)域。經(jīng)過(guò)幾十年的發(fā)展，許多元啟發(fā)式算法相繼被提出并已廣泛應(yīng)用于研究領(lǐng)域和實(shí)際場(chǎng)景，如粒子群算法（particle swarm optimization, PSO）[3]、蟻群算法（ant colony optimization, ACO）[4]、引力搜索算法（gravity search algorithm, GSA）[5]、蘑菇繁殖算法（mushroom reproduction optimization, MRO）[6]、人工蜂群算法（artificial bee colony algorithm, ABC）[7]、教學(xué)優(yōu)化算法（teaching-learning-based optimization,TLBO）[8]和哈里斯鷹算法（Harris hawks optimization,HHO）[9]等。但是，根據(jù)沒(méi)有免費(fèi)的午餐定理（NFL）[10]，沒(méi)有一種算法適用于求解所有問(wèn)題，因此，仍有必要對(duì)元啟發(fā)式算法進(jìn)行進(jìn)一步的研究。

在文獻(xiàn)[6]中，元啟發(fā)式算法可以根據(jù)不同的靈感來(lái)源而分為3 類：群智能算法、非群智能算法和物理/化學(xué)算法。2014 年Mirjalili 等[11]受自然界中灰狼領(lǐng)導(dǎo)等級(jí)制度和群體狩獵行為啟發(fā)，提出了一種新的群智能算法—灰狼算法（grey wolf optimization, GWO）?；依撬惴ㄊ且环N強(qiáng)大的優(yōu)化技術(shù)，已被成功應(yīng)用于許多領(lǐng)域。譚念等[12]提出了一種基于灰狼算法和支持向量機(jī)（GWO-SVM）的木材密度預(yù)測(cè)模型，為杉木性質(zhì)定量分析提供了理論依據(jù)；姚鵬等[13]提出了一種基于改進(jìn)的擾動(dòng)流體動(dòng)態(tài)系統(tǒng)與灰狼優(yōu)化理論的混合航路規(guī)劃算法，為復(fù)雜地形環(huán)境下的無(wú)人機(jī)三維航路規(guī)劃問(wèn)題提供了解決方案；王書(shū)芹等[14]提出了基于灰狼優(yōu)化算法的長(zhǎng)短期記憶模型（long short term memory, LSTM），解決了傳統(tǒng)LSTM 易于收斂到局部最優(yōu)的問(wèn)題。

然而，同其他元啟發(fā)式算法一樣，灰狼算法也存在著收斂速度較慢、求解精度不高、在一些復(fù)雜問(wèn)題上容易陷入局部最優(yōu)等缺點(diǎn)。為提升GWO 的性能，許多學(xué)者進(jìn)行了相關(guān)研究。Joshi等[15]提出了一種基于精英策略的增強(qiáng)灰狼算法（enhanced grey wolf optimization, EGWO）；Heidari等[16]提出了結(jié)合列維飛行和貪婪搜索機(jī)制的灰狼算法（levy-embedded GWO, LGWO）；文獻(xiàn)[17]提出了基于模糊層次算子的灰狼算法，并提出了兩種個(gè)體位置更新時(shí)的比例權(quán)重策略；Saremi 等[18]提出了基于進(jìn)化種群動(dòng)力學(xué)（evolutionary population dynamics, EPD）的灰狼算法（GWO-EPD），EPD 要求表現(xiàn)不佳的搜尋個(gè)體進(jìn)行動(dòng)態(tài)的位置調(diào)整；王敏等[19]提出了一種基于反向?qū)W習(xí)策略、非線性收斂因子和變異算法的新型灰狼算法（novel GWO,NGWO）。在GWO 中，收斂因子a在一定程度上決定了算法局部開(kāi)發(fā)過(guò)程和全局搜索過(guò)程之間的平衡。本研究首次提出了一種新的指數(shù)變化形式的收斂因子，而且與以前的所有針對(duì)a的改進(jìn)之處所不同的是：提出的收斂因子最終并不會(huì)減小為0。同時(shí)為了加快算法的收斂速度，重新定義了灰狼算法中的 α，β，δ狼，并提出了一種基于自適應(yīng)調(diào)整策略的灰狼位置更新公式。此外，針對(duì)文獻(xiàn)中提出的動(dòng)態(tài)權(quán)重策略進(jìn)行修訂，以獲得更好的求解性能和魯棒性。最終，使用了10 個(gè)典型的基準(zhǔn)測(cè)試函數(shù)驗(yàn)證了以上改進(jìn)策略的有效性。

1 灰狼優(yōu)化算法

灰狼優(yōu)化算法是一種較新的優(yōu)化技術(shù)，其主要思想是灰狼群體中的領(lǐng)導(dǎo)等級(jí)制度和群體狩獵模式。在自然界中，灰狼群體具有森嚴(yán)的社會(huì)等級(jí)制，如圖1 所示，狼群中所有個(gè)體劃分為4 個(gè)等 級(jí)： α，β，δ 及 ω。 α是 狼 群 的 領(lǐng) 袖； β僅 次 于α，其職責(zé)是輔助頭狼 α制定決策及其他群體性活動(dòng)；接下來(lái)是δ，狼群中的偵察兵、哨兵、長(zhǎng)老、狩獵者及看守者都屬于此類；等級(jí)最低的狼則為ω。除等級(jí)制度外，小組狩獵也是灰狼種群中一個(gè)有趣的社會(huì)行為。

圖1 灰狼種群社會(huì)等級(jí)結(jié)構(gòu)Fig.1 Social hierarchy structure of grey wolf population

在GWO 中，適應(yīng)度值最優(yōu)的個(gè)體被視為 α，適應(yīng)度值第二和第三的分別被定義為 β 和 δ，其余個(gè)體均為 ω。在算法實(shí)施中，以下等式用來(lái)描述灰狼種群包圍獵物的行為：

式中：t表示當(dāng)前迭代次數(shù)；X(t)表示灰狼個(gè)體；Xp(t)表示當(dāng)前的灰狼個(gè)體；D表示個(gè)體與獵物之間的距離；A和C是系數(shù)向量，分別用式（3）和式（4）來(lái)計(jì)算。

式中：r1，r2是0～1 之間的隨機(jī)向量；a表示收斂因子，可用等式a=2×(1-t/T)來(lái)計(jì)算，T為設(shè)定的最大迭代次數(shù)。

在狩獵時(shí)，GWO 假定 α，β，δ 對(duì)于獵物的潛在位置更具有感知性，所以，灰狼個(gè)體會(huì)分別判斷自身與α ，β，δ 之 間的距離Dα，Dβ，Dδ（式（5）），分別計(jì)算其朝向三者的移步距離X1，X2，X3（式（6）），并最終在向著三者的包圍圈之內(nèi)進(jìn)行移步，移動(dòng)公式為式（7）。

2 改進(jìn)的灰狼優(yōu)化算法

2.1 指數(shù)規(guī)律收斂因子調(diào)整策略

全局搜索能力和局部開(kāi)發(fā)能力是描述元啟發(fā)式算法性能的兩種通用屬性，如何有效平衡兩者也是元啟發(fā)式算法首要思考的問(wèn)題。全局搜索是指在整個(gè)搜索空間搜尋最優(yōu)解，強(qiáng)的全局搜索性可使算法有效避免陷入局部最優(yōu)，增強(qiáng)算法的尋優(yōu)能力；局部開(kāi)發(fā)能力是指在局部區(qū)域進(jìn)行精確搜索，強(qiáng)的局部開(kāi)發(fā)性能夠提升算法的求解精度、加快算法收斂速度。一般情況下，迭代的前期使用全局搜索策略，而在后期則使用局部開(kāi)發(fā)策略。在GWO 中，當(dāng)|A|>1 時(shí)，灰狼種群在整個(gè)搜索域搜尋潛在獵物；當(dāng)|A|<1 時(shí)，灰狼種群將逐漸包圍并捕獲獵物。而A的值取決于收斂因子a的變化情況。

在GWO 中，a隨迭代的更新方式采用一種線性遞減策略，可用等式a=2(1-t/T)進(jìn)行計(jì)算。在Chiu 等[20]的研究中已經(jīng)證明重要參量的不同更新策略會(huì)極大影響算法的性能，而且線性策略往往不是最有效的。因此，本文提出了一種新的基于指數(shù)規(guī)律變化的收斂因子更新方式，如式（8）所示。

如圖2 所示，線性表示GWO 中所使用的收斂因子更新方式，指數(shù)形式表示本文提出的新的收斂因子更新方式。

圖2 兩種收斂因子Fig.2 Two convergence factors

2.2 自適應(yīng)位置更新策略

在GWO 中，初始化 α，β，δ 三組解會(huì)被記錄并保留，直至在迭代的過(guò)程出現(xiàn)適應(yīng)值更好的個(gè)體代替它們?yōu)橹埂＞褪钦f(shuō)，如果在第t代，種群中未出現(xiàn)優(yōu)于記錄的 α，β，δ的解時(shí)，當(dāng)前的種群依然向著這三只狼更新位置。然而，當(dāng)這三者全部陷入局部最優(yōu)時(shí)，此時(shí)整個(gè)種群就難以尋求更好的解。可以理解為：當(dāng)狼群的決策者誤判了獵物出現(xiàn)的地方，那么，所有灰狼的包圍動(dòng)作都將是無(wú)效的，它們難以在錯(cuò)誤的地方找到獵物。

本文提出了一種新的 β 與 δ的定義方式以強(qiáng)化當(dāng)前代最優(yōu)個(gè)體的作用，從而提升算法的全局搜索能力。在算法實(shí)施中， α同GWO 一樣，但是，β與 δ被 定義為局部變量，它們是第t(t=1,2,···,T)代中的適應(yīng)度值最優(yōu)和第二的灰狼個(gè)體。

同時(shí)，為了更好地平衡算法的全局搜索和局部開(kāi)發(fā)過(guò)程，本文提出了一種新的自適應(yīng)的移步策略，其數(shù)學(xué)表達(dá)為

2.3 修訂的動(dòng)態(tài)權(quán)重策略

在文獻(xiàn)[17]中指出了GWO 位置更新公式的缺點(diǎn)，X1，X2，X3求平均值的方式無(wú)法凸顯α，β，δ三者之間的重要度，為此，作者提出了兩種新的比例權(quán)重策略，分別是加權(quán)平均策略和基于適應(yīng)度值的比例權(quán)重策略，并進(jìn)行了實(shí)驗(yàn)驗(yàn)證。文獻(xiàn)[21]提出了一種基于指導(dǎo)位置向量權(quán)值的動(dòng)態(tài)比例權(quán)重。文獻(xiàn)[22]對(duì)不同的權(quán)重策略進(jìn)行了分析與實(shí)驗(yàn)研究，并從理論上證明了動(dòng)態(tài)權(quán)重策略能高效尋優(yōu)的原因。

當(dāng)前灰狼個(gè)體到 α ，β，δ之間的距離權(quán)重

最終，灰狼位置更新公式可表示為

然而，在實(shí)際應(yīng)用中，式（10）～（12）很可能出現(xiàn)分母為0 的情況。因此，需添加一個(gè)很小的常數(shù)ε，取值為10-16。它們被修改為

結(jié)合前面的自適應(yīng)位置更新策略，最終的灰狼位置更新方式可表示為

2.4 MGWO 算法步驟

綜合使用以上3 種改進(jìn)策略，本文提出的改進(jìn)灰狼優(yōu)化算法（MGWO）步驟如下：

a.參數(shù)初始化。設(shè)定種群規(guī)模N，總迭代次數(shù)T，常數(shù)ε；初始化a，A，C。

b.在待解問(wèn)題的搜索空間內(nèi)隨機(jī)生成灰狼初始種群。

c.計(jì)算種群中所有灰狼個(gè)體的適應(yīng)值，選擇適應(yīng)度值最優(yōu)的個(gè)體，并記錄其位置為Xα。 令t=1，迭代開(kāi)始。

d.判斷算法是否滿足終止條件，若滿足，則輸出最優(yōu)解；否則，執(zhí)行步驟e。

e.根據(jù)式（8）計(jì)算收斂因子a。

f.灰狼種群按適應(yīng)度值排序，排在前兩位的位置分別被記錄為Xβ和Xδ。

g.根據(jù)式（14）～（16）計(jì)算動(dòng)態(tài)比例權(quán)重，灰狼個(gè)體根據(jù)式（17）更新位置。

h.將超出搜索空間的灰狼個(gè)體放回搜索空間。

i.更新Xα， 令t=t+1，返回步驟d。

3 實(shí)驗(yàn)與分析

為了實(shí)驗(yàn)驗(yàn)證提出的3 種改進(jìn)策略的有效性，選取了國(guó)際上通用的10 種基準(zhǔn)測(cè)試函數(shù)進(jìn)行仿真實(shí)驗(yàn)。實(shí)驗(yàn)運(yùn)行環(huán)境為Intel（R）CPU 3550M，主頻2.30 GHz，內(nèi)存4 GB，Windows 10 64 位操作系統(tǒng)，所有代碼通過(guò)Python 3.7 進(jìn)行編程。在此基礎(chǔ)上，進(jìn)行了3 組對(duì)比實(shí)驗(yàn)。為保證無(wú)偏性，所有算法的種群數(shù)N均設(shè)定為30，總迭代次數(shù)設(shè)定為500。為了排除隨機(jī)性的影響，所有實(shí)驗(yàn)獨(dú)立運(yùn)行30 次，取30 次的平均值和標(biāo)準(zhǔn)差作為算法性能的度量標(biāo)準(zhǔn)。

3.1 基準(zhǔn)測(cè)試函數(shù)

基準(zhǔn)函數(shù)的數(shù)學(xué)表達(dá)如表1 所示，包括7 個(gè)單峰測(cè)試函數(shù)（F1～F7）和3 個(gè)多峰測(cè)試函數(shù)（F8～F10）。單峰測(cè)試函數(shù)是算法求解精度和收斂速度的度量，而多峰測(cè)試函數(shù)可以很好地表示算法的全局搜索能力和避免局部最優(yōu)的能力。部分測(cè)試函數(shù)的3 維圖像如圖3 所示。

表1 基準(zhǔn)測(cè)試函數(shù)Tab.1 Benchmark function

圖3 部分基準(zhǔn)函數(shù)的二維搜尋空間Fig.3 Two-dimensional search space of a part of benchmark function

3.2 與GWO 和其他改進(jìn)灰狼算法的比較

為了測(cè)試算法的搜尋性能，將提出的策略與其他文獻(xiàn)的策略進(jìn)行比較。各比較算法及其描述如表2 所示。實(shí)驗(yàn)結(jié)果如表3 所示。表3 中列出了各算法在基準(zhǔn)測(cè)試函數(shù)上的測(cè)試結(jié)果。表格中粗體顯示的是最佳求解結(jié)果。表2 中從EGWO 到GWO-Fuzzy 的實(shí)驗(yàn)結(jié)果直接來(lái)源于文獻(xiàn)[15-19]。其中，EGWO 的實(shí)驗(yàn)環(huán)境為QT Creator 2.4.0；LGWO實(shí)驗(yàn)運(yùn)行環(huán)境為Intel（R）Core i2，主頻2 GHz，內(nèi)存4 GB，Windows 7 64 位操作系統(tǒng)，編程軟件為Matlab R2013a。

從表3 可以看出，提出的MGWO_1，MGWO_2算法在函數(shù)F1～F4 和F6～F9 上均優(yōu)于GWO 算法，而且在大多數(shù)問(wèn)題上略優(yōu)于其他文獻(xiàn)提出的改進(jìn)灰狼算法，這說(shuō)明了提出的改進(jìn)策略的有效性。同時(shí)，依據(jù)數(shù)值結(jié)果，MGWO_2 相比MGWO_1獲得更高的求解精度，這說(shuō)明灰狼個(gè)體的自適應(yīng)位置更新策略是一種更優(yōu)的改進(jìn)策略。

表2 對(duì)比算法及描述Tab.2 Comparison algorithms and its description

MGWO_3 綜合使用了指數(shù)規(guī)律收斂因子調(diào)整策略和自適應(yīng)位置更新策略，從結(jié)果可以看出，在函數(shù)F1～F4 上，MGWO_3 明顯優(yōu)于除MGWO_4外的其他對(duì)比算法。在函數(shù)F7 上，MGWO_3 排在MGWO_4 和NGWO 之后，勝過(guò)其他8 種對(duì)比算法。在函數(shù)F8 上，MGWO_3 僅次于MGWO_4，NGWO 和LGWO。在函數(shù)F9 上，MGWO_3 僅次于MGWO_4 和LGWO，優(yōu)于其他8 種算法。但是，在F6 和F10 上，MGWO_3 的結(jié)果并不理想。

MGWO_4 結(jié)合了本文提出的3 種改進(jìn)策略，從表3 可以看出，它在函數(shù)F1，F(xiàn)3，F(xiàn)8 和F10上均收斂到理論最優(yōu)值。在函數(shù)F2，F(xiàn)4，F(xiàn)7 和F9 上，其求解精度明顯優(yōu)于其他算法。然而，對(duì)于函數(shù)F5 和F6，MGWO_4 表現(xiàn)一般。

在所有測(cè)試函數(shù)上各算法30 次獨(dú)立運(yùn)行的運(yùn)行總時(shí)間列于表4。從運(yùn)行時(shí)間來(lái)看，在函數(shù)F1～F4，F(xiàn)7 和F8 上，GWO_EPD 所耗費(fèi)的時(shí)間成本最小。在函數(shù)F5 上，MGWO_2 的運(yùn)行時(shí)間要小于其他的對(duì)比算法。在函數(shù)F6，F(xiàn)9 和F10 上，GWO 的運(yùn)行時(shí)間最短，優(yōu)于其他算法。在函數(shù)F7 和F8 上，GWO_EPD 的運(yùn)行時(shí)間依然具有競(jìng)爭(zhēng)性。

表3 11 種對(duì)比算法在10 個(gè)測(cè)試函數(shù)上的實(shí)驗(yàn)結(jié)果Tab.3 Experimental results of 11 comparison algorithms on 10 test functions

表4 對(duì)比算法運(yùn)行時(shí)間Tab.4 Running time of the comparison algorithms

綜上所述，本文提出的3 種改進(jìn)策略都是有效的。綜合所有數(shù)值結(jié)果，MGWO_4 明顯優(yōu)于其他10 種對(duì)比算法，它具有很高的尋優(yōu)能力，獲得了較高的求解精度，同時(shí)從標(biāo)準(zhǔn)差來(lái)看，具有較高的魯棒性。MGWO_3，NGWO 和LGWO 僅次于MGWO_4，它們?cè)谔囟ǖ膯?wèn)題上表現(xiàn)突出，且都明顯改進(jìn)了GWO 算法的性能。但是，也應(yīng)該認(rèn)識(shí)到，MGWO_4 和MGWO_3 高效尋優(yōu)結(jié)果的代價(jià)是一定程度上時(shí)間成本的增加。而MGWO_1 和MGWO_2 可能是最佳選擇，它們對(duì)GWO 的性能提升幅度不如MGWO_4 和MGWO_3 那么大，但運(yùn)行時(shí)間無(wú)明顯變化。具體選擇哪種算法，使用者要視自己的需求情況而定，如果對(duì)待求解問(wèn)題的精度要求較高，MGWO_4 和MGWO_3 算法顯然是更好的選擇。

3.3 與其他元啟發(fā)式算法的比較

為進(jìn)一步驗(yàn)證本文改進(jìn)策略的有效性，將表現(xiàn)較佳的MGWO_3，MGWO_4 與GA，TLBO，HHO 進(jìn)行對(duì)比實(shí)驗(yàn)。其中，GA 和TLBO 是經(jīng)典的優(yōu)化算法，HHO 是該領(lǐng)域最新的研究成果。其中，對(duì)于GA，使用輪盤(pán)賭的染色體復(fù)制策略、均勻雜交策略和隨機(jī)選擇突變3 種策略來(lái)進(jìn)行編碼，染色體交配概率設(shè)置為0.8，突變概率設(shè)為0.3。其余算法均只有種群規(guī)模和總迭代次數(shù)2 個(gè)通用參數(shù)。實(shí)驗(yàn)結(jié)果如表5 所示。

從表5 中的數(shù)值結(jié)果可知，對(duì)于函數(shù)F1 和F3，MGWO_4 在500 次迭代內(nèi)收斂到理論最優(yōu)值0，遠(yuǎn)遠(yuǎn)優(yōu)于其他的優(yōu)化算法；對(duì)于函數(shù)F8 和F10，MGWO_4，TLBO 和HHO 這3 種算 法 均達(dá)到理論最優(yōu)值0，但是，從收斂曲線（圖4）可以看出，MGWO_4 的收斂速度明顯高于TLBO 和HHO，所以，它依然是最佳算法；對(duì)于函數(shù)F2 和F4，MGWO_4 分別得到了3.14E-205 和1.23E-181 的最優(yōu)解，其求解精度相當(dāng)高，遠(yuǎn)勝于其他的優(yōu)化算法；對(duì)于函數(shù)F5 和F6，HHO 算法取得了最好的尋優(yōu)結(jié)果，其次是MGWO_3 與MGWO_4；對(duì)于函數(shù)F7，MGWO_4 獲得最佳求解結(jié)果，其次是TLBO 與MGWO_3。值得注意的是，對(duì)于8 個(gè)基準(zhǔn)函數(shù)（F1，F(xiàn)2，F(xiàn)3，F(xiàn)4，F(xiàn)7，F(xiàn)8，F(xiàn)9，F(xiàn)10），MGWO_4 不僅取得了最優(yōu)的平均值，而且其標(biāo)準(zhǔn)差也是最小的，這說(shuō)明了MGWO_4 算法具有較高的魯棒性。

表6 列出了GA，TLBO，HHO，MGWO_3 和MGWO_4 在10 個(gè)測(cè)試函數(shù)上30 次運(yùn)行的總時(shí)間。從表6 中的數(shù)據(jù)可以看出，TLBO 和HHO 算法的求解時(shí)間明顯優(yōu)于其他的對(duì)比者，特別是對(duì)于單峰測(cè)試函數(shù)，這兩種算法的運(yùn)行時(shí)間很短。而MGWO_3 和MGWO_4 算法的運(yùn)行時(shí)間相對(duì)較長(zhǎng)。

6 種算法（GWO，GA，TLBO，HHO，MGWO_3 和MGWO_4）在部分測(cè)試函數(shù)上的收斂曲線如圖4 所示。由圖4 可知，MGWO_3 提升了GWO的收斂速度和求解精度，MGWO_4 具有相當(dāng)高的收斂速度和求解精度，遠(yuǎn)優(yōu)于GWO 和其他的對(duì)比優(yōu)化算法。

以上兩組對(duì)比實(shí)驗(yàn)驗(yàn)證了本文提出了的3 種改進(jìn)策略的實(shí)施效果，從運(yùn)行時(shí)間來(lái)看，MGWO_1和MGWO_2 在提升GWO 的性能的基礎(chǔ)上運(yùn)行時(shí)間無(wú)明顯變化，而MGWO_4 和MGWO_3 雖然在求解精度和收斂速度上有大幅度的提升，但帶來(lái)了程序運(yùn)行的時(shí)間成本的增加。綜上所述，可以看出本文提出的MGWO 算法依然是一種更好的優(yōu)化算法。

表6 5 種對(duì)比算法在10 個(gè)測(cè)試函數(shù)上的運(yùn)行時(shí)間Tab.6 Running time of the 5 comparison algorithms on 10 test functions

3.4 拉伸/壓縮彈簧設(shè)計(jì)優(yōu)化問(wèn)題

以上兩組對(duì)比實(shí)驗(yàn)所針對(duì)的問(wèn)題都是無(wú)約束優(yōu)化問(wèn)題，為進(jìn)一步驗(yàn)證本文改進(jìn)策略的有效性，使用約束優(yōu)化問(wèn)題繼續(xù)進(jìn)行實(shí)驗(yàn)。拉伸/壓縮彈簧設(shè)計(jì)優(yōu)化問(wèn)題是一個(gè)經(jīng)典的帶約束的工程優(yōu)化問(wèn)題，該問(wèn)題的目的是在滿足3 個(gè)設(shè)計(jì)變量的4 個(gè)約束條件的情況下求得彈簧重量的最小值。3 個(gè)設(shè)計(jì)變量分別是線徑d，平均線圈直徑D和有效線圈數(shù)N；4 個(gè)約束條件與剪應(yīng)力、彈簧顫動(dòng)頻率和最小撓度相關(guān)。圖5 描述了這一問(wèn)題和其設(shè)計(jì)變量，該問(wèn)題的數(shù)學(xué)表達(dá)如下：

設(shè)計(jì)變量

優(yōu)化函數(shù)

約束條件：

圖5 拉伸/壓縮彈簧設(shè)計(jì)優(yōu)化問(wèn)題Fig.5 Tension/compression spring design problem

設(shè)計(jì)變量范圍：

使用7 種對(duì)比算法來(lái)求解這一約束優(yōu)化問(wèn)題，算法在30 次實(shí)驗(yàn)中所獲得的最佳實(shí)驗(yàn)結(jié)果如表7所示。從表7 中的數(shù)據(jù)可以看出，MGWO_2 取得了最好的結(jié)果，優(yōu)于其他對(duì)比算法。MGWO_1 和MGWO_3 的數(shù)值結(jié)果也要優(yōu)于GWO 和其他算法。MGWO_4 雖然在無(wú)約束優(yōu)化問(wèn)題時(shí)表現(xiàn)突出，但在求解拉伸/壓縮彈簧設(shè)計(jì)優(yōu)化問(wèn)題時(shí)不如其他3 種MGWO 算法和GWO 算法。

表7 拉伸/壓縮彈簧問(wèn)題的結(jié)果比較Tab.7 Comparison of results for the tension/compression spring problem

4 結(jié)束語(yǔ)

針對(duì)基本GWO 求解精度不高、收斂速度較慢和易于陷入局部最優(yōu)的問(wèn)題，提出了一種改進(jìn)的灰狼優(yōu)化算法。首先，提出了一種收斂因子指數(shù)變化規(guī)律的修正策略，處于這種策略下，a在迭代末期并不會(huì)減小為0，實(shí)驗(yàn)結(jié)果驗(yàn)證了其有效性。其次，通過(guò)重新定義 β 和δ 而強(qiáng)化每代最佳搜尋個(gè)體的作用，進(jìn)而提出了一種灰狼個(gè)體自適應(yīng)位置更新策略。此外，引入動(dòng)態(tài)權(quán)重策略，并對(duì)其進(jìn)行了修訂。接著，針對(duì)無(wú)約束優(yōu)化問(wèn)題的對(duì)比測(cè)試實(shí)驗(yàn)證明了改進(jìn)策略的有效性，尤其是綜合使用3 種改進(jìn)策略的MGWO_4 具有較強(qiáng)的局部最優(yōu)避免性、較快的收斂速度、較高的求解精度，它明顯提升了基本GWO 的性能，更好地平衡了算法的全局搜索性和局部開(kāi)發(fā)性。最后，對(duì)于工程設(shè)計(jì)約束優(yōu)化問(wèn)題的實(shí)驗(yàn)結(jié)果進(jìn)一步證明了本文提出的MGWO 的有效性。

未來(lái)的研究思路是編寫(xiě)二進(jìn)制的MGWO，并將MGWO 運(yùn)用到更多的實(shí)際問(wèn)題中。