進(jìn)化論動態(tài)參量自適應(yīng)消噪算法

第一作者肖會芳女，博士，講師， 1984年12月生

進(jìn)化論動態(tài)參量自適應(yīng)消噪算法

肖會芳1，邵毅敏2，周曉君2

(1.北京科技大學(xué)國家板帶生產(chǎn)先進(jìn)裝備工程技術(shù)研究中心，北京100083；2.重慶大學(xué)機(jī)械傳動國家重點實驗室，重慶400044)

摘要：由于以輸出信號平均能量的倒數(shù)為適應(yīng)度函數(shù)、以及定值克隆系數(shù)和定值匹配系數(shù)的因素，制約了傳統(tǒng)進(jìn)化論自適應(yīng)消噪算法的收斂特性，影響了傳統(tǒng)進(jìn)化論自適應(yīng)算法降噪的有效性。針對該問題，提出了基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法。新算法中，匹配系數(shù)與進(jìn)化代數(shù)相關(guān)聯(lián)，克隆系數(shù)由進(jìn)化代數(shù)和適應(yīng)度值決定，并采用動態(tài)適應(yīng)度函數(shù)，改善了濾波器的收斂特性和噪聲抑制能力。模擬仿真分析表明，新算法較同類進(jìn)化論算法有較快的收斂速度和良好的消噪效果。物理臺架實驗數(shù)據(jù)驗證了該算法對故障信號的有效提取能力。

關(guān)鍵詞：進(jìn)化論濾波器；自適應(yīng)消噪；變進(jìn)化系數(shù)；動態(tài)適應(yīng)度函數(shù)

基金項目：國家自然科學(xué)基金重點項目(51035008)，國家自然科學(xué)基金青年科學(xué)

收稿日期：2014-01-20修改稿收到日期：2014-03-14

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

Variable evolutionary coefficients adaptive noise cancellation algorithm

XIAOHui-fang1,SHAOYi-mim2,ZHOUXiao-jun2(1. National Engineering Research Center of Flat Rolling Equipment, University of Science and Technology Beijing, Beijing 100083, China;2. State Key Laboratory of Mechanical Transmission, Chongqing University, Chongqing 400044, China)

Abstract:The fitness function using the reciprocal of average energy of output signal, and the constant cloning and mating coefficients all restrict the convergence and noise cancellation performances of the traditional evolutionary adaptive noise cancellation algorithm. Aiming at improving it, a variable coefficients evolutionary adaptive noise cancellation algorithm using dynamic fitness function was proposed. In the new algorithm, the mating coefficient was related to evolutionary generation, the cloning coefficient was determined by evolutionary generation and fitness value, and a dynamic fitness function was introduced. The results from a simulation example and test datasets show that the proposed algorithm exhibits better convergence and noise cancellation performances than other similar algorithms.

Key words:evolutionary filter; adaptive noise cancellation; variable evolutionary coefficients; dynamic fitness function

自適應(yīng)消噪算法可以依據(jù)某種預(yù)先的準(zhǔn)則，在迭代過程中自動調(diào)整自身的參數(shù)或結(jié)構(gòu)，以實現(xiàn)最優(yōu)準(zhǔn)則下的最優(yōu)噪聲消除，提取出能夠反映機(jī)械設(shè)備故障的特征信號[1-8]。以達(dá)爾文生物進(jìn)化論策略建立的進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法可實現(xiàn)全局最優(yōu)估計自適應(yīng)消噪濾波器系數(shù)，從而有效地抑制可加性噪聲[9-10]。

Abe等[11]論證了進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法的全局搜索收斂能力，且比最小均方(Least Mean Square, LMS)、遺傳算法(Genetic Algorithm, GA)等自適應(yīng)算法有更好的全局收斂性能[12]。邵毅敏等[10]研究了克隆系數(shù)和匹配系數(shù)對進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法的影響規(guī)律。在進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法中，適應(yīng)度函數(shù)是評價群體中個體好壞的標(biāo)準(zhǔn)，是模擬自然選擇的唯一依據(jù)，因此，適應(yīng)度函數(shù)的優(yōu)劣將直接影響算法的收斂速度。

本文首先研究了傳統(tǒng)的進(jìn)化論自適應(yīng)消噪算法中，以輸出信號平均能量的倒數(shù)為適應(yīng)度函數(shù)對算法的收斂特性和消噪性能的影響，結(jié)果表明以輸出信號平均能量的倒數(shù)為適應(yīng)度函數(shù)的評價方法，不能完全正確判斷進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法是否已完全收斂和實現(xiàn)最優(yōu)噪聲抑制功能；進(jìn)而提出了基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法，改善了算法的收斂特性和消噪性能。模擬仿真和物理臺架實驗數(shù)據(jù)分析表明，新算法比其他同類算法有更快的收斂速度且具有良好的消噪效果。

