軸承matlab處理程序

1．?dāng)?shù)據(jù)導(dǎo)入matlab啟動(dòng)Matlab軟件1.2點(diǎn)擊岀載入故障數(shù)據(jù)中的G2015,Workspace窗口出現(xiàn)：N：=uTie Value Class5G201S <100000x1double>double1.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ì)原始數(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)過不同方法做處理前后的變化。1.1零均值化處理(原理公式見報(bào)告P8)命令窗口輸入：G201I二G201-sum(G201)/20000;%G201l為零均值處理后的數(shù)據(jù)?！?0000”為采樣點(diǎn)數(shù)。采樣點(diǎn)數(shù) xw4從時(shí)域圖形上看，是波形整體在Y軸的平移。再看看頻域變化，命令窗口輸入:N=20000; %采樣點(diǎn)數(shù)fs=10000;%采樣頻率

f=(0:N-1)'*fs/N;%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換G201p二abs(fft(G201));%進(jìn)行fft變換,G201p為G201進(jìn)行fft變換后結(jié)果G201lp二abs(fft(G201l));%進(jìn)行fft變換，G201lp為G2011進(jìn)行fft變換后結(jié)果subplot(2,1,1),plot(f(1:N/2),G201p(1:N/2));subplot(2,1,2),plot(f(1:N/2),G201lp(1:N/2));%顯示G201與G201p的頻譜圖得到下面圖形：從頻域圖可以明顯看出，零均值后消除?二0處出現(xiàn)一個(gè)由直流分量產(chǎn)生的大譜峰(將近達(dá)到4.5x104)，處理后避免了其對(duì)周圍小峰值產(chǎn)生的負(fù)面影響，便于頻域分析。消除趨勢(shì)項(xiàng)(原理公式見報(bào)告P10)使用最小二乘法，命令窗口輸入：t=(0:1/fs:(N-1)/fs)';%離散時(shí)間列向量G201x二polyfit(t,G201,6); %計(jì)算多項(xiàng)式待定系數(shù)向量G201x二G201-polyval(G201x,t);%用G201減去多項(xiàng)式系數(shù)生成的趨勢(shì)項(xiàng),G201x即為消除趨勢(shì)項(xiàng)后的數(shù)據(jù)subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201x);%顯示G201與G201x得到以下圖形：G2011 1 1 1 1 r--□ 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2采樣點(diǎn)數(shù) x1Q4G201xI I I I I_3 I | | | | | | | | '"0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2采樣點(diǎn)數(shù) x1Q4與前面零均值化處理中做頻域圖的方法一樣，做出G201與G201x的頻譜圖G201p與G201xp，得到圖形如下：從時(shí)域圖形和頻域圖形上看，消除趨勢(shì)項(xiàng)與零均值化處理的功能相似。不過，需要注意的是，它更重要的消除趨勢(shì)項(xiàng)，因?yàn)楸緮?shù)據(jù)中的多項(xiàng)式趨勢(shì)項(xiàng)很小，所以沒有明顯的變化平滑處理(原理公式見報(bào)告P11)使用五點(diǎn)三次平滑，命令窗口輸入：a=G201';fork=1:2b(1)=(69*a(1)+4*(a(2)+a(4))-6*a(3)-a(5))/70;b(2)=(2*(a(1)+a(5))+27*a(2)+12*a(3)-8*a(4))/35;forj=3:N-2b(j)=(-3*(a(j-2)+a(j+2))+12*(a(j-1)+a(j+1))+17*a(j))/35;end

