EMD+相关系数+Hilbert+包络谱

clc
clear all
% z=load（‘C:\Users\wangkun\Desktop\轴承诊断\test2.mat‘）;
% x=z.X105_DE_time（1:1024）;
% load （‘IR007_0.mat‘）; 
load（‘F:\轴承数据\CaseWesternReserveUniversityData-master\12k_Drive_End_OR007@3_2_146.mat‘）
a=X146_DE_time;

x=a（10001:15000）;
N=5000;
fs=12000;
n=0:N-1;
f=n*fs/N;
lag=N;
n=0:N-1;
t=n/fs;      
imf=emd（x）;
[mn]=size（imf）; %imf为一m*n阶矩阵，m是imf分量，n为数据点
emd_visu（x1:length（x）imfm）; %实信号的信号重构及emd结果显示函数
for i=1:m
  a1（i）=kurtosis（imf（i:））;%峭度
  b1（i）=mean（imf（i:））; %均值;
  c1（i）=var（imf（i:））;  %方差;
  d1（i）=std（imf（i:））;  %均方值
  e1（i）=std（imf（i:））.^0.5; %均方根值
  f1（i）=s