張?jiān)茝?qiáng)，張培林，李　兵

(軍械工程學(xué)院車輛與電氣工程系，石家莊　050003)

第一作者張?jiān)茝?qiáng)男，博士生，1987年生

分?jǐn)?shù)階時(shí)頻局部二值模式譜在齒輪故障診斷中的應(yīng)用

張?jiān)茝?qiáng)，張培林，李兵

(軍械工程學(xué)院車輛與電氣工程系，石家莊050003)

摘要：針對齒輪故障信號分析，提出采用分?jǐn)?shù)階時(shí)頻LBP譜表達(dá)齒輪故障信號的時(shí)頻特征。首先為克服S變換對高頻信號的時(shí)頻分辨性能差的不足，基于分?jǐn)?shù)階Fourier變換良好的時(shí)頻旋轉(zhuǎn)特性設(shè)計(jì)了一種分?jǐn)?shù)階S變換，用于獲取齒輪信號的二維時(shí)頻表示；然后引入局部二值模式(LBP)算子，將LBP算子作用于分?jǐn)?shù)階S變換時(shí)頻圖，提取分?jǐn)?shù)階時(shí)頻LBP譜；最后結(jié)合“uniform”模式LBP的概念和類內(nèi)類間距準(zhǔn)則，對分?jǐn)?shù)階時(shí)頻局部二值模式譜進(jìn)行特征優(yōu)選，用于表達(dá)齒輪故障特征。對5種不同狀態(tài)的齒輪信號進(jìn)行了分析，結(jié)果表明優(yōu)選后的分?jǐn)?shù)階時(shí)頻LBP譜具有較強(qiáng)的特征描述能力，是齒輪故障信號的一類新的有效特征參數(shù)。

關(guān)鍵詞：故障診斷；齒輪；特征提??；分?jǐn)?shù)階S變換；局部二值模式

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目(E51205405,51305454)

收稿日期：2014-03-05修改稿收到日期：2014-06-06

中圖分類號:TH165.3；TH113

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2015.11.022

Abstract:A novel scheme utilizing fractional order time-frequency local binary pattern(LBP) spectra was proposed to express time-frequency characteristics of gear fault signals. Since the time-frequency resolution of S transformation was low to high-frequency signals, the fractional Fourier transformation with excellent time-frequency rotating character was integrated into S transformation and as a result, a fractional S transformation was designed. The fractional S transformation was employed to obtain 2D time-frequency representations of gear fault signals. Then, LBP operator was introduced and imposed on time-frequency images of the fractional S transformation to extract fractional order time-frequency LBP spectram. With the concept of uniform pattern of LBP and the separability evaluation criterion, the fractional order time-frequency LBP spectra were optimally selected for succinct description of gear fault features. Gear vibration signals under five different states were analyzed. The results indicated that the optimally selected fractional order time-frequency LBP spectra have an outstanding feature description capacity; they are a kind of new and effective feature parameters for gear fault signals.

Application of fractional order time-frequency local binary pattern spectra in gear fault diagnosis

ZHANGYun-qiang,ZHANGPei-lin,LIBing(Department of Vehicles and Electrical Engineering, Ordnance Engineering College, Shijiazhuang 050003, China)

Key words:fault diagnosis; gear; feature extraction; fractional order S transformation; local binary pattern (LBP)

齒輪是在機(jī)械設(shè)備中廣泛應(yīng)用的傳動(dòng)部件之一，當(dāng)發(fā)生故障時(shí)，其振動(dòng)信號屬于典型的非線性、非平穩(wěn)信號[1]。在非平穩(wěn)信號處理領(lǐng)域，時(shí)頻分析技術(shù)采用時(shí)域和頻域的二維聯(lián)合表示，可以對非平穩(wěn)信號的局部特性精確描述，具有時(shí)域和頻域等傳統(tǒng)技術(shù)無法比擬的優(yōu)勢[2-3]。

