耿健， 張雨飛， 范赫

(東南大學(xué) 能源與環(huán)境學(xué)院，江蘇 南京 210096)

基于TB-PC算法的MISO系統(tǒng)辨識(shí)的研究

耿健， 張雨飛， 范赫

(東南大學(xué) 能源與環(huán)境學(xué)院，江蘇 南京 210096)

如何保持MISO系統(tǒng)辨識(shí)的精確度、收斂度、耗時(shí)以及跳出局部最優(yōu)解，是當(dāng)今研究的熱點(diǎn)和難點(diǎn)。提出了一種基于協(xié)同進(jìn)化策略和禁忌搜索的蝙蝠和粒子群混合算法(TB-PC)，分析蝙蝠算法控制參數(shù)，提出了合理的脈沖頻度和音強(qiáng)初值，將蝙蝠算法與粒子群算法的優(yōu)勢(shì)用協(xié)同進(jìn)化結(jié)合起來(lái)，并引入了禁忌搜索。通過(guò)對(duì)四個(gè)測(cè)試函數(shù)和MISO系統(tǒng)實(shí)例的辨識(shí)仿真，驗(yàn)證了TB-PC算法具有穩(wěn)定性能好、收斂精度高等優(yōu)點(diǎn)，對(duì)MISO系統(tǒng)有優(yōu)良的辨識(shí)效果。

TB-PC算法；MISO系統(tǒng)辨識(shí)；BA算法；POS算法；協(xié)同進(jìn)化；禁忌搜索

0 引 言

在現(xiàn)實(shí)應(yīng)用中，很多的實(shí)際過(guò)程和應(yīng)用系統(tǒng)，比如熱工過(guò)程，生物化學(xué)反應(yīng)過(guò)程以及微型燃?xì)廨啓C(jī)冷熱電聯(lián)供系統(tǒng)(MGT-CCHP系統(tǒng))等都存在著時(shí)間延遲。所以，進(jìn)行這些系統(tǒng)、過(guò)程辨識(shí)建模，對(duì)控制、預(yù)測(cè)、綜合和分析具有重要的意義[1]。辨識(shí)帶有時(shí)間延遲的系統(tǒng)，現(xiàn)階段主流的辨識(shí)方法有兩種，一種是使用有理傳遞函數(shù)法來(lái)進(jìn)行逼近時(shí)間延遲；但是此方法需要用到大量的估計(jì)參數(shù)，對(duì)于具有很大時(shí)間延遲的系統(tǒng)過(guò)程來(lái)說(shuō)，會(huì)帶來(lái)很大的誤差，偏離實(shí)際系統(tǒng)。另一種方法是采用非線(xiàn)性?xún)?yōu)化法，利用梯度逐級(jí)遞推來(lái)搜索得到最優(yōu)解；但是搜索得到最優(yōu)解常常會(huì)收斂于局部最優(yōu)解，不能得到全局最優(yōu)解[2]。

對(duì)于未知時(shí)間延遲多輸入單輸出( MISO) 系統(tǒng)來(lái)說(shuō)，非線(xiàn)性?xún)?yōu)化法方法更容易收斂到局部最優(yōu)值。并且，如何保持辨識(shí)模型的精確度、收斂度以及所耗時(shí)長(zhǎng)，都是需要解決的問(wèn)題，是當(dāng)今研究的熱點(diǎn)和難點(diǎn)。本文的研究對(duì)象，即為含有多個(gè)未知時(shí)間延遲的MISO系統(tǒng)。首先，介紹了基本蝙蝠算法，其理論基礎(chǔ)易于理解。接著，提出了一種基于協(xié)同進(jìn)化策略和禁忌搜索的蝙蝠和粒子群混合算法(bat algorithm cooperative with particle swarm optimization based on tabu search strategy，TB-PC)，改進(jìn)了蝙蝠算法難以跳出局部最優(yōu)值、收斂速度不高、準(zhǔn)確度低的缺點(diǎn)。然后，對(duì)TB-PC算法做了性能測(cè)試，其具有全局尋優(yōu)搜索性能好、穩(wěn)定性能好和收斂精度高等優(yōu)點(diǎn)。最后利用仿真實(shí)驗(yàn)，驗(yàn)證這種算法對(duì)MISO系統(tǒng)辨識(shí)顯著有效性以及準(zhǔn)確性。

