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;

EMD分解为高频IMF、低频IMF和残余信号，对高频IMF进行小波阈值去噪，再对去噪后的IMF’、低频IMF和残余信号求和，重构信号。