Matlab实现EMD（经验模式分解）

Matlab实现EMD（经验模式分解）

背景：对8通道表面肌电信号进行EMD去噪

%%
%EMD分解
[m,n] = size(emg_data);
refactor = zeros(m,n);%重构信号储存于此
for i = 1:n%每一列单独进行emd分解
    c = emd(emg_data(:,i),'Interpolation','pchip');%c储存分解后的IMF1、IMF2到IMF9和残差
    for j = 1:4%取前4列IMF，求和-->重构信号
        refactor(:,i) = refactor(:,i)+c(:,j);
    end
end
emg_data = refactor;