1 蝙蝠算法基本原理

蝙蝠算法(bat algorithm ,BA)是將自己追捕獵物的動(dòng)作模仿為智能體隨機(jī)分布在搜索空間范圍內(nèi)，衡量個(gè)體與獵物空間距離采用的是目標(biāo)函數(shù)；蝙蝠自身搜索目標(biāo)和移動(dòng)滑行行動(dòng)作模擬成用更加接近目標(biāo)簡(jiǎn)單迅速動(dòng)作來(lái)代替尋優(yōu)過(guò)程可行解[3]。此時(shí)，對(duì)蝙蝠的回聲定位系統(tǒng)提出以下假設(shè)，用來(lái)模擬其追捕獵物的過(guò)程[4-5]：

(1) 感應(yīng)蝙蝠個(gè)體與獵物之間的距離，用的是自身的回聲定位；并且其能夠準(zhǔn)確的區(qū)分出獵物與其他障礙物之間的區(qū)別。

(2) 蝙蝠在固有位置xi上是利用速度vi任意飛行，音強(qiáng)A0、可改變波長(zhǎng)λ和頻率f大小值去尋找獵物的；并且其可以利用獵物以及蝙蝠個(gè)體之間距離差大小來(lái)改變波長(zhǎng)，在離獵物很近時(shí)及時(shí)改變蝙蝠個(gè)體發(fā)出波頻度r∈[0,1]。

(3) 音強(qiáng)A從最大值(正值)A0減小到固定最小值A(chǔ)min。

蝙蝠算法中并沒(méi)有采取光線(xiàn)追蹤來(lái)預(yù)算空間規(guī)模以及時(shí)間延遲大小，這就大大減少了多維問(wèn)題的計(jì)算量。由于波長(zhǎng)λ和頻率f符合v=λf規(guī)律，即在特定的頻率內(nèi)[fmin,fmax]必與特定波長(zhǎng)[λmin,λmax]相一致，波長(zhǎng)λ這個(gè)量就可以舍棄。但是在特定波長(zhǎng)λ，通過(guò)調(diào)整其頻率f的大小來(lái)完成尋找目標(biāo)，所以蝙蝠個(gè)體的飛行行為跟音強(qiáng)Ai以及頻度ri相關(guān)。如果Amin=0，就是蝙蝠個(gè)體馬上發(fā)現(xiàn)目標(biāo)之后，沒(méi)進(jìn)行其他行動(dòng)，音強(qiáng)Ai以及頻度ri迭代式如式(1)、式(2)所示。

(1)

(2)

(3)

(4)

式中ε∈[-1,1]是一個(gè)任意的數(shù)字；xbest是當(dāng)前全局最優(yōu)位置。

fi=fmin+(fmax-fmin)×rand

(5)

(6)

(7)

式中fi∈[fmin,fmax]為蝙蝠個(gè)體i尋找目標(biāo)所采取的頻率值，其掌握著粒子移動(dòng)節(jié)拍與界限。

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

蝙蝠算法有著參數(shù)較少、很簡(jiǎn)單就能夠?qū)崿F(xiàn)以及操作簡(jiǎn)單方便等優(yōu)點(diǎn)，但也存在著求解精度偏低、收斂速度慢與難以跳出局部最優(yōu)等缺點(diǎn)[6]。針對(duì)蝙蝠算法上述不足之處，本文一種提出基于協(xié)同進(jìn)化策略和禁忌搜索的蝙蝠和粒子群混合算法(bat algorithm cooperative with particle swarm optimization based on tabu search strategy，TB-PC)，鑒于禁忌策略能夠?qū)崿F(xiàn)算法逃離局部最優(yōu)解，進(jìn)而避免陷入到其中去，提高群體多樣性；同時(shí)結(jié)合粒子群算法，對(duì)算法搜索機(jī)制穩(wěn)定性能夠加以保障，提高系統(tǒng)辨識(shí)收斂精確度。

