魏永合　牛保國　劉雪麗　趙旭寧

沈陽理工大學(xué)，沈陽，110159

波形曲率延拓在局域均值分解中的應(yīng)用

魏永合牛保國劉雪麗趙旭寧

沈陽理工大學(xué)，沈陽，110159

摘要：局域均值分解(LMD)是一種能夠?qū)?fù)雜的調(diào)幅調(diào)頻信號(hào)自適應(yīng)地分解為一系列單分量的調(diào)幅調(diào)頻信號(hào)的處理方法，其分解過程存在端點(diǎn)效應(yīng)，分解結(jié)果有一定程度的失真。針對(duì)此問題，提出根據(jù)波形曲率特征對(duì)信號(hào)端點(diǎn)進(jìn)行極值延拓，通過特征波的曲率波動(dòng)來篩選與邊界波形最為相似的數(shù)據(jù)段，在此基礎(chǔ)上將波形匹配曲率估計(jì)應(yīng)用于LMD分解過程中, 并與鏡像延拓及自適應(yīng)波形匹配延拓方法相比較，驗(yàn)證了所提方法的優(yōu)點(diǎn)。使用仿真信號(hào)與實(shí)際的齒輪故障數(shù)據(jù)進(jìn)行試驗(yàn)與檢測(cè)，結(jié)果表明，所提方法可以有效改善LMD分解過程的端點(diǎn)效應(yīng)，提高分解精度。

關(guān)鍵詞：局域均值分解(LMD)；端點(diǎn)效應(yīng)；波形匹配；曲率特征

0引言

Smith[1]于2005年首次提出局域均值分解(local mean decomposition，LMD)方法。作為能夠處理信號(hào)數(shù)據(jù)的新的時(shí)域故障診斷方法，LMD可以自適應(yīng)地處理繁瑣的非線性、非平穩(wěn)、多分量的調(diào)幅調(diào)頻信號(hào)，將其分解為多個(gè)乘積函數(shù)(product function,PF)分量。每一個(gè)PF分量都由一個(gè)包絡(luò)信號(hào)和一個(gè)純調(diào)頻信號(hào)相乘組成，包絡(luò)信號(hào)就是該P(yáng)F分量的瞬時(shí)幅值，由純調(diào)頻信號(hào)可以直接求出該P(yáng)F分量的瞬時(shí)頻率。進(jìn)一步組合所有乘積函數(shù)的瞬時(shí)幅值和瞬時(shí)頻率，便可以獲得原始信號(hào)完整的時(shí)頻分布。

經(jīng)典LMD方法自身有一些缺陷，端點(diǎn)效應(yīng)就是其中之一。對(duì)于迭代篩選分解信號(hào)方法中存在的端點(diǎn)效應(yīng)問題，許多學(xué)者已提出了很多改進(jìn)方法。文獻(xiàn)[2-4]針對(duì)LMD分解過程，提出使用自適應(yīng)波形匹配方法來延拓原始信號(hào)，然后再進(jìn)行LMD分解，該方法抑制了端點(diǎn)效應(yīng)，并在工程實(shí)際中得到應(yīng)用；文獻(xiàn) [5]提出在信號(hào)中添加窗函數(shù)來抑制經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition, EMD)過程的端點(diǎn)效應(yīng)。文獻(xiàn)[6]對(duì)比分析了端點(diǎn)效應(yīng)在LMD和EMD分解過程中對(duì)分解結(jié)果的影響，提出了一種使用鏡像延拓方法抑制端點(diǎn)效應(yīng)的處理技術(shù)。文獻(xiàn)[7]提出了基于對(duì)稱差分能量算子解調(diào)的LMD端點(diǎn)效應(yīng)抑制方法。文獻(xiàn)[8]針對(duì)局域波分解端點(diǎn)效應(yīng)問題提出使用三角波形相似估計(jì)法對(duì)信號(hào)端點(diǎn)局部均值進(jìn)行估計(jì)。

