一種基于椋鳥(niǎo)群行為的改進(jìn)型蝙蝠算法

胡 飛, 孫自強(qiáng)

(華東理工大學(xué)化工過(guò)程先進(jìn)控制和優(yōu)化技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200237)

一種基于椋鳥(niǎo)群行為的改進(jìn)型蝙蝠算法

胡 飛, 孫自強(qiáng)

(華東理工大學(xué)化工過(guò)程先進(jìn)控制和優(yōu)化技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200237)

蝙蝠算法是一種新興的元啟發(fā)式算法,基本蝙蝠算法(BA)存在尋優(yōu)精度低、易陷入局部最優(yōu)等缺點(diǎn)。將椋鳥(niǎo)群的集體性行為引入到基本蝙蝠算法中,有效地提高了算法的搜索范圍;引入線性遞減權(quán)重,用于平衡全局搜索和局部搜索。通過(guò)一些測(cè)試函數(shù)對(duì)該算法進(jìn)行仿真研究,結(jié)果表明改進(jìn)的蝙蝠算法有效地避免了種群個(gè)體陷入局部最優(yōu),提高了算法的尋優(yōu)精度,優(yōu)化效果得到改善。

蝙蝠算法(BA); 椋鳥(niǎo)群行為; 權(quán)重; 局部最優(yōu)

蝙蝠算法是一種新興的元啟發(fā)優(yōu)化算法,該算法可以看作是和聲搜索算法和粒子群算法的混合[1]。目前,已有許多經(jīng)典的優(yōu)化算法應(yīng)用于諸多方面。文獻(xiàn)[2]將一種改進(jìn)的混沌蜂群算法應(yīng)用于流水車(chē)間調(diào)度,優(yōu)化調(diào)度過(guò)程中的最大完成時(shí)間,優(yōu)化效果較好。文獻(xiàn)[3]針對(duì)化工過(guò)程的動(dòng)態(tài)多目標(biāo)問(wèn)題,利用一種改進(jìn)型粒子群算法對(duì)其進(jìn)行優(yōu)化,測(cè)試結(jié)果表明該算法能較好地解決工程問(wèn)題。蝙蝠算法作為新興的優(yōu)化算法,同樣具有較好的應(yīng)用前景,已經(jīng)廣泛應(yīng)用于工程和自然科學(xué)領(lǐng)域,如函數(shù)優(yōu)化問(wèn)題、生產(chǎn)調(diào)度問(wèn)題以及多目標(biāo)優(yōu)化問(wèn)題等。文獻(xiàn)[4]將蝙蝠算法應(yīng)用于一些工程函數(shù)尋優(yōu),結(jié)果表明該算法尋優(yōu)精度較高。文獻(xiàn)[5]將蝙蝠算法應(yīng)用于混合流水車(chē)間調(diào)度問(wèn)題,能夠最小化車(chē)間調(diào)度的完工時(shí)間。文獻(xiàn)[6]針對(duì)多目標(biāo)優(yōu)化問(wèn)題提出了一種多目標(biāo)蝙蝠算法,應(yīng)用于一些多目標(biāo)函數(shù)優(yōu)化,尋優(yōu)效果較佳。然而蝙蝠算法亦具有粒子群等優(yōu)化算法的通病,即收斂速度慢、容易陷入局部最優(yōu)、搜索精度較低。為了提高蝙蝠算法的控制效果,許多研究者對(duì)蝙蝠算法開(kāi)展研究,并對(duì)其進(jìn)行改進(jìn)。劉長(zhǎng)平等[7]提出基于Lévy飛行特征的蝙蝠算法,從算法仿生原理入手,采用Lévy飛行搜索策略,更真實(shí)地模擬蝙蝠的捕食行為。Topal 等[8]提出一種動(dòng)態(tài)虛擬蝙蝠算法,將蝙蝠種群分成探索類(lèi)和開(kāi)拓類(lèi)子種群,當(dāng)探索類(lèi)蝙蝠進(jìn)行對(duì)外搜索最優(yōu)值時(shí),開(kāi)拓類(lèi)蝙蝠進(jìn)行局部搜索,尋找最優(yōu)值。李枝勇等[9]將量子計(jì)算引入到蝙蝠算法中,利用量子位對(duì)蝙蝠的位置進(jìn)行編碼,利用量子旋轉(zhuǎn)門(mén)進(jìn)行蝙蝠最優(yōu)位置的搜索。Gandomi 等[10]提出具有混沌搜索策略的蝙蝠算法,并利用不同的混沌函數(shù)產(chǎn)生的混沌序列對(duì)精英個(gè)體進(jìn)行混沌優(yōu)化,驗(yàn)證了混沌策略的有效性。Niknam 等[11]提出了一種自適應(yīng)蝙蝠算法,將自適應(yīng)方法與蝙蝠算法結(jié)合,提高算法的種群多樣性及搜索能力,并應(yīng)用于電力系統(tǒng)中UC問(wèn)題,能夠快速解決UC問(wèn)題。Wang等[12]提出將差分進(jìn)化算法與蝙蝠算法相結(jié)合,并用于解決UCVA三維路徑規(guī)劃問(wèn)題,結(jié)果表明控制效果較好。