2.1 BA與PSO協(xié)同進(jìn)化

協(xié)同進(jìn)化是一種模仿自然界生態(tài)系統(tǒng)中物種間的協(xié)作進(jìn)化機(jī)制而得到的策略，它借鑒了生態(tài)學(xué)的種群協(xié)同理論，運(yùn)用了種群間的自動(dòng)調(diào)節(jié)、自動(dòng)適應(yīng)原理。協(xié)同進(jìn)化的各個(gè)種群相互驅(qū)使，相互影響和制約，以提高各自和全局的性能?；镜腂A算法為了得到全局最優(yōu)解，使運(yùn)行速度得到提高并且降低xbest陷入到局部最優(yōu)的機(jī)率，會(huì)跟著脈沖頻度ri和音強(qiáng)Ai的進(jìn)化；但是種群的多樣性沒(méi)有被有效的保護(hù)，比如當(dāng)xbest是局部最優(yōu)值時(shí)，就會(huì)提前收斂，從而群體空間的多樣性就消失了?；镜腜SO算法，不須要大量的梯度信息，收斂精度高以及參數(shù)較少，便于實(shí)現(xiàn)；但是其尋優(yōu)所花費(fèi)時(shí)間較長(zhǎng)，算法運(yùn)行到中后期時(shí)極其容易陷入到局部最優(yōu)中去。綜上兩種算法，各有各的尋優(yōu)優(yōu)勢(shì)，尋優(yōu)的原理不盡相同，但是都存在著較容易陷入到局部最優(yōu)中的缺點(diǎn)[7]。

TB-PC算法綜合了兩種算法的優(yōu)點(diǎn)，以協(xié)同進(jìn)化為算法的基本框架，各個(gè)信息要素之間的交流方式被重新建立，信息要素在蝙蝠群體與粒子群體進(jìn)行轉(zhuǎn)達(dá)，兩個(gè)種群相互借鑒對(duì)方的優(yōu)良粒子要素，使得尋優(yōu)可以從局部尋優(yōu)跳出來(lái)，達(dá)到全局尋優(yōu)的最佳效果。因此，TB-PC算法對(duì)復(fù)雜優(yōu)化問(wèn)題以及大范圍搜索高維度系統(tǒng)辨識(shí)問(wèn)題特別有效。

“分而治之”思想在TB-PC算法中得以實(shí)現(xiàn)，將包含m個(gè)初始值的種群向量xi=[xi1,xi2,…,xid]分開(kāi)為兩個(gè)部分，其中mba個(gè)初始值進(jìn)行單獨(dú)BA算法搜索，mpso個(gè)初始值進(jìn)行單獨(dú)PSO算法搜索，這兩個(gè)算法一起尋找全局空間，每回運(yùn)算以后，蝙蝠、粒子群群體又從新合成為m行d維的向量，由適應(yīng)度更換以往最優(yōu)解當(dāng)作兩子群下一次更換的方向。擴(kuò)大了尋優(yōu)規(guī)模，信息要素交流傳遞得以加強(qiáng)，使算法收斂精確度和效率都得到了大幅度提高。

2.2 禁忌搜索

