emdc_fix.m

经验模式分解源码 学习希尔波特黄变换的参考程序 源于Internet网络

%EMDC_FIX  computes Empirical Mode Decomposition
%
%
%   Syntax
%
%
% [IMF,NB_ITERATIONS]=EMDC_FIX(T,X,NB_ITERATONS,MAX_IMFS);
%
%
%   Description
%
%
% computes EMD according to [1] with NB_ITERATONS sifting iterations for each IMF
%
%   mean of boolean array {(mean_amplitude)/(envelope_amplitude) > THRESHOLD} < TOLERANCE
%   &
%   |#zeros-#extrema|<=1
%
% inputs:	
%       - T: sampling times. If T=[], the signal is assumed uniformly sampled.
%       - X: analyzed signal
%       - NB_ITERATIONS: number of sifting iterations to be performed to
%         extract each IMF. If NB_ITERATIONS is empty or unspecified, 10 iterations 
%         are performed by default.
%         Note: The effective number of sifting iterations might be less 
%         than NB_ITERATIONS for the last modes if the sifting process has 
%         to be stopped because of a lack of extrema.
%       - MAX_IMFS: maximum number of IMFs to be extracted. If MAX_IMFS is
%         zero, empty or unspecified, the default behavior is to extract as
%         many IMFs as possible.
%         
% outputs: 
%		- IMF: intrinsic mode functions (IMFs) (last line = residual)
%		- NB_ITERATIONS: effective number of sifting iterations for each mode
%
%
%   Examples
%
%
% workspace: 
%  T: 1xN time instants
%  X: 1xN signal data 
%
%>>IMF = EMDC_FIX(T,X);
%>>[IMF,NB_IT] = EMDC_FIX([],X);
%>>IMF = EMDC_FIX(T,X,0.1);
%>>IMF = EMDC_FIX(T,X,[0.1,0.1]);
%>>[IMF,NB_IT] = EMDC_FIX([],X,[],4);
%
%
%   References
%
%
% [1] N. E. Huang et al., "The empirical mode decomposition and the
% Hilbert spectrum for non-linear and non stationary time series analysis",
% Proc. Royal Soc. London A, Vol. 454, pp. 903-995, 1998
%
%
% See also
%  emd_visu (visualization),
%  emd (slow but has many options),
%  emdc (fast implementation of EMD with different stopping criterion)
%  hhspectrum (Hilbert-Huang spectrum)
%
%
% G. Rilling, last modification: 3.2007
% [e-mail]
%
% code based on a student project by T. Boustane and G. Quellec, 11.03.2004
% supervised by P. Chainais (ISIMA - LIMOS - Universite Blaise Pascal - Clermont II
% email : [e-mail]).