衛(wèi)曉娟，丁旺才，李寧洲，周學(xué)舟

（蘭州交通大學(xué)機(jī)電工程學(xué)院，蘭州 730070）

基于改進(jìn)粒子群算法的Volterra模型參數(shù)辨識(shí)

衛(wèi)曉娟，丁旺才，李寧洲，周學(xué)舟

（蘭州交通大學(xué)機(jī)電工程學(xué)院，蘭州 730070）

針對(duì)非線性系統(tǒng)Volterra泛函級(jí)數(shù)模型，結(jié)合混沌優(yōu)化策略和種群多樣性控制思想，提出了一種改進(jìn)粒子群算法，并應(yīng)用于Volterra模型參數(shù)的辨識(shí)，將非線性系統(tǒng)的辨識(shí)問(wèn)題轉(zhuǎn)化為高維參數(shù)空間上的優(yōu)化問(wèn)題。利用混沌序列增加初始種群的多樣性，通過(guò)構(gòu)建動(dòng)態(tài)子群以進(jìn)行協(xié)作尋優(yōu)，且各子群采用不同的參數(shù)自適應(yīng)調(diào)整策略，并定義算法收斂性測(cè)度以對(duì)精英粒子進(jìn)行合理的混沌變異，避免了算法早熟收斂，提高了算法的尋優(yōu)速度和尋優(yōu)精度。仿真實(shí)驗(yàn)中，將該方法與基于標(biāo)準(zhǔn)粒子群算法、遺傳算法、量子粒子群算法的Volterra模型參數(shù)辨識(shí)方法相比較，驗(yàn)證了該辨識(shí)方法的有效性和魯棒性。

改進(jìn)粒子群算法；非線性系統(tǒng)；Volterra級(jí)數(shù)；系統(tǒng)辨識(shí)

Volterra級(jí)數(shù)模型由意大利科學(xué)家Volterra提出，它能完全刻畫(huà)系統(tǒng)非線性傳遞特性［1］，因此Volterra級(jí)數(shù)可以描述一大類(lèi)非線性現(xiàn)象。對(duì)于該模型參數(shù)辨識(shí)的研究是利用Volterra級(jí)數(shù)模型描述并研究非線性系統(tǒng)動(dòng)態(tài)特性、并由此監(jiān)測(cè)和評(píng)估非線性系統(tǒng)健康狀況的關(guān)鍵。

非線性系統(tǒng)Volterra核的辨識(shí)實(shí)質(zhì)上是一個(gè)優(yōu)化問(wèn)題。目前Volterra模型的辨識(shí)方法有傳統(tǒng)的最小二乘算法［2－5］、遺傳算法［6－7］、蟻群算法［8］、量子粒子群算法［9－10］等，但是，最小二乘算法辨識(shí)時(shí)要求目標(biāo)函數(shù)連續(xù)可導(dǎo)，且采用梯度信息進(jìn)行搜索，容易陷入局部最小，而遺傳算法等智能計(jì)算方法雖然克服了傳統(tǒng)辨識(shí)算法的相應(yīng)缺點(diǎn)，但在算法結(jié)構(gòu)和優(yōu)化效率上依然有改進(jìn)的余地。

粒子群算法（Particle Swarm Optimization，PSO）是基于生物種群模擬的群體智能優(yōu)化算法［11－12］，具有結(jié)構(gòu)簡(jiǎn)單、易于實(shí)現(xiàn)等特點(diǎn)［13］。因此，本文把粒子群算法引入到非線性Volterra級(jí)數(shù)模型的辨識(shí)中，且為了提高算法性能和避免陷入局部極值點(diǎn)，提出了一種改進(jìn)粒子群算法（Improved Particle Swarm Optimization Algo-rithm，IPSO）。IPSO算法在初始化階段以混沌序列初始化種群，增加了初始搜索的遍歷性，進(jìn)化過(guò)程中則依據(jù)粒子性能將粒子群劃分為不同子群，且各子群采用不同的參數(shù)自適應(yīng)調(diào)整策略，同時(shí)定義了算法收斂性測(cè)度以便采用混沌變異策略對(duì)部分粒子進(jìn)行合理變異，不僅增加了種群的多樣性、避免了算法陷入局部最優(yōu)，而且實(shí)現(xiàn)了全局探索和局部開(kāi)發(fā)能力的平衡，提高了算法的整體尋優(yōu)效率。將IPSO算法應(yīng)用于一個(gè)典型的Volterra級(jí)數(shù)模型并進(jìn)行仿真實(shí)驗(yàn)，驗(yàn)證了本文提出的辨識(shí)方法的有效性和魯棒性。

