EMD的MATLAB程序

用MATLAB编写的emd程序，有比较详细的注释，能得到结果。

function imf = emd（x）
% Empiricial Mode Decomposition （Hilbert-Huang Transform）
% EMD分解或HHT变换
% 返回值为cell类型，依次为一次IMF、二次IMF、...、最后残差
x   = transpose（x（:））;%TRANSPOSE函数可返回转置单元格区域，即将行单元格区域转置成列单元格区域
imf = [];

while ~ismonotonic（x） %当x不是单调函数，分解终止条件
    x1 = x;
sd = Inf;%均值
%直到x1满足IMF条件，得c1
    while （sd > 0.1） || ~isimf（x1） %当标准偏差系数sd大于0.1或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;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-04-01 17:48  EMD\
     文件         763  2018-03-31 21:06  EMD\emd.m
     文件         206  2018-03-31 21:08  EMD\FFTAnalysis.m
     文件         992  2018-03-31 22:02  EMD\findpeaks.m
     文件         125  2018-03-31 21:07  EMD\getspline.m
     文件         235  2018-03-31 21:08  EMD\HilbertAnalysis.m
     文件         294  2018-03-31 21:09  EMD\isimf.m
     文件         176  2018-03-31 21:07  EMD\ismonotonic.m
     文件        1688  2018-03-31 21:09  EMD\plot_hht.m
     文件        1784  2018-03-31 21:46  EMD\test.m