基于形態非抽樣融合與DCT高階奇異熵的液壓泵退化特征提取

基于形態非抽樣融合與DCT高階奇異熵的液壓泵退化特征提取

孫健，李洪儒，王衛國，許葆華

(軍械工程學院，石家莊050003)

摘要：針對軸向柱塞式液壓泵性能退化中振動信號非線性強、退化特征提取困難等問題，提出基于形態非抽樣融合與DCT(Discrete Cosine Transform)高階奇異熵的退化特征提取方法。在一般框架下提出形態非抽樣小波融合方法，通過構建特征能量因子篩選各分解層近似信號，據融合規則實現雙通道振動信號融合重構、改善重構信號的特征信息；并利用DCT高階譜分析法對融合信號進一步處理，通過奇異值分解分別計算Shannon、Tsallis奇異熵作為液壓泵性能退化特征向量；用仿真信號及液壓泵實測振動信號驗證該方法的有效性。

關鍵詞：退化特征提取；形態非抽樣小波融合；DCT；高階奇異熵

中圖分類號:Tp12；Tp13.3文獻標志碼：A

基金項目：國家自然科學基金資助項目(11302046，11172063)

收稿日期：2014-10-16修改稿收到日期：2014-11-20

基金項目：國家自然科學基金(11461069,51075346);國家重點基礎研究發展計劃973項目(2011CB706600)

收稿日期：2014-11-21修改稿收到日期：2015-07-14

Degradation feature extraction of hydraulic pump based on morphological undecimated decomposition fusion and DCT high order singular entropy

SUNJian,LIHong-ru,WANGWei-guo,XUBao-hua(Ordnance Engineering College, Shijiazhuang 050003, China)

Abstract:To solve the problem that vibration signals of a hydraulic pump as usual are strongly nonlinear and its degradation features are difficult to extract, a degradation feature extraction method based upon morphological undecimated wavelet decomposition fusion (MUWDF) and DCT high order singular entropy was proposed. The MUWDF algorithm was presented under the general framework of morphological undecimated decomposition. The approximate signals of all decomposition layers were selected by using the feature energy factor and dual-channel vibration signals were fused according to the presented fusion rules so as to increase the proportion of feature information. On this basis, a high order spectrum analysis algorithm modified by DCT was proposed for further dealing with the fused signal. Shannon and Tsallis singular entropies, which were considered as fault degradation features of hydraulic pump, were respectively achieved by singular value decomposition. Finally, the proposed method was verified by using simulation signals and real pump vibration signals in various working conditions.

Key words:degradation feature extraction; morphological undecimated wavelet decomposition fusion; DCT; high order singular entropy

液壓泵作為液壓系統關鍵部件之一，其性能好壞直接影響整個系統的可靠性。由于流體的壓縮性、泵源與伺服系統流固耦合作用及自身的大幅固有機械振動，液壓泵在性能退化過程中振動信號呈較強非線性、非平穩性，傳統的線性信號分析方法已不再適用[1]。需尋求有效的液壓泵故障特征提取方法，以提高故障預測的準確度。

信息融合技術因能對冗余、互補及協同信息進行綜合分析處理，提高系統的可靠性[2]，已廣泛用于信號分析[3]、故障診斷[4]等領域，且呈現極大優越性。考慮液壓泵自身結構特點，為提高故障預測準確度，提取更理想的退化特征信息，需采用合適的融合方法對多方向振動信號進行融合。常用方法有加權融合法[5-6]、卡爾曼濾波法[7-8]及小波分析法[9-10]。而加權融合法對權值系數的設置與調整存在較大主觀性，卡爾曼濾波法對非線性系統缺乏嚴格的濾波公式，小波分析法重構過程會導致特征部分信息遺漏。作為有效的非線性信號處理手段，形態非抽樣小波方法(Morphological Un-decimated Wavelet Decomposition, MUWD)省去傳統小波分解下抽樣及重構上抽樣，可避免因降噪引起的失真、信息遺漏問題[11]，在振動信號處理中得到一定程度應用[12-14]。傳統MUWD算法僅分析單個信號，以最高分解層近似信號作為處理結果，會致重構信號中特征信息比重不高。因此，引入融合思想對MUWD方法進行改進，以提高重構信號中的特征信息量。