當(dāng)前針對(duì)蝙蝠算法的改進(jìn)主要有兩個(gè)方面:(1)針對(duì)蝙蝠算法自身機(jī)理存在的變量進(jìn)行改進(jìn),包括參數(shù)自適應(yīng)、輔助參數(shù)的引入以及搜索策略等;(2)將其他的優(yōu)化算法與蝙蝠算法相結(jié)合,提高基本蝙蝠算法的控制效果[13]。本文提出的改進(jìn)策略屬于第1種,將椋鳥(niǎo)群行為引入到蝙蝠算法中,當(dāng)全局最優(yōu)值陷入局部最優(yōu)時(shí),根據(jù)椋鳥(niǎo)群的集體性行為對(duì)蝙蝠的位置和速度進(jìn)行更新。

1 基本蝙蝠算法

1.1基本蝙蝠算法的位置和速度更新

將蝙蝠類(lèi)比為在當(dāng)前可行域內(nèi)分布的搜索點(diǎn),用目標(biāo)函數(shù)值(適應(yīng)度函數(shù)值)的大小來(lái)判斷當(dāng)前蝙蝠所處位置的優(yōu)劣;將蝙蝠飛行移動(dòng)探測(cè)目標(biāo)的過(guò)程類(lèi)比為搜索目標(biāo)函數(shù)位置最優(yōu)解的過(guò)程,則可用以下流程和公式來(lái)模擬蝙蝠回聲定位的行為[14]。初始時(shí)刻,需對(duì)蝙蝠的頻率進(jìn)行初始化,同時(shí),蝙蝠在飛行過(guò)程中也需對(duì)其速度和位置進(jìn)行更新,兩者的更新公式如式(1)～ 式(3)所示。

(1)

(2)

(3)

在搜索過(guò)程中還需進(jìn)行局部搜索,公式如下:

(4)

其中:xnew為蝙蝠更新后的位置;xold為蝙蝠更新前的位置;ε∈[-1,1]為隨機(jī)數(shù);At為整個(gè)蝙蝠種群在t時(shí)刻的平均響度。

1.2響度和脈沖發(fā)射率

在蝙蝠搜尋目標(biāo)獵物時(shí),需對(duì)響度Ai和脈沖發(fā)射率ri進(jìn)行更新。初始階段,Ai較大,ri較小,搜尋到目標(biāo)獵物后,Ai減小,ri增大[15]。Ai和ri更新公式如下:

(5)

(6)

2 改進(jìn)的蝙蝠算法

2.1椋鳥(niǎo)群行為

椋鳥(niǎo)是一種群體性種群,在日常生活中會(huì)產(chǎn)生一些壯觀的景象。當(dāng)遇到天敵時(shí),其中一些個(gè)體會(huì)改變飛行方向,并影響周?chē)鷤€(gè)體,如此使得整個(gè)種群方向迅速改變,從而躲避天敵的捕食。

