余玲珍 楊雪梅

（貴陽(yáng)信息科技學(xué)院智能工程系，貴州 貴陽(yáng) 550025）

群體智能優(yōu)化算法是以自然界生物覓食、避障等行為方式為靈感創(chuàng)造的一類啟發(fā)式方法，其包括粒子群優(yōu)化算法[1]、蝴蝶優(yōu)化算法[2]和灰狼優(yōu)化算法[3]等。由于蝙蝠算法[4]（BA）具有結(jié)構(gòu)簡(jiǎn)單、收斂速度快和參數(shù)少等優(yōu)點(diǎn)，因此廣泛應(yīng)用于無線傳感器網(wǎng)絡(luò)定位[5]、圖像分割[6]和路徑規(guī)劃[7]等。作為一種有效解決復(fù)雜優(yōu)化問題的啟發(fā)式算法，雖然BA 算法廣泛應(yīng)用于實(shí)際問題中，但是該算法執(zhí)行后期存在尋優(yōu)精度不足、局部搜索能力較差的缺點(diǎn)。針對(duì)上述缺點(diǎn)，文獻(xiàn)[8]引入開關(guān)函數(shù)來控制蝙蝠個(gè)體有序發(fā)生變異操作，并將均勻變異和高斯變異加速算法定位到全局最優(yōu)解區(qū)域。文獻(xiàn)[9]將幾種邊界變異策略進(jìn)行比較，并提出利用越界重置策略對(duì)飛越解空間的蝙蝠位置進(jìn)行重新分配。文獻(xiàn)[10]針對(duì)BA 算法速度更新公式的不足，在速度更新公式中引入慣性權(quán)重因子來改變速度更新的方向，便于算法跳出局部最優(yōu)。

該文提出了一種基于混沌映射與高斯擾動(dòng)的蝙蝠優(yōu)化算法（TGBA），該算法采用Tent 映射和高斯擾動(dòng)策略對(duì)標(biāo)準(zhǔn)BA 算法進(jìn)行了改進(jìn)，并將測(cè)試結(jié)果與BA 算法、基于慣性權(quán)重的BA 算法（IWBA）[10]和新型BA 算法（IBA）[11]進(jìn)行了比較。

1 蝙蝠算法

蝙蝠算法是受微型蝙蝠利用回聲定位系統(tǒng)覓食行為的啟發(fā)提出的一種新型群智能優(yōu)化算法[4]。在搜索過程中，蝙蝠通過相互傳遞各自的信息來尋找群體的最優(yōu)解。在整個(gè)解空間中，蝙蝠種群初始化、蝙蝠頻率、速度和位置更新計(jì)算如公式（1）～公式（4）所示。

式中：yi，j為第i只蝙蝠在第j維搜索下的空間位置；i=1，2，…，n；j=1，2，…，d；ymin，j為j維搜索的下邊界；rand為隨機(jī)數(shù)，rand∈[0，1]；ymax，j為j維搜索的上邊界。

式中：fi為第i只蝙蝠脈沖的當(dāng)前頻率值；fmin為脈沖頻率最小值；fmax為脈沖頻率最大值；α為[0，1]中服從均勻分布的隨機(jī)數(shù)；vti和yti分別為第i只蝙蝠在t時(shí)刻的速度和位置；y*為t時(shí)刻全局搜索過程中的最優(yōu)位置。

當(dāng)蝙蝠個(gè)體逐漸靠近全局最優(yōu)解時(shí)，采用局部搜索機(jī)制，局部位置更新公式如公式（5）所示。

式中：ynew為局部解；η∈[-1，1]的隨機(jī)數(shù)；At為所有蝙蝠在t時(shí)刻發(fā)出聲波的平均響度。

響度和脈沖發(fā)射率隨迭代次數(shù)增加而變化。當(dāng)蝙蝠向獵物靠近時(shí)，響度A的值降低，脈沖發(fā)射率r增加。響度和脈沖發(fā)射率更新方式如公式（6）所示。

式中：為第i只蝙蝠在t+1時(shí)刻發(fā)出的響度；為第i只蝙蝠在t+1時(shí)刻發(fā)出的脈沖發(fā)射率；ri0為初始時(shí)刻蝙蝠個(gè)體發(fā)出的脈沖發(fā)射率；響度衰減系數(shù)β∈[0，1]；脈沖頻度增加系數(shù)γ＞0。

