西储大学轴承故障分析matlab代码

针对西储大学轴承故障数据进行了EMD分解，并对故障特征频率等信息进行了提取

% This is an EMD/EEMD program
%
%   function allmode=eemd（YNstdNE）
%
% INPUT:
%       Y: Inputted data;
%       Nstd: ratio of the standard deviation of the added noise and that of Y;
%       NE: Ensemble number for the EEMD
% OUTPUT:
%       A matrix of N*（m+1） matrix where N is the length of the input
%       data Y and m=fix（log2（N））-1. Column 1 is the original data columns 2 3 ...
%       m are the IMFs from high to low frequency and comlumn （m+1） is the
%       residual （over all trend）.
%
% NOTE:
%       It should be noted that when Nstd is set to zero and NE is set to 1 the
%       program degenerates to a EMD program.
%
% References can be found in the “Reference“ section.
%
% The code is prepared by Zhaohua Wu. For questions please read the “Q&A“ section or
% contact
%   zhwu@cola.iges.org
%

function allmode=eemd（YNstdNE）
xsize=length（Y）;
dd=1:1:xsize;
Ystd=std（Y）;
Y=Y/Ystd;

TNM=fix（log2（xsize））-1;
TNM2=TNM+2;
for kk=1:1:TNM2 
    for ii=1:1:xsize
        allmode（iikk）=0.0;
    end
end

for iii=1:1:NE
    for i=1:xsize
        temp=randn（11）*Nstd;
        X1（i）=Y（i）+temp;
    end

    for jj=1:1:xsize
        mode（jj1） = Y（jj）;
    end
    
    xorigin = X1;
    xend = xorigin;
    
    nmode = 1;
    while nmode <= TNM
        xstart = xend;
        iter = 1;
   
        while iter<=10
            [spmax spmin flag]=extrema（xstart）;
            upper= spline（spmax（:1）spmax（:2）dd）;
            lower= spline（spmin（:1）spmin（:2）dd）;
            mean_ul = （upper + lower）/2;
            xstart = xstart - mean_ul;
            iter = iter +1;
        end
        xend = xend - xstart;
   
   	    nmode=nmode+1;
        
        for jj=1:1:xsize
            mode（jjnmode） = xstart（jj）;
        end
    end
   
    for jj=1:1:xsize
        mode（jjnmode+1）=xend（jj）;
    end
   
    allmode=allmode+mode;
    
end

allmode=allmode/NE;
allmode=allmode*Ystd;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    2928192  2018-08-06 10:55  EMD\130.mat

     文件       2011  2018-10-16 09:26  EMD\eemd.m

     文件     838457  2018-09-28 08:57  EMD\EMD.docx

     文件      22275  2015-05-29 10:37  EMD\emd.m

     文件    1122115  2018-09-27 11:18  EMD\emd.png

     文件     905216  2018-09-27 10:00  EMD\emd.ppt

     文件        435  2015-05-31 11:20  EMD\exzample.m

     文件       2326  2015-05-31 15:55  EMD\hilbert.m

     文件        427  2015-05-31 16:04  EMD\hua_baoluo.m

     文件       1435  2015-05-31 11:24  EMD\hua_fft.m

     文件       2010  2018-10-16 16:01  EMD\innerguzhang.m

    ..A..H.       162  2018-09-27 10:30  EMD\~$EMD.docx

     目录          0  2018-10-16 15:38  EMD

----------- ---------  ---------- -----  ----

              5825061                    13