經(jīng)過人們的不斷研究，已經(jīng)出現(xiàn)了許多抑制端點(diǎn)效應(yīng)的技術(shù)方法,但是這些方法或多或少存在一些缺陷，如鏡像延拓方法通常需要將鏡像點(diǎn)置于端點(diǎn)極值處，但實(shí)際中往往無法確定端點(diǎn)是否為極值，因此通常會(huì)導(dǎo)致結(jié)果具有一定程度的失真。本文在目前經(jīng)典抑制端點(diǎn)效應(yīng)方法的基礎(chǔ)上，提出一種基于信號(hào)波形曲率特征的估計(jì)匹配法來抑制LMD分解過程中的端點(diǎn)效應(yīng)，并將該方法與自適應(yīng)波形曲率延拓以及鏡像延拓方法進(jìn)行對(duì)比，選取相關(guān)性數(shù)據(jù)以及能量評(píng)價(jià)指標(biāo)[5-6]綜合說明該方法的優(yōu)點(diǎn)。在此基礎(chǔ)上將其應(yīng)用到LMD分解中，通過仿真信號(hào)與工程實(shí)際案例分析，證明該方法能夠有效地抑制端點(diǎn)效應(yīng)對(duì)LMD分解結(jié)果的影響，顯著提高分解精度，從而為準(zhǔn)確提取故障特征提供參考。

1LMD方法及端點(diǎn)效應(yīng)

1.1LMD基本原理

LMD是一種能夠有效處理非線性、非平穩(wěn)信號(hào)的方法，其具體步驟如下：

(1)求出原始信號(hào)x(t)的全部局部極值點(diǎn)ni，計(jì)算相鄰極值點(diǎn)ni和ni+1之間的均值mi:

(1)

將所有的均值mi進(jìn)行線性直線延拓，獲取局部均值線段，對(duì)局部均值線段進(jìn)行平滑處理，獲取局部均值函數(shù)m11(t)。

(2)根據(jù)局部極值點(diǎn)及其時(shí)刻，計(jì)算局部幅值ai：

(2)

將獲得的局部幅值ai進(jìn)行線性直線延拓，獲取信號(hào)的局部幅值線段，對(duì)局部幅值線段進(jìn)行平滑處理，獲取信號(hào)的包絡(luò)估計(jì)函數(shù)a11(t)。

(3)從原始信號(hào)x(t)中分離出局部均值函數(shù)m11(t)，對(duì)分離后的信號(hào)按照下式進(jìn)行解調(diào)處理：

(3)

h11(t)=x(t)-m11(t)

(4)

假如s11(t)的包絡(luò)估計(jì)函數(shù)不等于1，則對(duì)s11(t)重復(fù)上述迭代過程，直到s1n(t)表現(xiàn)出純調(diào)頻信號(hào)的特征，即它的包絡(luò)函數(shù)值為1。所以有

(6)

(4)求包絡(luò)信號(hào)a1(t)，即

(7)

(5)計(jì)算乘積函數(shù)，即

PF1(t)=a1(t)s1n(t)

(8)

(6)從原始信號(hào)x(t)中分離PF分量，獲取新的信號(hào)u1(t)，將其作為對(duì)象數(shù)據(jù)重復(fù)上述迭代過程，迭代n次，直至un(t)呈現(xiàn)單調(diào)函數(shù)的特征為止，即

u1(t)=x(t)-PF1(t)u2(t) =u1(t)-PF2(t) ?un(t)=un-1(t)-PFn(t)üty????

(9)

原始信號(hào)x(t)被解調(diào)處理分解為 n個(gè)PF分量和一個(gè)余量un(t)，即

(10)

1.2端點(diǎn)效應(yīng)

LMD分解是以信號(hào)的局部極值點(diǎn)作為基礎(chǔ)，循環(huán)地多次對(duì)局部均值函數(shù)與局部包絡(luò)函數(shù)等進(jìn)行迭代處理。但是在迭代過程中，端點(diǎn)以外的極值點(diǎn)情況不能確定，這樣就會(huì)使分解結(jié)果在兩端處產(chǎn)生誤差，迭代次數(shù)越多，這種誤差就會(huì)越向信號(hào)內(nèi)部擴(kuò)散，對(duì)于低頻PF分量，由于時(shí)間尺度大，這種誤差影響就會(huì)更加明顯，這就是端點(diǎn)效應(yīng)。

