经验模态分解CEEMDAN

function [modes its]=ceemdan（xNstdNRMaxIter）

% WARNING: for this code works it is necessary to include in the same
%directoy the file emd.m developed by Rilling and Flandrin.
%This file is available at %http://perso.ens-lyon.fr/patrick.flandrin/emd.html
%We use the default stopping criterion.
%We use the last modification: 3.2007
% 
% This version was run on Matlab 7.10.0 （R2010a）
%----------------------------------------------------------------------
%   INPUTs
%   x: signal to decompose
%   Nstd: noise standard deviation
%   NR: number of realizations
%   MaxIter: maximum number of sifting iterations allowed.
%
%  OUTPUTs
%  modes: contain the obtained modes in a matrix with the rows being the modes        
%   its: contain the sifting iterations needed for each mode for each realization （one row for each realization）
% -------------------------------------------------------------------------
%  Syntax
%
%  modes=ceemdan（xNstdNRMaxIter）
%  [modes its]=ceemdan（xNstdNRMaxIter）
%
%--------------------------------------------------------------------------
% This algorithm was presented at ICASSP 2011 Prague Czech Republic
% Plese if you use this code in your work please cite the paper where the
% algorithm was first presented. 
% If you use this code please cite:
%
% M.E.TORRES M.A. COLOMINAS G. SCHLOTTHAUER P. FLANDRIN
%  “A complete Ensemble Empirical Mode decomposition with adaptive noise“ 
%  IEEE Int. Conf. on Acoust. Speech and Signal Proc. ICASSP-11 pp. 4144-4147 Prague （CZ）
%
% -------------------------------------------------------------------------
% Date: June 062011
% Authors:  Torres ME Colominas MA Schlotthauer G Flandrin P.
% For problems with the code please contact the a