matlab emd代码 详细中文注释

% EMD 计算经验模式分解
%
%
%   语法
%
%
% IMF = EMD（X）
% IMF = EMD（X...‘Option_name‘Option_value...）
% IMF = EMD（XOPTS）
% [IMFORTNB_ITERATIONS] = EMD（...）
%
%
%   描述
%
%
% IMF = EMD（X） X是一个实矢量，计算方法参考[1]，计算结果包含在IMF矩阵中，每一行包含一个IMF分量，
% 最后一行是残余分量，默认的停止条件如下[2]：
%
%   在每一个点 mean_amplitude < THRESHOLD2*envelope_amplitude （注：平均幅度与包络幅度的比值小于门限2）
%   &
%   mean of boolean array {（mean_amplitude）/（envelope_amplitude） > THRESHOLD} < TOLERANCE 
%  （注：平均幅度与包络幅度比值大于门限的点数占信号总点数中的比例小于容限）
%   &
%   |#zeros-#extrema|<=1 （注：过零点和极值点个数相等或者相差1）
%
% 这里 mean_amplitude = abs（envelope_max+envelope_min）/2 （注：平均幅度等于上下包络相互抵消后残差的一半的绝对值，理想情况等于0）
% 且 envelope_amplitude = abs（envelope_max-envelope_min）/2 （注：包络幅度等于上下包络相对距离的一半，理想情况等于上下包络本身的绝对值）
% 
% IMF = EMD（X） X是一个实矢量，计算方法参考[3]，计算结果包含在IMF矩阵中，每一行包含一个IMF分量，
% 最后一行是残余分量，默认的停止条件如