关于HHT应用的问题

svca_0

这是我用HHT变换做出来的,信号是声发射的泄露信号,效果为什么这么不好啊!!!是不是噪声干扰的问题,信号已经用小波除噪了

svca_0

[n,s]=size(imf);
for m=1:n

u(:,m)=imf(m,:)';
end
for k=1:n

subplot(n,1,k)

plot(u(:,k))
end



[imf,ort,nbits] = emd (x,t,[0.05,0.5,0.05]);
emd_visu(x,t,imf,1);
figure(1)
[iha ihf] = hhspectrum(imf);
[im hrtp] = toimage(iha,ihf);
disp_hhs(im,hrtp);
figure;
fs = 2048;
gfs = sum(im,2);
plot([length(gfs):-1:1]/length(gfs)*fs/2,gfs);
xlabel(['Frequency (KHz)'])
ylabel(['Sqrt(Power)'])
title('Hilbert-Huang Marginal Spectrum')

是应用这个程序运行的

kevin19821

出现这种情况是正常的,不知道你用的是什么信号,就从你的IMF来看,并没有把原来信号分解成为一系列的单组分分量,你分解的每个IMF分量都出现了模态混叠,你可以用FFT做下频谱看看.出现这种情况我想有两种原因:
(1)  你分析的信号各个频率接近
(2)  分析的信号的各个幅值造成

svca_0

先谢谢楼上给予的解答
  我用的是管道泄漏的声发射信号,采样频率是2M,楼上说的用FFT做频谱,是不是可以理解成先进行FFT预处理,再用HHT进行变换,刚接触HHT,很多东西还是一知半解,虚心求教.:@)

zhangnan3509

kevin19821 的意思是 用FFT看看，主要的频率成分是什么，这些成分是不是很接近。

并不是先用FFT处理，再用HHT

svca_0

好的,先谢谢各位.我先试试看

chenjh

或者对各个IMF分量做FFT

asdf229955

不知问题解决了没？