1進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法

1.1進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法簡介

進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法主要由濾波器系數(shù)的構(gòu)造、適應(yīng)度值的計算、濾波器參數(shù)進(jìn)化迭代、最優(yōu)輸出信號和算法收斂判斷等步驟組成[9-10]。算法結(jié)構(gòu)圖，如圖1所示。傳統(tǒng)的進(jìn)化論自適應(yīng)消噪算法中，適應(yīng)度函數(shù)采用輸出信號平均能量的倒數(shù)，且克隆系數(shù)和匹配系數(shù)為恒定值。

圖1　進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法結(jié)構(gòu)圖 [9-10] Fig.1 Diagram of the original Evolutionary Digital Filter [9-10]

為了研究進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法的收斂特性和消噪性能，以非線性軸承故障沖擊信號作為故障仿真特征信號D(k)。仿真信號數(shù)學(xué)模型[5，13]為：

(1)

式中：A為幅值，ζ為阻尼率，f為中心頻率，t為時間，τ為沖擊發(fā)生時刻τ=mod(t,1/fm)，fm為沖擊的間隔。選取各參數(shù)A=1，ζ=0.026，f=800 Hz，fm=10 Hz，采樣頻率為2 000 Hz，t=8.192 s，數(shù)據(jù)總長度為16 384點。

其他機(jī)械的振動信號X(k)采用頻率為27 Hz、幅值為1 g的余弦信號；隨機(jī)噪聲N(k)采用均值為0，方差為1，功率為10，服從高斯正態(tài)分布的白噪聲。

主輸入信號S(k)=D(k)+μ×(X(k)+N(k))是特征信號和噪聲信號以信噪比-5 dB相疊加的混合信號。0~1 s時間內(nèi)的時域波形，如圖2所示。此時特征信號完全被強(qiáng)噪聲所淹沒，無法直接從時域波形識別特征信號。參考輸入信號Z(k)與噪聲信號(X(k)+N(k))具有一定相關(guān)性，但具有不同幅值和相位[5]。

圖2　主輸入信號的時域波形 Fig.2 Waveform of the noisy signal

傳統(tǒng)的進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法以輸出信號平均能量的倒數(shù)作為適應(yīng)度函數(shù)[2-5]。不同進(jìn)化代數(shù)時最優(yōu)個體對應(yīng)的輸出信號適應(yīng)度值變化圖，如圖3所示。

圖3顯示，進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法在第40代時實現(xiàn)了收斂，此時進(jìn)化論自適應(yīng)消噪算法的輸出信號波形與模擬特征信號對比圖，如圖4所示。圖4(b)顯示，與特征信號相比，輸出信號中還存在著噪聲信號，表明噪聲信號并沒有得到最優(yōu)抑制，進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法并沒有完全收斂，即以輸出信號平均能量的倒數(shù)為適應(yīng)度函數(shù)的評價方法，不能完全正確判斷進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法是否已經(jīng)完全收斂和實現(xiàn)最優(yōu)噪聲抑制功能。這是由于，當(dāng)進(jìn)化論自適應(yīng)消噪算法以殘差信號平均能量的倒數(shù)作為適應(yīng)度函數(shù)時，適應(yīng)度函數(shù)只能表征信號的整體特性，不能體現(xiàn)出信號的局部特性，且對噪聲不敏感，導(dǎo)致適應(yīng)度函數(shù)已經(jīng)收斂時，進(jìn)化論數(shù)字濾波自適應(yīng)算法處理后的輸出信號與特征信號之間還存在著差別。

圖3 在不同進(jìn)化代數(shù)時最優(yōu)個體對應(yīng)的輸出信號適應(yīng)度值Fig.3Evolutionofthefitnesswithgeneration圖4 輸出信號與模擬特征信號對比圖Fig.4Comparisonoftheoutputsignalsandfeaturesignal

2基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法

2.1基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法

為了有效地解決適應(yīng)度函數(shù)問題，提高收斂效率，本文提出了基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法。在該算法中，采用動態(tài)適應(yīng)度函數(shù)實時度量，并將匹配系數(shù)與進(jìn)化代數(shù)相關(guān)聯(lián)，克隆系數(shù)由進(jìn)化代數(shù)和適應(yīng)度值決定，由此改善濾波器的收斂特性和噪聲抑制能力。算法框圖，如圖5所示。

