|
|
请教如何画HHT的三维时频图???
imf = emd(ss);
emd_visu(ss,imf)
for j=1:2
[A,f,tt]=hhspectrum(imf(j,:));
%axis([0 1024 0 0.1])
[im,tt]=toimage(A,f);
disp_hhs(im,[],fs);
mesh(tt,f,A);
%axis([0 512 0 6000])
%set(gca,'YDir','normal')
%colormap(flipud(hot));
end
程序运行出错 问题是 tt是1*512 点 f 也是1*512点 但是A是400*512 所以没能画出三维图 。 [A,f,tt]=hhspectrum(imf(j,:)); 这一句就是求出的瞬时频率,时间 幅值的函数啊 为什么A的维数 不对呢? |
|