第一作者李良敏女，副教授，1977年6月生

基于粒子群算法優(yōu)化的獨(dú)立分量分析算法

李良敏，任景巖

(1.長(zhǎng)安大學(xué)汽車學(xué)院汽車運(yùn)輸安全保障技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室，西安710064;2.西安格蒂電力有限公司，西安710075)

摘要：通過(guò)兩組模擬信號(hào)對(duì)三種主流獨(dú)立分量分析算法—JADE、FastICA、擴(kuò)展Infomax算法的性能進(jìn)行了對(duì)比分析，結(jié)果表明三種算法均無(wú)法完全分離超高斯源與亞高斯源形成的混合信號(hào)，F(xiàn)astICA算法對(duì)能量強(qiáng)弱差別大的混合信號(hào)失效?；谶@一現(xiàn)象，提出了一種新的獨(dú)立分量分析算法，以粒子群算法為優(yōu)化工具，以分離矩陣為優(yōu)化變量，最小化分離信號(hào)聯(lián)合概率與邊緣概率乘積的差值，并給出了具體的計(jì)算流程。仿真實(shí)驗(yàn)結(jié)果表明，該算法的性能顯著優(yōu)于上述三種獨(dú)立分量分析算法。同時(shí)，新提出算法實(shí)施過(guò)程中不需要任何先驗(yàn)知識(shí)，相比其他三種ICA算法，更適合解決工程實(shí)際問題。最后，將該算法應(yīng)用于對(duì)滾動(dòng)軸承實(shí)驗(yàn)臺(tái)實(shí)測(cè)信號(hào)的處理，通過(guò)對(duì)分離信號(hào)的分析實(shí)現(xiàn)了對(duì)滾動(dòng)軸承故障類型的準(zhǔn)確識(shí)別，進(jìn)一步證明了算法的有效性。

關(guān)鍵詞：獨(dú)立分量分析；FastICA；JADE；擴(kuò)展Infomax算法；粒子群算法；滾動(dòng)軸承

基金項(xiàng)目：國(guó)家自然科學(xué)

收稿日期：2013-11-08修改稿收到日期：2014-04-10

中圖分類號(hào):TN911.7文獻(xiàn)標(biāo)志碼：A

A universal particle swarm-optimized independent component analysis algorithm

LILiang-min1,RENJing-yan2(1. Key Laboratory of Automotive Transportation Safety Enhancement Technology, Ministry of Communication,School of Automobile, Chang’an University, Xi’an 710064, China;2. Grid Electric Power Science & Technology Co., Ltd., Xi’an 710075, China)

Abstract:Two sets of simulated signals were made to test the separation ability of three popular independent component analysis (ICA) algorithms including JADE, FastICA, and extended-Infomax. The results showed that the three ICA algorithms can’t recover source signals from mixtures of super-Gaussian sources and sub-Gaussian ones precisely; FastICA fails in solving the separation problem of strong sources mixed with weak sources. A particle swarm optimized ICA algorithm minimizing the difference between joint probabilities and products of marginal probabilities of separated signals was proposed. The computing procedure was derived. Simulation tests showed that compared with the above three ICA algorithms, the proposed algorithm is the best; furthermore, the implementation of the proposed algorithm needs no prior knowledge, thus it is more suitable for solving practical engineering problems. Finally, the proposed algorithm was used to process the actual signals sampled from a rolling bearing test rig. The separated signals were analyzed to indentify the fault types of rolling bearings, the effectiveness of the proposed algorithm was verified.

Key words:independent component analysis (ICA); FastICA; JADE; extended-Infomax; particle swarm optimization (PSO); rolling bearing

獨(dú)立分量分析(Independent Component Analysis，ICA)是一種應(yīng)用比較廣泛的用來(lái)解決盲源分離問題(Blind Source Separation, BSS)的方法，可以在源信號(hào)的數(shù)量、位置及傳輸通道參數(shù)等信息未知的情況下，僅利用源信號(hào)的一些統(tǒng)計(jì)特性，從觀測(cè)信號(hào)中恢復(fù)提取出源信號(hào)[1]。由于這一特性，獨(dú)立分量分析在許多領(lǐng)域得到了廣泛應(yīng)用[2-4]。