常用的時(shí)頻分析技術(shù)有小波變換[4]、S變換[5]和經(jīng)驗(yàn)?zāi)Ｊ椒纸鈁6]等。這些技術(shù)雖然在機(jī)械信號處理中具有廣泛應(yīng)用，但是都存在一些不足。在小波分析中，母小波的選擇對信號分析結(jié)果影響很大，而選擇合適的母小波并非易事，并且小波尺度與信號頻率沒有良好的對應(yīng)關(guān)系。在短時(shí)Fourier變換和小波變換的基礎(chǔ)上提出的S變換雖然克服了小波變換的一些不足，但是其高斯窗函數(shù)的標(biāo)準(zhǔn)差固定為頻率的倒數(shù)，導(dǎo)致對高頻信號的時(shí)頻分辨力較差。經(jīng)驗(yàn)?zāi)Ｊ椒纸鈱υ肼暶舾?，存在模式混疊現(xiàn)象，容易產(chǎn)生虛假分量。

分?jǐn)?shù)階Fourier變換作為Fourier變換的推廣，是一種統(tǒng)一的時(shí)頻分析方法，可以看成時(shí)頻面內(nèi)的旋轉(zhuǎn)因子，在不同階次的分?jǐn)?shù)階頻域具有不同的時(shí)頻聚集性[7-8]。為了提高信號時(shí)頻分析的靈活性和S變換的時(shí)頻分辨性能，本文基于分?jǐn)?shù)階Fourier變換良好的時(shí)頻旋轉(zhuǎn)特性，提出一種新的時(shí)頻分析技術(shù)，分?jǐn)?shù)階S變換，并應(yīng)用于齒輪故障信號處理。

分?jǐn)?shù)階S變換二維時(shí)頻圖雖能有效描述信號的時(shí)頻局部特性，但其維數(shù)巨大，不能直接作為信號的時(shí)頻特征參數(shù)，還需要進(jìn)一步提取低維特征。鑒于時(shí)頻圖本質(zhì)上是一種圖像，不同狀態(tài)齒輪信號的時(shí)頻圖會呈現(xiàn)出不同的紋理，引入常用于描述圖像紋理特征的局部二值模式[9-10]理論，在分?jǐn)?shù)階S變換對齒輪信號處理的基礎(chǔ)上，提出一種分?jǐn)?shù)階時(shí)頻局部二值模式譜，并結(jié)合“uniform”局部二值模式[11]概念和類內(nèi)類間距準(zhǔn)則[12],對分?jǐn)?shù)階時(shí)頻局部二值模式譜進(jìn)行特征優(yōu)選，用于表達(dá)齒輪故障信號特征。

1分?jǐn)?shù)階S變換

1.1S變換和分?jǐn)?shù)階Fourier變換

S變換結(jié)合了短時(shí)Fourier變換和小波變換的優(yōu)點(diǎn)，是一種較新的時(shí)頻分析方法。對于非平穩(wěn)時(shí)間信號x(t)，其S變換如下[13]

ST(τ,f)=

(1)

式中:w(t)為高斯窗函數(shù)，表達(dá)式為

(2)

由式(2)可知，S變換的高斯窗函數(shù)的標(biāo)準(zhǔn)差σ為頻率f的倒數(shù)，即σ=1/|f|，其時(shí)窗寬度隨著頻率f的增大而減小。由此可知，S變換對低頻信號具有較高的頻率分辨力，而對高頻信號具有較高的時(shí)間分辨力。

分?jǐn)?shù)階Fourier變換是Fourier變換的廣義形式，能將時(shí)域信號變換到分?jǐn)?shù)階頻域，并且得到信號新的表示。信號x(t)的分?jǐn)?shù)階Fourier變換定義為

(3)

式中:Ka(t,u)為分?jǐn)?shù)階Fourier變換核。

(4)

分?jǐn)?shù)階Fourier變換可以看成二維時(shí)頻面內(nèi)的旋轉(zhuǎn)因子，具有一定的時(shí)頻旋轉(zhuǎn)特性。a階分?jǐn)?shù)階Fourier變換可以將信號從時(shí)間軸逆時(shí)針旋轉(zhuǎn)任意角度α到分?jǐn)?shù)階頻率軸，并且獲得其在分?jǐn)?shù)階頻域中新的表示。當(dāng)a=1時(shí)，分?jǐn)?shù)階Fourier變換即為傳統(tǒng)的Fourier變換。

1.2分?jǐn)?shù)階S變換

S變換雖然具有比小波變換更好的時(shí)頻特性，但是對高頻信號的時(shí)頻分辨力不理想。為了提高S變換的時(shí)頻分辨力，本文利用分?jǐn)?shù)階Fourier變換良好的時(shí)頻旋轉(zhuǎn)特性，定義分?jǐn)?shù)階S變換如下