禁忌搜索(tabu search，TS)是對(duì)局部鄰域搜索的一種擴(kuò)展，是一種全局逐步尋優(yōu)的算法。禁忌搜索算法中充分體現(xiàn)了集中和擴(kuò)散兩個(gè)策略，它的集中策略體現(xiàn)在局部搜索，即從一點(diǎn)出發(fā)，在這點(diǎn)的鄰域內(nèi)尋求更好的解，以達(dá)到局部最優(yōu)解而結(jié)束;而為了跳出局部最優(yōu)解，擴(kuò)散策略通過(guò)禁忌表的功能來(lái)實(shí)現(xiàn)。禁忌表中記下已經(jīng)到達(dá)的某些信息，算法通過(guò)對(duì)禁忌表中點(diǎn)的禁忌，而達(dá)到一些沒(méi)有搜索的點(diǎn)，從而實(shí)現(xiàn)更大區(qū)域的搜索[8-9]。

(1) 適配值函數(shù)和初始值

本文算法適配值函數(shù)由目標(biāo)函數(shù)變化而來(lái)的，就是用來(lái)測(cè)評(píng)算法尋優(yōu)結(jié)果適應(yīng)度函數(shù)，如式(3)所示，尋優(yōu)性能完全依賴(lài)于鄰域結(jié)構(gòu)和初始值。TB-PC算法利用定義研究對(duì)象的穩(wěn)態(tài)變化大小gd，就是相鄰兩次迭代適應(yīng)度函數(shù)變化值，計(jì)算公式如式(8)所示，此為使禁忌觸發(fā)的前提條件。

gd(t)=|goodness(t)-goodness(t-1)|

(8)

式中t=2,3,…,tmax，即為迭代次數(shù)。算法在接連10次迭代運(yùn)算gd值大小都不大于0.01，就認(rèn)定此時(shí)的最優(yōu)值xbest就是最接近真實(shí)值的最優(yōu)值，運(yùn)算過(guò)程穩(wěn)定，應(yīng)當(dāng)引進(jìn)禁忌搜索且把xbest值作為初始解，記下當(dāng)前迭代運(yùn)算次數(shù)以及上次引進(jìn)禁忌搜索迭代運(yùn)算次數(shù)差值nLi，i就是為禁忌引入次數(shù)值。

(2) 禁忌表、禁忌長(zhǎng)度和禁忌對(duì)象

設(shè)置禁忌的目的是最大可能避免迂回反復(fù)尋優(yōu)，禁止重復(fù)前面工作，跳出局部最優(yōu)點(diǎn)。禁忌表指的是被禁對(duì)象具有定長(zhǎng)nJ且?guī)в忻鞔_順序進(jìn)出序列的隊(duì)列。禁忌長(zhǎng)度nJ指的是被禁對(duì)象不允許被選取的迭代步數(shù)，一般是給被禁對(duì)象aj一個(gè)數(shù)lj，要求aj在lj步迭代內(nèi)被禁，在禁忌表中采用Tabu(aj)=lj記憶；每迭代一步，該指標(biāo)做運(yùn)算Tabu(aj)=lj-1一次，直到Tabu(aj)=0時(shí)解禁。禁忌對(duì)象指的是禁忌表中被禁的那些變化元素。

TB-PC算法把蝙蝠個(gè)體群當(dāng)成空間范圍內(nèi)一組粒子，粒子所處具體位置當(dāng)作成算法一組遴選解[10]，然后由式(8)來(lái)判斷引進(jìn)禁忌之后，把當(dāng)前最優(yōu)解xbest跟禁忌表做比對(duì)，假如此全局極值和禁忌表不重合，那么xbest將會(huì)被以特定記憶順序記錄到禁忌表當(dāng)中去，并被視為禁忌對(duì)象；假如此全局極值跟禁忌表重合，此全局極值禁用，接著任意在粒子群體選取另一個(gè)粒子為全局極值，這個(gè)全局極值適應(yīng)度和禁忌表不重合，然后停止此次禁忌搜索，TB-PC算法正常繼續(xù)尋優(yōu)。最后重復(fù)以上算法過(guò)程，直到符合算法運(yùn)算結(jié)束規(guī)則為止。采取先進(jìn)先出的排列順序來(lái)進(jìn)行禁忌對(duì)象的替換，即是每一個(gè)禁忌對(duì)象在經(jīng)歷了禁忌長(zhǎng)度nJ次運(yùn)算后，又會(huì)被重新釋放出來(lái)。其中，禁忌長(zhǎng)度nJ大小由tmax和nLi所決定，多次試驗(yàn)之后可以選取最佳適合值。

