EMD和BEMD算法实现（matlab代码）

今天结合自己在网上找的一些实现代码，稍微修改进行测试，没有进行更多的实验，可能在一些问题的处理上还是比较毛糙的。
EMD

function imf = emd(x)
% Empiricial Mode Decomposition (Hilbert-Huang Transform)
% EMD分解或HHT变换
% 返回值为cell类型，依次为一次IMF、二次IMF、...、最后残差
x   = transpose(x(:));
imf = [];
while ~ismonotonic(x)
    x1 = x;
    sd = Inf;
    while (sd > 0.1) || ~isimf(x1)
        s1 = getspline(x1);         % 极大值点样条曲线
        s2 = -getspline(-x1);       % 极小值点样条曲线
        h = x1-(s1+s2)/2;