由于每個(gè)個(gè)體都會(huì)與周?chē)膫€(gè)體交流信息,因此由該個(gè)體發(fā)出的信息會(huì)很快傳播至整個(gè)種群,這種行為稱為椋鳥(niǎo)群行為。在椋鳥(niǎo)群中,椋鳥(niǎo)個(gè)體一般與其附近的7個(gè)個(gè)體進(jìn)行信息交流,這7個(gè)個(gè)體再分別與其附近的7個(gè)個(gè)體進(jìn)行信息交流,進(jìn)而遍布整個(gè)種群,達(dá)到整個(gè)種群的信息共享。一項(xiàng)分析表明,6個(gè)或7個(gè)個(gè)體組成一個(gè)群體交流網(wǎng)絡(luò)將有效平衡整個(gè)種群凝聚性及個(gè)體交流,因此個(gè)體與其附近的7個(gè)個(gè)體進(jìn)行信息交流足以影響整個(gè)種群的行為[16]。

2.2引入椋鳥(niǎo)群行為的蝙蝠算法[17](SFBA)

基本蝙蝠算法在搜索過(guò)程中易陷入局部最優(yōu),缺乏種群多樣性。將椋鳥(niǎo)群行為引入基本蝙蝠算法中,有助于增加種群多樣性,避免陷入局部最優(yōu)。當(dāng)種群迭代結(jié)果發(fā)生停滯或更差超過(guò)上限值(count_limit)時(shí),挑選出一定數(shù)目(max_num)適應(yīng)度值較差的個(gè)體,對(duì)這些個(gè)體的位置和速度進(jìn)行變換,挑選出全局最優(yōu)值,替換原先的種群全局最優(yōu)值。

針對(duì)位置變換,依據(jù)式(7)對(duì)粒子i的位置進(jìn)行更新。

(7)

針對(duì)速度變換,依據(jù)式(8)對(duì)粒子i的速度進(jìn)行更新。

(8)

2.3針對(duì)位置和速度引入權(quán)重

由于固定的慣性權(quán)重控制精度不高,控制效果較差,本文針對(duì)蝙蝠算法速度和位置更新公式,分別引入權(quán)重w1和w2,見(jiàn)式(9)和式(10)。將權(quán)重w1由固定權(quán)重變?yōu)榫€性遞減權(quán)重。起初權(quán)重較大,加強(qiáng)全局搜索,隨著迭代次數(shù)的增加,迭代后期,權(quán)重減小,有利于提高局部搜索能力。線性遞減的權(quán)重公式如下:

(9)

(10)

(11)

其中:wmax為w1的上限值;wmin為w1的下限值;Iter為當(dāng)前迭代次數(shù);Iter_max為最大迭代次數(shù)。

2.4SFBA算法流程

(2) 對(duì)所有蝙蝠個(gè)體計(jì)算適應(yīng)度函數(shù)值f(x),找出當(dāng)前位置最優(yōu)解x*。

(3) 按照式(1)隨機(jī)產(chǎn)生頻率的新解。

(4) 根據(jù)式(9)和式(10)對(duì)蝙蝠的速度和位置進(jìn)行更新,針對(duì)速度設(shè)置的權(quán)重w1使用線性遞減權(quán)重,見(jiàn)式(11)。

(5) 產(chǎn)生一個(gè)隨機(jī)值rand,若rand大于脈沖發(fā)射率ri,則根據(jù)式(4),從處在最佳位置的蝙蝠i個(gè)體中選取一個(gè)解,對(duì)該位置最優(yōu)解進(jìn)行隨機(jī)擾動(dòng),產(chǎn)生一個(gè)局部解xnew替代蝙蝠i的當(dāng)前位置。

(6) 計(jì)算適應(yīng)度函數(shù)值(目標(biāo)函數(shù)值);

