改进的谱峭度方法对于滚动轴承故障诊断的应用.docVIP

改进的谱峭度方法对于滚动轴承故障诊断的应用 一、论文步骤： 对信号做四层基于db10小波包分解 对每层分量信号做改进的峭度值计算 计算公式为： 做峭度分布图 选出最大峭度的分量，并做解调 对选出的分量做包络谱分析 诊断故障类型 二、Matlab代码 %轴承外圈 clear all Fs=15360; T=1/Fs; N=8192; n=0:N-1; t=n/Fs; f=n*Fs/N; X=load(C:\Users\loveh\Desktop\实验数据\s3.TXT);%载入数据 %数据波形 figure(1) subplot(211) plot(t,X) title(实验数据曲线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(信号频谱图);xlabel(f);ylabel(c); %Daubechies小波包变换，WPT nlevel=4; Dau=wpdec(X,nlevel,db10); cfs=zeros(31,8192); k=zeros(1,31); NN=zeros(1,31); for i=1:31 cfs(i,1:length(wpcoef(Dau,i-1)))=wpcoef(Dau,i-1); N