求大神 EMD分解后如何重构信号

挽一轮明月

这是程序

fs = 1000;
t = 0:1/fs:1;
f2= 4;
w2 = 0;
B2 = 0;
x=sin(2*pi*f2*t+w2)+B2;
x=awgn(x,10,0)
num = length(t);
figure(1)
plot(x)
%emd
imf =emd (x);
%plot emd 结果
[N_imf,n_imf] = size(imf);
% compute the correlation coefficients
for i= 1:N_imf
%     corr = corrcoef(imf(i,:),x);
    figure(2)
    subplot(N_imf,1,i);
    plot(t(1:num),imf(i,:))
% figure(2)
%  plot(t(1:num),imf(1,:))
end

529899778

回复 1 # 挽一轮明月 的帖子

不知道现在是否有人证明HHT是可逆的，这是重构信号的前提，建议先查查相关文献。

yghit08

原文写的很明白，可以重构，具体方法不知道。

粥糊了耶

楼主注意这一步“imf =emd (x)”，得到的imf应该是一个m*n阶矩阵，意味着m各imf分量和n个数据点。可以自编程序将你所感兴趣的imf分量进行叠加。