我在用EEMD处理凯斯西储大学的轴承故障数据,具体代码如下
close all;
clear all;
clc;
load('C:\Users\Administrator\Desktop\123\130.mat');%导入数据
fs=12000;%采样率
N=10240;%采样点数
n=0:N-1;
t=n/fs;%信号时间序列
X=X130_DE_time(1:N);%装载驱动计数端的内圈故障数据
y=X;%信号幅值序列
k_in=kurtosis(y);%峭度系数,正常轴承为3左右
figure;%画原始信号时域和频域图
subplot(211);plot(t,y);title('原始信号时域波形');
subplot(212);hua_fft(y,fs,1);title('原始信号频谱');
figure;%原始信号的包络谱
subplot(211);hua_baoluo(y,fs,1);title('原始信号包络谱');
subplot(212);hua_baoluo(y,fs,1);title('原始信号部分频段包络谱');
subplot(212);hua_baoluo(y,fs,1);axis([0 400 0 0.7]);title('原始信号部分频段包络谱');
NE=50;
Nstd=0.01;
imf=eemd(y,0.01,50);%经验模态分解
figure;%画出所有imf及残值
%画出每个IMF分量的辐频曲线
[a,b]=size(im