國(guó)內(nèi)外學(xué)者提出了許多ICA實(shí)現(xiàn)算法，其中較為典型的算法包括FastICA[5-6]、JADE[7-8]和擴(kuò)展Infomax算法[9]，大量文獻(xiàn)充分證明了這些算法處理常規(guī)BSS問題(僅存在超高斯源或亞高斯源，且各源信號(hào)能量強(qiáng)弱較為接近)的有效性。然而，在實(shí)際使用過(guò)程中，可能會(huì)遇到一些非常規(guī)情況，例如分離超高斯源與亞高斯源形成的混合信號(hào)、從強(qiáng)背景噪聲中檢測(cè)弱源信號(hào)等，這些情況下上述三種算法的應(yīng)用效果難以判定。本文通過(guò)模擬信號(hào)對(duì)三種ICA算法處理非常規(guī)信號(hào)的能力進(jìn)行了分析，在此基礎(chǔ)上提出了一種基于概率論獨(dú)立性判據(jù)最小化的ICA算法，建立了優(yōu)化數(shù)學(xué)模型，選擇粒子群算法進(jìn)行優(yōu)化，并給出了算法流程，最后通過(guò)仿真實(shí)驗(yàn)及滾動(dòng)軸承實(shí)驗(yàn)臺(tái)實(shí)測(cè)信號(hào)對(duì)新提出算法的性能進(jìn)行了檢驗(yàn)。

1典型ICA算法性能分析

本節(jié)通過(guò)模擬信號(hào)對(duì)FastICA、JADE、擴(kuò)展Infomax算法處理非常規(guī)混合信號(hào)的性能進(jìn)行對(duì)比分析，評(píng)價(jià)指標(biāo)采用分離信號(hào)信噪比RSN(yi)，其計(jì)算公式如下[10]：

(1)

式中，si為源信號(hào)，yi為對(duì)應(yīng)分離信號(hào)。RSN(yi)越大，算法分離效果越好。

1.1　超高斯源混合亞高斯源

圖1 源信號(hào)波形Fig.1Sourcesignals圖2 混合信號(hào)波形Fig.2Mixedsignals圖3 擴(kuò)展Infomax算法分離結(jié)果Fig.3Signalsrecoveredbyextended-Infomaxalgorithm

表1　三種算法分離信號(hào)信噪比(y 1—源信號(hào)1，y 2—源信號(hào)2)

通過(guò)對(duì)比可以看出，三種算法中擴(kuò)展Infomax算法的分離效果最好；FastICA算法在根據(jù)源信號(hào)的峭度特征選擇了恰當(dāng)?shù)腉(u)函數(shù)后，分離效果略優(yōu)于JADE算法。由圖3可以看出，三種算法中表現(xiàn)最佳的擴(kuò)展Infomax算法其分離結(jié)果仍不夠理想，分離信號(hào)存在輕度混疊。

1.2　強(qiáng)源信號(hào)混合弱源信號(hào)

在實(shí)際應(yīng)用中，經(jīng)常會(huì)遇到強(qiáng)源信號(hào)混合弱源信號(hào)，例如對(duì)機(jī)械設(shè)備進(jìn)行早期故障診斷時(shí)，由于設(shè)備故障狀態(tài)不明顯，故障信號(hào)極為微弱，完全淹沒在背景噪聲信號(hào)中，屬于典型的強(qiáng)源信號(hào)混合弱源信號(hào)。

圖4　源信號(hào)波形 Fig.4 Source signals

圖5　混合信號(hào)波形 Fig.5 Mixed signals

表2給出了各種算法分離信號(hào)信噪比，其中FastICA算法認(rèn)為源信號(hào)只有一組，分離失敗。對(duì)比結(jié)果表明，JADE算法和擴(kuò)展Infomax算法均成功的將源信號(hào)分離了出來(lái)，其中JADE算法的效果略優(yōu)于擴(kuò)展Infomax算法。

表2　三種算法分離信號(hào)信噪比(y 1—源信號(hào)1，y 2—源信號(hào)2)

結(jié)論：仿真實(shí)驗(yàn)表明，三種ICA算法均無(wú)法完全分離由超高斯源與亞高斯源形成的混合信號(hào)，F(xiàn)astICA算法對(duì)強(qiáng)源信號(hào)混合弱源信號(hào)分離失敗。

