%EMD  computes Empirical Mode Decomposition%%%   Syntax%%% IMF = EMD(X)% IMF = EMD(X,...,'Option_name',Option_value,...)% IMF = EMD(X,OPTS)% [IMF,ORT,NB_ITERATIONS] = EMD(...)%%%   Description%%% IMF = EMD(X) where X is a real vector computes the Empirical Mode% Decomposition [1] of X, resulting in a matrix IMF containing 1 IMF per row, the% last one being the residue. The default stopping criterion is the one proposed% in [2]:%%   at each point, mean_amplitude < THRESHOLD2*envelope_amplitude%   &%   mean of boolean array {(mean_amplitude)/(envelope_amplitude) > THRESHOLD} < TOLERANCE%   &%   |#zeros-#extrema|<=1%% where mean_amplitude = abs(envelope_max+envelope_min)/2% and envelope_amplitude = abs(envelope_max-envelope_min)/2% % IMF = EMD(X) where X is a complex vector computes Bivariate Empirical Mode% Decomposition [3] of X, resulting in a matrix IMF containing 1 IMF per row, the% last one being the residue. The default stopping criterion is similar to the% one proposed in [2]:%%   at each point, mean_amplitude < THRESHOLD2*envelope_amplitude%   &%   mean of boolean array {(mean_amplitude)/(envelope_amplitude) > THRESHOLD} < TOLERANCE%% where mean_amplitude and envelope_amplitude have definitions similar to the% real case%% IMF = EMD(X,...,'Option_name',Option_value,...) sets options Option_name to% the specified Option_value (see Options)%% IMF = EMD(X,OPTS) is equivalent to the above syntax provided OPTS is a struct % object with field names corresponding to option names and field values being the % associated values %% [IMF,ORT,NB_ITERATIONS] = EMD(...) returns an index of orthogonality%                       ________%         _  |IMF(i,:).*IMF(j,:)|%   ORT = \ _____________________%         /%