吳占濤，程軍圣，曾　鳴，鄭近德

(1.湖南大學(xué) 汽車車身先進設(shè)計制造國家重點實驗室，長沙　410082； 2.湖南大學(xué) 現(xiàn)代工程訓(xùn)練中心，長沙　4100823.安徽工業(yè)大學(xué) 機械工程學(xué)院，安徽 馬鞍山　243002)

?

基于分段多項式的局部特征尺度分解方法及應(yīng)用

吳占濤1,2，程軍圣1，曾鳴1，鄭近德3

(1.湖南大學(xué) 汽車車身先進設(shè)計制造國家重點實驗室，長沙410082； 2.湖南大學(xué) 現(xiàn)代工程訓(xùn)練中心，長沙4100823.安徽工業(yè)大學(xué) 機械工程學(xué)院，安徽 馬鞍山243002)

摘要：針對局部特征尺度分解(Local Characteristic-scale Decomposition,LCD)方法中兩極值點連線在極值點處一階微分不連續(xù)，引起分解精度降低問題，提出新非平穩(wěn)信號分析方法-基于分段多項式的局部特征尺度分解(Piecewise Polynomial based Local Characteristic-scale Decomposition, PPLCD)。用分段多項式取代LCD中直線連接，且均值曲線插值點由相鄰3個同類極值點構(gòu)成的多項式計算產(chǎn)生。通過仿真信號將PPLCD與LCD對比，結(jié)果表明，PPLCD在提高分量正交性、精確性等具有一定優(yōu)越性；由轉(zhuǎn)子碰摩故障診斷表明該方法的有效性。

關(guān)鍵詞：局部特征尺度分解；分段多項式；故障診斷；轉(zhuǎn)子；內(nèi)稟尺度分量

局部特征尺度分解(Local Characteristic-scale Decomposition，LCD)為新的時頻分析方法[1]，能自適應(yīng)將1個復(fù)雜信號分解為若干瞬時頻率具有物理意義內(nèi)稟尺度分量(Intrinsic Scale Component，ISC)之和[2]，獲得原始信號完整的時頻分布。與經(jīng)驗?zāi)B(tài)分解[3](Empirical Mode Decomposition，EMD)相比， LCD在減小擬合誤差、提高計算速度、抑制模態(tài)混淆等有一定優(yōu)勢[4]，已用于信號分析及機械故障診斷等，并取得較好效果[5-7]。EMD與LCD基于篩分的自適應(yīng)時頻分析方法關(guān)鍵為定義合理的均值曲線，該曲線定義優(yōu)劣直接決定方法的有效性及精確性[8]。

LCD均值曲線插值點由連接兩相鄰?fù)悩O值點的連線計算產(chǎn)生，但存在連線在各極值點處一階微分不連續(xù)會致分解精度降低。為此，本文用分段多項式取代LCD中直線連接，且均值曲線插值點由相鄰3個同類極值點構(gòu)成的多項式計算產(chǎn)生，對LCD均值曲線進行改進，提出基于分段多項式的局部特征尺度分解(Piecewise Polynomial based Local Characteristic-scale Decomposition,PPLCD)，以解決LCD中的不連續(xù)問題。

1PPLCD

LCD用線性插值方法計算所有相鄰極大(或極小)值點確定的直線及與二者間極值點對應(yīng)的函數(shù)值，進而求出均值曲線插值點，存在連接直線在兩極值點處一階微分不連續(xù)問題。采用低次分段多項式擬合相鄰3個同類極值點，不僅在與均值曲線各插值點相鄰兩極值點處一階微分連續(xù)，擬合曲線能較好體現(xiàn)相關(guān)極值點數(shù)據(jù)的變化趨勢，且能規(guī)避用高次多項式插值可能產(chǎn)生的龍格現(xiàn)象。為此，通過研究PPLCD原理，本文提出基于分段多項式的局部特征尺度分解方法。

1.1PPLCD算法