目前，對液壓泵故障特征提取研究多集中于故障模式特征提取，如時、頻域特征等，而對性能退化特征提取的研究較少，導致現有特征無法有效反映退化過程，影響故障預測準確度。雙譜分析法含二階統計量沒有的大量信息，對信號中非對稱、非線性信息具有高度敏感性，雖一定程度上反映性能退化過程[15]，但僅以基于傅里葉系數的雙相干值為特征量，導致對退化程度的敏感性不理想。而離散余弦變換(Discrete Cosine Transform，DCT)為傅里葉變換的推廣[16]，具有能量聚集性質，其系數能敏感反應能量變化[17]。因此，將DCT引入高階譜分析中代替傅里葉變換過程，改善對退化程度的敏感性。

為此，本文提出基于形態非抽樣融合與DCT高階奇異熵的液壓泵退化特征提取方法。在形態非抽樣一般框架下提出形態非抽樣小波融合方法，利用特征能量因子據融合規則對雙通道振動信號進行融合；并用DCT代替傅里葉變換提出DCT高階譜分析方法，將Shannon、Tsallis奇異熵作為液壓泵退化特征向量表征性能退化；利用仿真信號及不同松靴程度下液壓泵實測數據對本文所提退化特征提取方法進行驗證。

1形態非抽樣小波融合方法

1.1MUWDF基礎算法

設集合Vi，Wi分別為第i層信號空間和第i層細節空間，T()為形態算子，傳統MUWD方法的基礎運算框架[18]為

(1)

(2)

(3)

分析MWUD基礎運算框架知，其關鍵運算點在于形態算子T()的選擇[21]。常用有平均組合算子、形態差值算子、形態梯度算子及混合算子[22-23]。考慮液壓泵自身結構特點及退化特征提取需求，本文選用形態差值算子。該算子含黑Top-Hat及白Top-Hat變換，可同時提取信號中的正負脈沖[24]。因此，MUWDF基礎運算可描述為

(4)

f(xi)°(i+1)g]

(5)

(6)

式中：f(x)為原始信號；°，?分別為形態開、閉運算；g為扁平結構元素，(i+1)g為對結構元素進行i次膨脹操作。

1.2基于MUWDF的振動信號融合

令X1表示振動傳感器1采集的通道1信號，X2表示通道2信號。設分解層數均為N。為綜合利用X1與X2信號各分解層所含特征冗余及互補信息，結合MUWDF基礎運算，給出基于MUWDF的振動信號融合方法。

FEF=(E1+E2+…+Eh)/E

(7)

式中：Eh為特征頻率第h倍頻處對應的能量值(本文h=3)；E為頻域總能量。

ki=wi/sum(wi)

(8)

據所得權重值對集合C中各層近似信號進行融合實現信號重構，即

xFinal=k1y1+k2y2+…+kNyN

(9)

不同分解層的融合規則即規則2。對各層篩選出的近似信號集合C，據FEF計算各自融合權重值，利用式(9)進行信號融合重構。

規則2為在規則1基礎上對多傳感器信號特征冗余及互補信息的綜合處理，利用近似信號實現融合重構，完成數據層決策，為特征提取奠定基礎。

2基于DCT高階奇異熵的液壓泵退化特征提取

2.1DCT高階奇異譜分析方法

作為高階統計量方法雙譜分析能有效抑制信號中高斯噪聲干擾，提取二階統計量無法表征的特征信息。設零均值平穩隨機過程{x(n)}的3階累積量[25]為

c3x(i,j)=E{x(n)x(n+i)x(n+j)}

(10)

(11)

作為傅里葉變換的推廣，DCT系數能更敏感反映能量變化。因此可用DCT變換代替雙譜分析的傅里葉變換，通過DCT高階譜分析提取反映性能退化的敏感參量。奇異熵作為常用譜分析參數恰好能表征該變化趨勢。考慮性能退化實際過程體現的系統復雜特性，本文分別選Shannon熵(Shannon Entropy，SE)及Tsallis熵(Tsallis Entropy，TE)作為譜參數特征對退化狀態進行度量，計算式[26]分別為

(12)

(13)

