emd分解程序代码，里面有详细介绍，请下载使用(function imf = emd(x)
% Empiricial Mode Decomposition (Hilbert-Huang Transform)
% imf = emd(x)
% Func : findpeaks

x　 = transpose(x(:));
imf = [];
while ~ismonotonic(x)
　 x1 = x;
　 sd = Inf;
　 while (sd > 0.1) || ~isimf(x1)
　　　s1 = getspline(x1);
　　　s2 = -getspline(-x1);
　　　x2 = x1-(s1+s2)/2;
　　　
　　　sd = sum((x1-x2).^2)/sum(x1.^2);
　　　x1 = x2;
　 end
　
　 imf{end+1} = x1;
　 x　　　　　= x-x1;
end
imf{end+1} = x;)