1 改進(jìn)粒子群算法

1.1 標(biāo)準(zhǔn)粒子群算法

標(biāo)準(zhǔn)粒子群算法（Standard Particle Swarm Optimization，SPSO）的核心思想是通過(guò)跟蹤粒子當(dāng)前的局部最優(yōu)解和全局最優(yōu)解來(lái)更新粒子速度和位置，達(dá)到中止條件時(shí)，當(dāng)前全局最優(yōu)解即作為優(yōu)化問(wèn)題的最優(yōu)解。

假設(shè)粒子群規(guī)模為n，粒子j的速度和位置分別用M維向量Vj和Xj表示，則SPSO算法的進(jìn)化方程為［14］

式中，Vj（t）為粒子j在第t代的速度，Xj（t）為粒子j在第t代的位置，Pj（t）為粒子j在第t代的個(gè)體歷史最優(yōu)位置，Pg（t）為粒子群在第t代的歷史最優(yōu)位置，w為慣性權(quán)重，c1為認(rèn)知系數(shù)，c2為社會(huì)系數(shù)，r1、r2為各分量在［0，1］之間的M維隨機(jī)向量。Pj（t）和Pg（t）的更新規(guī)則為［14］

式中，f（Pj（t））表示Pj（t）的適應(yīng)值，對(duì)于最小化問(wèn)題而言，適應(yīng)值越小，對(duì)應(yīng)的適應(yīng)度越優(yōu)；argmin ｛f（Pj（t））｜j＝1，2，…，n｝表示使f（Pj（t））取最小值的Pj（t），該P(yáng)j（t）則作為第t代群體最優(yōu)位置Pg（t）。

1.2 改進(jìn)粒子群算法的改進(jìn)策略

目前大多數(shù)文獻(xiàn)并未考慮按照粒子的性能進(jìn)行參數(shù)調(diào)整［17－19］，從而在一定程度上影響了算法的尋優(yōu)效率；同時(shí)，搜索過(guò)程中由于粒子群的聚集性，群體多樣性的丟失不可避免，可能陷入局部極值而造成早熟收斂。

基于此，本文提出的IPSO算法的改進(jìn)策略主要體現(xiàn)在三個(gè)方面：①粒子群的混沌序列初始化；②根據(jù)粒子性能，在進(jìn)化過(guò)程中動(dòng)態(tài)構(gòu)建精英子群和普通子群，并對(duì)各子群采用不同的參數(shù)自適應(yīng)調(diào)整策略，以增加種群多樣性、同時(shí)實(shí)現(xiàn)全局探索和局部開(kāi)發(fā)能力的平衡；③對(duì)部分精英粒子進(jìn)行混沌變異操作。

1.2.1 粒子群的混沌序列初始化

采用混沌序列初始化粒子群，是在不改變PSO算法初始化時(shí)所具有的隨機(jī)性本質(zhì)的前提下，利用混沌思想提高種群的多樣性和粒子搜索的遍歷性。文獻(xiàn)［14－16］通過(guò)分析指出式（3）所示Tent模型比Logisitic映射具有更好的遍歷性和更高的搜索效率，并驗(yàn)證了Tent模型在增加種群多樣性及改善算法尋優(yōu)性能方面的有效性。

