kurtosis 滚动轴承故障诊断

滚动轴承故障诊断， 小波包分解 峭度最大尺度重构

%外圈故障特征频率 理论值 107.36Hz    
clc
clear all
load OR007@12_0_drive;   %outer-race fault（drive）
data1=X156_DE_time;  
load normal_0_1797;     %Normal baseline Data（drive）
data2=X097_DE_time;
N=4096;                %采样点数
fs=12000; 
t=1:N;
t=t/fs;
xdata1=data1（1:4096）;xdata_outer=（xdata1-mean（xdata1）/std（xdata11））;
xdata2=data2（1:4096）;xdata_normal=（xdata2-mean（xdata2）/std（xdata21））;%零均值处理
 %% ====================故障信号小波包=====================%%
wpt1=wpdec（xdata_outer5‘db1‘‘shannon‘）;
%% =================故障信号峭度值========
for i=1:32
    c1=wpcoef（wpt1[5i-1]）;
    K1（i）=kurtosis（c1）;
end
figure;
subplot（311）;
bar（K1）;
title（‘外圈故障信号小波包分解各节点峭度值‘）;
%% ==============正常信号小波包===============%%
wpt2=wpdec（xdata_normal5‘db1‘‘shannon‘）;
%% ===============正常信号峭度值 =============

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1540  2015-04-23 16:35  outer_normal_kurtosis.m
     文件     2935152  2015-01-02 15:28  OR007@12_0_drive.mat