b(N-1)=(2*(a(N)+a(N-4))+27*a(N-1)+12*a(N-2)-8*a(N-3))/35;b(N)=(69*a(N)+4*(a(N-1)+a(N-3))-6*a(N-2)-a(N-4))/70;a=b;endG201ph=a'; %G201ph為五點(diǎn)三次平滑法處理的數(shù)據(jù)subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201ph);%顯示G201與G201ph得到以下圖形：與前面零均值化處理中做頻域圖的方法一樣，做出G201與G201ph的頻譜圖G201p與G201php，得到圖形如下：從時(shí)域圖形上看，平滑處理使圖形變得平滑，去除毛刺，從頻域圖形上看，高頻部分明顯變少變小，而低頻部分基本無變化。因?yàn)楣收系念l率主要集中在低中頻部分，這樣處理后不僅對(duì)故障的分析無影響，而且去除部分噪音，減少干擾。濾波處理(原理公式見報(bào)告P13)%使用巴特沃斯濾波器進(jìn)行濾波，命令窗口輸入：wp=2400;%通帶截至頻率2400hzws=2800;%阻帶截至頻率2800hzrp=2;%通帶波動(dòng)系數(shù)rs=60;%阻帶波動(dòng)系數(shù)[N,wn]二buttord(wp/(fs/2),ws/(fs/2),rp,rs,'z');%建立巴特沃斯濾波器[num,den]二butter(N,wn);%建立數(shù)字濾波器[H,W]二freqz(num,den);%分析濾波器的幅頻特性plot(W*fs/(2*pi),abs(H));grid;%巴特沃斯濾波器頻率響應(yīng)圖得到巴特沃斯濾波器頻率響應(yīng)圖：繼續(xù)輸入：G201lb=filtfilt(num,den,G201);%G201lb為G201濾波后的數(shù)據(jù)subplot(2,1,1),plot(G201);subplot(2,1,2),plot(G201lb);%顯示G201與G201lb得到以下圖形：G201與前面零均值化處理中做頻域圖的方法一樣，做出G201與G201lb的頻譜圖G201p與G201lbp，得到圖形如下:在時(shí)域內(nèi)不能明顯的看出處理前后的區(qū)別。但從頻域圖可以看出，2500Hz后的頻率幾乎不存在。因?yàn)榈屯V波器的通帶截至頻率為2400hz，阻帶截至頻率為2800hz?？梢姙V波效果是很好的。以上介紹了一些數(shù)據(jù)預(yù)處理的方法，鑒于本文采集的原始信號(hào)數(shù)據(jù)較好，故只做零均值化這一項(xiàng)處理。3．時(shí)域特征值提?。ㄔ砉揭奝15）命令窗口輸入：

G201m二sum(G201l)/20000; %G201m為均值，G201I為零均值化處理后結(jié)果，下同G201仁sum((G201l-G201m)?八2); %G201f為方差G201rms=sqrt(sum(G201I?八2)/20000); %G201rms均方根值G201peak=(max(G201I)—min(G201I))/2; %G201peak為峰值G201c=G201peak/G201rms; %G201c為峰值因子G201k二sum(G201I?八4)/((G201rms.八4)*20000); %G201k為峭度系數(shù)G201s=(G201rms*20000)/sum(abs(G201l)); %G201s為波形因子G201cI=G201peak/(sum(sqrt(abs(G201I)))/20000).八2; %G201cI裕度因子G201i=(G201peak*20000)/sum(abs(G201I)); %G201i脈沖因子由此得到G201的時(shí)域特征值根據(jù)前述方法一次得到G202~G2010,Z201~Z2010的時(shí)域特征值，建立表格狀態(tài)樣本時(shí)域特征值均值（107）方差均方根值RMS峰值peak峭度系數(shù)K峰值因子C裕度因子CL脈沖因子1波形因子SG2016.85222340.800.34212.270113.3236.635718.22512.4641.8785499G20222.34522605.740.36102.354914.2176.524218.82112.6271.9355027G203—32.38542902.360.38092.488613.5326.532618.93212.6361.9343030G20415.32902630.680.36272.480313.8756.838819.02812.9641.8957684故G20515.89442510.690.35432.337913.2636.598518.57412.5271.8985障201軸G206—3.73632647.010.36382.393613.5386.579318.24012.4291.8892承285G207—2.06752379.660.34492.287112.3886.630517.56212.1191.8278230G208—4.51022548.620.35702.451214.0456.866619.81913.2831.9346759G2098.08802496.800.35332.337112.6306.614617.37512.0821.8266473G20104.40802871.860.37892.316711.7416.113816.91311.3861.8625

