軸承matlab處理程序

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

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

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

4、消除趨勢項與零均值化處理的功能相似。不過，需要注意的是，它更重要的消除趨勢項，因為本數據中的多項式趨勢項很小，所以沒有明顯的變化。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(