希尔伯特黄变换

%非主函数，被调用
function n = findpeaks(x)%用于寻找极值点，该函数只会求极大值
%   Find peaks.
%   n = findpeaks(x)
    n = find(diff(diff(x)>0)<0);%一阶导数大于0二阶导数小于0的点
    u = find(x(n+1)>x(n));
    n(u) = n(u) + 1;
end
1
2
3
4
5
6
7
8
9
%非主函数，被调用<br data-filtered="filtered">%判断x是否单调，返回0代表不是单调，返回1代表是单调
function u = ismonotonic(x)
    u1 = length(findpeaks(x))*length(findpeaks(-x));%如果最大/最小值有一个为0即可判断程序满足退出条件了
    if u1 > 0
        u = 0;
    else
        u = 1;
    end
end　

%非主函数，被调用。判断当前x是不是真IMF
function u = isimf(x)
    N  = length(x);