(7) 若隨機(jī)值rand小于響度Ai且更新位置后的目標(biāo)函數(shù)值小于當(dāng)前目標(biāo)函數(shù)最優(yōu)值,則接受這個(gè)位置新解,并增大脈沖發(fā)射率ri和響度Ai,見(jiàn)式(5)、式(6)。

(8) 對(duì)位置變動(dòng)后的蝙蝠重新排列,找出當(dāng)前位置最優(yōu)解和適應(yīng)度函數(shù)最優(yōu)值。

(9) 判斷當(dāng)前適應(yīng)度最優(yōu)值是否大于或等于上一代適應(yīng)度最優(yōu)值,若是,則計(jì)數(shù)器count+1(count初始值設(shè)為0)。

(10) 判斷count是否大于上限值count_limit,若是,則挑選出適應(yīng)度較差的max_num個(gè)個(gè)體,利用式(7)與式(8)更新該部分個(gè)體位置和速度值,替換原先個(gè)體最優(yōu)值,并從中找出全局最優(yōu)值,作為當(dāng)前種群全局最優(yōu)值,并置count=0。

(11) 當(dāng)?shù)螖?shù)小于最大迭代次數(shù)時(shí),回到步驟(2)繼續(xù)迭代搜索最優(yōu)解。

(12) 算法迭代結(jié)束,輸出位置全局最優(yōu)解。

3 算法的參數(shù)選擇

SFBA算法參數(shù)主要包含w1,w2,count_limit,max_num,其中w1主要由wmax和wmin決定,因此需要選擇的參數(shù)主要為wmax,wmin,w2,count_limit,max_num,這5個(gè)參數(shù)可以通過(guò)實(shí)驗(yàn)確定。以Sphere函數(shù)為基準(zhǔn)進(jìn)行實(shí)驗(yàn),種群規(guī)模設(shè)為100,迭代次數(shù)Iter_max為500次,分別運(yùn)行20次,取20次優(yōu)化結(jié)果的平均值作為參考值,確定權(quán)值wmax,wmin以及w2的最優(yōu)值。

表1示出了wmax為0.45、0.46、0.47、0.48、0.49、0.50、0.51、0.52、0.53、0.54、0.55,wmin=0.1,w2=0.50時(shí)的運(yùn)行結(jié)果,可以看出,wmax=0.50時(shí)算法控制效果較好。

表2示出了wmin為0.05、0.06、0.07、0.08、0.09、0.10、0.11、0.12、0.13、0.14、0.15,wmax=0.50,w2=0.50時(shí)的運(yùn)行結(jié)果,可以看出,wmin=0.10時(shí)算法控制效果較好。

表3示出了w2為0.45、0.46、0.47、0.48、0.49、0.50、0.51、0.52、0.53、0.54、0.55,wmax=0.50,wmin=0.10時(shí)的運(yùn)行結(jié)果,可以看出,w2=0.50時(shí)算法控制效果較好。

表1 wmax參數(shù)取值表Table 1 Parameter value table of wmax

表2 wmin參數(shù)取值表Table 2 Parameter value table of wmin

表3 w2參數(shù)取值表Table 3 Parameter value table of w2

表4示出了count_limit為1、2、3、4、5時(shí)的運(yùn)行結(jié)果,可以看出,count_limit=3時(shí)算法控制效果較好。

表5示出了max_num為13、17、19、22、23時(shí)的運(yùn)行結(jié)果,可以看出,max_num=19時(shí)算法控制效果較好。

表4 count_limit參數(shù)取值表Table 4 Parameter value table of count_limit

表5 max_num參數(shù)取值表Table 5 Parameter value table of max_num

綜上所述,對(duì)于改進(jìn)型蝙蝠算法,取wmax=0.50,wmin=0.10,w2=0.50,count_limit=3,max_num=19,采用這組參數(shù)進(jìn)行算法測(cè)試,可以進(jìn)一步驗(yàn)證改進(jìn)型蝙蝠算法的可行性。

4 測(cè)試結(jié)果

