ceemd 一段eemd分解的程序

这是一段eemd分解的程序，可以直接使用，可用于故障信号的分解，提取故障特征，进行故障判断的一段程序。

% Y: Inputted data;
% Nstd: ratio of the standard deviation of the added noise and that of Y;
% NE: Ensemble member being used
% TNM: total number of modes （not including the trend）
%
function allmode=ceemd（YNstdNETNM）
% find data length
xsize=length（Y）;
dd=1:1:xsize;
% Nornaliz data
Ystd=std（Y）;
Y=Y/Ystd;
% Initialize saved data
TNM2=TNM+2;
for kk=1:1:TNM2
    for ii=1:1:xsize
        allmode（iikk）=0.0;
    end
end

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

    % sifting X1
    xorigin = X1;
    xend = xorigin;
% save the initial data into the first column
    for jj=1:1:xsize
        mode（jj1） = xorigin（jj）;
    end
    nmode = 1;
    while nmode <= T

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2603  2014-04-01 21:11  ceemd.m