圖5　基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法 Fig.5 Flow chart of the proposed evolutionary digital filter

圖5中，基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法核心步驟為：

①適應(yīng)度值的計算新準(zhǔn)則

采用Goladberg的線性拉伸適應(yīng)度函數(shù)[14]作為新算法的適應(yīng)度計算函數(shù)，用以克服傳統(tǒng)進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法因采用各濾波器輸出信號與參考信號的殘差信號的平均能量的倒數(shù)為適應(yīng)度函數(shù)，而不能有效地描述個體性能特性的缺點。第j代的第i個個體濾波器Fi的輸出信號的適應(yīng)度值為：

(2)

式中，ei是原來的適應(yīng)度計算值，ej_ave是第j代所有個體的平均原有適應(yīng)度值，ej_max是第j代所有個體中最大原有適應(yīng)度值，c為常數(shù)，一般取值1~2。

②進(jìn)化規(guī)則

進(jìn)化論規(guī)則中克隆法：第j代第i個個體所產(chǎn)生的Nij個子代個體的特征向量可由式(3)確定：

Wc,i,j=Wp,i,j+r*·nl

(3)

式中，i=1,2…Nap，l=1,2…Nac，Nap和Nac分別是參與克隆的父代和子代個體個數(shù)，r*是改進(jìn)的克隆系數(shù)，nl是均值為零的高斯分布的隨機(jī)數(shù)(小范圍隨機(jī)波動)，Wp,i,j和Wc,i,j分別是參與克隆的父代特征矢量和子代特征矢量。最后，從由一個父代和Nij個子代組成的新種群中選取最優(yōu)秀的一個個體作為第j+1代第i個個體。

進(jìn)化論規(guī)則中的匹配法：第j+1代第i個個體的特征向量可由式(4)確定：

(4)

式中，i=1,2…Nsp，Nsp是參與匹配的個體個數(shù)，q*是改進(jìn)的匹配系數(shù)，k(i)，s(i){1,2…Nap}，nl是平均值為零的高斯隨機(jī)數(shù)，Wm,i,j，Wc,k(i),j是參與匹配法的父代特征矢量和子代特征矢量。

為了提高算法的搜索效率，提出進(jìn)化代數(shù)與各個體適應(yīng)度均相關(guān)的進(jìn)化系數(shù)計算新準(zhǔn)則，以克服定值進(jìn)化系數(shù)算法引起的即使到了運算最后的迭代期間，仍然存在較大的隨機(jī)波動性的缺點[10]。新的克隆系數(shù)和匹配系數(shù)計算準(zhǔn)則，分別如式(5)和式(6)所示：

(5)

(6)

2.2各進(jìn)化論自適應(yīng)消噪算法的定義

文獻(xiàn)[12]已經(jīng)論證了進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法比LMS、GA等自適應(yīng)算法具有更好的全局收斂性能，因此本文重點研究具有不同適應(yīng)度函數(shù)和不同進(jìn)化系數(shù)計算方法時的進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法。各算法的定義，如下表1所示。

表1　各種進(jìn)化論數(shù)字濾波自適應(yīng)消噪算法列表

3基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法特性研究

3.1新DC-EDF算法消噪性能和收斂特性分析

以式(1)的仿真信號作為分析信號，研究DC-EDF算法的消噪性能和收斂特性。各進(jìn)化參數(shù)的選取準(zhǔn)則為[11]：父代種群個數(shù)為64，子代種群個數(shù)為16，進(jìn)化系數(shù)恒定值為0.9，總預(yù)期進(jìn)化迭代數(shù)J為100。動態(tài)適應(yīng)度函數(shù)中c=1.5。

3.1.1新DC-EDF算法消噪性能分析

DC-EDF算法最后最優(yōu)輸出信號的時域波形圖，如圖6所示。與圖2所示的被噪聲嚴(yán)重污染的主輸入信號相比，圖6(a)所示的消噪后的輸出信號顯示，DC-EDF算法能有效地抑制主輸入信號中與參考信號相關(guān)的噪聲信號，凸顯特征信號，以進(jìn)行準(zhǔn)確的故障診斷。

3.1.2新DC-EDF算法收斂特性分析

每代中最優(yōu)個體所對應(yīng)的適應(yīng)度值與進(jìn)化代數(shù)的關(guān)系，如圖7所示。圖7顯示，最初最優(yōu)個體所對應(yīng)的輸出信號的峭度值為3，經(jīng)過約20代的迭代進(jìn)化后，快速達(dá)到了峭度最大值18，且隨后穩(wěn)定收斂，表明新DC-EDF算法具有穩(wěn)定和快速的收斂特性。