746Z2015.24191940.060.31151.58504.32035.08898.18286.73971.3244Z20227.71791805.180.30041.50304.46845.00278.06126.63511.3263Z203-23.98241698.730.29141.37644.62554.72287.70336.31551.3372Z2042.58211677.680.28961.73994.88596.00739.68087.97701.3279Z2053.42741890.520.30751.52314.50354.95397.94036.55221.3226Z20628.42331688.290.29051.32473.92824.55947.21485.96471.3082Z20716.67021629.540.28541.46184.58805.12138.23386.79211.3262Z208-17.69651605.220.28331.34904.43664.76187.67706.32011.3272Z20920.48481714.370.29281.45734.66104.97748.09836.64661.3354Z2010-4.03201790.060.29921.68775.19085.64129.37957.64671.3555正常軸承列出時(shí)域參數(shù)的數(shù)字表后可以簡(jiǎn)單分析，故障軸承和正常軸承在方差，峰值，峭度系數(shù)，裕度因子，脈沖因子，波形因子差別較為明顯，而在均值，均方根值，峰值因子差別不明顯。4．頻域特征值提?。ㄔ砉揭奝18）頻域參數(shù)命令窗口輸入：fori=2:20000G201g(i)=(G201l(i)-G201l(i-1))/(1/10000);endfori=2:20000G201gg(i)=G201g(i)*G201l(i);endG201ms仁(sum((G201g)?八2))/(4*(pL2)*sum(G201l?八2)); %G201msf為均方頻率G201fc二(sum(G201gg))/(2*pi*sum(G201l.八2)); %G201fc重心頻率G201vf二G201msf-G201fc.八2; %G201vf為頻率方差由此得到G201I的頻域參數(shù)。根據(jù)前述方法一次得到G202~G2010,Z201~Z2010的時(shí)域特征值，建立表格狀態(tài)樣本頻域參數(shù)重心頻率頻率方差均方頻率故G201727.4390768683.09001297850.5751障G202732.7359755254.24021292156.1458軸G203783.7554830793.71491445066.2883承G204772.6559858228.74341455225.8325G205708.5887743009.38451245107.3676

G206807.02478766520240G207776.64488252423696G208752.5714842258.30101408621.9914G209839.0520952037.65141656045.8589G2010774.6435842014.28501442086.8989Z2011947.72046177234.38169970849.1539Z2022063.44766752047.589611009863.5570Z2032124.78896969666.443511484394.4028Z2042129.98197038297.661611575120.4785Z2052049.52726681185.600010881747.3107Z2062047.97916659325.170810853543.6293Z2072151.69967045509.332011675320.4950Z2082181.39907213367.433511971868.8764Z2092171.50927046245.891511761698.1452Z20102088.47246637350.777710999067.8211正常軸承從上表可以看出，頻域參數(shù)的特征值重復(fù)性和差異性都是比較良好的。4.2傅里葉變換（原理公式見P20）將G201I和Z201I(Z201I為Z201零均值化后數(shù)據(jù))的fft變換后的G201lp與Z201lp做出，程序如下：G201lp=abs(fft(G201l,16384));G201lp=G201lp(1:8192,1);Z201lp=abs(fft(Z201l,16384));Z201lp=Z201lp(1:8192,1);subplot(2,1,1),plot(G201lp);subplot(2,1,2),plot(Z201lp);如下圖所示：

能夠區(qū)分兩個(gè)狀態(tài)且能代表自己頻譜的區(qū)域有：點(diǎn)(326,1).區(qū)域(2560?3000)、點(diǎn)(3278,1)、區(qū)域(6310~6646)、區(qū)域(6850~7300)用屹標(biāo)記。對(duì)故障軸承數(shù)據(jù)隨機(jī)抽取G202fft、G206fft、G207fft、G209fft數(shù)據(jù)對(duì)比圖形如下:從故障軸承抽樣數(shù)據(jù)對(duì)比圖形可以看出,各個(gè)特征值的性重復(fù)較好。-…比T X：326 -止?uJiufc對(duì)故障軸承數(shù)據(jù)隨機(jī)抽取G202fft、G206fft、G207fft、G209fft數(shù)據(jù)對(duì)比圖形如下:從故障軸承抽樣數(shù)據(jù)對(duì)比圖形可以看出,各個(gè)特征值的性重復(fù)較好。-…比T X：326 -止?uJiufcX：3278￥：179.8 …5一亠—1…X:326 : :1136.13 :It-1對(duì)正常軸承數(shù)據(jù)隨機(jī)抽取Z203fft、Z204fft、Z206fft、Z208fft數(shù)據(jù)對(duì)比圖形如下:X：326Y：187.6500唾圖3-12正常軸承重復(fù)性FFT譜X：6461Y:252.8.一i3■■UX:6549Y:242.9X：3278Y：6.227—從正常軸承抽樣數(shù)據(jù)對(duì)比圖形可以看出,各個(gè)特征值的性重復(fù)較好。利用以下程序?qū)?/p>