由LMD原理及過程可知，產(chǎn)生端點(diǎn)誤差的主要原因是端點(diǎn)可能不是極值點(diǎn)。假設(shè)在求局部均值函數(shù)以及局部包絡(luò)函數(shù)時(shí)，所使用的端點(diǎn)均是原信號(hào)的極值點(diǎn)，這樣端點(diǎn)處的誤差就會(huì)得到抑制，端點(diǎn)效應(yīng)的影響就會(huì)減小。進(jìn)行LMD分解處理之前，先對(duì)對(duì)象信號(hào)數(shù)據(jù)進(jìn)行端點(diǎn)極值延拓，以求更加真實(shí)地表現(xiàn)對(duì)象信號(hào)端點(diǎn)以外的極值分布情況，便于進(jìn)行后續(xù)迭代處理。分解完成以后，按照左右延拓的點(diǎn)數(shù)，再對(duì)分解結(jié)果進(jìn)行適當(dāng)?shù)慕厝?。這樣就可以抑制端點(diǎn)。

2波形曲率估計(jì)延拓及仿真

2.1波形曲率估計(jì)延拓

波形曲率估計(jì)匹配延拓是一種基于曲線采樣點(diǎn)曲率的端點(diǎn)極值預(yù)測(cè)方法。對(duì)信號(hào)端部進(jìn)行極值延拓，在一定程度上可以減輕端點(diǎn)效應(yīng)，其實(shí)現(xiàn)方法關(guān)鍵步驟如下：

(1)截取匹配特征波形，計(jì)算特征曲率。定義e(i)為待處理信號(hào)x(t)的極值點(diǎn)位置矩陣，選取x(1)-x(e(1))-x(e(2))波段為特征波段，采用多項(xiàng)式擬合方法擬合特征波段離散數(shù)據(jù)，按照下式計(jì)算波段內(nèi)采樣點(diǎn)附近曲率，將其作為特征曲率C：

(11)

(2)劃分采樣信號(hào)。以極值點(diǎn)為基準(zhǔn)、特征波段長(zhǎng)為步長(zhǎng)劃分整個(gè)采樣信號(hào)。對(duì)模擬信號(hào)進(jìn)行左端延拓，進(jìn)一步說明該延拓方法步驟。定義采樣信號(hào)x(t)的極大值、極小值位置的矩陣數(shù)組分別為max(i)、min(j)。其中i是從1至m的整數(shù)，m是極大值個(gè)數(shù)；j為從1至n的整數(shù)，n是極小值個(gè)數(shù)。假設(shè)信號(hào)左端第一個(gè)極值點(diǎn)是極大值點(diǎn)，則特征波段就是x(1)-x(max(1))-x(min(1))，特征波形長(zhǎng)度為l，圖1中波形段a為特征波形段。再以剩余信號(hào)內(nèi)部極小值序列為匹配波形段終點(diǎn)，向左取步長(zhǎng)為l的點(diǎn)為起點(diǎn)，截取波形段為匹配波形段，圖1中a1即為匹配的第一個(gè)波形段，以此類推。

圖1　波形段劃分示意圖

(3)選取指標(biāo)，尋找最優(yōu)波形段。選擇相關(guān)系數(shù)ρ及曲率匹配誤差ε作為指標(biāo)進(jìn)行波形段篩選。采用多項(xiàng)式擬合方法擬合每個(gè)波形段采樣數(shù)據(jù)，然后利用式(11)計(jì)算n個(gè)波形段內(nèi)采樣點(diǎn)附近的曲率Cn(i)。根據(jù)下面的公式分別計(jì)算特征波段曲率與匹配波段曲率之間的相關(guān)系數(shù)和匹配誤差：

(12)

(13)

式中，Cov(C，Ci)表示特征曲率C與第i個(gè)波形段曲率Ci的協(xié)方差；D(C)為特征波形曲率的標(biāo)準(zhǔn)差；D(Ci)為第i個(gè)匹配波形段內(nèi)曲率的標(biāo)準(zhǔn)差。

分別設(shè)定閾值ρ1、ε1。在匹配波形段集合里搜索出閾值ρC-Ci>ρ1且ε(i)<ε1的波形段為目標(biāo)波形段。為了更加準(zhǔn)確地找出符合信號(hào)內(nèi)部規(guī)律的波形段，閾值ρ1、ε1可分別選定為0.999、0.0001，具體閾值設(shè)定要視實(shí)際情況而定。

(4)以步驟(3)中的目標(biāo)波形段為基準(zhǔn)，截取兩個(gè)極值長(zhǎng)度數(shù)據(jù)段，延拓到信號(hào)一端，按照同樣的方法延拓?cái)?shù)據(jù)另一端。

