geoyang 发表于 2013-3-25 14:20

G Rilling 程序中Hilbert谱问题

各位大侠！我想得到Hilbert时频图中振幅大小的图标，但是用Ｇ Rilling程序中disp_hhs.m程序时，显示的振幅数值为啥为负数啊？

M=max(max(im));

warning off
im = 10*log10(im/M);
warning on

figure

if fs == 0
%imagesc(t,,im,);
imagesc(t,,im,);
ylabel('normalized frequency')
else
imagesc(t,,im,);
% imagesc(t,,im,);
ylabel('frequency')
end
set(gca,'YDir','normal')
xlabel('time')
title('Hilbert-Huang spectrum') 这段程序咋改，才能得到振幅大小啊？

yghit08 发表于 2013-3-25 14:54

嵴提取试试看，应该可以！

geoyang 发表于 2013-3-25 15:41

yghit08 发表于 2013-3-25 14:54 static/image/common/back.gif
嵴提取试试看，应该可以！

M=max(max(im));

warning off
im = 10*log10(im/M);
warning on

哦，我的意思是如何把真实的幅值画出来，这个程序中im好像是dB值啊？

yghit08 发表于 2013-3-25 15:51

geoyang 发表于 2013-3-25 15:41 static/image/common/back.gif
M=max(max(im));

warning off


这就不清楚了！

qqcd 发表于 2013-3-26 18:24

yghit08 发表于 2013-3-25 15:51 static/image/common/back.gif
这就不清楚了！

M=max(max(im)); 提取幅值当中最大值
im = 10*log10(im/M); 这个式子算出来的都是幅值，最大值是0，就是db单位了。

要想把真实的幅值画出来，注释删除上门两句即可了吧...
页: [1]
查看完整版本: G Rilling 程序中Hilbert谱问题