(3)特赦準(zhǔn)則和終止準(zhǔn)則

在TB-PC算法中，設(shè)置特赦準(zhǔn)則是為了使算法能夠保持最佳的局部搜索效果，進(jìn)而保證了較為優(yōu)良的全局尋優(yōu)。特赦準(zhǔn)則的設(shè)定，包括以下三個(gè)方面：

①基于適配值的規(guī)則。若出現(xiàn)一個(gè)解的目標(biāo)值好于前面任何一個(gè)最佳候選解，可特赦，用來(lái)提高收斂精度和速度；

②基于最小錯(cuò)誤的規(guī)則，若所有對(duì)象都被禁忌，特赦一個(gè)適配值最小的解；

③基于影響力的規(guī)則，特赦對(duì)目標(biāo)值影響最大的對(duì)象。

考慮到算法的實(shí)際操作情況，設(shè)計(jì)算法的時(shí)候不可能讓禁忌長(zhǎng)度nJ無(wú)限大，因此用終止準(zhǔn)則來(lái)使得禁忌搜索算法的進(jìn)程結(jié)束終止，且收斂近似準(zhǔn)則就是最為常用的方法之一，包含最高要求的收斂精度、最大的迭代次數(shù)、某個(gè)最大禁忌頻率和最大迭代步數(shù)等因素。本文在設(shè)計(jì)改進(jìn)算法時(shí)，當(dāng)?shù)螖?shù)達(dá)到最大迭代次數(shù)tmax90%時(shí)候，就會(huì)當(dāng)作算法終止準(zhǔn)則條件，提高了算法穩(wěn)定性和收斂精確度，同時(shí)最大程度保留了禁忌搜索的優(yōu)良效果。

2.3 TB-PC算法流程

TB-PC算法的流程是由九個(gè)步驟構(gòu)成，算法流程圖如圖1所示。

圖1 TB-PC算法流程圖

3 TB-PC算法性能測(cè)試

為了測(cè)試TB-PC算法全局尋優(yōu)效果、避免過(guò)早收斂的性能以及群體多樣性的保持能力，在此使用Shaffer、Rastrigin、Schwefel以及Rosenbrock這四個(gè)非線(xiàn)性的多峰值函數(shù)，用適應(yīng)度、標(biāo)準(zhǔn)差和耗時(shí)這三個(gè)方面與基本的蝙蝠算法和基本粒子群算法來(lái)進(jìn)行算法性能好壞的分析比較。Shaffer這個(gè)測(cè)試函數(shù)全局極值是-1，其余三個(gè)測(cè)試函數(shù)為0。

由前面的理論分析可知，TB-PC算法對(duì)初始值等參數(shù)大小設(shè)置很敏感，直接影響到算法的穩(wěn)定性和精確度。此時(shí)，函數(shù)基本參數(shù)的設(shè)置如下：解群體數(shù)目m=50，蝙蝠子群mba=25，粒子群子群mpso=25，最大迭代次數(shù)tmax=1 000，搜索脈沖頻率范圍[-1.0,1.0]，最大脈沖頻度r0=0.35，頻度增大系數(shù)γ=0.07，最大脈沖音強(qiáng)A0=1.95，音強(qiáng)衰減系數(shù)α=0.9，慣性權(quán)重采用線(xiàn)性遞減策略，從0.9遞減到0.4，學(xué)習(xí)因子c1=c2=2.0，禁忌表長(zhǎng)度nJ=10。算法重復(fù)運(yùn)行10次并記錄最佳適應(yīng)度、標(biāo)準(zhǔn)差和算法耗時(shí)，結(jié)果如表1[6]和圖2至圖5所示，其中，縱坐標(biāo)適應(yīng)度無(wú)單位量綱，用來(lái)評(píng)價(jià)算法搜索結(jié)果。

