集成经验模态分解EEMD代码

EEMD是Ensemble Empirical Mode Decomposition的缩写，中文是集合经验模态分解，是针对EMD方法的不足，提出了一种噪声辅助数据分析方法。EEMD分解原理是当附加的白噪声均匀分布在整个时频空间时，该时频空间就由滤波器组分割成的不同尺度成分组成。

% 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
%   zwu@fsu.edu
%
%



function allmode=eemd（YNstdNE）
xsize=length（Y）;
dd=1:1:xsize;
Ystd=std（Y）;
Y=Y/Ystd;Y=Y/Ystd;
%标准差归一化，这样设定加入噪声额标准差才有意义，最后还要结果还要乘上 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;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2156  2012-04-11 08:35  Copy_of_eemd.m
     文件         727  2012-04-17 20:10  Data_Partition.m
     文件         854  2008-03-27 22:46  dist_value.m
     文件     8205843  2012-06-13 22:33  doppler_data.mat
     文件        2172  2012-04-11 08:25  eemd.asv
     文件        2011  2008-03-27 21:10  eemd.m
     文件        8980  2012-04-26 09:25  EEMD_ANN.m
     文件         374  2011-11-25 15:35  eemd_emd_quzao.asv
     文件         403  2011-11-25 15:36  eemd_emd_quzao.m
     文件        2187  2008-03-27 21:11  extrema.m
     文件        1020  2012-04-17 22:19  Hilbert_FFT_Plot.m
     文件        1825  2008-03-27 21:11  ifndq.m
     文件        1047  2012-04-15 16:24  Node_Index_generate.m
     文件        1094  2012-04-10 21:54  read me.txt
     文件        2945  2008-03-27 22:46  significance.m
     文件        1398  2012-04-18 22:50  zhu.m
     文件     1216738  2011-10-29 07:22  大型旋转机械非平稳振动信号的EEMD降噪方法.pdf
     文件      886664  2011-04-09 15:35  集合经验模式分解在柴油机机械故障诊断中的应用.pdf
     文件     1544371  2011-04-09 15:34  集合经验模式分解在旋转机械故障诊断中的应用.pdf