2 基于混沌序列與高斯擾動(dòng)的蝙蝠優(yōu)化算法

2.1 基于混沌序列的種群初始化

種群的全局最優(yōu)解是一個(gè)未知量，并且無法提前確定具體位置。由于蝙蝠算法的種群初始化是通過隨機(jī)函數(shù)完成的，因此導(dǎo)致蝙蝠的覆蓋率具有一定的隨機(jī)性，影響蝙蝠算法的整體求解性能。采用混沌映射函數(shù)生成混沌序列的方法對(duì)種群初始化進(jìn)行改進(jìn)，有利于提高蝙蝠個(gè)體在初始解空間的覆蓋率。

分別采用Logistic 映射（LBA）、Chebyshev 映射（CBA）以及Tent 映射（TBA）3 種混沌映射函數(shù)生成不同的混沌序列，對(duì)蝙蝠種群進(jìn)行初始化操作，并利用Sphere函數(shù)測(cè)試3 種映射方法優(yōu)化蝙蝠算法后的尋優(yōu)效果，仿真結(jié)果如圖1所示。

圖1 基于混沌映射的蝙蝠算法尋優(yōu)結(jié)果

根據(jù)圖1 可知，利用混沌序列對(duì)蝙蝠種群進(jìn)行優(yōu)化能夠明顯提高算法的尋優(yōu)精度和收斂速度。與Chebyshev 映射、Logistic 映射相比，Tent 映射具有更快的收斂速度和更高的尋優(yōu)精度。因此，該文采用Tent 映射方法來進(jìn)行種群的初始化。

2.2 高斯擾動(dòng)

標(biāo)準(zhǔn)蝙蝠算法的局部搜索機(jī)制是在當(dāng)前的解空間中搜索最優(yōu)解，采用隨機(jī)游走的方式產(chǎn)生新的局部解。該搜索機(jī)制主要通過對(duì)平均響度進(jìn)行調(diào)整，從而產(chǎn)生局部新解，但是這種搜索方式對(duì)當(dāng)前全局最優(yōu)個(gè)體攜帶的信息利用率并不高，在蝙蝠種群中某只蝙蝠一旦找到局部最優(yōu)值，就會(huì)吸引其他蝙蝠向其靠攏，從而導(dǎo)致整個(gè)算法進(jìn)入早熟收斂，優(yōu)化性能下降。為了克服標(biāo)準(zhǔn)蝙蝠算法中局部搜索的缺點(diǎn)，該文采用高斯擾動(dòng)策略對(duì)當(dāng)前全局最優(yōu)位置進(jìn)行擾動(dòng)，改進(jìn)后的局部搜索機(jī)制如公式（7）所示。

式中：Gauss為數(shù)學(xué)期望為0、方差為1 的高斯分布。

2.3 TGBA 算法實(shí)現(xiàn)流程

基于混沌映射與高斯擾動(dòng)的蝙蝠優(yōu)化算法（TGBA）的具體流程如下：1）初始化種群大小N、頻率fi、脈沖發(fā)射率ri和響度Ai等參數(shù)，引入混沌映射函數(shù)初始化蝙蝠種群位置yi（i=1，2，...d），并求解當(dāng)前空間最優(yōu)值fmin和最優(yōu)解y*。2）蝙蝠個(gè)體根據(jù)公式（2）～公式（4）更新速度和位置。3）利用rand生成1 個(gè)隨機(jī)數(shù)，如果rand＞ri，那么利用公式（7）的機(jī)制對(duì)當(dāng)前最優(yōu)解進(jìn)行擾動(dòng)，產(chǎn)生一個(gè)新的局部解ynew。4）將局部解ynew代入目標(biāo)函數(shù)，計(jì)算其適應(yīng)度值Fnew。5）生成1 個(gè)隨機(jī)數(shù)，如果rand＜Ai且Fnew＜f（y*），那么將步驟3 產(chǎn)生的局部新解記作當(dāng)前最優(yōu)解y*，按照公式（6）更新脈沖發(fā)射率ri和響度Ai。6）如果Fnew＜fmin，更新最優(yōu)值fmin和最優(yōu)解y*。7）重復(fù)步驟2～步驟6，直到算法達(dá)到設(shè)定的最大迭代次數(shù)Tmax。8）輸出最優(yōu)目標(biāo)函數(shù)適應(yīng)度值fmin和全局最優(yōu)解y*。