假設(shè)粒子群規(guī)模為n，變量的維數(shù)為M維。粒子群初始化時(shí)，首先隨機(jī)產(chǎn)生一個(gè)M維，且每個(gè)分量數(shù)值在［0，1］之間的向量T1＝（T11，T12，…，T1M），然后將其代入式（3），經(jīng)過(guò)N次（N＞n）迭代產(chǎn)生一個(gè)包含N個(gè)向量T1，T2，…，Tm，…，TN的混沌序列，再將Tm（m＝1，2，…，N）的各個(gè)分量經(jīng)映射xmk＝xkmin+（xkmaxxkmin）Tmk（k＝1，2，…，M）轉(zhuǎn)化到對(duì)應(yīng)變量的取值區(qū)間，構(gòu)成包含N個(gè)初始粒子的群體，再通過(guò)比較這N個(gè)粒子的適應(yīng)值，從中選出適應(yīng)值相對(duì)較小的n個(gè)解作為初始解，此即為粒子群的初始位置。

1.2.2 參數(shù)自適應(yīng)調(diào)整策略

（1）子群構(gòu)建

設(shè)X（t）＝｛X1（t），X2（t），…，Xn（t）｝為第t代的群體，f（Xj（t））為Xj（t）的適應(yīng)值，f（Xj（t））越小，粒子適應(yīng)度越優(yōu)（即其性能越優(yōu)）；為了依據(jù)粒子性能構(gòu)建子群，定義群體平均適應(yīng)度為

定義群體最優(yōu)適應(yīng)度為

定義優(yōu)秀子群為

式中，q為優(yōu)秀粒子數(shù)。

定義優(yōu)秀粒子平均適應(yīng)度為

則可構(gòu)建精英子群和普通子群如下：

定義1精英子群 精英子群定義為

式中，s為精英粒子數(shù)。

定義2普通子群 普通子群定義為

式中，p為普通粒子數(shù)。

（2）慣性權(quán)重的自適應(yīng)調(diào)整策略

①精英子群慣性權(quán)重調(diào)整策略

設(shè)wj（t）為粒子j在第t代的慣性權(quán)重。為了說(shuō)明IPSO中精英子群的wj（t）調(diào)整策略，引入下列定義：

定義3 精英子群最差適應(yīng)度 精英子群最差適應(yīng)度定義為

定義4 精英粒子性能指標(biāo) 精英粒子性能指標(biāo)定義為

式中，Perforj（t）∈［0，1］。精英粒子j在第t代的適應(yīng)度越優(yōu)，Perforj（t）越大，反之，則越?。沪恋娜≈蹈鶕?jù)實(shí)際問(wèn)題確定。

將Perforj（t）作為精英粒子j在第t代對(duì)慣性權(quán)重的影響因子，則精英粒子j在第t代的慣性權(quán)重為

式（7）表明，Perforj（t）越大的精英粒子，全局最優(yōu)出現(xiàn)在其周?chē)母怕试酱?，而較小的慣性權(quán)重能增強(qiáng)算法的局部開(kāi)發(fā)能力，所以wj（t）應(yīng)較小，以便保證這類(lèi)粒子進(jìn)行局部搜索；反之，Perforj（t）越小，全局最優(yōu)出現(xiàn)在其周?chē)母怕试叫?，因較大的慣性權(quán)重能增強(qiáng)算法的全局探索能力，所以wj（t）應(yīng)較大，以便保證這類(lèi)粒子進(jìn)行全局搜索。

②普通子群慣性權(quán)重調(diào)整策略

IPSO算法中，普通子群按照式（8）調(diào)節(jié)wj（t）。

式中，iter為迭代步數(shù)，max iter為最大迭代步數(shù)。此時(shí)，wj（t）隨迭代進(jìn)行按線性規(guī)律遞減。

（3）c1、c2的自適應(yīng)調(diào)整策略

c1和c2的調(diào)整會(huì)影響到算法的收斂性能和收斂速度。因此，IPSO算法中，采用式（9）對(duì)c1、c2進(jìn)行迭代更新。

式中，j，k表示粒子j的第k維，λ1，λ2為收縮因子，λ1、 λ2∈（0，3），α1、α2為［0，1］之間的隨機(jī)數(shù)。