式中：M為微觀狀態系統數；k為Boltzmann常數；pi為第i個狀態發生概率；q為非廣度參數(一般q=3)。

從SE定義知，其描述的系統通常是廣度的，反應系統總體變化。而實際系統或多或少具有時間、空間相關性。TE為SE的廣義形式，能反映Shannon熵無法揭示的信息。基于此，在譜分析基礎上提取Shannon高階奇異熵(High-order Singular Shannon Entropy, HSSE)及Tsallis高階奇異熵(High-order Singular Tsallis Entropy, HSTE)作為退化特征。具體步驟為：

步驟1對MUWDF融合信號進行標準化處理，將所得{x(t),t=1,2,…,N}分為K段，每段含M個觀測值，即N=KM；

步驟2計算DCT系數，即

(14)

式中：xi(k)為第i段數據；k=0,1,…,M-1；i=1,2…,K；λ=1,…,M/2。

步驟3計算DCT系數的三重相關，即

k2)Y(i)(λ1+λ2+k1+k2)

(15)

式中：μ=fs/N；s1，s2為三階累積區間；λ1，λ2為正數且二者之和不超過采樣頻率fs的一半。

步驟4計算DCT高階譜估計矩陣，即

(16)

步驟5對估計矩陣進行SVD分解，計算各奇異值權重，即

B=PDQ

(17)

式中：D=diag(x1,x2,…xL)；L為B的階次。

(18)

步驟6分別計算Shannon、Tsallis定義的DCT高階奇異熵HSSE及HSTE，即

(19)

(20)

2.2基于DCT高階奇異熵退化特征提取

由以上分析，基于DCT高階奇異熵的液壓泵退化特征提取方法流程見圖1，具體為：①采集液壓泵雙通道振動信號X1、X2。②利用MUWDF方法對X1、X2進行融合，即確定分解層數N與結構元素初始長度L；據融合規則1篩選出每層近似信號，用式(7)計算相應特征能量因子值FEF；據融合規則2，用式(8)計算融合權重按式(9)進行融合，獲得重構信號xFinal。③以xFinal為基礎，利用DCT高階奇異熵提取液壓泵退化特征：即對xFinal進行標準化及分段處理；用式(14)分別計算每段數據的DCT系數；用式(15)計算DCT系數的三重相關，由式(16)～式(18)獲得估計矩陣各奇異值權重；用式(19)、(20)計算高階奇異熵HSSE及HSTE，作為液壓泵退化特征向量。

圖1　基于DCT高階奇異熵的退化特征提取流程 Fig.1 Flow chart of degradation feature extraction based on DCT high order singular entropy

3仿真信號分析

為驗證本文退化特征提取方法的有效性，用仿真信號進行研究。取采樣頻率fs=1 024 Hz，采樣時間t=10 s，采樣點數N=10 240，采集仿真信號為

y1(t)=0.1a1(t)+0.3a2(t)+0.5n(t)

y2(t)=0.2a1(t)+0.5a2(t)+0.3n(t)

可以看出，y1(t)，y2(t)主要由三部分組成，即性能退化信號a1(t)、諧波信號a2(t)及白噪聲n(t)。a1(t)為模擬滾動軸承性能退化中產生的周期性沖擊信號，沖擊頻率為f0=16 Hz，每周期內沖擊函數為t2e-200t×sin(512πt)，共振頻率為256 Hz；a2(t)模擬諧波信號cos(80πt)+cos(100πt)，含40 Hz、50 Hz兩個頻率成分，差頻10 Hz。通過設置故障信號、諧波信號及白噪聲信號的不同系數，組成仿真信號y1(t)及y2(t)模擬不同位置傳感器采集的雙通道振動信號。y1(t)，y2(t)的時、頻域圖分別見圖2、圖3。由圖3看出，仿真雙通道信號在共振頻率256 Hz處均有明顯的調制現象，且故障特征頻率16 Hz完全湮滅在噪聲中。

圖2 雙通道仿真信號時域圖Fig.2Timedomainofdual-channelsimulationsignals圖3 雙通道仿真信號頻域圖Fig.3Frequencydomainofdual-channelsimulationsignals圖4 雙通道仿真信號第一分解層近似信號Fig.4Approximatesignalsinthefirstdecompositionlayerofdual-channelsignals

