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和残余信号求和,重构信号。
在这里插入图片描述

©️2020 CSDN 皮肤主题: 数字20 设计师:CSDN官方博客 返回首页