對實信號x(t)(t>0)，將所有極值點記為(ti，Xi)(i=1,2,…,K)。其中(tk，Xk)、(tk+2，Xk+2)及(tk+4，Xk+4)為相鄰極大(小)值點，(tk+1，Xk+1)及(tk+3，Xk+3)為相鄰極小(大)值點，(k=1,2,…,K-4)。記函數(shù)y(t)在區(qū)間[tk,tk+4]部分為yk(t)，(k=1,2,…,K-4)。若yk(t)滿足條件

(1)yk(t)≥x(t)或yk(t)≤x(t)；

(1)

(2)

設(shè)Pk為條件(2)方程組系數(shù)矩陣，即

(3)

則該方程組可表示為

MkPk=Nk

(4)

由Mk反除Nk可求得系數(shù)矩陣Pk，即

PK=MkNk

(5)

令

(6)

得信號x(t)唯一確定的局部上(下)包絡(luò)線yk(t)

yk(t)=QkPk，t∈[tk,tk+4]

(7)

可求出yk(t)在兩相鄰極值點(tk+2，Xk+2)、(tk+4，Xk+4)間極值點對應(yīng)時刻tk+3處Ak+3點的函數(shù)值yk(tk+3)，見圖1。

圖1　PPLCD方法示意圖Fig.1 Schematic of PPLCD method

對應(yīng)的PPLCD均值曲線插值點Lk+3值為

Lk+3=(Xk+3+yk(tk+3))/2，

(k=1,2,…,K-4)

(8)

由于Lk值下標k從4～K-1，用端點延拓法[9-10]求得L1、L2、L3及LK值?；诖_定信號實信號x(t)所有Li值(i=1,2,…,K)，按LCD計算步驟[11]，將原始信號分解為n個ISC分量與1個趨勢項un(t)之和，即

(9)

1.2仿真信號分析

為說明PPLCD方法的優(yōu)越性，考慮式(9)混合信號

x(t)=x1(t)+x2(t)，(t∈[0,1])

(10)

式中：x1(t)=(1+0.35sin(25πt))sin(200πt+60πt2)，x2(t)=sin(85πt)，其時域波形見圖2。

圖2　仿真信號x(t)及各成分時域波形Fig.2 The time domain waveforms of simulation signal x(t) and its componments

分別采用LCD、PPLCD對仿真信號x(t)進行分解，結(jié)果見圖3、圖4，兩種方法分解絕對誤差見圖5。分解絕對誤差定義為分解所得ISC分量與真實分量之差絕對值。圖3中i1、i2、r2分別表示第1、2個ISC分量及剩余項，圖4中I1、I2、R2分別表示第1、2個ISC分量及剩余項，圖5中E1、E2分別表示兩種方法第1、2個ISC分量分解絕對誤差。由圖3～圖5看出，PPLCD分解結(jié)果較理想，分解分量與真實分量非常接近，分解絕對誤差較小。LCD分解結(jié)果雖亦較理想，但分量i1出現(xiàn)局部波形失真(圖3實線圈部分)，分量i2幅值與真實分量幅值相差較大，且兩分量絕對誤差略大于PPLCD。

圖3　仿真信號x(t)的LCD分解結(jié)果Fig.3 The LCD decomposition results of simulation signal x(t)

圖4　仿真信號x(t)的PPLCD分解結(jié)果Fig.4 The PPLCD decomposition results of simulation signal x(t)

圖5　兩種方法分解仿真信號x(t)絕對誤差Fig.5 The absolute error of the LCD and PPLCD

為比較兩種方法的分解效果，考察其分解正交性指標(IO)[12]及所得前兩分量與真實分量均方根誤差(RMSE)及相關(guān)系數(shù)(CC)[12]。IO值越小表示所有分解分量間正交性越好；RMSE值越小表示分解誤差越小；CC值越大表示分解準確性越高。各評價指標值見表1，其中RMSEi、CCi分別表示第i個分解分量與對應(yīng)的真實分量均方根誤差及相關(guān)系數(shù)，i=1,2。由表1看出，與LCD方法相比，PPLCD方法正交性指標及均方根誤差指標值均最小，相關(guān)系數(shù)指標值均最大，說明PPLCD的正交性、精確性等具有一定優(yōu)越性。

