jiangxiaolong 发表于 2010-5-23 10:39

一个例子的hht图错误

t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.4*sin(2*pi*40*t+140);
z=x+y;
plot(t,z);
imf=emd(z);
emd_visu(z,1:length(z),imf,1);
=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray));
这是我在论坛上找的例子，但是运行了一下得到figure4与他给的图不一样，我的图的频率在上边（应该在下边的），请各位大侠指点一二，小弟不胜感激。

Minnie 发表于 2010-5-23 10:51

我知道这个，你先帮我看看我发的“调频信号”帖子

Minnie 发表于 2010-5-23 10:55

是不是需要这条语句，或者去掉这条语句E=flipud(E);

jiangxiaolong 发表于 2010-5-23 19:15

回复 板凳 Minnie 的帖子

谢谢你的提议，但是我弄不出来，希望能给一个详细点的解答。“ 调频信号”的帖子我没有找到，能否给个链接，我看看会不会。

xiangyu537 发表于 2010-5-23 19:42

在程序中加入这条语句im=flipud(im);    然后运行试试。
t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.4*sin(2*pi*40*t+140);
z=x+y;
plot(t,z);
imf=emd(z);
emd_visu(z,1:length(z),imf,1);
=hhspectrum(imf);
=toimage(A,f);
im=flipud(im);
disp_hhs(im);
colormap(flipud(gray));
页: [1]
查看完整版本: 一个例子的hht图错误