FrST(τ,u)=

(5)

式中:高斯窗函數(shù)w(t,u)是時(shí)間t和分?jǐn)?shù)階頻率u的函數(shù)，如式(6)所示。

(6)

式中:p為窗函數(shù)調(diào)整參數(shù)，p∈(0,1]。

由式(5)和式(6)可知，隨著a和p取值的變化，分?jǐn)?shù)階S變換能將信號變換到不同的分?jǐn)?shù)階頻域分析，具有更大的靈活性。同時(shí)，分?jǐn)?shù)階S變換的時(shí)窗寬度隨著分?jǐn)?shù)階頻率的增大而減小，所以它對分?jǐn)?shù)階低頻信號具有較高的頻率分辨力，而對分?jǐn)?shù)階高頻信號具有較高的時(shí)間分辨率，繼承了S變換良好的時(shí)頻特性。通過引入調(diào)整參數(shù)p，可以有效調(diào)節(jié)信號頻率對時(shí)窗寬度的影響，使低頻和高頻信號的分辨力達(dá)到較好的折中，從而提高信號的整體時(shí)頻分辨性能。需要注意的是，當(dāng)a=1且p=1時(shí)，分?jǐn)?shù)階頻域即傳統(tǒng)的Fourier頻域，此時(shí)分?jǐn)?shù)階S變換退化為S變換。

1.3分?jǐn)?shù)階S變換快速算法

為了快速實(shí)現(xiàn)分?jǐn)?shù)階S變換，根據(jù)卷積的定義，首先將分?jǐn)?shù)階S變換的定義式改寫為

FrST(τ,u)=[x(t)Ka(t,u)]*w(t,u)

(7)

令F(x(t))和F-1(x(t))分別表示x(t)的Fourier變換及其逆變換，根據(jù)時(shí)域卷積定理有

F[FrST(τ,u)]=F[x(t)Ka(t,u)]F[w(t,u)]

(8)

對上式兩邊同時(shí)進(jìn)行Fourier逆變換得

FrST(τ,u)=

F-1{F[x(t)Ka(t,u)]F[w(t,u)]}

(9)

由式(9)可知，信號x(t)的分?jǐn)?shù)階S變換可以通過Fourier變換及其逆變換得到。由于Fourier變換及其逆變換已有成熟的快速算法，因此可以借助它們實(shí)現(xiàn)分?jǐn)?shù)階S變換的快速計(jì)算，具體步驟不再贅述。

1.4仿真信號分析

為了驗(yàn)證分?jǐn)?shù)階S變換的有效性，構(gòu)造仿真信號x(t)如下

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

x1(t)=sin(30πt)

x2(t)=sin(30+50πt·4t)

x3(t)=sin[200πt+200πt2+20πtcos(2πt)]

(10)

由式(10)可知，x(t)由1個(gè)正弦分量和2個(gè)非線性調(diào)頻分量組成。設(shè)置采樣頻率和采樣時(shí)間分別為1 024 Hz和1 s，其時(shí)域波形見圖1。

圖1　仿真信號Fig.1 The simulated signal

對仿真信號x(t)分別進(jìn)行短時(shí)Fourier變換、Morlet連續(xù)小波變換、S變換和分?jǐn)?shù)階S變換(見圖2)。

對比圖2中子圖可以看出，由于短時(shí)Fourier變換的時(shí)窗寬度固定，其低頻和高頻的時(shí)頻分辨率都很差；連續(xù)小波變換由于尺度的大小與信號的頻率沒有較好的對應(yīng)關(guān)系，導(dǎo)致時(shí)頻圖的頻率分辨力較差，不能很好地表達(dá)信號各分量的頻率隨時(shí)間的變化情況；S變換雖然對低頻信號具有較好的時(shí)頻分辨力，但在高頻的時(shí)頻分辨力較差；分?jǐn)?shù)階S變換由于具有類似于分?jǐn)?shù)階Fourier變換的時(shí)頻旋轉(zhuǎn)特性，并引入了窗函數(shù)調(diào)整參數(shù)，時(shí)頻分辨性能明顯提高，尤其是在信號高頻部分。因此，分?jǐn)?shù)階S變換具有更好的時(shí)頻特性，本文選其獲取齒輪信號的二維時(shí)頻圖。