圖6 輸出信號與模擬特征信號對比圖Fig.6Comparisonoftheoutputsignalsandfeaturesignal圖7 每代中最優(yōu)個體對應(yīng)的適應(yīng)度值Fig.7Fitnessevolutionofthebestoutputsignal

圖8 輸出信號的時域波形對比Fig.8Comparisonoftheoutputsignalsfromdifferentalgorithms圖9 不同算法,最優(yōu)個體對應(yīng)的輸出信號峭度值與進(jìn)化迭代的關(guān)系Fig.9KurtosisevolutioncomparisonamongthealgorithmsofOR-EDF,CI-EDF,DF-EDFandDC-EDF

3.2不同算法消噪性能和收斂特性對比分析

3.2.1不同算法消噪性能分析

表1中所列不同算法最優(yōu)輸出信號的時域波形對比圖，如圖8(a)所示。為了清晰對比4種算法消噪效果，放大0.32~0.4 s區(qū)間的波形，如圖8(b)所示。其中，OR-SIG表示仿真特征信號。

圖8(a)顯示，不同算法均能有效地消除噪聲，從被噪聲污染的信號中提取出故障特征信號。圖8(b)的放大區(qū)間波形對比顯示，與仿真特征信號相比，最優(yōu)化消噪之后的輸出信號仍存在局部噪聲，但噪聲幅值已遠(yuǎn)小于特征信號的幅值數(shù)量級，且DC-EDF和CI-EDF算法較OR-EDF和DF-EDF算法有更好表現(xiàn)。

3.2.2不同算法收斂特性對比分析

對表1所列的不同算法，每代最優(yōu)個體對應(yīng)輸出信號的峭度值隨進(jìn)化迭代的關(guān)系，如圖9所示。

圖9顯示，經(jīng)各算法處理后的最優(yōu)輸出信號的峭度值均趨于恒定值，表明不同的進(jìn)化論消噪算法均收斂，即各算法均能有效地抑制噪聲，從噪聲污染的信號中提取特征信號。圖9同時顯示，算法DC-EDF較其他算法具有更快的收斂速度：算法DC-EDF對應(yīng)輸出信號的峭度值在進(jìn)化迭代15代時已經(jīng)收斂，而算法OR-EDF，CI-EDF和DF-EDF的收斂迭代數(shù)分別為37代，30代和35代。該結(jié)果表明，本文提出的基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法(DC-EDF)取得了最優(yōu)的收斂特性，較其他3種同類算法減少了約一半的迭代代數(shù)。

4新DC-EDF算法在軸承故障信號檢測中的應(yīng)用

采用美國西儲大學(xué)軸承故障模擬實驗臺獲取的滾動軸承故障數(shù)據(jù)[15]，驗證本文提出的DC-EDF算法的有效性。選取故障直徑為0.007英寸，電機(jī)轉(zhuǎn)速為1 797 r/min，編號為105的數(shù)據(jù)進(jìn)行分析，此時采樣頻率為12000Hz，內(nèi)圈的故障特征頻率為fi=162.2 Hz。將驅(qū)動端采集的加速度振動信號作為主輸入信號，靠風(fēng)扇端的振動信號作為參考輸入信號。主輸入信號和參考信號波形，分別如圖10和圖11所示。主輸入信號經(jīng)DC-EDF算法處理后的時域波形圖，如圖12所示。主輸入信號和經(jīng)DC-EDF算法處理后的輸出信號的包絡(luò)解調(diào)譜，分別如圖13和圖14 所示。圖13和圖14的結(jié)果顯示，經(jīng)DC-EDF算法處理后信號的故障特征頻率及其倍頻更明顯，能更有效地識別故障特征。

圖10　主輸入信號時域波形圖 Fig.10 Waveform of the primary input signal

圖11　參考信號時域波形圖 Fig.11 Waveform of the reference signal

圖12　經(jīng)DC-EDF算法處理后的輸出信號時域波形圖 Fig.12 Output signal using the proposed DC-EDF algorithm

圖13　主輸入信號包絡(luò)解調(diào)譜圖 Fig.13 Demodulation spectrum of the primary input signal

圖14　輸出信號包絡(luò)解調(diào)譜圖 Fig.14 Demodulation spectrum of the output signal

5結(jié)論

