eemd的r语言序列_EMD与EEMD程序

%%%%%%%%%%%

载入信号

x=load('1.txt');%

产生信号

N=le ngth(x);%

采样点数

fs=2000;%

采样频率

dt=1/fs;%

采样时间间隔

t=(0:N-1)*dt;%

产生时间序列

%%%%%%%%EMD

imf=emd(x);

%EMD.M(EMDW)

%G.Rilling,July2002

%

%computesEMD(EmpiricalModeDecomposition)accordingto:

%

%N.E.Huangetal.,"Theempiricalmodedecompositionandthe %Hilbertspectrumfornon-

linearandnonstationarytimeseriesanalysis," %Proc.RoyalSoc.LondonA,Vol.454,pp.903-995,1998

%

%withvariationsreportedin:

%

%G.Rilli ng,P.Fla ndr inan dP.Go n?alv

es

%"OnEmpiricalModeDecompositionanditsalgorithms"

%IEEE-EURASIPWorkshoponNonlinearSignalandImageProcessing %NSIP-03,Grado(I),June2003

%

%stoppingcriterionforsifting:

%ateachpoint:me