利用本文基于MUWDF與DCT高階奇異熵方法提取退化特征。

3.1基于MUWDF的振動信號融合

用MUWDF方法對信號y1(t)，y2(t)融合，提高特征信息在重構信號中的比重。分解層數N及結構元素初始長度L不宜過大，否則會增加算法復雜性，降低冗余信息利用率。一般N取2～5，L取3～6。綜合考慮算法復雜性及特征信息利用率，本文取N=4，L=5。用式(4)～式(6)分別對信號y1(t)及y2(t)進行分解，由式(7)分別計算獲得二者各分解層近似信號的FEF，結果分別為{0.006 8, 0.008 6, 0.010 6, 0.028 4}及{0.001 6, 0.012 9, 0.012 2, 0.008 7}。

據融合規則1，對近似信號進行篩選。以信號第一層為例，見圖4。由圖4看出，二者特征信息均不明顯，特征頻率16 Hz及其2、3倍頻成分仍湮滅在干擾成分中。僅y1近似信號特征頻率稍顯突出。此現象與兩信號對應的FEF值相符合，y1近似信號FEF=0.006 8，y2的FEF= 0.001 6，兩值相對偏小，特征信息比重不高，且y1近似信號所含特征信息成分稍優于y2。因此，據融合規則1，在第一分解層中，保留y1的近似信號。同理，對各分解層，經規則1篩選后保留的近似信號為：y1的第一、四層，y2的第二、三層，對應FEF值為{0.006 8, 0.012 9, 0.012 2, 0.028 4}。

利用式(8)可得各層近似信號的融合權值為k1= 0.113，k2=0.214，k3=0.202，k4=0.471。據融合規則2，利用式(9)對所選4個近似信號進行融合，重構信號為

xFinal=0.113y11+0.214y22+0.202y23+0.471y14

式中：yij為信號yi第j層近似信號。

MUWDF對性能退化仿真雙通道信號的融合結果見圖5。由圖5看出，通過對y1、y2各分解層冗余、互補信息的綜合利用，特征頻率16 Hz及其2、3倍頻信息成分突出，而直頻干擾成分40 Hz、50 Hz得到有效抑制，且256 Hz處的調制現象也得到改善。計算xFinal的FEF=0.033 7，與各層近似信號相比，特征信息比重得到有效提高。

圖5　MUWDF對雙通道仿真信號融合效果 Fig.5 Fusion performance of dual-channel simulation signals based on MUWDF

3.2基于DCT高階奇異熵的退化特征提取

為模擬性能退化過程，將融合重構信號xFinal等分10段，形成10個信號{xF1,xF2,…,xF10}，表示退化程度由輕到重不同階段的振動信號，每個信號含1 024個數據點。利用DCT高階奇異熵分別提取xF1～xF10的退化特征。以xF1為例，對數據標準化處理后進行分段，數據點總數為1 024，分8段，每段含數據點128，即K=8，M=128。利用式(14)計算每段的DCT系數。利用式(15)對所得DCT系數進行三階累積，據式(16)～式(18)計算估計矩陣并進行SVD分解，得到奇異值矩陣D1。計算各奇異值所占比重，分別利用式(19)、(20)得到xF1的特征向量為{HSSE=4.654, HSTE=0.977}。按以上方法計算獲得其余xF2～xF10信號高階奇異熵構成特征向量，見表1。為對比說明本文所提特征向量的有效性，分別計算xF1～xF10信號奇異熵FSE(Fusion Singular Entropy)、雙譜奇異熵BFSE(Bispectral Fusion Singular Entropy)及源信號y1直接計算的奇異熵SSE(Single Singular Entropy)，結果見表1。

表1　不同退化階段特征向量值

選取xF1～xF10相應的奇異值矩陣前5個值作為主奇異值Q1～Q5，隨退化程度加劇，其變化見圖6。表1中各特征向量隨退化過程變化見圖7。

由圖6看出，隨退化程度加深，主奇異值曲線大體呈增加趨勢，尤以Q1最顯著，但從曲線細微變化可知，其斜率并未持續增加。在機械設備實際退化過程中，隨退化程度加劇直至失效前，其退化曲線應類似指數型增函數曲線，斜率呈持續增加趨勢。因此，主奇異值只能反映退化的基本規律，不能作為特征量準確衡量退化狀態。