為了評(píng)估本文提出的改進(jìn)型蝙蝠算法的尋優(yōu)效果,選擇5個(gè)標(biāo)準(zhǔn)測(cè)試函數(shù)(如表6所示)進(jìn)行測(cè)試。

表6 標(biāo)準(zhǔn)測(cè)試函數(shù)Table 6 Benchmark functions

表7 標(biāo)準(zhǔn)測(cè)試函數(shù)尋優(yōu)結(jié)果對(duì)比Table 7 Comparison of Benchmark functions optimization results

測(cè)試函數(shù)f1,f2,f3均為多峰值測(cè)試函數(shù),BA尋優(yōu)精度較低,控制效果較差,而SFBA均能取得較高的尋優(yōu)精度,控制效果較好。測(cè)試函數(shù)f4是一種典型的病態(tài)非凸單模函數(shù),難以極小化,從表7中數(shù)據(jù)可以看出BA和SFBA尋優(yōu)效果均不理想。測(cè)試函數(shù)f5為簡(jiǎn)單的單峰值函數(shù),在自變量取值范圍內(nèi)僅有一個(gè)全局極小值點(diǎn),因而B(niǎo)A的尋優(yōu)效果還可以,但是相比之下,SFBA的尋優(yōu)精度更高,控制效果更好。限于篇幅,本文僅給出部分維度的尋優(yōu)效果,測(cè)試函數(shù)迭代曲線對(duì)比如圖1～圖5所示。

圖1 測(cè)試函數(shù)f1迭代曲線對(duì)比(d=20)Fig.1 Comparison of Benchmark functionf1′s iterative curves (d=20)

圖2 測(cè)試函數(shù)f2迭代曲線對(duì)比(d=20)Fig.2 Comparison of Benchmark functionf2′s iterative curves (d=20)

為了進(jìn)一步驗(yàn)證本文算法的可行性,將SFBA與文獻(xiàn)[18]基于微分算子與Lévy飛行策略的蝙蝠算法(DLBA)及文獻(xiàn)[19]基于遺傳交叉因子的蝙蝠算法(GHBA)進(jìn)行對(duì)比,結(jié)果如表8所示。從表8中數(shù)據(jù)可知,SFBA的尋優(yōu)效果優(yōu)于GHBA、DLBA。

圖3 測(cè)試函數(shù)f3迭代曲線對(duì)比(d=10)Fig.3 Comparison of Benchmark functionf3′s iterative curves (d=10)

圖4 測(cè)試函數(shù)f4迭代曲線對(duì)比(d=10)Fig.4 Comparison of Benchmark functionf4′s iterative curves (d=10)

圖5 測(cè)試函數(shù)f5迭代曲線對(duì)比(d=10)Fig.5 Comparison of Benchmark functionf5′s iterative curves (d=10)表8 改進(jìn)的BA算法在5個(gè)測(cè)試函數(shù)尋優(yōu)結(jié)果對(duì)比Table 8 Comparison of improved BA algorithms in 5 Benchmark function optimization results

函數(shù)SFBAGHBADLBA最優(yōu)值平均值最優(yōu)值平均值最優(yōu)值平均值f12.8040×10-121.1585×10-91.31×10-73.41×10-24.7174×10-76.4097×10-6f2003.65×10-71.41×10-51.3818×10-75.0103×10-6f300005.9998×10-84.4705×10-6f41.8935×101.8991×107.64×10-51.515.66727.3903f54.2686×10-222.7087×10-188.08×10-195.65×10-184.3057×10-84.8036×10-6

5 改進(jìn)型蝙蝠算法在PID參數(shù)整定中的應(yīng)用

為了進(jìn)一步驗(yàn)證本文算法的有效性,將改進(jìn)型蝙蝠算法應(yīng)用于PID控制器參數(shù)優(yōu)化中[20]。選取二階被控對(duì)象,傳遞函數(shù)如下:

(12)

