emdcmp 对心电信号进行EMD分解

对心电信号进行EMD分解，比较EMD、EEMD、CEEMD、CEEMDAN的分解差别

% 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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2931  2015-08-22 14:57  emd_cmp.m
     文件        3552  2015-08-05 21:12  ceemdan.m
     文件        2295  2015-08-05 10:54  ecg.mat
     文件        1616  2015-08-05 20:42  eemd_n.m
     文件        2187  2015-08-05 17:11  extrema.m
     文件       22275  2015-07-23 20:33  emd.m
     文件        2603  2015-08-05 16:26  ceemd.m