軸承matlab處理程序

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

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

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

4、消除趨勢(shì)項(xiàng)與零均值化處理的功能相似。不過(guò)，需要注意的是，它更重要的消除趨勢(shì)項(xiàng)，因?yàn)楸緮?shù)據(jù)中的多項(xiàng)式趨勢(shì)項(xiàng)很小，所以沒(méi)有明顯的變化。1.3平滑處理（原理公式見(jiàn)報(bào)告P11）使用五點(diǎn)三次平滑，命令窗口輸入：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(