大神,麻烦帮我看看,这是我画的三维时频幅值图,我的样本是西储大学滚动轴承数据样本,那些彩色的块状的不知道是什么东西,也不知道是不是代码有问题,代码如下,
s=gdt(:,1);
imf=emd(s);
[A,fa,tt]=hhspectrum(imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for i=1:size(E,1);
faa=fa(i,:);
[FA,TT1]=meshgrid(faa,tt1);%三维图显示HHT时频图
figure(11);
surf(FA,TT1,E);
title('HHT时频谱三维显示');
hold on
end