提出了基于動態(tài)適應(yīng)度函數(shù)的進(jìn)化論變系數(shù)自適應(yīng)消噪算法(DC-EDF)。動態(tài)適應(yīng)度函數(shù)能更好地描述個體的生存狀態(tài)，而與適應(yīng)度函數(shù)值和迭代代數(shù)相關(guān)的變進(jìn)化系數(shù)則能有效地提高進(jìn)化論算法消噪性能和收斂性能。模擬仿真分析研究結(jié)果表明，本算法可有效地提高收斂特性和保持消噪性能，較其他同類算法有較優(yōu)的表現(xiàn)。物理臺架實驗數(shù)據(jù)驗證了新DC-EDF算法對故障信號的有效提取能力。

參考文獻(xiàn)

[1]Lin J. Gearbox fault diagnosis using adaptive wavelet filter [J]. Mechanical Systems and Signal Processing, 2003, 17 (6): 1259-1269.

[2]Chiementin X, Bolaers F, Dron J P. Early detection of fatigue damage on rolling element bearings using adapted wavelet [J]. ASME Journal of Vibration and Acoustics, 2007, 129: 495-506.

[3]Panda G, Pradhan P M, Majhi B. IIR system identification using cat swarm optimization [J]. Expert Systems with Applications, 2011, 38: 12671-12683.

[4]Ibrahim G, Albarbar A, Abouhnik A, et al. Adaptive filtering based system for extracting gearbox condition feature from the measured vibrations [J]. Measurement, 2013, 46: 2029-2034.

[5]Zhou X J, Shao Y M. Adaptive noise cancellation based on beehive pattern evolutionary digital filter [J]. Mechanical Systems and Signal Processing, 2014, 42: 225-235.

[6]嚴(yán)保康,周鳳星. 一種基于形態(tài)提升的自適應(yīng)軸承微沖擊提取方法[J]. 振動與沖擊, 2013, 32(24): 198-203.

YAN Bao-kang, ZHOU Feng-xing. Adaptive weak impulse extraction method of rolling bearings based on morphological lifting wavelet[J]. Journal of Vibration and Shock, 2013, 32(24):198-203.

[7]陽子婧,蔡力鋼,高立新.自適應(yīng)冗余提升小波降噪分析及軸承故障識別應(yīng)用[J]. 振動與沖擊,2013,32(7): 54-57.

YANG Zi-jing, CAI Li-gang, GAO Li-xin. Adaptive redundant lifting wavelet denoising analysis and its application in bearing fault identification [J].Journal of Vibration and Shock, 2013, 32(7): 54-57.

[8]楊輝躍, 涂亞慶, 張海濤,等. 振動信號頻率跟蹤的反饋修正自適應(yīng)陷波器法[J]. 振動與沖擊, 2014, 33(3):145-149.

YANG Hui-yue, TU Ya-qing, ZHANG Hai-tao, et al. Feedback corrected adaptive notch filter for vibration signal frequency tracking [J]. Journal of Vibration and Shock, 2014, 33(3): 145-149.

[9]Kawamata M, Abe M, Higuchi T. Evolutionary digital filters-adaptive digital filters based on evolutionary strategies of organisms[C]. Proceedings of International Workshop on Intelligent Signal Processing and Communication Systems, 1994: 263-268.

[10]邵毅敏,周曉君,馬祥華,等.進(jìn)化論自適應(yīng)濾波算法及降噪特性研究[J].振動與沖擊, 2007, 26(12): 157-159.

SHAO Yi-min, ZHOU Xiao-jun, MA Xiang-hua, et al. Noise cancellation using an evolutionary adaptive algorithm [J]. Journal of Vibration and Shock, 2007, 26(12): 157-159.

[11]Abe M, Kawamata M. Convergence behavior of evolutionary digital filters on a multi-peak surface[C]. In Proceedings of IEEE International Symposium on Circuits and Systems, 1996, 2: 185-188.

[12]Abe M, Kawamata M. Comparison of the convergence of IIR evolutionary digital filters and other adaptive digital filters on a multiple-peak surface [C]. Conference Record of the Thirty-First Asilomar Conference on Signal, Systems, and Computers, 1997, 2: 1674-1678.

[13]Sheen Y T. A complex filter for vibration signal demodulation in bearing defect diagnosis [J]. Journal of Sound and Vibration, 2004, 276: 105-119.

[14]Goldberg D E. Genetic algorithms in search, optimization and machine learning [M]. Massachusetts: Addison-Wesley Professional, 1989.

[15]http://csegroups.case.edu/bearingdatacenter/pages.