用希爾伯特變換法提取兩種分解法分解的前兩分量瞬時頻率及瞬時幅值，結(jié)果見圖6、圖7。由兩圖看出，無論瞬時頻率或幅值，I1均較i1更準確，波動性更?。籭2的瞬時頻率與I2結(jié)果接近，效果均較好，但I2的瞬時幅值較i2更準確。進一步說明PPLCD分解方法有效、可行。

表1　兩分解方法評價指標

圖6　由LCD所得分量i1、i2瞬時頻率及幅值Fig.6 IFs and IAs of ISCi calculated by LCD

圖7　由PPLCD所得分量I1、I2瞬時頻率及幅值Fig.7 IFs and IAs of ISCi calculated by PPLCD

2應(yīng)用實例

為說明PPLCD方法的有效性、實用性，將其用于具有碰摩故障的轉(zhuǎn)子振動位移實驗信號分析。實驗采樣頻率fs=2 048 Hz，采樣時長0.5 s，轉(zhuǎn)速3 000 r/min，轉(zhuǎn)頻fr=50 Hz。實驗數(shù)據(jù)時域波形見圖8，幅值譜見圖9。轉(zhuǎn)子產(chǎn)生碰摩故障時其振動位移信號會產(chǎn)生1個以轉(zhuǎn)頻為幅值調(diào)制頻率的高頻碰摩分量[13-14]，由圖8可看出高頻碰摩信號淹沒于強大背景信號中。而圖9中主要為與轉(zhuǎn)頻fr相關(guān)的基頻分量，看不出碰摩故障信息。

圖8　實驗數(shù)據(jù)時域波形Fig.8 The time domain waveforms ofexperimental data

圖9　實驗數(shù)據(jù)幅值譜Fig.9 The amplitude spectrum of experimental data

圖10　實驗數(shù)據(jù)的PPLCD分解結(jié)果Fig.10 The PPLCD decomposition results of experimental data

為提取高頻碰摩信息，分別采用PPLCD及LCD對實驗數(shù)據(jù)進行分解，結(jié)果見圖10、圖11。PPLCD、LCD分解所得第1個分量I1、i1的包絡(luò)譜見圖12、圖13。由圖10看出，PPLCD方法對實驗數(shù)據(jù)分解的第1分量I1調(diào)制特征明顯。由圖12可看到明顯的轉(zhuǎn)頻及其倍頻成分，說明I1的調(diào)制波頻率與轉(zhuǎn)頻fr相同，此由轉(zhuǎn)子每旋轉(zhuǎn)一周動件均與靜件摩擦所致[15-16]；I1分量主要成分為碰摩信號，亦會包含噪聲。由圖11看出，LCD方法對實驗數(shù)據(jù)分解的第1分量i1亦有調(diào)制特征，而圖13包絡(luò)譜雖亦能看出轉(zhuǎn)頻及其倍頻成分，但幅值較小，無PPLCD方法明顯，且出現(xiàn)兩個較明顯的不明低頻分量。

圖11　實驗數(shù)據(jù)的LCD分解結(jié)果Fig.11 The LCD decomposition results of experimental data

圖12　PPLCD第1分量包絡(luò)譜Fig.12 The envelope spectrum of calculated by PPLCD

圖13　LCD第1分量包絡(luò)譜Fig.13 The envelope spectrum of calculated by LCD

PPLCD的分量I2、I3為與轉(zhuǎn)頻有關(guān)的背景信號，其幅值譜見圖14。由圖14看出，I2、I3所含主要頻率成分為轉(zhuǎn)頻的1倍頻及1/2倍頻。剩余信號R3為低頻噪聲。圖11中LCD分量i2、i3、r3出現(xiàn)嚴重的波形失真(實線所圈部分)，其幅值譜見圖15。由圖15看出，分量i2的主要頻率為轉(zhuǎn)頻的1倍頻，同時也出現(xiàn)轉(zhuǎn)頻的1/2倍頻；分量i3主要頻率成分為轉(zhuǎn)頻的1/2倍頻、1倍頻。由于LCD對實驗數(shù)據(jù)分解時分量i2、i3出現(xiàn)波形失真，分別混入轉(zhuǎn)頻1/2倍頻、1倍頻的頻率成分，可見與圖11分解結(jié)果一致。