采樣時(shí)間為2 ms,輸入一個(gè)階躍信號(hào),采用誤差絕對(duì)值時(shí)間積分性能指標(biāo)作為參數(shù)選擇的最小目標(biāo)函數(shù),為防止能量過(guò)大,在目標(biāo)函數(shù)中加入控制輸入平方項(xiàng),最優(yōu)指標(biāo)如式(13)所示。

(13)

式中:J表示目標(biāo)函數(shù);e(t)表示系統(tǒng)誤差;u(t)為控制器輸出;tu為上升時(shí)間;w1,w2,w3為權(quán)重值。

為避免超調(diào),采取懲罰功能,一旦出現(xiàn)超調(diào),就將超調(diào)量作為最優(yōu)指標(biāo)中的一項(xiàng),此時(shí)最優(yōu)指標(biāo)如式(14)所示。

ifey(t)<0

(14)

其中:w4為權(quán)重值,且w4>>w1;ey(t)=y(t)-y(t-1),y(t)為被控對(duì)象輸出。

SFBA算法中,種群數(shù)為100,wmax=0.50,wmin=0.10,w2=0.5,count_limit=3,max_num=19,PID參數(shù)Kp,Ki,Kd取值范圍分別為[0,40],[0,4],[0,4],取w1=0.999,w2=0.001,w3=100,w4=2.0。迭代100次后,用BA算法優(yōu)化PID三參數(shù)優(yōu)化后Kp=31.711 6,Ki=3.354 7,Kd=0.147 2,性能指標(biāo)J=79.189 6,適應(yīng)度函數(shù)J、Kp、Ki、Kd的優(yōu)化過(guò)程如圖6所示。用SFBA算法優(yōu)化PID三參數(shù)優(yōu)化后Kp=35.053 4,Ki=3.527 6,Kd=0.003 1,性能指標(biāo)J=78.352 7,適應(yīng)度函數(shù)J,Kp,Ki,Kd的優(yōu)化過(guò)程如圖7所示。BA和SFBA算法優(yōu)化后,PID控制階躍響應(yīng)如圖8所示。

由圖8可知,SFBA算法優(yōu)化效果優(yōu)于BA算法,SFBA基本無(wú)超調(diào),調(diào)節(jié)時(shí)間較短,響應(yīng)速度較快。

圖6 BA算法迭代曲線Fig.6 Iterative curves of BA algorithm

圖7 SFBA算法的迭代曲線Fig.7 Iterative curves of SFBA algorithm

圖8 BA算法與SFBA算法的階躍響應(yīng)曲線Fig.8 Unit-step response of BA algorithm and SFBA algorithm

6 結(jié) 論

本文針對(duì)BA算法收斂速度慢、容易陷入局部最優(yōu)、搜索精度不高等缺點(diǎn),提出了一種基于椋鳥(niǎo)群行為的蝙蝠算法(SFBA),同時(shí)引入線性遞減權(quán)重平衡全局搜索與局部搜索。通過(guò)5個(gè)不同的單峰值和多峰值測(cè)試函數(shù)進(jìn)行驗(yàn)證,與BA算法相比,SFBA算法優(yōu)化效果較佳。此外,將SFBA算法應(yīng)用于PID參數(shù)優(yōu)化,結(jié)果表明SFBA算法優(yōu)化效果優(yōu)于BA算法。

[1] YANG X S.A new metaheuristic bat-inspired algorithm[J].Physics,2010,284:65-74.

[2] 劉華,顧幸生.改進(jìn)的混沌蜂群算法在流水線調(diào)度中的應(yīng)用[J].華東理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,39(3):345-350.

[3] 王珊珊,杜文莉,陳旭,等.基于約束骨干粒子群算法的化工過(guò)程動(dòng)態(tài)多目標(biāo)優(yōu)化[J].華東理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,40(4):449-457.

[4] YANG X S,GANDOMI A H.Bat algorithm:A novel approach for global engineering optimization[J].Engineering Computation,2012,29(5):267-289.

