toniling310 发表于 2008-5-26 12:53

还是边际谱问题 请高手帮助!

我是名正在做毕业设计的大四学生 ,题目是hht在功率密度谱估计的应用,以前完全没有接触过hht,所以很多很多问题都不懂,现在做到边际谱了 程序很简单 就是把破凰的程序改了一下clear;
fs=1000;   
N=1000;   
t=1/fs:1/fs:1;
xn = sin(2*pi*80*t)+3*sin(2*pi*140*t) ;
imf=emd(xn);
=hhspectrum(imf(1:end-1,:));
=toimage(A,fa,tt,length(tt));
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
但是出来的图却是这样:

这是怎么回事,我觉得谱峰应该很窄而且比较平滑才对,
还有边际谱最后能不能得到功率谱呢 ?

Rovis 发表于 2008-5-27 00:45

回复 楼主 的帖子

主要的原因不是因为边际谱本身出的问题,而是EMD分解出来的IMF的质量有问题,现在正在尝试解决中

toniling310 发表于 2008-5-27 09:27

我的直接就用了emd(xn)这个语句 是不是这样不妥呢?
页: [1]
查看完整版本: 还是边际谱问题 请高手帮助!