圖2　仿真信號二維時(shí)頻圖Fig.2 2D time-frequency images of the simulated signal

2分?jǐn)?shù)階時(shí)頻局部二值模式譜特征

2.1局部二值模式算子

局部二值模式是由Ojala等為了度量圖像局部對比度而提出的，其基本思想是根據(jù)圖像局部區(qū)域的中心像素與鄰域像素的灰度差異進(jìn)行二進(jìn)制編碼，從而對圖像紋理進(jìn)行描述[14]。

LBP算子最初定義在3×3的矩形鄰域上，編碼過程見圖3，編碼規(guī)則為：當(dāng)fi≥fc時(shí)，對應(yīng)位置編碼為1，否則編碼為零，然后按順時(shí)針方向構(gòu)成8位二進(jìn)制數(shù)，即為中心像素fc的LBP。由于定義在3×3矩形鄰域上的LBP存在不能提取大尺寸紋理特征的局限性，Ojala等又將LBP算子的鄰域從3×3的矩形擴(kuò)展到了任意半徑R和任意鄰域點(diǎn)數(shù)P的圓環(huán)形，對應(yīng)的算子記為LBPP,R。圖4為三種常見的圓環(huán)形鄰域。

圖3　LBP編碼過程Fig.3 The coding process of LBP

圖4　三種常見的圓環(huán)形鄰域Fig.4 Three circular neighborhoods for LBP

LBPP,R的計(jì)算公式如下[11]

(11)

(12)

2.2分?jǐn)?shù)階時(shí)頻局部二值模式譜

齒輪信號的時(shí)頻圖經(jīng)過局部二值模式算子處理后仍然是一個(gè)維數(shù)巨大的矩陣。為了有效描述齒輪信號時(shí)頻圖的紋理特征，在分?jǐn)?shù)階S變換時(shí)頻圖經(jīng)局部二值模式算子處理的基礎(chǔ)上，定義分?jǐn)?shù)階時(shí)頻局部二值模式譜H(h)，數(shù)學(xué)描述如下

(13)

式中，n=2P；h=0,1,…，n-1；I是一個(gè)語句判斷函數(shù)，當(dāng)LBPP,R(x,y)=h時(shí)，I{LBPP,R(x,y)=h}=1，否則I{LBPP,R(x,y)=h}=0；M、N分別為時(shí)頻圖矩陣的行數(shù)和列數(shù)。

2.3分?jǐn)?shù)階局部二值模式譜特征優(yōu)選

通過研究發(fā)現(xiàn)，LBP只有少部分模式屬于描述圖像紋理的重要模式，其出現(xiàn)概率達(dá)到90%以上，這些模式稱為“uniform”模式。如果把二進(jìn)制串連成一個(gè)圓，“uniform”模式就是串中從0到1和1到0的轉(zhuǎn)換次數(shù)不超過2次的模式。以P=8，R=1為例，由式(11)可知LBPP,R的取值可達(dá)256種之多，而“uniform”模式LBPP,R的模式數(shù)目最多只有59種。雖然“uniform”模式僅占LBP的小部分，但是它反映了絕大部分紋理信息。另外， “uniform”模式LBP譜中不同參數(shù)的可分性差異很大，因此對齒輪信號特征描述和分類的貢獻(xiàn)差異也很大。

鑒于特征優(yōu)選的目的是尋求類間散度最大而類內(nèi)散度最小，本文首先利用“uniform”模式的概念對分?jǐn)?shù)階時(shí)頻局部二值模式譜進(jìn)行初選，即將分?jǐn)?shù)階時(shí)頻局部二值模式譜中非“uniform”模式對應(yīng)的數(shù)值舍棄，僅保留“uniform”模式對應(yīng)的數(shù)值。在此基礎(chǔ)上，進(jìn)一步通過類內(nèi)類間判據(jù)對特征進(jìn)一步優(yōu)選，其定量評價(jià)指標(biāo)為

(14)

式中:Sb和Sw分別為類內(nèi)散度和類間散度。

3齒輪故障信號分析