按式（9）調(diào)整c1、c2時(shí)，對(duì)于優(yōu)秀子群，因其比較接近全局最優(yōu)，所以c1應(yīng)較小、c2應(yīng)較大，以加速收斂到全局最優(yōu)，從而提高算法的尋優(yōu)效率和尋優(yōu)精度；而對(duì)于普通子群，因其具有良好的全局探索能力，所以應(yīng)使c1較大、c2較小，以擴(kuò)大算法的搜索范圍，防止陷入局部最優(yōu)。

1.2.3 混沌變異策略

（1）算法收斂性測(cè)度

為了定量描述粒子群的狀態(tài)，避免盲目變異引起的大量運(yùn)算以及對(duì)算法收斂速度的影響，采用粒子群平均點(diǎn)距作為算法收斂性測(cè)度來(lái)衡量算法的收斂情況，以便合理進(jìn)行變異操作以擴(kuò)大搜索空間。

定義5 粒子群平均點(diǎn)距 設(shè)群體的平均中心為

（2）混沌變異策略

對(duì)粒子速度進(jìn)行變異以使粒子在相應(yīng)進(jìn)化代數(shù)條件下對(duì)搜索空間進(jìn)行大范圍的遍歷，當(dāng)算法陷入局部最優(yōu)時(shí)，使之能夠跳出局部極值點(diǎn)，進(jìn)入解空間的其它區(qū)域繼續(xù)搜索。

根據(jù)式（3）所示Tent模型迭代產(chǎn)生一個(gè)包含若干個(gè)向量T1，T2，…的混沌序列，再將Tk的各分量經(jīng)過(guò)映射變換至原粒子速度的解空間，得到若干個(gè)向量V1，V2，…，以此替代精英子群中按比例選取的精英粒子的速度，即得到變異后該部分粒子的當(dāng)前速度。

2 非線性系統(tǒng)Volterra模型辨識(shí)

2.1 Volterra模型描述

對(duì)于具有輸入x（k）和輸出y（k）的離散時(shí)間非線性動(dòng)態(tài)系統(tǒng)，可用Volterra級(jí)數(shù)將其描述為［9］：

