一个例子的hht图错误

jiangxiaolong

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);
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray));
这是我在论坛上找的例子，但是运行了一下得到figure4与他给的图不一样，我的图的频率在上边（应该在下边的），请各位大侠指点一二，小弟不胜感激。

Minnie

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

Minnie

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

jiangxiaolong

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

xiangyu537

在程序中加入  这条语句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);
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
im=flipud(im);  
disp_hhs(im);
colormap(flipud(gray));