由圖7看出，隨退化程度加劇設備趨于某種特定故障，熵值應呈負指數型曲線，斜率應持續增加。對比5條曲線可知，未經融合處理直接利用源信號計算的奇異熵SSE曲線，因受噪聲干擾熵值最高，且對退化過程變化敏感度不高；經融合處理后的奇異熵FSE曲線，因利用MUWDF融合方法能提高重構信號的特征信息比重，較SSE曲線熵值有所下降，但曲線變化不明顯；雙譜奇異熵BFSE由于采用雙譜分析法提取特征，一定程度上能提高對退化過程的敏感性，但退化后期仍存在斜率減小情況。本文所提高階奇異熵HSSE及HSTE，由于采用DCT系數進行三階累積，對能量變化具有良好敏感性，其曲線變化符合退化規律，斜率持續增加，對退化過程跟隨性良好，而HSSE對退化狀態變化的整體反映情況更好。

圖6 主奇異值隨退化程度加劇變化Fig.6Changesofmainsingularvaluesalongwithdegradationaggravating圖7 各特征量隨退化過程的變化Fig.7Changesofvariousfeaturevectorsalongwithdegradationaggravating圖8 歸一化處理后各特征量隨退化過程的變化Fig.8Changesofvariousfeaturevectorsalongwithdegradationaggravatingafternormalization

為從細節上對各特征量進一步對比，將表1數據歸一化處理，對退化過程細節的描述見圖8。由圖8看出，SSE曲線、FSE曲線、BFSE曲線變化趨勢不顯著； HSSE及HSTE呈斜率持續增加的變化趨勢，而HSTE敏感性更顯著。圖8直觀反映出各特征量對退化過程的敏感性。

仿真信號分析知，本文的退化特征提取方法通過對雙通道信號融合，能提高重構信號中特征信息比重，所提HSSE及HSTE特征向量對退化過程敏感性良好，能全面定量反映性能退化情況。

4實驗數據驗證

為驗證本文方法的有效性、實用性，將其用于液壓泵實測雙通道振動信號中。以松靴故障退化過程為例。通常，隨松靴程度加劇柱塞球頭與滑靴間隙不斷增加，因此可用間隙距離描述退化程度。選正常、松靴間隙0.15 mm、0.24 mm、0.38 mm、0.57 mm五種柱塞分別實驗，觀測松靴退化過程不同階段，見圖9。

圖9　實驗柱塞 Fig.9 Experimental pistons

圖10　振動傳感器安裝 Fig.10 Installation of vibration sensors

實驗用軸向柱塞式液壓泵型號為SY-10MCY14- 1EL，共7個柱塞(每次更換其中1個)，驅動電機型號為Y132M-4，額定轉速1 480 r/min，周期0.041 s。將振動加速度傳感器S-1、S-2安裝于泵端蓋正面、側面采集雙通道振動信號，見圖10。采樣頻率12 kHz，采樣時間1 s，用 DH-5920動態信號測試分析系統將采集的振動信號存入電腦，見圖11。

傳感器S-1采集的通道1信號為X1，傳感器S-2采集的通道2信號為X2。本文以0.24 mm松靴為例具體分析。X1、X2時域信號采集過程已顯示在圖11中，頻域見圖12。由圖11看出，在時域上雙通道信號X1、X2的每個周期內均有7個沖擊，主要由液壓泵結構決定，與其旋轉1周7個柱塞均沖擊斜盤情況相符。由圖12看出，在頻域上X1、X2均存在調制現象。由于液壓泵軸的實際轉速為1 480 r/min，轉軸頻率1 480/60=24.6 Hz，而柱塞泵有7個柱塞，液壓泵振動信號固有沖擊頻率為24.6×7=172.2 Hz，該頻率應同時存在于X1、X2中，且沖擊幅度應較大；而對松靴故障，其特征頻率理論上應等于或接近轉軸頻率24.6 Hz；但無論固有頻率信息或故障特征頻率信息均湮滅在噪聲信息中，僅在517 Hz處有1個較明顯的沖擊成分(517 Hz為固有頻率172.2 Hz的3倍頻)。

