求大神 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 结果
= 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
回复 1 # 挽一轮明月 的帖子
不知道现在是否有人证明HHT是可逆的,这是重构信号的前提,建议先查查相关文献。 原文写的很明白,可以重构,具体方法不知道。 楼主注意这一步“imf =emd (x)”,得到的imf应该是一个m*n阶矩阵,意味着m各imf分量和n个数据点。可以自编程序将你所感兴趣的imf分量进行叠加。
页:
[1]