表1 TB-PC算法測(cè)試結(jié)果

從表1 TB-PC算法測(cè)試結(jié)果和圖2到圖5各個(gè)函數(shù)適應(yīng)度測(cè)試比較情況能夠看出，對(duì)于Shaffer函數(shù)來(lái)說(shuō)，在維度較低的情況下，這三種算法的適應(yīng)度都能夠很好的收斂到全局最優(yōu)值1，尋優(yōu)效果普遍較好；對(duì)于Rastrigin、Schwefel和Rosenbrock測(cè)試函數(shù)來(lái)說(shuō)，對(duì)于維度較高情況下，PSO算法和BA算法適應(yīng)度偏離全局最優(yōu)值0很大，收斂誤差大；對(duì)于維度較低Rastrigin、Schwefel函數(shù)測(cè)試下，PSO算法收斂精確度明顯要比BA算法高，但是對(duì)于高維度Rosenbrock函數(shù)測(cè)試下，PSO算法收斂精度要稍微比BA算法低；TB-PC算法以協(xié)同進(jìn)化和禁忌搜索為依托，使得尋優(yōu)能夠從局部尋優(yōu)中跳出來(lái)，達(dá)到全局尋優(yōu)，以上四個(gè)函數(shù)適應(yīng)度普遍比較小，其收斂精度較高。

從標(biāo)準(zhǔn)差來(lái)比較分析，以上四個(gè)測(cè)試函數(shù)的標(biāo)準(zhǔn)差都是TB-PC算法最低，從中得出在10次實(shí)驗(yàn)測(cè)試中，TB-PC算法的成功收斂率和穩(wěn)定性也都是最高的。從耗時(shí)長(zhǎng)短來(lái)比較分析，TB-PC算法耗時(shí)比BA算法長(zhǎng)比PSO算法短，介于二者之間。綜合以上三個(gè)方面的數(shù)據(jù)分析，TB-PC算法具有以下幾個(gè)優(yōu)點(diǎn)：

(1) 全局尋優(yōu)搜索性能好。TB-PC算法以協(xié)同進(jìn)化和禁忌搜索為依托，結(jié)合禁忌觸發(fā)準(zhǔn)則、特赦準(zhǔn)則和終止準(zhǔn)則，使得尋優(yōu)能夠從局部尋優(yōu)中跳出來(lái)，十分有效的避免局部尋優(yōu)，能夠達(dá)到最佳全局尋優(yōu)效果。

(2)穩(wěn)定性能好。TB-PC算法對(duì)以上四個(gè)測(cè)試函數(shù)的收斂精度和標(biāo)準(zhǔn)差都最低，測(cè)試結(jié)果穩(wěn)定性好，對(duì)于不同維度的系統(tǒng)辨識(shí)問(wèn)題，都能以最高的機(jī)率取得最低的標(biāo)準(zhǔn)差，系統(tǒng)辨識(shí)結(jié)果穩(wěn)定。

(3) 收斂精度高。從表1以及圖2到圖5中能夠看出，TB-PC算法對(duì)以上四個(gè)測(cè)試函數(shù)都有較好的收斂精度，比PSO算法、BA算法要好，并且很適合大范圍高維度系統(tǒng)辨識(shí)問(wèn)題。

圖2 Schaffer函數(shù)適應(yīng)度測(cè)試比較

圖3 Rastrigin函數(shù)適應(yīng)度測(cè)試比較

圖4 Schwefel函數(shù)適應(yīng)度測(cè)試比較

圖5 Rosenbrock函數(shù)適應(yīng)度測(cè)試比較

4 MISO系統(tǒng)辨識(shí)仿真