圖14　PPLCD分量I2、I3幅值譜Fig.14 The amplitude spectrum of I2 and I3 calculated by PPLCD

圖15　LCD分量i2、i3幅值譜Fig.15 The amplitude spectrum of i2 and i3 calculated by LCD

由此可認為，PPLCD能更有效分離出實驗數(shù)據(jù)中與碰摩相關(guān)的高頻分量，便于提取故障信息，實現(xiàn)轉(zhuǎn)子系統(tǒng)故障診斷。而該方法亦有不足之處，也存在模態(tài)混淆等問題，尚待進一步完善。

3結(jié)論

通過基于分段多項式局部特征尺度分解(PPLCD)對仿真信號及實驗數(shù)據(jù)分析，結(jié)論如下：

(1)PPLCD所得分量更具有物理意義，分量與真實值間誤差更小。

(2)PPLCD無需預(yù)設(shè)運行參數(shù)，為自適應(yīng)的信號分析方法。

(3)PPLCD能有效地提取振動信號中與故障特征有關(guān)的分量，且效果優(yōu)于LCD。

參 考 文 獻

[ 1 ] 程軍圣,鄭近德,楊宇. 一種新的非平穩(wěn)信號分析方法-局部特征尺度分解法[J]. 振動工程學(xué)報, 2012, 25(2): 215-220.

CHENG Jun-sheng,ZHENG Jin-de,YANG Yu. A nonstationary signal analysis approach-the local characteristic-scale decomposition method[J]. Journal of Vibration Engineering, 2012, 25(2): 215-220.

[ 2 ] 楊宇,曾鳴,程軍圣. 局部特征尺度分解方法及其分量判據(jù)研究[J]. 中國機械工程, 2013, 24(2):195-201.

YANG Yu, ZENG Ming, CHENG Jun-sheng. Research on local characteristic-scale decomposition and its stopping criteria[J]. China Mechanical Engineering, 2013, 24(2):195-201.

[ 3 ] Lei Y, Lin J, He Z, et al. A review on empirical mode decomposition in fault diagnosis of rotating machinery[J]. Mechanical Systems and Signal Processing,2013,35(1):108-126.

[ 4 ] 楊宇,曾鳴,程軍圣. 局部特征尺度分解方法及其分解能力研究[J]. 振動工程學(xué)報, 2012, 25(5):602-609.

YANG Yu, ZENG Ming, CHENG Jun-sheng. Research on local characteristic-scale decomposition and its capacities[J]. Journal of Vibration Engineering, 2012, 25(5):602-609.

[ 5 ] Zheng Jin-de, Cheng Jun-sheng, Yang Yu. A rolling bearing fault diagnosis approach based on LCD and fuzzy entropy[J]. Mechanism and Machine Theory, 2013, 70: 441-453.

[ 6 ] 程軍圣,楊怡,楊宇. 局部特征尺度分解方法及其在齒輪故障診斷中的應(yīng)用[J]. 中國機械工程學(xué)報,2012,48(9):64-71.

CHENG Jun-sheng, YANG Yi, YANG Yu. Local characteristic-scale decomposition method and its application to gear fault diagnosis[J]. Journal of Mechanical Engineering, 2012, 48(9): 64-71.

[ 7 ] 鄭超,郭奇,郭麗杰,等. 基于局部特征尺度分解的旋轉(zhuǎn)機械故障欠定盲源分離方法研究[J]. 燕山大學(xué)學(xué)報,2014, 38(2): 168-174.

ZHENG Chao, GUO Qi, GUO Li-jie, et al. Underdetermined blind source separation method of rotating machinery faults based on local characteristic-scale decomposition[J]. Journal of Yanshan University, 2014, 38(2):168-174.

[ 8 ] Zheng Jin-de, Cheng Jun-sheng, Yang Yu. Generalized empirical mode decomposition and its applications to rolling element bearing fault diagnosis[J]. Mechanical Systems and Signal Processing, 2013, 40(1): 136-153.