3 仿真試驗(yàn)

為驗(yàn)證改進(jìn)算法（TGBA）的尋優(yōu)性能，選取10 個(gè)基準(zhǔn)測(cè)試函數(shù)進(jìn)行仿真測(cè)試，并與BA 算法、IBA 算法和IWBA算法進(jìn)行比較。

3.1 試驗(yàn)環(huán)境與參數(shù)設(shè)置

該文仿真測(cè)試環(huán)境的硬件配置為Windows 10 操作系統(tǒng)，Intel（R）Core（TM）i52.50GHz 處理器，內(nèi)存為8G，軟件配置為MATLAB 2017b。參數(shù)統(tǒng)一設(shè)置如下：種群大小N=40，響度A=1，脈沖發(fā)射率r=0.75，最大頻率fmax=1，最小頻率fmin=-1，最大迭代次數(shù)Tmax=1 000。

3.2 測(cè)試函數(shù)

該文選取10 個(gè)測(cè)試函數(shù)，測(cè)試維度為30/100，分別給出測(cè)試函數(shù)名、對(duì)應(yīng)的函數(shù)表達(dá)式、搜索范圍以及理論最優(yōu)值，其中f1（x）、f2（x）、f6（x）、f7（x）、f8（x）和f9（x）為單峰函數(shù)，可用于檢測(cè)算法的收斂能力。f3（x）、f4（x）、f5（x）和f10（x）為復(fù)雜的非線性多峰函數(shù)，存在大量的局部極值點(diǎn)，用來檢驗(yàn)算法的種群多樣性、全局搜索能力以及逃離局部最優(yōu)等性能。

第一，Sphere 函數(shù)，如公式（8）所示。

式中：xi為第i維分量；n為維度。

f1函數(shù)搜索范圍為[-100，100]，理論最優(yōu)值為0。第二，Schwefel2.22 函數(shù)，如公式（9）所示。

f2函數(shù)搜索范圍為[-10，10]，理論最優(yōu)值為0。

第三，Ackley函數(shù)，如公式（10）所示。

f3函數(shù)搜索范圍為[-32，32]，理論最優(yōu)值為0。

第四，Rastrigin函數(shù)，如公式（11）所示。

f4函數(shù)搜索范圍為[-5.12，5.12]，理論最優(yōu)值為0。

第五，Griewank函數(shù)，如公式（12）所示。

f5函數(shù)搜索范圍為[-600，600]，理論最優(yōu)值為0。

第六，Rosenbrock函數(shù)，如公式（13）所示。

f6函數(shù)搜索范圍為[-2.048，2.048]，理論最優(yōu)值為0。

第七，Noise函數(shù)，如公式（14）所示。

f7函數(shù)搜索范圍為[-1.28，1.28]，理論最優(yōu)值為0。

第八，Quartic函數(shù)，如公式（15）所示。

f8函數(shù)搜索范圍為[-1.28，1.28]，理論最優(yōu)值為0。

第九，Sumsquares函數(shù)，如公式（16）所示。

f9函數(shù)搜索范圍為[-10，10]，理論最優(yōu)值為0。

第十，Alpine函數(shù)，如公式（17）所示。

f10函數(shù)搜索范圍為[-10，10]，理論最優(yōu)值為0。

3.3 尋優(yōu)性能分析

為了減少隨機(jī)因素對(duì)算法尋優(yōu)性能的影響，對(duì)目標(biāo)函數(shù)進(jìn)行求解時(shí)分別將每種算法獨(dú)立運(yùn)行30 次，記錄最優(yōu)值、最差值、平均值和方差，并對(duì)求解精度最高的解加粗。搜索維度為30 時(shí)的仿真結(jié)果比較見表1。

表1 4 種算法對(duì)目標(biāo)函數(shù)的測(cè)試結(jié)果

