WJzhengdong 发表于 2010-11-14 10:24

HHT频谱图为什么会这样

求组:函数为z=sin(8*pi*t)+sin(16*pi*t)+sin(32*pi*t),为什么我自己的做的图,和原图有那么大的差异,而且右边的彩条刻度还是负的。请问哪位大哥大姐能给说说这个问题出在哪里,

WJzhengdong 发表于 2010-11-14 10:28

nt=1:1000;
fs=2000;
t=(nt-1)*1/fs;
z1=sin(8*pi*t)+sin(16*pi*t)+sin(32*pi*t)
imf=emd(z1);
emd_visu(z1,1:length(z1),imf);
=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im,[],fs);
ylim()
这是程序,能给看看问题出在哪里了吗,另外提一句采样频率是我自己设置的,不知道原图用的是多少。

hahaer 发表于 2010-11-14 14:49

个人认为你的也就是第二个图得到的结果是正确的,但是你给出的程序反而觉得非常奇怪,利用你的程序是得不到你所贴出的结果的。
两个图不一样结果的原因:第一个给出的幅度值
                                          第二个给出的dB值

WJzhengdong 发表于 2010-11-14 20:15

可是我运行出来就是这样啊,请问你运行后出来什么样的结果啊?和原图一样吗?怎么从DB图得到幅值图呢?不是很懂这个DB

shawnbon 发表于 2010-11-14 20:30

改一下程式碼,把10*log10這段消掉,直接輸出就行啦!

WJzhengdong 发表于 2010-11-14 21:28

我试过的,不行,出现的图一片深红色,什么都没有,你对这方面懂得多吗,加我QQ,请教一下445876674
页: [1]
查看完整版本: HHT频谱图为什么会这样