[ 9 ] Cheng J, Yu D, Yang Y. Application of support vector regression machines to the processing of end effects of Hilbert-Huang transform[J]. Mechanical Systems and Signal Processing, 2007, 21(3): 1197-1211.

[10] Rilling G, Flandrin P, Goncalves P. On empirical mode decomposition and its algorithms[C]//IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing NSIP. Grodo, Italy,2003.

[11] 程軍圣,鄭近德,楊宇. 基于局部特征尺度分解的經(jīng)驗包絡(luò)解調(diào)方法及其在機械故障診斷中的應(yīng)用[J]. 機械工程學(xué)報, 2012, 48(19): 87-99.

CHENG Jun-sheng, ZHENG Jin-de, YANG Yu. Empirical envelope demodulation approach based on local characteristic-scale decomposition and its applications to mechanical fault diagnosis[J]. Journal of Mechanical Engineering, 2012, 48(19): 87-99.

[12] Huang N E, Shen Z, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 1998, 454(1971): 903-995.

[13] 馬輝,太興宇,牛和強,等. 轉(zhuǎn)定子靜態(tài)不對中條件下碰摩故障仿真[J]. 振動工程學(xué)報, 2013, 26(2):260-268.

MA Hui, TAI Xing-yu, NIU He-qiang, et al. Rub-impact fault simulation of a rotor system under static misalignment conditions[J]. Journal of Vibration Engineering,2013,26(2):260-268.

[14] Yang Y, Cheng J S, Zhang K. An ensemble local means decomposition method and its application to local rub-impact fault diagnosis of the rotor systems[J]. Measurement, 2012, 45(3): 561-570.

[15] 韓清凱,于濤,王德友,等. 故障轉(zhuǎn)子系統(tǒng)的非線性振動分析與診斷方法[M]. 北京:科學(xué)出版社, 2010:37-65.

[16] 聞邦椿,武新華,丁千,等. 故障旋轉(zhuǎn)機械非線性動力學(xué)的理論與試驗[M]. 北京:科學(xué)出版社, 2004:112-135.

Piecewise polynomial based local characteristic-scale decomposition and its application

WUZhan-tao1,2,CHENGJun-sheng1,ZENGMing1,ZHENGJin-de3

(1. State Key Laboratory of Advanced Design and Manufacture for Vehicle Body, Hunan University, Changsha 410082, China;2.Modern Engineering Training Center, Hunan University, Changsha 410082, China;3. School of Mechanical Engineering, Anhui University of Technology, Ma’anshan 243002, China)

Abstract:A novel nonstationary signal analysis method, piecewise polynomial based local characteristic-scale decomposition (PPLCD), was proposed for improving the shortcomings of conventional local characteristic-scale decomposition (LCD), in which the first derivative at extremum points on the connecting line of two extreme value points is discontinuous, thereby the decomposition precision is lowered. Piecewise polynomial was used in PPLCD to replace line connection in LCD, and the interpolation points on the mean curve were computed by the polynomial which was generated according to three adjacent similar extremum points. The theory of PPLCD was studied and then simulation experiments were used to compare the PPLCD with the LCD. The results indicate that the PPLCD is more efficient in improving the orthogonality and veracity of components than the LCD. The proposed method was applied to diagnose the rotor with rub-impact fault successfully which indicates the effectiveness of PPLCD.

Key words:local characteristic-scale decomposition; piecewise polynomial; fault diagnosis; rotor; intrinsic scale component

中圖分類號:TH165+.3；TH911.7

文獻標志碼：A

DOI:10.13465/j.cnki.jvs.2016.04.037

通信作者程軍圣 男，博士，教授，博士生導(dǎo)師，1968年10月生

收稿日期：2015-01-09修改稿收到日期：2015-03-31

基金項目：國家自然科學(xué)基金(51375152；51575168)；湖南省科技計劃(2014WK3005)；智能型新能源汽車國家2011協(xié)同創(chuàng)新中心；湖南省綠色汽車2011協(xié)同創(chuàng)新中心資助項目

第一作者 吳占濤 男，博士，工程師，1982年6月生