二维经验模态分解的matlab代码

以下是二维经验模态分解 (2D-EMD) 的 MATLAB 代码示例：

function [IMFs, res] = emd2(IMF)
% 2D-EMD
% IMF：输入的二维数据
% IMFs：输出的二维 IMFs
% res：剩余项

% 初始化
IMFs = {};
res = IMF;
n = 1;

% 处理剩余项
while ~isempty(res)
    % 计算当前剩余项的极值点
    max_idx = imregionalmax(res);
    min_idx = imregionalmin(res);
    extr_idx = max_idx | min_idx;
    
    % 如果没有极值点，则将剩余项作为最后一个 IMF
    if ~any(extr_idx(:))
        IMFs{n} = res;
        break;
    end
    
    % 极值点插值
    extr_points = find(extr_idx);
    [extr_y, extr_x] = ind2sub(size(res), extr_points);
    extr_val = res(extr_points);
    [xq, yq] = meshgrid(1:size(res, 2), 1:size(res, 1));
    Vq_max = griddata(extr_x, extr_y, extr_val, xq, yq, 'v4');
    Vq_min = griddata(extr_x, extr_y, extr_val, xq, yq, 'v4');
    
    % 计算当前 IMF
    imf = res - (Vq_max + Vq_min) / 2;
    IMFs{n} = imf;
    n = n + 1;
    
    % 更新剩余项
    res = res - imf;
end
end

这里的 `IMF` 是一个二维矩阵，`IMFs` 是一个二维矩阵数组，包含所有的 IMFs，`res` 是一个二维矩阵，表示分解后的剩余项。代码中使用了 MATLAB 自带的函数 `imregionalmax` 和 `imregionalmin` 来计算极值点，使用了 `griddata` 函数来进行插值。