圖11　DH-5920動態測試分析系統采集雙通道信號 Fig.11 Acquisition of dual-channel signals by DH-5920 dynamic testing analysis system

圖12　雙通道松靴信號頻域圖 Fig.12 Frequency domain of dual-channel signals of loose slipper

利用本文的基于MUWDF與DCT高階奇異熵的方法提取退化特征。用MUWDF方法對信號X1、X2進行融合。取N=4，L=5，用式(4)～式(6)分別對信號X1、X2進行分解，據式(7)分別計算得X1、X2各分解層近似信號的FEF為{0.028 1, 0.047 7, 0.030 9, 0.031 8}及{0.021 1, 0.031 6, 0.054 4, 0.035 6}。據融合規則1，對近似信號進行篩選。保留近似信號為X1的第一、二層，X2的第三、四層，對應的FEF值為{0.028 1, 0.047 7, 0.054 4, 0.035 6}。利用式(8)得各層近似信號的融合權值為k1=0.169，k2=0.288，k3=0.328，k4=0.215。據融合規則2，利用式(9)對所選4個近似信號進行融合，重構信號為

xFinal=0.169X11+0.288X12+0.328X23+0.215X24

式中：Xij為信號Xi第j層近似信號。

利用MUWDF所得正常及不同松靴程度下的融合結果，見圖13。由圖13看出，固有沖擊頻率170.9 Hz(近似為172.2 Hz)及其倍頻得到較好體現；隨松靴退化程度加劇特征頻率24.28 Hz(近似為24.6 Hz)及其倍頻越明顯，特征信息所占比重不斷增大。

圖13　基于MUWDF的液壓泵雙通道振動信號融合效果 Fig.13 Fusion performance of dual-channel vibration signals of hydraulic pump based upon MUWDF

利用DCT高階奇異熵分別提取以上5種情況特征向量。按退化程度由輕到重對融合結果編號：{xF1,xF2,…,xF5}，每個信號含12 000個數據點。對每個信號經標準化處理后進行分段，數據點總數12 000，分60段，每段含數據點200，即K=60，M=200。利用式(14)計算每段的DCT系數。利用式(15)對所得DCT系數進行三階累積，據式(16)～式(18)計算估計矩陣并進行SVD分解，得到各奇異值權重。分別利用式(19)、(20)得到xF1～xF5信號的高階奇異熵，構成相應特征向量，見表2。為對比說明本文所提特征向量的有效性，分別計算xF1～xF5信號的FSE、FBSE及用單一源信號的SSE，結果亦匯于表2。

表2　不同松靴程度下特征向量值

xF1～xF5對應的主奇異值Q1～Q5隨退化程度加劇的變化見圖14，表2中各特征向量隨退化過程的變化見圖15，歸一化后各特征向量隨退化過程的變化見圖16。

由圖14看出，隨退化程度加深主奇異值曲線大體呈增加趨勢，尤以Q1最顯著，但從曲線的細微變化知，斜率未呈持續增加趨勢，尤其在松靴0.38 mm處，Q4、Q5甚至出現降低。因此，主奇異值對退化過程的跟隨能力有限，與仿真分析結論一致。

由圖15看出，隨松靴退化過程熵值大體呈降低趨勢，表明狀態趨于某種特定故障模式，與客觀規律一致。對比5條曲線知，隨退化程度加劇SSE曲線、FSE曲線、BFSE曲線變化不顯著，與其算法固有缺陷有關，導致對退化狀態的跟隨性、敏感性受限；本文的高階奇異熵HSSE及HSTE，由于綜合MUWDF與DCT高階譜分析優點，曲線斜率隨退化加深而呈持續增加趨勢，表明對退化過程較好的跟隨性。

圖14 主奇異值隨退化程度加劇的變化Fig.14Changesofmainsingularvaluesalongwithdegradationaggravating圖15 各特征量隨退化過程的變化Fig.15Changesofvariousfeaturevectorsalongwithdegradationaggravating圖16 歸一化后各特征量隨退化過程的變化Fig.16Changesofvariousfeaturevectorsalongwithdegradationaggravatingafternormalization

