matlab写一段EMD滤波的程序

以下是用MATLAB编写的EMD滤波程序：

function [IMF,residue] = emd_filter(signal,threshold,max_imf)
% signal: 输入信号
% threshold: 停止分解的阈值
% max_imf: 最大的IMF分量个数

% 初始化
IMF = zeros(max_imf,length(signal));
residue = signal;
imf_num = 0;

% 循环迭代
while imf_num < max_imf
    imf_num = imf_num + 1;
    c = residue;
    n = 0;
    while true
        n = n + 1;
        [upper_env,lower_env] = envelope(c);
        m = (upper_env + lower_env) / 2;
        h = c - m;
        if norm(h) / norm(c) < threshold || n > 100
            break;
        end
        c = h;
    end
    IMF(imf_num,:) = m;
    residue = residue - m;
    if norm(residue) / norm(signal) < threshold
        break;
    end
end

% 剩余项即为趋势项
IMF(imf_num+1,:) = residue;
IMF = IMF(1:imf_num+1,:);
end

使用方法：将需要滤波的信号作为参数传入函数 `emd_filter` 中即可，函数会返回分解得到的所有 IMFs 和剩余项 residue。其中，`threshold` 是一个停止分解的阈值，当分解得到的剩余项的幅值相对于原始信号幅值的比例小于该值时停止分解；`max_imf` 是最大的IMF分量个数，当分解得到的分量个数达到该值时停止分解。