2基于粒子群優(yōu)化的ICA算法

簡(jiǎn)而言之，ICA是一種多變量分析方法，要求分離信號(hào)間盡可能的統(tǒng)計(jì)獨(dú)立，因此在ICA算法設(shè)計(jì)中，首先需要建立度量輸出變量獨(dú)立性的目標(biāo)函數(shù)，然后選擇恰當(dāng)?shù)膬?yōu)化算法加以解決。

2.1　優(yōu)化目標(biāo)函數(shù)

在概率論中，獨(dú)立性通過(guò)下式加以定義：

(2)

式中P(y1,y2,…,ym)為變量y1,y2,…,ym的聯(lián)合概率，P(yi)為變量yi的邊緣概率。因此可以通過(guò)信號(hào)聯(lián)合概率與邊緣概率乘積間的差值來(lái)衡量信號(hào)間的獨(dú)立性，建立如下所示的優(yōu)化數(shù)學(xué)模型：

s.t.:y(t)=Wx(t)

(3)

式中，x(t)為n維觀測(cè)信號(hào)，y(t)為m維分離信號(hào)，W為分離矩陣。

2.2　優(yōu)化過(guò)程

由于優(yōu)化目標(biāo)函數(shù)的梯度信息難以獲取，不適合采用傳統(tǒng)優(yōu)化算法，本文選擇采用粒子群優(yōu)化算法(Particle Swarm Optimization，PSO)[12-13]來(lái)解決式(3)所示優(yōu)化問題。

基于粒子群優(yōu)化的ICA算法流程如下：

步驟1確定粒子群規(guī)模M，并初始化粒子群。

步驟2將各粒子位置代入優(yōu)化目標(biāo)函數(shù)，評(píng)價(jià)粒子性能好壞。

(1)將粒子位置與混合信號(hào)相乘，得到呈離散狀態(tài)的分離信號(hào)Yj(k),j=1,…,m;k=1,…,N。

(2)采用直方圖法估計(jì)分離信號(hào)的聯(lián)合概率及邊緣概率。

① 對(duì)于第j個(gè)分離信號(hào)Yj(k),k=1,…,N，參照以下經(jīng)驗(yàn)公式確定直方圖的組距寬度hj[14]：

(4)

③統(tǒng)計(jì)分離信號(hào)落入各個(gè)區(qū)間的離散觀測(cè)點(diǎn)數(shù)Dr1r2,…,rm,r1=1,…,N1,r2=1,…,N2,…,rm=1,…,Nm;

④信號(hào)Yj(k),j=1,…,m,k=1,…,N的聯(lián)合概率可用下式近似估計(jì)：

(5)

第j個(gè)分離信號(hào)Yj(k),k=1,…,N的邊緣概率估計(jì)公式如下：

(6)

(3)將上述計(jì)算結(jié)果代入式(3)，即可得到該粒子所對(duì)應(yīng)的目標(biāo)函數(shù)值，該值越小，該粒子性能越好。

步驟5根據(jù)下列公式更新每個(gè)粒子的速度與位置：

(7)

(8)

(9)

式中,wmax為最大權(quán)重，wmin為最小權(quán)重，T為算法迭代總次數(shù)。

步驟6如未達(dá)到結(jié)束條件(通常為達(dá)到預(yù)先設(shè)定的最大迭代次數(shù)T)，返回步驟2。

3基于PSO優(yōu)化的ICA算法性能檢驗(yàn)

采用相同的模擬信號(hào)對(duì)基于PSO優(yōu)化的ICA算法性能進(jìn)行檢驗(yàn)，為簡(jiǎn)化描述，以下將該算法記為PSO-ICA算法。

3.1　超高斯源混合亞高斯源

本例中PSO各控制參數(shù)取值分別為：M=30，T=800，c1=1.4，c2=1.5，wmax=0.95，wmin=0.4，優(yōu)化變量取值范圍[-1, 1]。