G201~G2010,Z201~Z2010的傅里葉變換特征值與特征區(qū)域提取出來：G201ffttz=[G201lp(326,1),sum(G201lp(2560:3000,1)),G201lp(3278,1),sum(G201lp(6310:6646,1)),sum(G201lp(6850:7300,1))];%G201ffttz為G2011進(jìn)行fft變換后提取的特征值建立表格：狀態(tài)樣本FFT頻域特征值(326,1)(2560：3000,1)(3278,1)(6310：6646,1)(6850：7300,1)故G20124.037016457.4591223.41115511.60235070.6419

G20236.131815337.1272216.35055283.45645046.5340G20328.267118178.6398213.62015642.92055046.5340G20411.124618178.6398163.98695662.49845687.7720G205114.039916154.9008168.20395562.42985151.0799G20636.935119195.7369179.82815528.26136344.3349G20711.462918551.4539162.56425491.62415443.5653G208141.138515024.3193148.25905700.04945710.0924G20937.458620519.3069148.70935955.78156390.3742G201092.740220333.6468138.53706092.24495809.4190Z201247.171610356.38888.211015539.769111057.7255Z202179.159810083.433627.494314772.372010632.8833Z203187.585210018.71816.226915396.491310692.6573Z204212.92359502.587717.208415035.716910475.6605Z205132.327710125.071418.669915691.313911033.8520Z206210.835610029.603417.639215133.589610612.5788Z207205.152910220.055918.711914764.755710577.8874Z208227.172010414.902020.944014186.354211097.2232Z209173.456810161.368810.417015232.020110877.8867Z2010198.617310366.74376.862415144.132911205.6632障軸承正常軸承4.3功率譜處理((原理公式見P22)采用Welch平均周期法，采樣頻率為10000Hz,長(zhǎng)度為16384點(diǎn)，分段時(shí)每段長(zhǎng)度為4096,相鄰兩段重疊的點(diǎn)數(shù)為2048，因此分成了7段，窗函數(shù)為缺省。命令窗口輸入以下程序：[p,f]=spectrum(G201l,4096,2048,[],fs);G201gl二p(:,1);%采用Welch平均周期法，G201功率譜處理結(jié)果將G201gl和Z201gl顯示出來，得到：G201glZ201gl能夠區(qū)分兩個(gè)狀態(tài)且能代表自己頻譜的區(qū)域有：點(diǎn)（82,1）、區(qū)域（660?739）、點(diǎn)（820,1）、點(diǎn)（1473,1）、點(diǎn)（1616,1）、點(diǎn)（1639，1）點(diǎn)（1778，1）。對(duì)故障軸承數(shù)據(jù)隨機(jī)抽取G202gl、G204gl、G206gl、G208gl數(shù)據(jù)對(duì)比圖形如下:111J.Qgf一一…iL. i........i111J.Qgf一一…iL. i........iX：820:Y:0.5068■—X:1473X:1616X:1639X:1778Y:0.1865Y:0.02457Y:0.04409Y:0.01281140016001800X:1473X:1616X:1639X:1778Y:0.07676Y:0.01558Y:0.06597Y:0.0080171400 "16001800迺j圖3-14故障軸承重復(fù)性功率譜從故障軸承抽樣數(shù)據(jù)對(duì)比圖形可以看出,各個(gè)特征值的性重復(fù)較好。對(duì)正常軸承數(shù)據(jù)隨機(jī)抽取Z203gl、Z204gl、Z209gl、Z2010gl數(shù)據(jù)對(duì)比圖形如下:一…LX：1616Y:1.276X:1639X：82X：820Y:0.00682Y：1443 x：1778]Y：0.4774..■口 uv...L.■■Ajl,,z-. ■.f:0.2682圖3-15正常軸承重復(fù)性功率譜從正常軸承抽樣數(shù)據(jù)對(duì)比圖形可以看出，各個(gè)特征值的性重復(fù)較好。將G201~G2010Z201~Z2010的傅里葉變換特征值與特征區(qū)域提取出來，建立表格：狀態(tài)樣本W(wǎng)elch平均周期法功率譜特征值（82，1）（820，1）（1473,1）（1616,1）（1639，1）（1778，1）（660：739，1）G2010.36430.56650.21450.01790.04310.00878.8665G2020.23510.50680.18650.02460.04410.012810.8937G2030.26840.52910.12290.01150.03530.011012.0233故G2040.27000.34990.07680.01560.06600.008010.2235障G2050.42580.36920.13480.01790.04520.00919.7017軸G2060.22950.41230.11890.01260.08690.014714.1895承G2070.15160.32130.02270.01630.04390.006511.9705G2080.29760.28780.04740.01300.05170.00747.6765G2090.21950.32900.10010.02170.03190.009415.106G20100.36300.28870.15710.02510.04680.011413.565Z2010.65380.00910.02171.22961.31050.24273.8739Z2020.36680.01870.01171.43032.31200.31313.5024Z2030.47740.00680.00971.27641.44310.26823.4486正Z2040.56580.00970.01601.44681.97070.20482.9251常Z2050.46120.01110.00950.42011.42860.22223.3279軸Z2060.59870.00920.02771.23900.94610.2747承Z2070.43350.01320.01340.89882.08890.1968Z2080.55980.01360.01630.54911.89510.1672Z2090.41640.00950.02321.02071.47640.1930Z20100.39260.01640.01631.28561.42980.18733.29203.39163.57903.53333.8456從上表可以看出，故障軸承和正常軸承功率譜的特征值重復(fù)性和差異性都是比較良好的。5．時(shí)頻分析法——小波包絡(luò)(原理公式見P24)(0^17(25J00V5000)Hz+J(OT5)H(0^500)Hz+JHz*1(0^17(25J00V5000)Hz+J(OT5)H(0^500)Hz+JHz*1(17E-ri"2500)H5000)Hz+J在這里選擇(30)頻段的信號(hào)為例進(jìn)行重構(gòu)，再包絡(luò)解調(diào)，程序如下：wpt二wpdec(G201l,3,'db4');%小波包進(jìn)行3層分解，分解用的小波基為db4cz二wprcoef(wpt,[30]);%對(duì)(30)節(jié)點(diǎn)即(0~875Hz)進(jìn)行重構(gòu)G201xb二abs(hilbert(cz));%G201xb為G201I的包絡(luò)分析G201xbp=abs(fft(G201xb,4096));%G201xbp為G201xb進(jìn)行FFT變換后數(shù)據(jù)故障軸承G202與正常軸承Z202的(30)節(jié)點(diǎn)的包絡(luò)譜：