2.2仿真實(shí)驗(yàn)

將波形曲率延拓方法應(yīng)用于LMD分解中，在進(jìn)行LMD分解之前進(jìn)行一次波形曲率延拓。進(jìn)行仿真實(shí)驗(yàn)檢驗(yàn)。選取仿真信號(hào)x(t)如下：

x(t)=x1(t)+x2(t)

(14)

x1(t)=(1+0.9(sin(0.04πt))cos(0.1πt+sin(0.01πt))

(15)

x2(t)=sin(0.03πt)cos(0.08πt)

(16)

x(t)為兩個(gè)分量信號(hào)的疊加，其中x1(t)為調(diào)幅調(diào)頻信號(hào)；x2(t)為調(diào)幅信號(hào)。t取值范圍為[0,500 s]。首先對(duì)原始信號(hào)x(t)進(jìn)行波形曲率延拓，延拓左右端兩個(gè)極值，由于信號(hào)的規(guī)律性特別強(qiáng)，故設(shè)定誤差閾值ρ1=0.9999，ε1=0.0001。

延拓前后的信號(hào)如圖2所示，可以看出，左右端點(diǎn)延拓后的信號(hào)趨勢(shì)非常符合原始信號(hào)內(nèi)部規(guī)律。

圖2　波形延拓前后信號(hào)圖像對(duì)比

求取延拓前后信號(hào)的局域均值線段以及局部包絡(luò)線段，如圖3所示。從圖3中可以看出，延拓前后的局部均值線段及局部包絡(luò)線段端部有所不同。

(a) 延拓前后信號(hào)局部均值線段

(b) 延拓前后信號(hào)局部包絡(luò)線段圖3　延拓前后信號(hào)線段對(duì)比

分別求取原始信號(hào)及延拓信號(hào)的局部均值函數(shù)和局部包絡(luò)函數(shù)，如圖4所示。從圖4中可以看出，原始信號(hào)的局部均值函數(shù)和局部包絡(luò)函數(shù)在端部有一定程度的失真。

(a)延拓前后信號(hào)局部均值函數(shù)

(b)延拓前后信號(hào)局部包絡(luò)函數(shù)圖4　延拓前后信號(hào)局部均值函數(shù)和包絡(luò)函數(shù)對(duì)比

曲率波形延拓后進(jìn)行LMD分解，得到4階PF分量如圖5所示。由圖5可以看出，信號(hào)的主要信息集中在前兩階PF分量PF1(t)、PF2(t)中，第3、4階PF分量PF3(t)、PF4(t)幅值明顯減小，不是信息集中區(qū)域，為L(zhǎng)MD分解的“虛假”成分；u4(t)是殘余分量信號(hào)。理論上，分解結(jié)果只有兩階，出現(xiàn)這些“虛假”成分的主要原因是求局部均值函數(shù)和局部包絡(luò)函數(shù)產(chǎn)生的誤差以及估計(jì)端點(diǎn)極值所產(chǎn)生的延拓誤差，這些誤差不可避免。

圖5　波形曲率延拓后再進(jìn)行LMD分解的結(jié)果

2.3與其他延拓方法對(duì)比

將波形曲率匹配延拓法與鏡像延拓法及自適應(yīng)波形匹配方法相比較。取式(14)的仿真信號(hào)x(t)為分析研究對(duì)象，分析參數(shù)的設(shè)置也相同。分別比較波形曲率估計(jì)匹配法、鏡像延拓法及自適應(yīng)波形匹配方法在LMD方法中的應(yīng)用及其對(duì)端點(diǎn)效應(yīng)的改進(jìn)程度。

2.3.1延拓處理及LMD分解對(duì)比

采用波形曲率匹配法、自適應(yīng)波形匹配法、鏡像方法對(duì)模擬采樣信號(hào)進(jìn)行延拓，延拓結(jié)果如圖6所示。由圖6可得出：鏡像延拓缺乏自適應(yīng)性，延拓結(jié)果和其他兩種方法相比，不是特別滿足信號(hào)內(nèi)部規(guī)律；波形曲率延拓和自適應(yīng)波形匹配延拓都能夠自適應(yīng)地選取最優(yōu)波形段進(jìn)行極值延拓，但是相比較來說，波形曲率延拓選取的波形段更加符合信號(hào)內(nèi)部規(guī)律。

圖6　各種延拓方法延拓后的信號(hào)

信號(hào)經(jīng)過波形曲率匹配延拓之后的LMD分解結(jié)果如圖5所示；信號(hào)經(jīng)過自適應(yīng)波形匹配法、鏡像延拓法延拓后的LMD分解結(jié)果分別如圖7、圖8所示；信號(hào)未經(jīng)延拓處理的LMD分解結(jié)果如圖9所示。

圖8　鏡像延拓后信號(hào)的LMD分解結(jié)果

圖9　未經(jīng)延拓信號(hào)的LMD分解結(jié)果

2.3.2端點(diǎn)效應(yīng)檢測(cè)

根據(jù)LMD分解結(jié)果分別選取各延拓情況下所得的PF分量，將PF分量分別與原始模擬信號(hào)x(t)之間進(jìn)行相關(guān)性分析，結(jié)果如表1所示。

根據(jù)LMD分解理論，原始模擬信號(hào)x(t)延拓之后通過LMD分解處理，其主要特征信息集中在前兩階PF分量PF1(t)、PF2(t)中。理論上，前兩階PF分量相關(guān)性越大則后兩階PF分量相關(guān)性越小，說明延拓效果越好則端點(diǎn)效應(yīng)影響越小。

表1　各PF分量與原始模擬信號(hào)間的相關(guān)系數(shù)

從表1中的相關(guān)性數(shù)值能明顯看出，經(jīng)過波形曲率延拓后再進(jìn)行LMD分解所得的第一階PF分量PF1(t)與x(t)相關(guān)性最大，高達(dá)0.963；后兩階PF分量與x(t)相關(guān)性非常小，幾乎可以忽略，說明波形曲率特征匹配法對(duì)端點(diǎn)效應(yīng)的抑制效果非常顯著。

采用文獻(xiàn)[5-6]提出的基于能量評(píng)價(jià)指標(biāo)檢驗(yàn)曲率延拓前后端點(diǎn)效應(yīng)對(duì)乘積函數(shù)的影響程度。信號(hào)有效值計(jì)算公式為

(17)

式中，R為信號(hào)的有效值；i為信號(hào)序列編號(hào)；n為信號(hào)長(zhǎng)度。

端點(diǎn)效應(yīng)影響因子θ定義為

(18)

根據(jù)θ定義可以得出，端點(diǎn)效應(yīng)對(duì)PF分量影響越小，θ值就越?。蝗绻麤]有產(chǎn)生端點(diǎn)效應(yīng)，則θ值為0。相應(yīng)的θ值結(jié)果如表2所示。