通過(guò)PSO-ICA算法得到的分離信號(hào)信噪比分別為：RSN(y1)=142.482 3，RSN(y2)=183.759 6?？梢钥闯?，PSO-ICA算法得到的分離信號(hào)波形極為干凈，信噪比遠(yuǎn)大于其他三種算法。對(duì)于超高斯源混合亞高斯源，PSO-ICA算法表現(xiàn)出了優(yōu)異的分離能力。

3.2　強(qiáng)源信號(hào)混合弱源信號(hào)

本例中優(yōu)化變量取值范圍設(shè)定為[-100,100]，其余各控制參數(shù)取值與上例相同。

通過(guò)PSO-ICA算法得到的分離信號(hào)信噪比分別為：RSN(y1)=256.566 7，RSN(y2)=167.476 5。同樣，PSO-ICA算法得到的分離信號(hào)波形極為干凈，信噪比遠(yuǎn)大于其他三種算法。對(duì)于強(qiáng)源信號(hào)混合弱源信號(hào)，PSO-ICA算法表現(xiàn)出了同樣優(yōu)異的分離能力。

4應(yīng)用實(shí)例

本節(jié)利用所提出算法處理滾動(dòng)軸承實(shí)驗(yàn)臺(tái)實(shí)測(cè)聲音信號(hào)，發(fā)聲源包括6308滾動(dòng)軸承和電機(jī)，其中已知滾動(dòng)軸承存在故障，但故障類型未知；電機(jī)工作正常。采用兩組聲級(jí)計(jì)同時(shí)采集，聲級(jí)計(jì)位于電機(jī)和軸承的中間位置，得到聲音信號(hào)見圖6。

圖6　聲級(jí)計(jì)采集得到信號(hào) Fig.6 Sound signals picked up by sound level meters

圖7　分離信號(hào) Fig.7 Sound signals recovered by PSO-ICA algorithm

PSO-ICA算法分離結(jié)果見圖7。當(dāng)滾動(dòng)軸承出現(xiàn)疲勞點(diǎn)蝕、剝落等局部損傷類故障時(shí)，振動(dòng)信號(hào)中會(huì)出現(xiàn)周期性的沖擊成分，正常電機(jī)主要表現(xiàn)為低頻振動(dòng)，由此可判斷出分離信號(hào)1為軸承聲音信號(hào)，信號(hào)2為電機(jī)聲音信號(hào)。

根據(jù)6308滾動(dòng)軸承的結(jié)構(gòu)參數(shù)及安轉(zhuǎn)軸轉(zhuǎn)頻52.3 Hz，計(jì)算可知軸承內(nèi)圈、外圈或滾動(dòng)體剝落時(shí)故障特征頻率分別為257.6 Hz，161.2 Hz，107.2 Hz。由圖6可以清楚的看出，軸承聲音信號(hào)中包含有周期為0.062 s的周期性沖擊，對(duì)應(yīng)于外圈剝落故障特征頻率(161.2Hz)，說(shuō)明該軸承外圈出現(xiàn)了剝落類故障。對(duì)該軸承各部件進(jìn)行檢查，發(fā)現(xiàn)外圈滾道中央存在面積為7 mm2、深度為0.2 mm的剝落，證明了分析的正確性及所提出算法的有效性。

5結(jié)論

通過(guò)模擬信號(hào)對(duì)FastICA、JADE、擴(kuò)展Infomax算法處理非常規(guī)信號(hào)的能力進(jìn)行了對(duì)比分析，結(jié)果表明上述三種算法均無(wú)法完全分離超高斯與亞高斯的混合信號(hào)，F(xiàn)astICA算法對(duì)強(qiáng)源信號(hào)混合弱源信號(hào)分離失敗。針對(duì)這一現(xiàn)象，提出了一種基于粒子群優(yōu)化的獨(dú)立分量分析算法，建立了優(yōu)化數(shù)學(xué)模型，并給出了算法流程。采用模擬信號(hào)對(duì)新提出算法的性能進(jìn)行了分析，仿真實(shí)驗(yàn)表明，新提出算法對(duì)非常規(guī)盲源分離問題表現(xiàn)出了極為出色的分離能力，分離信號(hào)極為干凈，信噪比高于其他三種ICA算法，說(shuō)明其盲源分離能力優(yōu)于上述三種ICA算法。同時(shí)，新提出算法實(shí)施過(guò)程中不需要任何先驗(yàn)知識(shí)，相比其他三種ICA算法，更適合解決工程實(shí)際問題。最后，利用新提出算法對(duì)滾動(dòng)軸承實(shí)驗(yàn)臺(tái)實(shí)測(cè)信號(hào)進(jìn)行了分析處理，成功分離出軸承聲音信號(hào)和電機(jī)聲音信號(hào)，并通過(guò)對(duì)軸承聲音信號(hào)的分析實(shí)現(xiàn)了對(duì)軸承故障原因的準(zhǔn)確識(shí)別，進(jìn)一步驗(yàn)證了所提出算法的有效性和實(shí)用性。

