基于MATLAB最简洁的EMD程序，适合初学者

% EMD:  Emprical mode decomposition
%
% imf = emd（xn）
%
% x   - input signal （must be a column vector）
% n   - number of intrinsic mode functions
%
% imf - Matrix of intrinsic mode functions （each as a row）
%
% See:  Huang et al Royal Society Proceedings on Math Physical 
%       and Engineering Sciences vol. 454 no. 1971 pp. 903-995 
%       8 March 1998
%
% Author: Ivan Magrin-Chagnolleau  
% 
function imf = emd（xn）;%%最好把函数名改为emd1之类的，以免和Grilling的emd冲突                        
                                  %%n为你想得到的IMF的个数
c = x（:）‘; % copy of the input signal （as a row vector）
N = length（x）;-
% loop to decompose the input signal into n successive IMFs
imf = []; % Matrix which will contain the successive IMF and the residuefor t=1:n 

while