圖3-16節(jié)點(diǎn)（30）G202與Z202小波包包絡(luò)譜對(duì)比發(fā)現(xiàn)兩者并沒有明顯的區(qū)別，很難從此節(jié)點(diǎn)提取出特征值，說明故障頻率不在）?875Hz頻率段。再對(duì)（31）節(jié)點(diǎn)的樣本1數(shù)據(jù)進(jìn)行分析：圖3-17節(jié)點(diǎn)（31）G202與Z202小波包包絡(luò)譜從圖中可以得出在點(diǎn)（2，1）、點(diǎn)（513,1）和點(diǎn)（1025，1）處的值個(gè)狀態(tài)的差異性比較好，因此初步確定這三點(diǎn)為其特征值。接下來對(duì)（32）、（33）幾個(gè)節(jié)點(diǎn)進(jìn)行了包絡(luò)譜分析：圖3-18節(jié)點(diǎn)(32)G202與Z202小波包包絡(luò)譜圖3-18節(jié)點(diǎn)(33)G202與Z202小波包包絡(luò)譜在眾多的包絡(luò)譜分析中發(fā)現(xiàn)其規(guī)律是：他們都是在包絡(luò)譜的(2,1)、(513,1)、(1025,1)

點(diǎn)出現(xiàn)峰值并且各個(gè)狀態(tài)有區(qū)別，即差別性較好，并且通過重復(fù)性檢驗(yàn)，其可以作為特征值。因此將這三點(diǎn)作為特征值，并取(31)、(32)、(33)三個(gè)節(jié)點(diǎn)作為此種方法分析的特征值提取節(jié)點(diǎn)。特征值提取的程序如下：wpt=wpdec(G201l,3,'db4');fork=1:3cz=wprcoef(wpt,[3k]);hom=abs(hilbert(cz));czf=abs(fft(hom,4096));

