軸承matlab處理程序

1、1數(shù)據(jù)導入matlab1.1啟動Matlab軟件1.2點擊載入故障數(shù)據(jù)中的G2015，Wkspace窗口出現(xiàn)：1.3取第一組數(shù)據(jù)G201，命令窗口輸入：G201=G2015(1:1:20000)2.數(shù)據(jù)預處理在測試中由數(shù)據(jù)采集所得的原始信號，在分析前需要進行預處理，以提高數(shù)據(jù)的可靠性和真實性，并檢查信號的隨機性，以便正確地選擇分析處理方法。預處理工作主要包括三個方面：一是除去信號中的外界干擾信號和剔除異常數(shù)據(jù)，如趨勢項和異點；二是對原始

2、數(shù)據(jù)進行適當?shù)钠交驍M合；三是對原始信號的特性進行檢驗。當然這些處理工作不是全部必需的，可以選—項或兩項內容，當認為原始信號獲取工作十分可靠或原始數(shù)據(jù)簡單可以直接判斷的情況下，也可以不進行這些預處理工作。以下所做數(shù)據(jù)預處理，故障軸承以G201為例，正常軸承以Z201為例，觀察原始數(shù)據(jù)經過不同方法做處理前后的變化。1.1零均值化處理（原理公式見報告P8）命令窗口輸入：G201l=G201sum(G201)20000%G201l為零均值處理

3、后的數(shù)據(jù)?！?0000”為采樣點數(shù)。sum為求和語句subplot(211)plot(G201)subplot(212)plot(G201l)%顯示G201與G201l得到下面圖形：從時域圖形上看，是波形整體在Y軸的平移。再看看頻域變化，命令窗口輸入：N=20000%采樣點數(shù)fs=10000%采樣頻率與前面零均值化處理中做頻域圖的方法一樣，做出G201與G201x的頻譜圖G201p與G201xp，得到圖形如下：從時域圖形和頻域圖形上看，

4、消除趨勢項與零均值化處理的功能相似。不過，需要注意的是，它更重要的消除趨勢項，因為本數(shù)據(jù)中的多項式趨勢項很小，所以沒有明顯的變化。1.3平滑處理（原理公式見報告P11）使用五點三次平滑，命令窗口輸入：a=G201fk=1:2b(1)=(69a(1)4(a(2)a(4))6a(3)a(5))70b(2)=(2(a(1)a(5))27a(2)12a(3)8a(4))35fj=3:N2b(j)=(3(a(j2)a(j2))12(a(j1)a(