[發明專利]一種基于內秉時間尺度分解和譜峭度的包絡分析方法

有效

• 申請號: 201610492062.8
• 申請日: 2016-06-29
• 公開（公告）號: CN105954030B
• 公開（公告）日: 2018-03-23
• 發明（設計）人: 林近山;竇春紅
• 申請（專利權）人: 濰坊學院
• 主分類號: G01M13/02
• 分類號: G01M13/02;G01M13/04
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種基于內秉時間尺度分解和譜峭度的包絡分析方法，其特征在于，包括以下步驟：

步驟1：利用加速度傳感器以采樣頻率fs測取旋轉機械的振動信號x（k）, k=1, 2, …,N，N為采樣信號的長度；

步驟2：采用內秉時間尺度分解算法將信號x（k）分解成n個分量和一個趨勢項之和，即，其中，c_i（k）代表由內秉時間尺度分解算法得到的第i個分量，r_n（k）代表由內秉時間尺度分解算法得到的趨勢項；

步驟3：對c_i（k）執行重排操作和替代操作，經重排操作得到的數據用c_i^shuffle（k）表示，替代操作后得到數據用c_i^FTran（k）表示；

步驟4：對c_i（k）、c_i^shuffle（k）和c_i^FTran（k）分別執行多重分形去趨勢波動分析（Multifractal Detrended Fluctuation Analysis, MFDFA），得到廣義Hurst指數曲線，c_i（k）的廣義Hurst指數曲線用H_i（q）表示；c_i^shuffle（k）的廣義Hurst指數曲線用H_i^shuffle（q）表示；c_i^FTran（k）的廣義Hurst指數曲線用H_i^FTran（q）表示；

步驟5：如果H_i（q） 與H_i^shuffle（q）或H_i（q） 與H_i^FTran（q）之間的相對誤差小于5%，或者H_i（q） 、H_i^shuffle（q） 和H_i^FTran（q）三者都不隨q而變化，則拋棄對應的c_i（k）分量；

步驟6：對剩余的c_i（k）分量求和，將該和記為信號經重排和替代濾波后的結果x_f1（k）；

步驟7：對x_f1（k）執行譜峭度分析，求出信號峭度最大處所對應的中心頻率f₀和帶寬B；

步驟8： 根據中心頻率f₀和帶寬B對x_f1（k）進行帶通濾波，得到x_f2（k）；

步驟9：對信號x_f2（k）執行三次樣條迭代平滑包絡分析，得到信號包絡eov（k）；

步驟10：對得到的信號包絡eov（k）執行離散傅里葉變換得到包絡譜，根據包絡譜特征頻率判斷機器的故障類型。