G201xbp(1,3*(k)-2)=czf(2,1);G201xbp(1,3*(k)-1)=czf(513,1);G201xbp(1,3*(k))=czf(1025,1);end得到的G201xbp為4096*9的矩陣，第一行就是所求的各節(jié)點(diǎn)3個(gè)點(diǎn)的特征值。12345 |T891140.82375.985510.479257.291519.4564U.915652.201524.80815.8094依此建立如下表格：狀態(tài)樣本小波包包絡(luò)譜特征值（31）（32）（33）（2，1）(513,1)(1025,1)（2，1）(513,1)(1025,1)（2，1）(513,1)(1025,1)G201140.8237146.3217168.2915.985510.479257.291519.45640.915652.201524.80815.8094G2026.16529.797766.867420.45920.796290.958228.26856.2661G2039.481512.902764.289221.36120.6669103.759931.21457.3365G204187.2597.203512.036855.834720.15920.754997.767229.56787.1333故1G205120.1574.97749.117938.364721.64820.663064.718724.17255.8821障7軸G206126.1719.998613.977838.176219.99110.629763.416833.85788.0195承7G207172.1008.891213.615048.372116.89800.834875.363429.52816.89480G208129.8303116.4576.317313.272052.231019.16160.765483.543128.15307.9258G2094.113611.304852.985220.01131.237453.141328.90417.58049G2010135.1579.374115.721574.642928.54920.9465107.163344.01199.88002

Z20146.381212.393721.06194.710110.46040.295412.422318.86958.1549Z20223.954111.761418.22852.190710.23340.379916.762117.05978.3838Z20313.48868.445916.57215.43278.06850.345113.014517.45928.2029正Z20493.863514.109119.01434.954610.50310.452517.713519.01018.4352常Z20539.649412.643519.611511.75429.97510.350420.583019.69528.4686軸Z20638.831911.601818.92456.05379.57180.346915.048718.10158.2054承Z20727.221112.947618.80511.98989.88270.462919.416816.95098.4433Z20813.718510.653916.82794.08689.21350.36089.594918.36949.2244Z20963.891213.132218.80909.11319.85580.39784.604518.66667.9271Z201019.222713.034419.74573.828710.06120.353624.758616.12317.91866．特征值歸一化(原理公式見報(bào)告P29)程序如下：fori=1:33forj=1:20gy(i,j)=(tz(i,j)-min(tz(i,:)))/(max(tz(i,:))-min(tz(i,:)));%tz為原特征值矩陣gy為歸一化后的特征值矩陣。33為原始特征值個(gè)數(shù)，20為數(shù)據(jù)樣本個(gè)數(shù)。endend得到下列表格：故障軸承G201~G20100.64530.90000.78470.7940.47110.49860.45840.6650.60500610.56710.771310.79060.69810.80310.59700.72730.68730.97650.60250.796110.81350.72750.82480.63110.75510.71720.97950.81230.88510.99290.8700.91840.82690.96790.8690.85215830.913110.93340.96680.9070.93400.82230.98340.8450.7593840.89990.8510.85520.98800.8830.87550.897710.8900.67368870.87360.9200.92960.93730.9010.87480.820910.8060.76982140.88810.9100.91150.95630.8960.88330.840810.8350.740