表2　延拓前后θ能量表

由表2可以得出結(jié)論：相對(duì)于其他兩種方法，波形曲率匹配延拓對(duì)LMD分解中端點(diǎn)效應(yīng)的抑制效果最好，可以有效提高分解精度。

3工程實(shí)驗(yàn)

為進(jìn)一步驗(yàn)證波形曲率匹配延拓對(duì)LMD分解端點(diǎn)效應(yīng)的改善情況，進(jìn)行了故障模擬檢測(cè)實(shí)驗(yàn)。試驗(yàn)選取的齒輪箱故障狀態(tài)是大齒輪點(diǎn)蝕。小齒輪連接電機(jī)輸入軸，電機(jī)軸轉(zhuǎn)動(dòng)速度為880 r/mim，實(shí)驗(yàn)所用減速器大小齒輪齒數(shù)分別為75和55。采樣頻率為51 200 Hz，采樣點(diǎn)數(shù)為256 000。根據(jù)輸入軸轉(zhuǎn)速及齒輪參數(shù)，計(jì)算齒輪理論嚙合頻率為806.67 Hz，大小齒輪轉(zhuǎn)頻分別為10.76 Hz 、14.67 Hz。

對(duì)原始數(shù)據(jù)進(jìn)行中心化、異常點(diǎn)剔除、小波降噪等一系列預(yù)處理，以降低干擾項(xiàng)的干擾，然后進(jìn)行波形曲率延拓。工程信號(hào)數(shù)據(jù)摻雜其他無關(guān)分量較多，雖然周期性明顯，但還是會(huì)影響延拓結(jié)果，所以經(jīng)過多次試驗(yàn)總結(jié)，分別將誤差閾值ρ1、ε1設(shè)置為0.95、0.01。由于數(shù)據(jù)較長(zhǎng)，故截取左右端部局部信號(hào)放大，如圖10所示，圖中實(shí)線為原始信號(hào)，虛線為端點(diǎn)延拓部分。

