轴承故障诊断matlab代码四种方法

数据分析中主要针对内外圈故障，应用4种故障分析方法分别进行分析.通过对同一转动频率的信号采用不同处理分析方法及同一分析方法对不同转动频率的的比较，可以发现时域无量纲参数分析、FFT分析这两方法可以有效的发现故障，但是不容易得出故障信号频率，因此不能分析出具体是哪一种故障。Hilbert包络谱分析和Haar小波分析和数学形态学分析相对于前面的几种方法，可以更好的观测出故障信号频率，从而可以分析故障类型。

% ********** 内圈故障300r/min **********
clear;
clc;
t=0:1:16383;       %采样步长
fs=10000;          %采样频率
y=importdata（‘neiquan_10k_32768_5.dat‘）;   %读取数据
N=length（t）;                 
z=y-mean（y）;      %零均值处理
figure（1）;
subplot（211）;
plot（tz）;        %原始数据图形
title（‘内圈300r/min‘）;xlabel（‘采样数‘）;ylabel（‘幅值‘）;
subplot（212）;
f=fft（zN）;   
plot（（0:N/2-1）/N*fsabs（f（1:N/2）））;  
title（‘FFT变换‘）;xlabel（‘频率/Hz‘）;ylabel（‘幅值‘）;

% ********** 内圈故障450r/min **********
clear;
clc;
t=0:1:16383;       %采样步长
fs=10000;          %采样频率
y=importdata（‘neiquan_10k_32768_7.5.dat‘）;   %读取数据
N=length（t）;                 
z=y-mean（y）;      %零均值处理
figure（2）;
subplot（211）;
plot（tz）;        %原始数据图形
title（‘内圈450r/min‘）;xlabel（‘采样数‘）;ylabel（‘幅值‘）;
subplot（212）;
f=fft（zN）;   
plot（（0:N/2-1）/N*fsabs（f（1:N/2）））;    
title（‘FFT变换‘）;xlabel（‘频率/Hz‘）;ylabel（‘幅值‘）;

% ********** 外圈故障300r/min **********
clear;
clc;
t=0:1:16383;       %采样步长
fs=10000;          %采样频率
y=importdata（‘waiquan_10k_32768_5.dat‘）;   %读取数据
N=length（t）;                 
z=y-mean（y）;      %零均值处理
figure（3）;
subplot（211）;
plot（tz）;        %原始数据图形
title（‘外圈300r/min‘）;xlabel（‘采样数‘）;ylabel（‘幅值‘）;
subplot（212）;
f=fft（zN）;   
plot（（0:N/2-1）/N*fsabs（f（1:N/2）））;    %包络谱
title（‘FFT变换‘）;xlabel（‘频率/Hz‘）;ylabel（‘幅值‘）;

% ********** 外圈故障450r/min **********
clear;
clc;
t=0:1:16383;       %采样步长
fs=10000;          %采样频率
y=importdata（‘waiquan_10k_32768_7.5.dat‘）;   %读取数据
N=length（t）;                 
z=y-mean（y）;      %零均值处理
figure（4）;
subplot（211）;
plot（tz）;        %原始数据图形
title（‘外圈450r/min‘）;xlabel（‘采样数‘）;ylabel（‘幅值‘）;
subplot（212）;
f=fft（zN）;   
plot（（0:N/2-1）/N*fsabs（f（1:N/2）））;    
title（‘FFT变换‘）;xlabel（‘频率/Hz‘）;ylabel（‘幅值‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-05-24 10:47  程序\
     目录           0  2019-05-24 10:47  程序\FFT\
     文件        1905  2019-01-04 17:03  程序\FFT\fft_.m
     文件      287147  2018-12-28 20:12  程序\FFT\neiquan_10k_32768_5.dat
     文件      286814  2018-12-28 20:13  程序\FFT\neiquan_10k_32768_7.5.dat
     文件      288601  2018-12-28 20:27  程序\FFT\waiquan_10k_32768_5.dat
     文件      286772  2018-12-28 20:27  程序\FFT\waiquan_10k_32768_7.5.dat
     文件       87334  2019-01-04 16:45  程序\FFT\内300.jpg
     文件       86840  2019-01-04 16:51  程序\FFT\内450.jpg
     文件       86188  2019-01-04 17:04  程序\FFT\外300.jpg
     文件       88292  2019-01-04 17:05  程序\FFT\外450.jpg
     目录           0  2019-05-24 10:47  程序\包络\
     文件         793  2019-01-04 12:02  程序\包络\chengxu.m
     文件      286814  2018-12-28 20:13  程序\包络\neiquan_10k_32768_7.5.dat
     文件      286772  2018-12-28 20:27  程序\包络\waiquan_10k_32768_7.5.dat
     文件       29234  2019-01-04 12:02  程序\包络\内1.jpg
     文件       19564  2019-01-04 12:02  程序\包络\内2.jpg
     文件       29450  2019-01-04 12:04  程序\包络\外1.jpg
     文件       24179  2019-01-04 12:04  程序\包络\外2.jpg
     目录           0  2019-05-24 10:47  程序\小波外圈300\
     文件      288601  2018-12-28 20:27  程序\小波外圈300\10k_32768_5.dat
     文件        3334  2019-01-04 12:51  程序\小波外圈300\test_haar.m
     文件       63440  2019-01-04 12:56  程序\小波外圈300\外1.jpg
     文件      101175  2019-01-04 12:56  程序\小波外圈300\外2.jpg
     文件       92706  2019-01-04 12:57  程序\小波外圈300\外3.jpg
     文件       47091  2019-01-04 12:53  程序\小波外圈300\外4.jpg
     文件       46093  2019-01-04 12:52  程序\小波外圈300\外5.jpg
     文件       33861  2019-01-04 12:52  程序\小波外圈300\外6.jpg
     目录           0  2019-05-24 10:47  程序\无纲量内外\
     文件      287147  2018-12-28 20:12  程序\无纲量内外\neiquan_10k_32768_5.dat
     文件      286814  2018-12-28 20:13  程序\无纲量内外\neiquan_10k_32768_7.5.dat
............此处省略4个文件信息