G Rilling 程序中Hilbert谱问题

geoyang

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

M=max(max(im));

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

figure

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

yghit08

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

geoyang

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

                               
登录/注册后可看大图

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

M=max(max(im));

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

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

yghit08

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

                               
登录/注册后可看大图

M=max(max(im));

warning off

这就不清楚了！

qqcd

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

                               
登录/注册后可看大图

这就不清楚了！

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

要想把真实的幅值画出来，注释删除上门两句即可了吧...