圖16從細節上反映出各特征量隨退化過程的變化。由圖16看出，SSE曲線、FSE曲線、BFSE曲線變化趨勢仍不顯著，而HSSE及HSTE呈斜率持續增加變化趨勢，對退化過程敏感性較好，與仿真分析結論一致。

5結論

(1)退化特征提取方法結合MUWDF與DCT高階譜分析優點，不僅能提高重構信號特征信息比重，且能利用DCT系數對能量的敏感能力及奇異熵對不同狀態的跟隨能力。

(2)利用DCT高階譜分析法提取Shannon、Tsallis奇異熵，對退化過程具有較高的敏感性。所提特征向量能有效反映液壓泵性能退化過程。對實現基于狀態的維修具有一定促進意義。

參考文獻

[1]唐宏賓,吳運新,滑廣軍,等. 基于EMD包絡譜分析的液壓泵故障診斷方法[J]. 振動與沖擊, 2012, 31(9):44-48.

TANG Hong-bin, WU Yun-xin, HUA Guang-jun,et al. Fault diagnosis of pump using EMD and envelope spectrum analysis [J]. Journal of Vibration and Shock, 2012, 31(9):44-48.

[2]Ciuonzo D, Papa G, Romano G, et al.One-bit decentralized detection with a rao test for multi-sensor fusion [J].IEEE Signal Processing Letters, 2013, 20(9):861-864.

[3]Safizadeh M S, Latifi S K.Using multi-sensor data fusion for vibration fault diagnosis of rolling element bearings by accelerometer and load cell[J].Information Fusion, 2014, 18(1): 1-8.

[4]Higger M, Akcakaya M, Erdogmus D.A robust fusion algorithm for sensor failure[J].IEEE Signal Processing Letters, 2013, 20(8):755-758.

[5]Angelov P, Yager R. Density-based averaging-a new operator for data fusion [J].Information Science, 2013, 222(10):164-174.

[6]Wei C M, Blum R S.Theoretical analysis of correlation-based quality measures for weighted averaging image fusion[J]. Information Fusion, 2010, 11(1):301-309.

[7]Rodger J A.Toward reducing failure risk in an integrated vehicle health maintenance system: a fuzzy multi-sensor data fusion kalman filter approach for IVHMS[J].Expert Systems with Applications, 2012, 39(6):9821-9835.

[8]Yang W B, Li S Y.A switch-mode information fusion filter based on ISRUKF for autonomous navigation of spacecraft[J]. Information Fusion, 2014, 18(3):33-42.

[9]任亞飛,柯熙政.多小波基多尺度多傳感器數據融合[J]. 傳感器與微系統,2010, 29(9):77-79.

REN Ya-fei, KE Xi-zheng. Multi-wavelet-basis multi-scale multi-sensor data fusion[J].Transducer and Microsystem Technologies, 2010, 29(9):77-79.

[10]呂新艷,顧曉輝.多傳聲器小波多尺度信息融合濾波算法[J].儀器儀表學報, 2013, 33(4):787-793.

Lü Xin-yan, GU Xiao-hui. Filtering algorithm for multi-microphones based on wavelet multi-scale information fusion [J].Chinese Journal of Science Instrument, 2012,33(4):787-793.

[11]Zhang J F, Smith J S, Wu Q H. Morphological un-decimated wavelet decomposition for fault location on power transmission lines[J].IEEE Transactions on Circuits and Systems,2006, 53(6):1395-1402.

[12]林勇,楊友東,劉健. 形態非抽樣小波在主減速器振動特征提取中的應用[J].農業機械學報, 2010, 41(2):209-214.

LIN Yong, YANG You-dong, LIU Jian. Feature extraction methods of vibration signal in automabile main reducer based on morphological un-decimated wavelet[J].Journal of Agricultural Engineering, 2010, 41(2):209-214.

[13]黃兵峰,沈路,周曉軍,等. 基于形態非抽樣小波分解的滾動軸承故障特征提取[J].農業機械學報, 2010, 41(2):204-207.

HUANG Bing-feng, SHEN Lu, ZHOU Xiao-jun, et al. Fault feature extraction of rolling element bearing based on morphological undecimated wavelet decompostion[J].Journal of Agricultural Engineering, 2010, 41(2):209-214.