實(shí)驗(yàn)齒輪故障信號選自某型二級傳動(dòng)齒輪箱振動(dòng)試驗(yàn)。該試驗(yàn)系統(tǒng)中的兩對齒輪副的齒數(shù)分別為25/50和18/91，輸入軸由電機(jī)帶動(dòng)，其轉(zhuǎn)速為1 491 r/min，輸出軸連接一個(gè)磁性阻尼器作為負(fù)載。試驗(yàn)中模擬了齒輪在正常、中間軸齒根裂紋和齒面磨損、輸出軸齒根裂紋和齒面磨損5種狀態(tài)。選擇的采樣頻率和采樣點(diǎn)數(shù)為6 400 Hz和1 024。圖5為5種不同狀態(tài)的齒輪故障信號。

圖5　5種不同狀態(tài)的齒輪故障信號Fig.5 Gear fault signals from five different states

3.1齒輪故障信號的分?jǐn)?shù)階時(shí)頻表示

采用分?jǐn)?shù)階S變換分別對上述5種不同狀態(tài)下的齒輪故障信號進(jìn)行處理，獲取其二維時(shí)頻圖。對于同一種狀態(tài)下的齒輪信號，設(shè)定參數(shù)a和p的取值相同，并通過多次實(shí)驗(yàn)進(jìn)行確定。最終，5種信號的時(shí)頻表示見圖6。

圖6　齒輪信號的分?jǐn)?shù)階S變換時(shí)頻圖Fig.6 Time-frequency images of gear signals from fractional S transform

由圖6可知，不同狀態(tài)的齒輪信號的時(shí)頻圖具有明顯差異，呈現(xiàn)出不同的紋理特征。因此，齒輪故障信號的分?jǐn)?shù)階S變換二維時(shí)頻圖就紋理角度而言具有可分性。

3.2分?jǐn)?shù)階局部二值模式譜特征提取

在分?jǐn)?shù)階S變換的基礎(chǔ)上，分別采用圖4中三種常見的環(huán)形鄰域提取分?jǐn)?shù)階局部二值模式譜。

圖7為不同信號的分?jǐn)?shù)階時(shí)頻LBP8,1譜、LBP8,2譜和LBP16,2譜，其中每種信號包含5個(gè)樣本，維數(shù)分別為256、256和65 536?；凇皍niform”模式的概念對圖7中分?jǐn)?shù)階局部二值模式譜特征進(jìn)行初選的結(jié)果見圖8，其維數(shù)分別降為59、59和243。

對比圖7和圖8可知，由于“uniform”模式是局部二值模式中少量且重要的模式，初選后的“uniform”模式分?jǐn)?shù)階時(shí)頻LBP譜特征維數(shù)大大降低，并且包含了原分?jǐn)?shù)階LBP譜的主要信息。

根據(jù)類內(nèi)類間判據(jù)準(zhǔn)則易知，當(dāng)式(14)中定量評價(jià)指標(biāo)較大時(shí)，對應(yīng)的特征具有較好的可分性，但是評價(jià)指標(biāo)閾值不能選擇太大，太大會造成選出的特征數(shù)量太少，不能很好地描述信號。因此，通過多次實(shí)驗(yàn)，最終選擇定量評價(jià)指標(biāo)閾值為0.9，對上述初選后的分?jǐn)?shù)階時(shí)頻LBP譜特征進(jìn)一步優(yōu)選。優(yōu)選結(jié)果見圖9，特征維數(shù)分別為13、15和97。由圖9中優(yōu)選結(jié)果可以看出，雖然“uniform”模式是局部二值模式中比較重要的模式，但是就可分性而言，它們?nèi)源嬖诿黠@差異，表現(xiàn)出不同的類內(nèi)聚合性和類間分散性。

圖7～圖9中H(h)和h的含義同式(13)，且均無量綱。

綜上可知本文特征優(yōu)選方法具有合理性和有效性。

圖7　齒輪信號的分?jǐn)?shù)階時(shí)頻局部二值模式譜Fig.7 Fractional time-frequency LBP spectrums of gear signals

圖8　“uniform”模式分?jǐn)?shù)階時(shí)頻局部二值模式譜Fig.8 Uniform pattern fractional time-frequency LBP spectrums of gear signals

圖9　優(yōu)選后的分?jǐn)?shù)階時(shí)頻局部二值模式譜Fig.9 Optimally selected fractional time-frequency LBP spectrums of gear signals

3.3齒輪故障診斷