(a)左端點(diǎn)延拓

(b)右端點(diǎn)延拓圖10　原始信號(hào)左右延拓情況

對(duì)延拓后的信號(hào)及未經(jīng)波形曲率延拓的信號(hào)進(jìn)行LMD分解，結(jié)果如圖11、圖12所示。

圖11　波形曲率延拓后工程信號(hào)的LMD分解結(jié)果

圖12　未經(jīng)波形曲率延拓工程信號(hào)的LMD分解結(jié)果

使用相關(guān)系數(shù)理論方法來計(jì)算各個(gè)乘積函數(shù)與原始信號(hào)數(shù)據(jù)之間的相關(guān)情況，結(jié)果如表3所示。

表3　各PF分量與原始工程信號(hào)之間的相關(guān)系數(shù)表

由LMD分解原理可知，其分解過程實(shí)際上是從高頻到低頻的“篩選”過程，分解得到的PF分量的頻率是由高到低。在所有理論頻率中嚙合頻率及其倍頻是最大的，故第一個(gè)PF分量PF1(t)對(duì)應(yīng)的是嚙合頻率及其倍頻。由相關(guān)系數(shù)表可知，信號(hào)的主要特征信息集中在第一個(gè)PF分量PF1(t)中，故分別選取延拓之前信號(hào)及延拓之后信號(hào)經(jīng)LMD分解的第一個(gè)PF分量，對(duì)其作包絡(luò)譜分析，結(jié)果如圖13所示。

(a)延拓之前

(b)延拓之后圖13　延拓前后第一個(gè)PF分量PF1(t)的功率譜

由圖13可知，對(duì)未經(jīng)延拓的信號(hào)和經(jīng)過延拓之后的信號(hào)進(jìn)行特征信息提取，得到的特征故障頻率都能夠接近大齒輪的轉(zhuǎn)頻10.76 Hz，由此可判斷，大齒輪出現(xiàn)故障。但是未經(jīng)延拓得到的故障調(diào)制頻率是10.5 Hz，而經(jīng)過延拓之后得到的故障調(diào)制頻率是10.75 Hz，相比較而言，10.75 Hz與大齒輪的轉(zhuǎn)頻10.76 Hz更加接近，故可以得出結(jié)論：經(jīng)過波形曲率匹配延拓之后的信號(hào)能夠改善LMD分解過程中的端點(diǎn)效應(yīng)，提高分解精度，在后續(xù)的故障診斷中能夠更加精確地提取故障特征。

4總結(jié)

根據(jù)局域均值分解端點(diǎn)效應(yīng)產(chǎn)生的原因，提出使用波形曲率匹配延拓方法來從根源上抑制端點(diǎn)效應(yīng)。波形曲率匹配延拓法基于信號(hào)內(nèi)部的曲率波動(dòng)來匹配特征波形進(jìn)行端點(diǎn)延拓，同時(shí)設(shè)定閾值，提高了匹配精度以及估計(jì)精度。通過仿真信號(hào)與工程故障齒輪的實(shí)驗(yàn)，驗(yàn)證了信號(hào)經(jīng)波形曲率匹配延拓后再進(jìn)行局域均值分解，能夠?qū)Χ它c(diǎn)效應(yīng)進(jìn)行有效的抑制。

參考文獻(xiàn)：

[1]Smith J S. The Local Mean Decomposition and Its Application to EEG Perception Data[J]. Journal of the Royal Society Interface，2005，2(5)： 443-454．

[2]張亢，程軍圣，楊宇．基于自適應(yīng)波形匹配延拓的局部均值分解端點(diǎn)效應(yīng)處理方法[J]．中國機(jī)械工程，2010，21(4)：457-462．

Zhang Kang, Cheng Junsheng, Yang Yu. Processing Method for End Effects of Local Mean Decomposition Based on Self-adaptive Waveform Matching Extending[J]. China Mechanical Engineering, 2010,21(4): 457-462.

[3]張亢，程軍圣，楊宇．局部均值分解方法及其在軸承故障診斷中的應(yīng)用[J]．中國機(jī)械工程，2009，20(22)：2711-2717．