參考文獻(xiàn)

[1]Comon P， Jutten C. Handbook of blind source separation: Independent component analysis and applications[M]. Oxford UK: Academic Press, 2010: 155-178.

[2]李欣，梅德慶，陳子辰．基于ICA的鏜削過(guò)程顫振征兆信號(hào)分離方法研究[J]．振動(dòng)與沖擊，2013，32(9)：5-9．

LI Xin, MEI De-qing, CHEN Zi-chen. ICA based separation of chatter symptom signals for precision hole boring processing[J].Journal of Vibration and Shock,2013,32(9):5-9.

[3]Lei X, Caldes Sosa P A, Yao D. EGG/fMRI fusion based on independent component analysis: integration of data-driven and model-driven methods[J]. Journal of Integrative Neuroscience, 2012, 11(3): 313-337.

[4]Yu J, Chen J Y, Rashid M M. Multiway independent component analysis mixture model and mutual information based fault detection and diagnosis approach of multiphase batch processes[J].AiChE Journal, 2013, 59(8):2761-2779.

[5]Hyvarien A. Fast and roubust fixed-point algorithms for independent component analysis[J]. IEEE Transactions on Neural Networks, 1999, 10(3): 626-634.

[6]Reyhani N, Ylipaavalniemi J, Vigario R, et al. Consistency and asymptotic normality of FastICA and bootstrap FastICA[J]. Signal processing, 2012, 92(8):1767-1778.

[7]Cardoso J F. High-order contrasts for independent component analysis[J]. Neural Computation, 1999, 11(1):157-192.

[8]虞海強(qiáng)，王平．基于JADE算法的變壓器振動(dòng)信號(hào)分離的研究[J]．現(xiàn)代電力，2012,29(1)：42-46.

YU Hai-qiang, WANG Ping. Research on vibration signal separation of transformer based on JADE algorithm[J]. Modern Electric Power, 2012, 29(1): 42-46.

[9]Lee T, Girolami M, Sejnowski T. Independent component analysis using an extended infomax algorithm for mixed sub-Gaussian and super-Gaussian sources[J]. Neural Computation, 1999, 11(2): 417-441.

[10]徐明彪，朱維彰．關(guān)于信號(hào)盲分離分離效果評(píng)判指標(biāo)的分析[J]．杭州電子工業(yè)學(xué)院學(xué)報(bào)，2002，2(3)：63-66．

XU Ming-biao, ZHU Wei-zhang. Analyses of the index used for evaluating the effect of BSS[J]. Journal of Hangzhou Institute of Electronic Engineering, 2002, 2(3): 63-66.

[11]張賢達(dá)，保錚. 盲信號(hào)分離[J]. 電子學(xué)報(bào)，2001，29(S1)：1766-1771.

ZHANG Xian-da, BAO Zheng. Blind source separation[J]. Acta Electornica Sinica, 2001, 29(S1): 1766-1771.

[12]李士勇，李研. 智能優(yōu)化算法與應(yīng)用[M]. 哈爾濱：哈爾濱工業(yè)大學(xué)出版社，2012.

[13]莫思敏，曾建潮，徐衛(wèi)濱. 具有自組織種群結(jié)構(gòu)的微粒群算法[J]. 系統(tǒng)仿真學(xué)報(bào)，2013，25(3)：445-450.

MO Si-min, ZENG Jian-chao, XU Wei-bin. Particle swarm optimization based on self-organization topology[J]. Journal of System Simulation, 2013, 25(3): 445-450.

[14]肖云魁. 汽車故障診斷學(xué)[M].北京：北京理工大學(xué)出版社，2006: 05.