將分?jǐn)?shù)階時(shí)頻局部二值模式譜應(yīng)用于齒輪故障診斷。從每類狀態(tài)的齒輪故障信號中分別隨機(jī)抽取40個(gè)樣本，其中20個(gè)組成訓(xùn)練樣本，其余20個(gè)組成測試樣本，選擇樸素貝葉斯分類器進(jìn)行識別。同時(shí)，引入基于S變換的時(shí)頻局部二值模式譜(時(shí)頻局部二值模式譜)作為對比實(shí)驗(yàn)。實(shí)驗(yàn)一共重復(fù)10次，平均結(jié)果見圖10。

由圖10柱狀圖可知，由于分?jǐn)?shù)階S變換比S變換更靈活、時(shí)頻分辨性能更好，表達(dá)齒輪故障信號時(shí)頻特性更加精確，基于分?jǐn)?shù)階S變換提取的特征分類效果優(yōu)于基于S變換提取的特征，并且隨著特征優(yōu)選過程的進(jìn)行，這種優(yōu)勢更加明顯。

無論基于S變換還是分?jǐn)?shù)階S變換，所提取的局部二值模式譜特征在未進(jìn)行特征優(yōu)選前的分類精度較低，說明原始局部二值模式譜中包含的大量不重要信息不但對齒輪故障診斷沒有貢獻(xiàn)，反而嚴(yán)重影響了診斷精度。通過“uniform”模式概念優(yōu)選后的特征由于舍棄了不重要模式，特征維數(shù)大幅度降低，上述不重要信息對故障診斷的影響減小，因此診斷精度有所提高。經(jīng)過兩次優(yōu)選后的局部二值模式譜特征通過類內(nèi)類間判據(jù)準(zhǔn)則篩選出的特征參數(shù)具有較大的類間散度和較小的類內(nèi)散度，使對齒輪故障信號的描述能力和區(qū)分能力增強(qiáng)，因而樸素貝葉斯分類器獲得了更好的分類精度，進(jìn)一步提高了齒輪故障診斷精度。另外需要指出的是，隨著特征優(yōu)選過程的深入，特征維數(shù)迅速減少，故優(yōu)選后的特征具有更快的診斷速度。

總之，分?jǐn)?shù)階S變換改善了S變換的時(shí)頻特性，分析齒輪故障信號的靈活性更大、精度更高。結(jié)合“uniform”局部二值模的概念和類內(nèi)類間距準(zhǔn)則優(yōu)選后的分?jǐn)?shù)階時(shí)頻局部二值模式譜具有特征維數(shù)低、描述能力強(qiáng)和鑒別性能好的優(yōu)點(diǎn)，適用于齒輪故障診斷。

圖10　齒輪故障診斷精度Fig.10 Gear fault diagnosis accuracy

4結(jié)論

針對齒輪故障信號時(shí)頻特征提取問題，本文借鑒分?jǐn)?shù)階Fourier變換的時(shí)頻旋轉(zhuǎn)特性，改進(jìn)了S變換。通過在S變換中引入分?jǐn)?shù)階Fourier變換核和增加窗函數(shù)調(diào)整參數(shù)，提出了一種分?jǐn)?shù)階S變換。仿真信號分析表明，分?jǐn)?shù)階S變換具有更大的靈活性和更好的時(shí)頻分辨力。基于分?jǐn)?shù)階S變換，定義了分?jǐn)?shù)階時(shí)頻局部二值模式譜，并結(jié)合“uniform”局部二值模式的概念和類內(nèi)類間距準(zhǔn)則對分?jǐn)?shù)階時(shí)頻局部二值模式譜特征優(yōu)選方法進(jìn)行了研究。通過齒輪故障信號分析，驗(yàn)證了分?jǐn)?shù)階時(shí)頻局部二值模式譜用于表達(dá)齒輪故障信號時(shí)頻特征的可行性和有效性。與優(yōu)選前相比，優(yōu)選后的分?jǐn)?shù)階時(shí)頻局部二值模式譜具有特征維數(shù)低、可分性好的優(yōu)點(diǎn)，可以有效提高齒輪故障診斷精度。

參考文獻(xiàn)

[1]李兵, 米雙山, 劉鵬遠(yuǎn), 等. 二維非負(fù)矩陣分解在齒輪故障診斷中的應(yīng)用[J]. 振動(dòng)、測試與診斷, 2012, 32(5): 836-840.