式中，e（k）為截?cái)嗾`差。hn（m1，m2，…，mn）為非線性系統(tǒng)的n階Volterra時(shí)域核，又稱(chēng)為廣義脈沖響應(yīng)函數(shù)（GIRF）。GIRF具有對(duì)稱(chēng)性，且對(duì)稱(chēng)性是唯一的［1］。利用Volterra核的對(duì)稱(chēng)性，可以大大減小利用Volterra級(jí)數(shù)模型描述非線性系統(tǒng)所需參數(shù)的數(shù)目，從而可大大減小計(jì)算量。

利用前三階Volterra級(jí)數(shù)模型即可近似描述很多的實(shí)際非線性系統(tǒng)，因此，非線性系統(tǒng)Volterra模型的辨識(shí)通常是對(duì)其Volterra級(jí)數(shù)模型的一階、二階和三階核進(jìn)行辨識(shí)［1］。

考慮GIRF的對(duì)稱(chēng)性，根據(jù)式（11），則系統(tǒng)輸出可表述為：

式中，Mi為記憶長(zhǎng)度；Mi取適當(dāng)值時(shí)，e（k）可充分小，在誤差允許范圍內(nèi)，e（k）可忽略不計(jì)；s1（m1）、s2（m1，m2）、s3（m1，m2，m3）由式（13）確定。

若假設(shè)截?cái)郪olterra級(jí)數(shù)模型的最高階次為N，定義系統(tǒng)的輸入矩陣為：

式中，L為數(shù)據(jù)長(zhǎng)度。則k時(shí)刻的輸入（列）向量X（k）可定義為：

式中，e為誤差。

由式（15）可知，基于截?cái)郪olterra級(jí)數(shù)模型的非線性系統(tǒng)辨識(shí)，就是在已知系統(tǒng)的輸入、輸出序列｛x（k）｝、｛y（k）｝的基礎(chǔ)上求解核向量H。其辨識(shí)實(shí)質(zhì)上是一個(gè)最優(yōu)參數(shù)估計(jì)問(wèn)題。

假設(shè)粒子群規(guī)模為n，將Volterra模型的待估計(jì)參數(shù)（核矢量H）作為IPSO算法中的粒子位置，則粒子j的位置用M維向量表示，粒子群的位置可以用矩陣表示。對(duì)粒子j而言，算法中控制辨識(shí)誤差的適應(yīng)度函數(shù)定義為

式中，L為數(shù)據(jù)長(zhǎng)度，j＝1，2，…，n。

通過(guò)搜索使適應(yīng)度函數(shù)值（即適應(yīng)值）最小的全局最優(yōu)解Pg，即是Volterra核矢量H的最優(yōu)解H＾。

2.2 Volterra模型辨識(shí)流程

利用IPSO算法對(duì)形如式（15）所示有限階非線性離散系統(tǒng)的Volterra模型參數(shù)進(jìn)行辨識(shí)的算法流程如圖1所示。

圖1 基于IPSO的Volterra模型辨識(shí)算法流程圖Fig.1 The identification algorithm of Volterra kernel based on the IPSO

3 Volterra模型辨識(shí)測(cè)試

作為辨識(shí)目標(biāo)的Volterra模型如式（17）所示。

式（17）所示模型的一階、二階和三階核的記憶長(zhǎng)度均為3，核向量為：H＝［1.52，－0.55，0.24，0.78，0，0，1.36，－1.18，0，1.25，0，0，－0.88，0，0，0，0，1.06，0］T。

利用本文提出的IPSO算法對(duì)式（17）所示三階非線性系統(tǒng)的Volterra模型進(jìn)行辨識(shí)，系統(tǒng)的輸入信號(hào)為［0，1］之間的均勻白噪聲信號(hào)，種群規(guī)模為80，維數(shù)為19，最大進(jìn)化代數(shù)為800，系統(tǒng)辨識(shí)的對(duì)比測(cè)試重復(fù)進(jìn)行20次，結(jié)果取運(yùn)行20次結(jié)果的平均值。

表1為無(wú)噪聲情況下，IPSO算法、SPSO算法、遺傳算法（簡(jiǎn)稱(chēng)GA算法）、量子粒子群算法（簡(jiǎn)稱(chēng)QPSO算法）所得到的Volterra模型參數(shù)的估計(jì)值。

為了考察存在噪聲干擾時(shí)本文算法的辨識(shí)精度和抗噪性能，在輸入、輸出端疊加獨(dú)立平穩(wěn)的白噪聲，且輸入輸出信號(hào)的信噪比均為20 dB［9］。表2為四種算法測(cè)試得到的辨識(shí)結(jié)果。

表1 無(wú)噪聲干擾下Volterra模型辨識(shí)結(jié)果Tab.1 Estimation results of the Volterra model under the free－noise environment

由表1可知，在無(wú)噪聲干擾情況下，四種算法均可得到真值的近似值，說(shuō)明四種算法均可用于Volterra模型參數(shù)的辨識(shí)。其中應(yīng)用IPSO算法得到的估計(jì)值的平均值更加接近真實(shí)值，其辨識(shí)精度明顯較高。

由表2亦可看到，在加入噪聲干擾時(shí)，IPSO算法得到的估計(jì)值的誤差相對(duì)較小，顯然噪聲干擾并未對(duì)其辨識(shí)精度產(chǎn)生影響，體現(xiàn)了IPSO算法具有較強(qiáng)的魯棒抗噪性能。

圖2給出了無(wú)噪聲情況下，四種辨識(shí)方法分別得到的一階時(shí)域核參數(shù)h1（1）、二階時(shí)域核參數(shù)h2（1，2）和三階時(shí)域核參數(shù)h3（1，2，2）的進(jìn)化曲線；圖3給出了無(wú)噪聲情況下，四種辨識(shí)方法的平均適應(yīng)度值變化曲線。

從圖中可以看出，IPSO算法的收斂速度明顯快于其他三種算法，雖然其他三種算法最終也能較為準(zhǔn)確地估計(jì)出Volterra核參數(shù)的值，但是收斂速度相對(duì)較慢；且由圖3可看出，IPSO算法相較SPSO算法而言，其全局探索能力和局部開(kāi)發(fā)能力都得到了改善，在整個(gè)進(jìn)化過(guò)程中IPSO算法表現(xiàn)出了較強(qiáng)的持續(xù)尋優(yōu)能力，未陷入局部最優(yōu)，且辨識(shí)精度明顯優(yōu)于其他三種算法。

表2 信噪比20 dB時(shí)Volterra模型辨識(shí)結(jié)果Tab.2 Estimation results of the Volterra m odel（SNR＝20 dB）

圖2 無(wú)噪聲時(shí)h1（1），h2（1，2），h3（1，2，2）的收斂曲線Fig.2 The evolution curves of the Volterra kernel vector h1（1），h2（1，2），h3（1，2，2）under the free-noise environment

圖3 無(wú)噪聲時(shí)適應(yīng)度函數(shù)變化曲線Fig.3 The curves of the fitness function under the free-noise environment

圖4給出了信噪比20 dB情況下，四種辨識(shí)方法分別得到的一階時(shí)域核參數(shù)h1（1）、二階時(shí)域核參數(shù)h2（1，2）和三階時(shí)域核參數(shù)h3（1，2，2）的進(jìn)化曲線。從圖可以看出，IPSO算法的收斂速度明顯快于其他算法，雖然其他算法最終也能較為準(zhǔn)確地估計(jì)出Volterra核參數(shù)的值，但是由于噪聲干擾，收斂速度相對(duì)較慢，尤其是SPSO算法的收斂過(guò)程波動(dòng)也較大。

圖4 SNR＝20 dB時(shí)h1（1），h2（1，2），h3（1，2，2）的收斂曲線Fig.4 The evolution curves of the Volterra kernel vector h1（1），h2（1，2），h3（1，2，2）（SNR＝20 dB）

綜上所述，IPSO算法的尋優(yōu)性能和效率明顯優(yōu)于SPSO算法和GA算法，且相對(duì)于QPSO算法其性能也有一定提高。這主要是因?yàn)镮PSO算法是在SPSO算法的基礎(chǔ)上經(jīng)過(guò)優(yōu)化而來(lái)。參數(shù)的自適應(yīng)調(diào)整加快了算法的計(jì)算效率，組合混沌策略的引入有效避免了SPSO算法局部收斂與早熟的問(wèn)題，而SPSO算法早期收斂速度較快，但到尋優(yōu)后期其收斂速度則比較緩慢，這主要是因?yàn)樗惴ㄊ諗康骄植繕O小時(shí)缺乏有效機(jī)制使算法逃離極小點(diǎn)，從而導(dǎo)致算法在接近或進(jìn)入最優(yōu)點(diǎn)區(qū)域時(shí)的收斂速度比較緩慢。因此，相對(duì)而言，本文提出的方法在全局優(yōu)化能力和快速收斂能力上都有較大的提高，很好地克服了基于SPSO算法的Volterra模型辨識(shí)方法的不足。

4 結(jié) 論

本文將改進(jìn)粒子群算法（IPSO算法）引入到非線性系統(tǒng)Volterra模型參數(shù)的辨識(shí)中，并在仿真研究中，將該辨識(shí)方法與相應(yīng)的基于SPSO算法、GA算法和QPSO算法的辨識(shí)方法進(jìn)行了對(duì)比分析。結(jié)果表明，即使在有噪聲干擾的情況下，提出的方法也能得到很好的參數(shù)辨識(shí)精度，且收斂過(guò)程依然平穩(wěn)，表現(xiàn)出很好的魯棒抗噪性能。與其他三種算法相比，由于本文提出的方法將種群劃分為不同的子群，且不同子群采用相應(yīng)的參數(shù)自適應(yīng)調(diào)整策略，并利用混沌變異算子對(duì)一定比例的精英個(gè)體進(jìn)行變異操作，所以既增加了種群多樣性、又很好地利用了不同粒子分別進(jìn)行全局探索和局部開(kāi)發(fā)，提高了算法的計(jì)算效率。本文提出的IPSO算法及其在非線性系統(tǒng)Volterra模型辨識(shí)中的應(yīng)用研究為非線性系統(tǒng)辨識(shí)提供了一種新的有效方法，具有重要的理論價(jià)值和實(shí)際應(yīng)用價(jià)值。

［1］曹建福，韓崇昭，方洋旺.非線性系統(tǒng)理論及應(yīng)用［M］.西安：西安交通大學(xué)出版社，2006.

［2］張家良，曹建福，高峰.大型裝備傳動(dòng)系統(tǒng)非線性頻譜特征提取與故障診斷［J］.控制與決策，2012，27（1）：135 －138.

ZHANG Jia-liang，CAO Jian-fu，GAO Feng.Feature extraction and fault diagnosis of large-scale equipment transmission system based on nonlinear frequency spectrum ［J］.Control and Decision，2012，27（1）：135－138.

［3］張華君.基于遞推批量最小二乘的Volterra級(jí)數(shù)辨識(shí)方法［J］.小型微型計(jì)算機(jī)，2004，25（12）：2282－2285.

ZHANG Hua-jun.Volterra series identif icationmethod based on recursive least square algorithm［J］.Mini-Micro Systems，2004，25（12）：2282－2285.

［4］孔祥玉，韓崇昭，馬紅光，等.基于Volterra級(jí)數(shù)的全解耦RLS自適應(yīng)辨識(shí)算法［J］.系統(tǒng)仿真學(xué)報(bào)，2004，16（4）：807－809.

KONG Xiang-yu，HAN Chong-zhao，MA Hong-guang，et al.Fully decoupled RLS adaptive identification algorithm based on volterra series［J］.Journal of System Simulation，2004，16 （4）：807－809.

［5］孔祥玉，韓崇昭，馬紅光，等.一種總體最小二乘算法及在Volterra濾波器中的應(yīng)用［J］.西安交通大學(xué)學(xué)報(bào)，2004，38（4）：339－342.

KONG Xiang-yu，HAN Chong-zhao，MA Hong-guang，et al.Total least square algorithm and its application to volterra filter［J］.Journal of Xi'an Jiao Tong University，2004，38 （4）：339－342.

［6］唐浩，屈梁生，溫廣瑞.基于Volterra級(jí)數(shù)的轉(zhuǎn)子故障診斷研究［J］.中國(guó)機(jī)械工程，2009，20（4）：447－454.

TANG Hao，QU Liang-sheng，WEN Guang-rui.Fault diagnosis for rotor system based on volterra series［J］.China Mechanical Engineering，2009，20（4）：447－454.

［7］Abbas H M，Bayoumi M M.Volterra system identification using adaptive genetic algorithms［J］.Applied Soft Computing，5（2004）：75－86.

［8］李志農(nóng)，唐高松，肖堯先，等.基于自適應(yīng)蟻群優(yōu)化的Volterra核辨識(shí)算法研究［J］.振動(dòng)與沖擊，2011，30 （10）：35－38.

LIZhi-nong，TANGGao-song，XIAOYao-xian，etal.Volterra series identification method based on adaptive ant colony optimization［J］.Journal of Vibration and Hock，2011，30 （10）：35－38.

［9］李志農(nóng)，蔣靜，陳金剛，等.基于量子粒子群優(yōu)化的Volterra核辨識(shí)算法研究［J］.振動(dòng)與沖擊，2013，32（3）：60－63.

LI Zhi-nong，JIANG Jing，CHEN Jin-gang，et al.Volterra series identification method based on quantum particle swarm optimization［J］.Journal of Vibration and Hock，2013，32 （3）：60－63.

［10］李志農(nóng)，蔣靜，馮輔周，等.基于量子粒子群優(yōu)化Volterra時(shí)域核辨識(shí)的隱Markov模型識(shí)別方法［J］.儀器儀表學(xué)報(bào)，2011，32（12）：2693－2698.

LI Zhi-nong，JIANG Jing，F(xiàn)eng Fu-zhou，et al.Hidden markov model recognition method based on Volterra kernel identified with particle swarm optimization［J］.Chinese Journal of Scientific Instrument，2011，32（12）：2693－2698.

［11］Kennedy J，Eberhart R C.Particle swarm optimization［C］／／Proceedings of ICNN'95-IEEE International Conference on Neural Networks.IEEE，Piscataway，NJ，USA，1995：1942－1948.

［12］Eberhart R，Kennedy J.A new optimizer using particle swarm theory［C］／／Proc of 6th International Symposium on Micro Machine and Human Science，Nagoya，Japan.Piscataway NJ：IEEE Service Center，1995：39－43.

［13］毛鵬軒，肖揚(yáng).粒子群算法在穩(wěn)定時(shí)滯系統(tǒng)設(shè)計(jì)中的應(yīng)用［J］.計(jì)算機(jī)應(yīng)用研究，2012，29（1）：214－216.

MAO Peng-xuan，XIAO Yang.Application of PSO algorithm for stable time-delay system design［J］.Application Research of Computers，2012，29（1）：214－216.

［14］崔志華，曾建潮.微粒群優(yōu)化算法［M］.北京：科學(xué)出版社，2011.

［15］張浩，張鐵男，沈繼紅，等.Tent混沌粒子群算法及其在結(jié)構(gòu)優(yōu)化決策中的應(yīng)用［J］.控制與決策，2008，23（8）：857－862.

ZHANG Hao，ZHANG Tie-nan，SHEN Ji-hong，et al. Research on decision-makings of structure optimization based on improved Tent PSO［J］.Control and Decision，2008，23 （8）：857－862.

［16］魏玉琴，戴永壽，張亞南，等.基于Tent映射的自適應(yīng)混沌嵌入式粒子群算法［J］.計(jì)算機(jī)工程與應(yīng)用，2013，49 （10）：45－49.

WEIYu-qin，DAIYong-shou，ZHAN Ya-nan，et al.Adaptive chaotic embedded particle swarm optimization algorithm based on Tent mapping［J］.Computer Engineering and Applications，2013，49（10）：45－49.

［17］Cui Z H，Cai X J，Zeng J C，et al.Particle swarm optimization with FUSS and RWS for high dimensional functions［J］.Applied Mathematics and Computation，2008，205（1）：98－108.

［18］Arumugam M S，Rao M V C.On the improved performances of the particle swarm optimization algorithms with adaptive parameters，cross-over operators and rootmean square（RMS）variants for computing optimal control of a class ov hybrid system［J］.Applied Soft Computing，2008，8（1）：324－336.

［19］Bajpai P，Singh S N.Fuzzy adaptive particle swarm optimization for bidding strategy in uniform price spotmarket ［J］.IEEE Transactions on Power System，2007，22（4）：2152－2160.

Parametric identification of nonlinear volterra model based on improved PSO algorithm

WEIXiao-juan，DINGWang-cai，LINing-zhou，ZHOU Xue-zhou

（School of Mechatronic Engineering，Lanzhou Jiaotong University，Lanzhou 730070，China）

By combining the particle swarm optimization（PSO）with the chaotic optimization strategy and the control idea of population diversity，an improved particle swarm optimization（IPSO）algorithm was proposed for parametric identification of nonlinear Volterra seriesmodel.The basic idea of themethod was that the problem of nonlinear system identification was converted into an optimization problem in a high-dimensional parameter space.The chaotic optimization strategy was employed to increase the diversity of the initial population.By building dynamic subgroups，the optimization was realized through the collaboration of sub groups.The different adaptive adjustment strategies for control parameters of IPSO were used in subgroups.The convergencemeasure of IPSO was defined to perform the chaoticmutation operation reasonably.So the premature convergence was avoided，the speed and accuracy of IPSO were improved.In simulation tests，the proposed method was compared with Volterra model identification methods based on standard PSO algorithm，the genetic algorithm and the quantum-behaved PSO，respectively.Its effectiveness and robustness were verified.

improved particle swarm optimization algorithm（IPSO）；nonlinear system；Volterra series；system identification

TP18

A

10.13465／j.cnki.jvs.2015.21.019

國(guó)家自然科學(xué)基金項(xiàng)目（11162007，11462011）；甘肅省自然科學(xué)基金項(xiàng)目（1308RJZA149）

2014－05－23 修改稿收到日期：2014－08－22

衛(wèi)曉娟女，博士生，副教授，1973年生

丁旺才男，教授，博士生導(dǎo)師，1964年生