構(gòu)造一個(gè)帶有時(shí)間延遲多輸入單輸出系統(tǒng)(MISO)傳遞函數(shù)，如式(9)所示：

y=G1(s)x1+G2(s)x2+G3(s)x3

(9)

從式(9)構(gòu)造出MISO系統(tǒng)中可以看出，這個(gè)系統(tǒng)由慣性環(huán)節(jié)以及比例環(huán)節(jié)組成的三輸入一輸出有時(shí)間延遲的多維度低階次系統(tǒng)。需要辨識(shí)的結(jié)構(gòu)參數(shù)是{K1,K2,K3,T1,T2,T3}，過(guò)程參數(shù)為{τ1,τ2,τ3,n1,n2,n3}。在 MATLAB中，由頻域離散相似化法，把上面的傳遞函數(shù)模型進(jìn)行離散化處理，其采樣周期是1 s，共采樣200組數(shù)據(jù)；將采樣得到的200組數(shù)據(jù)按照傳遞函數(shù)模型分解的結(jié)果計(jì)算出樣本的辨識(shí)數(shù)據(jù)[x1,x2,x3,y]，在[1,3]之間使用遍歷法來(lái)辨識(shí)系統(tǒng)的各個(gè)階次，用TB-PC算法對(duì)于不同的階次辨識(shí){K1,K2,K3,T1,T2,T3,τ1,τ2,τ3}的適應(yīng)度大小，來(lái)比較確定待辨識(shí)系統(tǒng)具體各個(gè)階次。

圖6 辨識(shí)輸出與實(shí)際輸出比較

圖7 TB-PC算法適應(yīng)度隨迭代次數(shù)變化曲線(xiàn)

TB-PC算法參數(shù)設(shè)置：解群體數(shù)目m=50，蝙蝠子群mba=25，粒子群子群mpso=25，最大迭代次數(shù)tmax=1000，搜索脈沖頻率范圍[-1.0,1.0]，最大脈沖頻度r0=0.35，頻度增大系數(shù)γ=0.07，最大脈沖音強(qiáng)A0=1.95，音強(qiáng)衰減系數(shù)α=0.9，慣性權(quán)重采用線(xiàn)性遞減策略，從0.9遞減到0.4，學(xué)習(xí)因子c1=c2=2.0，禁忌表長(zhǎng)度nJ=10，反復(fù)實(shí)行10次試驗(yàn)得到辨識(shí)效果最好的，如式(10)所示。圖6是由式(10)算法辨識(shí)結(jié)果，離散采樣辨識(shí)輸出數(shù)據(jù)和實(shí)際輸出數(shù)據(jù)比較圖，縱坐標(biāo)輸出數(shù)據(jù)無(wú)單位量綱，僅表示數(shù)值大??；圖7為T(mén)B-PC算法適應(yīng)度隨迭代次數(shù)變化曲線(xiàn)圖，縱坐標(biāo)適應(yīng)度無(wú)單位量綱，用來(lái)評(píng)價(jià)算法搜索結(jié)果。

(10)

TB-PC算法辨識(shí)結(jié)果適應(yīng)度為0.004 990 45，平均耗時(shí)13.884 1s，通過(guò)對(duì)真實(shí)參數(shù)與辨識(shí)結(jié)果參數(shù)比較可知，TB-PC算法對(duì)有時(shí)間延遲高維度多變量MISO系統(tǒng)的結(jié)構(gòu)參數(shù)、過(guò)程參數(shù)以及階次參數(shù)辨識(shí)準(zhǔn)確度較高，全局搜索尋優(yōu)性能好，辨識(shí)的系統(tǒng)也具有很高的收斂精度和穩(wěn)定性。因此，TB-PC算法對(duì)大范圍搜索高維度帶時(shí)間延遲MISO系統(tǒng)，有優(yōu)良的辨識(shí)效果。

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

