改進(jìn)的譜峭度方法對(duì)于滾動(dòng)軸承故障診斷的應(yīng)用

1、改進(jìn)的譜峭度方法對(duì)于滾動(dòng)軸承故障診斷的應(yīng)用一、論文步驟：1、對(duì)信號(hào)做四層基于 db10 小波包分解2、對(duì)每層分量信號(hào)做改進(jìn)的峭度值計(jì)算計(jì)算公式為：bnN 1 bij n b 4K ij 4ij N 1 b43、做峭度分布圖4、選出最大峭度的分量，并做解調(diào)5、對(duì)選出的分量做包絡(luò)譜分析6、診斷故障類型二、 Matlab 代碼%軸承外圈 clear allFs=15360;T=1/Fs;N=8192; n=0:N-1; t=n/Fs; f=n*Fs/N;實(shí)驗(yàn)數(shù)據(jù) s3.TXT');%載入數(shù)據(jù)%數(shù)據(jù)波形figure(1) subplot(211)plot(t,X) title(' 實(shí)驗(yàn)

2、數(shù)據(jù)曲線 S');xlabel('t');ylabel('c');Y=fft(X,N); magY=abs(Y); subplot(212);plot(f(1:N/ 2),magY(1:N/2)*2/N) title(' 信號(hào)頻譜圖 ');xlabel('f');ylabel('c');%Daubechies小波包變換， WPT nlevel=4;Dau=wpdec(X,nlevel,'db10'); cfs=zeros(31,8192);k=zeros(1,31);NN=zeros(1,3

3、1);for i=1:31cfs(i,1:length(wpcoef(Dau,i-1)=wpcoef(Dau,i-1);NN(1,i)=length(wpcoef(Dau,i-1);ka=sum(cfs(i,:)-mean(cfs(i,:).4);kb=(NN(1,i)-1)*std(cfs(i,:)4;k(1,i)=ka/kb;end%求出最大峭度值 kur=zeros(5,16);kur(1,:)=k(1,1)*ones(1,16);kur(2,1:8)=k(1,2)*ones(1,8);kur(2,9:16)=k(1,3)*ones(1,8);kur(3,1:4)=k(1,4)*ones

4、(1,4);kur(3,5:8)=k(1,5)*ones(1,4);kur(3,9:12)=k(1,6)*ones(1,4);kur(3,13:16)=k(1,7)*ones(1,4);kur(4,1:2)=k(1,8)*ones(1,2);kur(4,3:4)=k(1,9)*ones(1,2);kur(4,5:6)=k(1,10)*ones(1,2);kur(4,7:8)=k(1,11)*ones(1,2);kur(4,9:10)=k(1,12)*ones(1,2);kur(4,11:12)=k(1,13)*ones(1,2);kur(4,13:14)=k(1,14)*ones(1,2);k

5、ur(4,15:16)=k(1,15)*ones(1,2);kur(5,1:16)=k(1,16:31);Level_x=1:nlevel;Level_x=Level_x;Level_x+log2(3)-1;Level_x=Level_x(:);Level_x=0 Level_x(1:2*nlevel-1)'%圖形縱坐標(biāo)Level_y=Fs*(0:3*2nlevel-1)/(3*2(nlevel+1)+1/(3*2(2+nlevel);figure imagesc(Level_y,Level_x,kur),colorbar;title('峭度分布圖 '); xlabel

6、('頻率 /f'); ylabel('尺度/k');%找出最大峭度值位置temp,tempI=max(kur);%求每列的最大值，并記錄每列最大值的行位 置M,J=max(temp);%求矩陣最大值，并記錄列位置I=tempI(J);%矩陣最大值的行位置%最大峭度值頻段信號(hào)重構(gòu)Cfs=wprcoef(Dau,2(I-1)+J-2);figuresubplot(211)plot(t,Cfs)title(' 峭度最大的頻段信號(hào)重構(gòu)圖像 ');xlabel('時(shí)間/t');ylabel('幅值 /c');%做包絡(luò)譜分析CFs=hilbert(Cfs-mean(Cfs);%hilber對(duì)t 信號(hào)進(jìn)行包絡(luò)分析 mag=abs(CFs);%取模subplot(212);plot(t,mag) title（'hilbert 方法包絡(luò)圖像 '）; xlabel（'時(shí)間/t'）;ylabel（'幅值 /c'）;CFS=fft（mag,N）;Mag=abs（CFS）;%取模 figure;plot（f（1:N/ 2）,Mag（1:N/2）*2/N） title（' 包絡(luò)解調(diào)譜圖 '）; xlabel（'頻率 /f'