由表1 可以看出，在同一搜索維度下，改進(jìn)的TGBA 算法對(duì)單峰測(cè)試函數(shù)的求解精度比BA 算法、IBA 算法和IWBA 算法更高。對(duì)于復(fù)雜的多峰函數(shù)來說，TGBA 算法能夠準(zhǔn)確找出函數(shù)適應(yīng)度最優(yōu)值。隨著測(cè)試函數(shù)維度增加，BA 算法、IBA 算法和IWBA 算法的最優(yōu)值、平均值和方差都大幅增加。而TGBA 算法的最優(yōu)值、平均值和方差基本保持不變。其中，對(duì)于函數(shù)f1（x）、f2（x）、f4（x）、f5（x）、f8（x）、f9（x）和f10（x），TGBA 算法的求解精度并沒有隨測(cè)試函數(shù)維度增加而降低，始終能夠找到目標(biāo)函數(shù)的全局最優(yōu)解。綜上所述，隨著測(cè)試函數(shù)維度增加，與BA、IBA 和IWBA 算法相比，TGBA 算法仍能保持較高的求解性能，改進(jìn)后的TGBA算法在尋優(yōu)精度和穩(wěn)定性方面具有顯著的優(yōu)越性。

3.4 收斂速度分析

比較BA、IBA、IWBA 以及TGBA 優(yōu)化算法求解不同測(cè)試函數(shù)時(shí)的適應(yīng)度值收斂曲線，并給出4 種算法求解部分測(cè)試函數(shù)的收斂曲線，如圖2～圖3 所示。

圖2 Rastrigin函數(shù)的收斂曲線（30D）

圖3 Quartic 函數(shù)的收斂曲線（30D）

根據(jù)優(yōu)化算法求解測(cè)試函數(shù)的收斂情況可知，對(duì)于不同的測(cè)試函數(shù)，BA 算法與IWBA 算法收斂速度緩慢，算法的迭代曲線在早期便呈平滑狀態(tài)，算法進(jìn)入早熟收斂，無法找出測(cè)試函數(shù)的全局最優(yōu)解。由圖2～圖3 可知，當(dāng)函數(shù)維度為30時(shí)，IBA 算法在求解f4（x）和f8（x）過程中具有較快的收斂速度，并且能夠收斂到函數(shù)f4（x）的最小值。當(dāng)函數(shù)維度增至100 時(shí)，求解維度增加，解的復(fù)雜性隨之增加，IBA 算法的尋優(yōu)性能驟降。隨著求解維數(shù)增加，IBA 算法的尋優(yōu)精度降低，收斂速度緩慢，在迭代過程中容易陷入局部最優(yōu)。而對(duì)于大多數(shù)單峰和多峰測(cè)試函數(shù)，TGBA 算法都能以接近直線的趨勢(shì)下降，收斂速度非?？?，尋優(yōu)精度也遠(yuǎn)超IBA 算法。最后,通過比較各函數(shù)的收斂曲線可知，與BA 算法和IWBA算法相比，TGBA 算法對(duì)單峰和多峰函數(shù)均能夠以較快的速度收斂到全局最優(yōu)解；與IBA 算法相比，TGBA 尋優(yōu)能力更強(qiáng)，具有更高的尋優(yōu)精度和更好的穩(wěn)定性。

4 結(jié)語(yǔ)

針對(duì)BA 算法在迭代后期收斂精度降低、易出現(xiàn)早熟收斂的問題，該文提出了基于混沌映射與高斯擾動(dòng)的蝙蝠優(yōu)化算法。利用Tent 映射產(chǎn)生混沌序列，對(duì)蝙蝠種群進(jìn)行初始化操作，提高了蝙蝠在解空間的覆蓋率，獲得了更豐富的種群。在產(chǎn)生局部新解的更新公式中對(duì)當(dāng)前最優(yōu)解添加高斯擾動(dòng)能夠提升算法的搜索性能，有效避免算法陷入局部最優(yōu)。仿真結(jié)果表明，與BA、IBA、IWBA 這3 種對(duì)照算法相比，TGBA 算法在尋優(yōu)精度、收斂速度和穩(wěn)定性方面表現(xiàn)更突出。