LI Bing, MI Shuang-shan, LIU Peng-yuan, et al.Application of two-dimensional non-negative factorization for gear fault diagnosis[J]. Journal of Vibration, Measurement & Diagnosis, 2012, 32(5): 836-840.

[2]姜鳴, 陳進(jìn), 汪慰軍. 幾種Cohen 類時(shí)頻分布的比較及應(yīng)用 [J]. 機(jī)械工程學(xué)報(bào), 2003, 39(8): 129-134.

JIANG Ming, CHEN Jin, WANG Wei-jun. Comparison and application of some time-frequency distributions belonging to Cohen class[J]. Chinese Journal of Mechanical Engineering, 2003, 39(8): 129-134.

[3]Meltzer G, Ivanov Y Y. Fault detection in gear drives with non-stationary rotational speed—part I: the time-frequency approach[J]. Mechanical Systems and Signal Processing, 2003, 17(5): 1033-1047.

[4]馬倫, 康建設(shè), 孟妍, 等. 基于Morlet小波變換的滾動(dòng)軸承早期故障特征提取研究[J]. 儀器儀表學(xué)報(bào), 2013, 34(4): 920-926.

MA Lun, KANG Jian-she, MENG Yan, et al. Research on feature extraction of rolling bearing incipient fault based on Morlet wavelet transform[J]. Chinese Journal of Scientific Instrument, 2013, 34(4): 920-926.

[5]趙淑紅, 朱光明. S變換時(shí)頻濾波去噪方法[J]. 石油地球物理勘探, 2007, 42(4): 402-408.

ZHAO Shu-hong, ZHU Guang-ming. Time-frequency filtering to denoise by S transform[J]. Oil Geophysical Prospecting, 2007, 42(4): 402-408.

[6]邵忍平, 曹精明, 李永龍. 基于EMD 小波閾值去噪和時(shí)頻分析的齒輪故障模式識別與診斷[J]. 振動(dòng)與沖擊, 2012, 31(8): 96-101.

SHAO Ren-ping, CAO Jing-ming, LI Yong-long. Gear fault pattern identification and diagnosis using time-frequency analysis and wavelet threshold de-noising based on EMD[J]. Journal of Vibration and Shock, 2012, 31(8): 96-101.

[7]陶然, 鄧兵, 王越. 分?jǐn)?shù)階傅里葉變換及其應(yīng)用[M]. 北京:清華大學(xué)出版社, 2009.

[8]Ervin S, Igor D, Ljubisa S. Fractional fourier transform as a signal processing tool: An overview of recent development[J]. Signal Processing, 2011, 91: 1351-1369.

[9]Loris N, Alessandra L, Sheryl B. Survey on LBP based texture descriptors for image classification[J]. Expert Systems with Applications, 2012, 39: 3634-3641.

[10]SONG Tie-cheng, LI Hong-liang. WaveLBP based hierarchical features for image classification[J]. Pattern Recognition Letters 2013, 34: 1323-1328.

[11]Ojala T, Pietikainen M, Maeenpaa T. Multiresolution gray-scale and rotation invariant texture classification with local binary patterns[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2002, 24(7): 971-987.

[12]周曉英, 李巍華, 丁康.基于關(guān)聯(lián)向量機(jī)的齒輪故障檢測方法研究[J]. 振動(dòng)與沖擊, 2008, 27(6): 51-54.

ZHOU Xiao-ying, LI Wei-hua, DING Kang. Gear fault detection based on relevance vector machine[J]. Journal of Vibration and Shock, 2008, 27(6): 51-54.

[13]Djurovi I, Sejdi E, Jiang J. Frequency-based window width optimization for S-transform[J]. AEU-International Journal of Electronics and Communications, 2008, 62(4): 245-250.

[14]白航, 趙擁軍, 胡德秀. 時(shí)頻圖像局部二值模式特征在雷達(dá)信號分類識別中的應(yīng)用[J]. 宇航學(xué)報(bào), 2013, 34(1): 139-146.

BAI Hang, ZHAO Yong-jun, HU De-xiu. Radar signal recognition based on the local binary pattern feature of time-frequency image[J]. Journal of Astronautics, 2013, 34(1): 139-146.