本文提出了一種適用于復(fù)雜大規(guī)模帶有時(shí)間延遲MISO系統(tǒng)模型辨識(shí)的TB-PC算法，通過(guò)對(duì)四個(gè)測(cè)試函數(shù)的測(cè)試比較，驗(yàn)證了TB-PC算法具有穩(wěn)定性能好、收斂精度高和全局尋優(yōu)搜索性能好等優(yōu)點(diǎn)。對(duì)MISO系統(tǒng)進(jìn)行的實(shí)例辨識(shí)仿真結(jié)果表明，TB-PC算法對(duì)其結(jié)構(gòu)、過(guò)程以及階次參數(shù)辨識(shí)準(zhǔn)確度和適應(yīng)度高，全局搜索尋優(yōu)性能和穩(wěn)定性好，辨識(shí)效果優(yōu)良。但是其平均耗時(shí)相對(duì)較長(zhǎng)，這是以后要加以改進(jìn)和進(jìn)一步研究的地方。

[ 1 ] 王建宏,王道波,王志勝.多個(gè)未知時(shí)延的MISO系統(tǒng)的遞推辨識(shí)[J].控制與決策,2010,25(1):93-98.

[ 2 ] GAWTHROP P J, NIHTILA M T, RAD A B. Recursive parameter estimation of continuous-time systems with unknown time delay [J] .Control Theory and Advanced Technology, 2014, 5(1):227-248.

[ 3 ] BORA T C, COELHO L S, LEBENSZTAJN L. Bat-Inspired optimization approach for the brushless DC wheel motor problem[J]. IEEE Transactions on Magentics, 2015, 48(2):947-950.

[ 4 ] 李煜.新型全局優(yōu)化蝙蝠算法[J].計(jì)算機(jī)科學(xué),2013,9(40):225-229.

[ 5 ] 李枝勇,馬良,張惠珍.遺傳變異蝙蝠算法在背包問(wèn)題上的應(yīng)用[J].計(jì)算機(jī)工程與應(yīng)用,2012,10(1):67-72.

[ 6 ] 張宇楠,劉付永.一種改進(jìn)的變步長(zhǎng)自適應(yīng)蝙蝠算法及其應(yīng)用[J].廣西民族大學(xué)學(xué)報(bào),2013,6(19):51-55.

[ 7 ] 黃宇,韓璞,劉長(zhǎng)良.改進(jìn)粒子群算法及其在系統(tǒng)辨識(shí)中的應(yīng)用[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,51(31): 114-120.

[ 8 ] YANG XINSHE. A new metaheuristic bat-inspired algorithm[J]. Nature Inspired Coperative Strategies for Optimization,2015,7(2):65-74.

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

[10] 賀一.禁忌搜索及其并行化研究[D].重慶:西南大學(xué),2006.

Research on MISO System Identification Based on TB-PC Algorithm

Geng Jian， Zhang Yufei， Fan He

(School of Energy and Environmental，Southeast University，Nanjing Jiangsu 210096，China)

How to maintain the accuracy, convergence, time consuming and jump out of the local optimal solution of MISO system identification is a hot and difficult problem in the research. It presented a bat algorithm cooperative with particle swarm optimization based on tabu search strategy (TB-PC), analyzed control parameters, put forward reasonable frequency and initial pulse intensity, combined bat algorithm and particle swarm optimization with co-evolution together and the introduced tabu search. Through the identification and simulation of four test functions and MISO system, the results showed that the TB-PC algorithm has the advantages of good stability, high convergence precision and good identification effect on the MISO system.

TB-PC algorithm; MISO system identification; BA algorithm; PSO algorithm; co-evolution; tabu search

10.3969/j.issn.1000-3886.2017.04.006

TP273

A

1000-3886(2017)04-0018-04

定稿日期: 2016-10-26

耿健( 1991-) ，男，江蘇徐州人，碩士生，研究方向?yàn)橄冗M(jìn)算法在復(fù)雜系統(tǒng)中建模與辨識(shí)。