EEMD的快速实现FEEMD

国立中央大学的最新研究成果，FEEMD是EEMD的一种快速实现，较为权威

% The example shows how to decomposed the signal by EEMD
function [] = example_eemd（）   

inputCase = 2;

if（inputCase == 1） % Input Signal
 n = 1024;
 t = 1:n;
 x =  sin（（2*pi*t/16）） + sin（2*pi*t/256）;
elseif（inputCase == 2） % Blood Flow Velocity
 % example in the paper [22]
 % K. Hu M. T. Lo C. K. Peng Y. Liu and V. Novak
 %A Nonlinear Dynamic Approach Reveals a Long-Term Stroke Effect on Cerebral Blood Flow Regulation at Multiple Time Scales
 % PLoS Computational Biology 8 （2012） e1002601.
     load BFVL.mat;
     x = BFVL;
end

t = 1:length（x）; figure; plot（x‘k‘‘LineWIdth‘1）; title（‘input Signal‘）;

fprintf（‘************ RUN RCADA EEMD ***********\n‘）;

Nstd = 0.4;
NE = 200; % # of ensemble
numImf = 8; % # of imfs
runCEEMD = 0;
maxSift = 10;
typeSpline = 2;
toModifyBC = 2;
randType = 2;
seedNo = 0;
checksignal = 1;

tic;

[imf] = rcada_eemd（xNstdNEnumImf）; % run EEMD [3]
   %[imf] = rcada_eemd（x01numImf）; % run EMD [2]
 %[imf] = rcada_eemd（xNstdNEnumImfrunCEEMDmaxSifttypeSplinetoModifyBCrandTypeseedNochecksignal）; % EEMD [3]
 % [imf] = rcada_eemd（xNstdNEnumImf1maxSifttypeSplinetoModifyBCrandTypeseedNochecksignal）; % CEEMD [4]
toc;

%imf = imf‘;

figure;  
nimf = size（imf1）; 
for （m=1:nimf）
  subplot（nimf1m）; plot（timf（m:）‘k‘‘LineWIdth‘1.5）;
end
subplot（nimf11）; title（‘IMFs‘）;

[1];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     756829  2014-04-01 17:01  ref.pdf

     文件       1411  2014-04-01 17:02  example_eemd.m

     文件     112697  2014-04-01 17:02  BFVL.mat

     文件       7900  2014-04-01 17:02  rcada_eemd.m

     文件      21504  2014-04-01 17:03  rcada_emd.mexw64

     文件      17408  2014-04-01 21:39  rcada_emd.mexw32

----------- ---------  ---------- -----  ----

               917749                    6