Zhang Kang, Cheng Junsheng, Yang Yu. Local Mean Decomposition Method and Its Application to Roller Bearing Fault Diagnosis[J].China Mechanical Engineering, 2009，20(22)：2711-2717．

[4]程軍圣，楊宇，于德介．局部均值分解方法及其在齒輪故障診斷中的應(yīng)用[J]．振動(dòng)工程學(xué)報(bào)，2009，22(1)：76-83．

Cheng Junsheng, Yang Yu, Yu Dejie. Local Mean Decomposition Method and Its Application in Gear Fault Diagnosis Fault Diagnosis[J]. Journal of Vibration Engineering, 2009，22(1)：76-83．

[5]任達(dá)千，吳昭同，嚴(yán)拱標(biāo)．EMD 端點(diǎn)效應(yīng)的評(píng)價(jià)指標(biāo)及抑制端點(diǎn)效應(yīng)的窗函數(shù)法[J]．制造業(yè)自動(dòng)化，2007，29(1)：21-24．

Ren Daqian, Wu Zhaotong, Yan Gongbiao. Evaluation of the EMD End Effect and Its Window Based Method[J]. Manufacturing Automation, 2007，29(1)：21-24．

[6]任達(dá)千，楊世錫，吳昭同，等．LMD 時(shí)頻分析方法的端點(diǎn)效應(yīng)在旋轉(zhuǎn)機(jī)械故障診斷中的影響[J]．中國機(jī)械工程，2012，23(8)：951-956．

Ren Daqian, Yang Shixi, Wu Zhaotong,et al. Research on End Effect of LMD Based Time-frequency Analysis in rotating Machinery Fault Diagnosis[J]. China Mechanical Engineering,2012,23(8):951-956．

[7]孟宗，李珊珊，季艷．基于對(duì)稱差分能量算子解調(diào)的局部均值分解端點(diǎn)效應(yīng)抑制方法[J]．機(jī)械工程學(xué)報(bào)，2014，50(13)：80-87．

Meng Zong, Li Shanshan, Ji Yan. Restraining Method for End Effect of Local Mean Decomposition Based on Energy Operator Demodulation of Symmetrical Differencing[J]. China Mechanical Engineering, 2014，50(13)：80-87．

[8]寧寧，張俊，秦文娟．抑制局域波分解端點(diǎn)效應(yīng)方法的研究[J]．計(jì)算機(jī)仿真，2008，25(11)：312-315．

Ning Ning, Zhang Jun, Qin Wenjuan. A Method for Restraining the Ending Effect in Local Wave Method[J]. Computer Simulation, 2008，25(11)：312-315．

(編輯蘇衛(wèi)國)

Waveform Curvature Extension in Applications of Local Mean Decomposition

Wei YongheNiu BaoguoLiu XueliZhao Xuning

Shenyang Ligong University,Shenyang,110159

Abstract:LMD was a kind of signal processing method to self-adaptively decompose the complex amplitude modulation-frequency modulation signals into a series of single-component amplitude modulation-frequency modulation. To some degrees, the decomposition results were of distortions because of the influences of end effects. In order to improve the end effects from LMD, the data segments which were most similar to wave boundary data segments were searched according to the curvature fluctuations of characteristic waves. The waveform curvature was applied in LMD. According to the comparison of mirror extension and self-adaptive waveform match, the advantages of waveform curvature extension were verified. The experiments of simulation signals and actual gear failure were made, the results show that this method can suppress the end effect effectively and improve the decomposition precision.

Key words：local mean decomposition(LMD); end effect; waveform matching;curvature characteristic

收稿日期：2015-08-31

基金項(xiàng)目：國家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2012AA041303)

中圖分類號(hào)：TH132.4；TH165.3

DOI：10.3969/j.issn.1004-132X.2016.03.011

作者簡(jiǎn)介：魏永合，男，1971年生。沈陽理工大學(xué)機(jī)械工程學(xué)院教授、博士。研究方向?yàn)闄C(jī)械設(shè)備的故障診斷、預(yù)測(cè)和健康管理。牛保國，男，1989年生。沈陽理工大學(xué)機(jī)械工程學(xué)院碩士研究生。劉雪麗，女，1990年生。沈陽理工大學(xué)機(jī)械工程學(xué)院碩士研究生。趙旭寧，女，1990年生。沈陽理工大學(xué)機(jī)械工程學(xué)院碩士研究生。