[14]孫志輝,呂文泉. 基于形態非抽樣小波和S變換的軋機振動信號分析[J].北京科技大學學報, 2013, 35(3):366-370.

SUN Zhi-hui, Lü Wen-quan. Signal analysis of rolling mill vibration based on morphological un-decimated wavelets and s-transform[J].Journal of University of Sceince and Technology Beijing, 2013, 35(3):366-370.

[15]李富強,劉國華,吳志根. 基于雙譜和奇異值分解的結構損傷試驗[J].浙江大學學報, 2012, 46(10):1872-1879.

LI Fu-qiang, LIU Guo-hua, WU Zhi-gen. Experimental study of structural damage based on bispectral analysis and singular value decomposition [J].Journal of Zhejiang University, 2012, 46(10):1872-1879.

[16]陳彥龍,張培林,伍曉亮. DCT和相關分析在軸承故障診斷中的應用[J].機械強度, 2013, 35(2):138-141.

CHEN Yan-long, ZHANG Pei-lin, WU Xiao-liang. Research on DCT based correlation analysis and its use in the bearing fault diagnosis[J].Journal of Mechanical Strength, 2013,35(2):138-141.

[17]范新橋,朱永利. 基于高階多分辨率奇異熵的高壓輸電線路故障選相[J].電力自動化設備, 2011, 31(4):50-54.

FAN Xin-qiao, ZHU Yong-li. Faulty phase selection based on high-order multi-resolution singular entropy for high-voltage transmission lines[J].Electric Power Automation Equipment, 2011, 31(4):50-54.

[18]章立軍,陽建宏,徐金梧,等. 形態非抽樣小波及其在沖擊信號特征提取中的應用[J].振動與沖擊, 2007, 26(10):56-59.

ZHANG Li-jun, YANG Jian-hong, XU Jin-wu, et al. Morphological undecimated wavelet and its application to feature extraction of impulsive signal[J].Journal of Vibration and Shock, 2007, 26(10):56-59.

[19]Morenilla A J, Carmona R M, Romero J L S. Mathematical morphology for design and manufacturing[J].Mathematical and Computer Modeling, 2011, 54(7):1753-1759.

[20]Cousty J, Najman L,Dias F,et al. Morphological filtering on graphs[J].Computer Vision and Image Understanding, 2013, 117(4): 370-385.

[21]Li Bing, Zhang Pei-lin, Wang Zheng-jun. Gear fault detection using multi-scale morphological filters[J].Measurement, 2011, 44(10):2078-2089.

[22]He Wei, Jiang Zhi-nong, Qin Qiang. A joint adaptive wavelet filter and morphological signal processing method for weak mechanical impulse extraction[J].Journal of Mechanical Science and Technology, 2011, 24(8):1709-1716.

[23]Hao R J, Chu F L. Morphological un-decimated wavelet decomposition for fault diagnostics of rolling element bearings [J].Journal of Sound and Vibration, 2009, 4(1):1164-1177.

[24]王冰,李洪儒,許葆華. 一種多元素多尺度形態非抽樣小波分解方法[J].軸承, 2013, 11(5):43-48.

WANG Bing, LI Hong-ru, XU Bao-hua. A kind of mult-element and multi-scale morphological undecimated wavelet decomposition method[J].Bearing, 2013, 11(5):43-48.

[25]錢祖平,許淵,邵尉,等. 基于高階譜和時域分析的電臺穩態特征提取算法[J].電子與信息學報, 2013, 35(7):1599-1606.

QIAN Zu-ping, XU Yuan, SHAO Wei,et al. Extraction algorithm of radio steady state characteristics based on high order spectrum and time-domain[J].Journal of Electronics & Information Technology, 2013, 35(7):1599-1606.

[26]Fiori E R, Plastino A. A shannon-tsallis transformation [J]. Physica A : Statistical Mechanics and its Applications, 2013, 392(8): 1742-1749.

第一作者張宇飛女，碩士，副教授，1973年4月生

通信作者王延慶男，博士，副教授，1982年2月生

第一作者艾孜海爾·哈力克男，講師，博士生，1978年生

通信作者買買提明·艾尼男，博士，教授，博士生導師，1958年生

郵箱:mgheni@263.net