36880.909110.99810.93660.94100.92620.82830.99860.82640.8836頻域0.01280.0160.05100.043500.06680.04620.02990.0880.044參468數(shù)0.00400.0010.01360.017800.02070.01270.01530.0320.0159330.00490.0040.01860.019600.02640.01710.01520.0380.018434FFT譜0.05470.1050.072600.4360.10930.00140.55080.1110.345特9068征0.63130.5290.78750.78750.6030.87980.82140.501210.983值68110.9670.95490.72640.7450.79930.71980.65400.6560.60958020.021900.03450.03640.0260.02350.02000.04000.0640.0778670.0039000.10410.0170.21070.06450.10770.2180.1230290.42350.1660.23260.23580.5460.155100.29070.1350.420功3029率10.8930.93320.61300.6470.72450.56190.50210.5750.503譜3577特10.8630.55320.32830.6110.53370.06440.18490.4420.72征420值0.00450.00900.00290.0040.00080.00330.00100.0070.0091510.00490.0050.00150.01500.0050.02410.00530.008700.0064850.00720.0200.01470.00490.0080.026700.00290.0090.01655500.48780.6540.74690.59920.5560.92480.74260.390110.8732350.73280.7640.890810.6130.64850.91280.66950.5920.700小4962波0.18730.2050.53700.30910.0860.58880.47800.220500.526包342絡(luò)0.11400.0560.31690.244400.40690.37650.34780.1830.552譜919特0.76120.8930.85750.74110.5000.49810.63840.69150.7011征079值0.55600.6050.64900.59030.6630.58210.43110.54160.58310010.65840.5310.39440.48780.3900.35490.57270.499010.6916220.46410.8420.96680.90840.5860.57340.68990.76970.47310130.31140.4350.54110.48210.2880.63590.48070.43140.458156300.1120.37500.32520.0170.54290.26660.51990.4351291正常軸承Z201~Z20100.6180.98840.1380.5750.58810.8060.2410.86950.46638209760.2580.15420.0720.0550.2190.06400.01800.08410.1425119970.2880.17520.0830.0640.2480.07380.02100.09730.162990505值0.2230.15320.0440.3560.17000.1170.0200.11390.31196475890.0380.05250.0670.0930.05500.0640.0490.07120.1227181940.2290.19210.0700.6270.17100.2430.0870.18120.46895860570.0760.06710.0380.1950.05700.0800.0360.070101050.09160.0470.2740.08000.1130.0480.09320.22989993060.0250.02890.0460.0310.02300.0280.0300.04340.0754824073頻域0.8410.91990.9610.9650.9100.90940.97910.99330.9369參36158數(shù)0.8390.92870.9620.9720.9170.91440.97410.97420.9110939810.8130.91030.9540.9630.8980.89570.97210.98040.909356044FFT譜10.71190.7470.8540.5130.84610.8220.9150.68770.7943特69503征0.0770.05270.04600.0560.04780.0650.0820.05980.0784值585180.0090.097900.0500.0570.05250.0570.0670.01930.0029163580.9850.91170.9710.93710.94640.9110.8550.95570.947447004

0.97600.90700.91670.88150.97210.90370.89810.98240.94681功10.42850.6480.8240.6160.89030.5610.8120.52730.4799率78538譜0.0040.021300.0050.0070.00430.0110.0120.00480.0172特12741征0.0590.01070.0010.03100.08880.0190.0330.06680.0332值507020.8480.9890.88110.2840.85520.6180.3740.70310.8877737260.56010.6180.8500.6120.40090.9020.8170.63350.61318936220.77010.8530.6460.7030.87480.6200.5240.60830.58974685710.0770.0470.04300.0330.03010.0380.0530.04990.0756901370.1890.060200.4620.1500.14580.0790.0010.29010.0330小35503波0.8280.76510.43310.8530.74910.8830.6540.90230.8925包44483絡(luò)10.76280.6240.8280.8780.82100.8110.6450.81140.8898譜16615特0.0370.00280.0470.0400.1340.055900.0280.09800.0253征44849值0.1160.105700.1180.0930.07340.0880.0550.08730.09738916900.08970.0520.1660.0580.05470.1770.0690.10870.0618884840.0760.11850.0820.1270.1550.10180.1440.04800.19652088470.0980.03360.0470.1030.1280.07090.0290.0800.091205951750.5760.63240