[發(fā)明專利]基于自適應(yīng)STFT的旋轉(zhuǎn)機(jī)械階次跟蹤方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201310256814.7
• 申請(qǐng)日: 2013-06-25
• 公開（公告）號(hào): CN103353344A
• 公開（公告）日: 2013-10-16
• 發(fā)明（設(shè)計(jì)）人: 陶學(xué)軍;初永剛;田勇;李佳;桓保軍;張舉良;趙洋;蘇靜
• 申請(qǐng)（專利權(quán)）人: 國家電網(wǎng)公司;許繼集團(tuán)有限公司;許昌許繼風(fēng)電科技有限公司
• 主分類號(hào): G01H17/00
• 分類號(hào): G01H17/00;G01M99/00
• 代理公司: 鄭州睿信知識(shí)產(chǎn)權(quán)代理有限公司 41119
• 代理人: 胡泳棋

1.基于自適應(yīng)STFT的旋轉(zhuǎn)機(jī)械階次跟蹤方法，其特征在于，該方法的步驟如下：

（1）采集旋轉(zhuǎn)機(jī)械機(jī)組的振動(dòng)信號(hào)x(n)，x(n),n=0,1,…,N-1；

（2）根據(jù)機(jī)組轉(zhuǎn)速的變化趨勢(shì)，依據(jù)轉(zhuǎn)速越小窗函數(shù)長度越大來確定窗函數(shù)長度的變化公式b(n)；

（3）用自適應(yīng)STFT計(jì)算振動(dòng)信號(hào)的時(shí)頻譜，得到變換后的時(shí)頻譜二維序列X[i,k],i,k=0,1,2,…,N-1；

（4）搜索時(shí)頻譜二維序列X[i,k]每一列中的極大值并記錄其行號(hào)，將要提取的階次分量之外的元素置為0，得到需要提取的階次分量的時(shí)頻域成分X'[i,k]；

（5）將提取的階次分量的時(shí)頻域成分帶入下式得提取的階次分量，實(shí)現(xiàn)階次跟蹤，

x(m)=1DNΣi=-∞+∞Σk=0N-1X′(i,k)b(k)-(m-k)Δb(k)(b(k))32gi(m-i)exp(j2πmNk),]]>

m=0,1,2,…,N-1。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟（2）中窗函數(shù)長度的變化公式為時(shí)間的線性函數(shù)。

3.根據(jù)權(quán)利要求2所述的方法，其特征在于，所述步驟（2）中窗函數(shù)長度的變化公式為：b(n)=M×(1-a×n/N),其中M為滿足最小頻率分辨率的窗長度，0<a<1是窗函數(shù)長度變化率，N為數(shù)據(jù)長度。