[5] MARICHELVAM M K,PRABAHARAN T,YANG X S,etal.Solving hybrid flow shop scheduling problems using bat algorithm[J].International Journal of Logistics Economics & Globalisation,2013,5(1):15-29.

[6] YANG X S.Bat algorithm for multi-objective optimization[J].International Journal of Bio-Inspired Computation,2011,3(5):267-274.

[7] 劉長(zhǎng)平,葉春明.具有 Lévy飛行特征的蝙蝠算法[J].智能系統(tǒng)學(xué)報(bào),2013(3):240-246.

[8] TOPAL A O,ALTUN O.A novel meta-heuristic algorithm:Dynamic virtual bats algorithm[J].Information Sciences,2016,354:222-235.

[9] 李枝勇,馬良,張惠珍.函數(shù)優(yōu)化的量子蝙蝠算法[J].系統(tǒng)管理學(xué)報(bào),2014(5):717-722.

[10] GANDOMI A H,YANG X S.Chaotic bat algorithm[J].Journal of Computational Science,2014,5(2):224-232.

[11] NIKNAM T,BAVAFA F,AZIZIPANAH-ABARGHOOEE R.New self-adaptive bat-inspired algorithm for unit comi--tment problem[J].IET Science Measurement Technology,2014,8(6):505-517.

[12] WANG G G,CHU H C E,MIRJALILI S.Three-dimensional path planning for UCAV using an improved bat algorithm[J].Aerospace Science & Technology,2016,49:231-238.

[13] AFRABANDPEY H,GHAFFARI M,MIRZAEI A,etal.A novel bat algorithm based on chaos for optimization tasks[C]//Iranian Conference on Intelligent Systems.USA:IEEE,2014:1-6.

[14] RAGHAVAN S,MARIMUTHU C,SARWESH P,etal.Bat algorithm for scheduling workflow applications in cloud[C]//International Conference on Electronic Design,Computer Networks & Automated Verification.USA:IEEE,2015:139-144.

[16] YOUNG G F,SCARDOVI L,CAVAGNA A,etal.Starling flock networks manage uncertainty in consensus at low cost[J].PloS Computational Biology,2013,9(1):e1002894.

[17] NETJINDA N,ACHALAKUL T,SIRINAOVAKUL B.Particle swarm optimization inspired by starling flock behavior[J].Applied Soft Computing,2015,35(C):411-422.

[18] XIE J,ZHOU Y,CHEN H.A novel bat algorithm based on differential operator and Lévy flights trajectory[J].Computational Intelligence & Neuroscience,2013(2013):453-812.

[19] 彭泓,丁玉成.基于遺傳交叉因子的蝙蝠算法的改進(jìn)[J].激光雜志,2015(2):23-26.

[20] 劉金琨.先進(jìn)PID控制MATLAB仿真[M].北京:電子工業(yè)出版社,2004.

AnImprovedBatAlgorithmBasedonStarlingFlockBehavior

HUFei,SUNZi-qiang

(KeyLaboratoryofAdvancedChemicalProcessControlandOptimizationTechnology,MinistryofEducation,EastChinaUniversityofScienceandTechnology,Shanghai200237,China)

Bat algorithm (BA) is a new metaheuristic algorithm.However,the standard BA has some shortcomings,e.g.,low convergence precision and easily relapsing into the local optima.In this work,by introducing the collective behavior of the starling group into BA algorithm,the searching range of the standard BA algorithm can be effectively improved.Besides,a linear decreasing weight is introduced to balance the global search and the local search.Simulation results from Benchmark functions show that the improved algorithm can effectively avoid the local optimum and attain higher convergence precision.

bat algorithm (BA); starling group behavior; weight; local optima

1006-3080(2017)04-0525-08

10.14135/j.cnki.1006-3080.2017.04.011

2016-11-01

胡 飛(1992-),男,安徽人,碩士生,研究方向?yàn)橹悄軆?yōu)化算法。 E-mail:xueshanfeihu1992@163.com

孫自強(qiáng),E-mail:sunziqiang@ecust.edu.cn

TP301.6

A