请教关于disp_hhs作图的问题,急

laughing

搞了两天,实在没弄明白,没办法,又上来求助了,实在不好意思.
以下是小弟练习作hht谱图用的小程序:
clear
fs=1000.1;%采样频率
t=0:10/1000:10;&共1001个数据点,一共采了10秒
x=exp(-0.3*t).*cos(pi*t.^2+2*pi*t);
plot(t,x);%原始数据序列
imf=emd(x);
emd_visu(x,1:length(x),imf);
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im,[],fs);

这是我画出的图,纵坐标表示的是真实频率(不是规则化频率),横坐标就不对了,我想让横坐标显示真实的时间,即从0到10.该怎么做?
我把disp_hhs(im) ; disp_hhs(im,t) ; disp_hhs(im,inf)
      disp_hhs(im,t,inf) ; disp_hhs(im,inf,fs) ; disp_hhs(im,[],fs); disp_hhs(im,t,inf,fs)
这些命令全用了,还是得不到横坐标为真实时间,纵坐标为真实频率的图.我感觉应该使用disp_hhs(im,[],fs)这条命令.
但是disp_hhs(im,t,[],fs) ; disp_hhs(im,t,inf,fs)这两条命令,运行后提示:Too many input arguments.
请问我该如何达到目的?问题可能很弱弱,烦请达人指点一下.

zhlong

直接到Figure的Edit菜单中修改axes properties！

laughing

原帖由 zhlong 于 2007-6-20 11:54 发表

                               
登录/注册后可看大图

直接到Figure的Edit菜单中修改axes properties！

不才刚才按您的意思改了一下，觉得还是不对.
设置x轴的范围是0:10,结果出来的图为原图的局部放大.

laughing

还有个问题,就是
运行disp_hhs(im,t,[],fs) ; disp_hhs(im,t,inf,fs)这两条命令,运行后提示:Too many input arguments
怎么回事?

w89986581

set(gca,'xtick',[0:100:1000],'xticklabel',[0:1:10])

laughing

多谢了:)
可是,谁能告诉我,运行disp_hhs(im,t,[],fs) ; disp_hhs(im,t,inf,fs)这两条命令,运行后提示:Too many input arguments
怎么回事?

zhlong

输入宗量太多，help一下disp_hhs看看它的函数使用说明！

laughing

disp_hhs(im,t,[],fs) ; disp_hhs(im,t,inf,fs)这两条命令就是disp_hhs的使用说明举的例子!

cboboc

楼主这个是怎样得到真实频率的呢，我按你的算法来运算，得出来的还是瞬时频率 啊 ？

bcyangbc

原帖由 laughing 于 2007-6-20 11:32 发表

                               
登录/注册后可看大图

搞了两天,实在没弄明白,没办法,又上来求助了,实在不好意思.
以下是小弟练习作hht谱图用的小程序:
clear
fs=1000.1;%采样频率
t=0:10/1000:10;&共1001个数据点,一共采了10秒
x=exp(-0.3*t).*cos(pi*t.^2+2*pi*t) ...

LZ，注意这一句？？？
我直接用楼主的代码，画出来的是这样的


我记得disp_hhs的代码关于参数这一块好像有点bug，我以前修改过但忘记是哪个位置了，现贴出来，请参考。

1. 
   
2. %DISP_HHS  display Hilbert-Huang spectrum
   
3. % DISP_HHS(im,t,inf)
   
4. % displays in a new figure the spectrum contained in matrix "im"
   
5. % (amplitudes in dB).
   
6. % inputs:  - im: image matrix (e.g., output of "toimage")
   
7. %          - t (optional): time instants (e.g., output of "toimage")
   
8. %          - inf (optional): -dynamic range in dB (wrt max)
   
9. %            default: inf = -20
   
10. %          - fs: sampling frequency
    
11. % use:  disp_hhs(im) ; disp_hhs(im,t) ; disp_hhs(im,inf)
    
12. %       disp_hhs(im,t,inf) ; disp_hhs(im,inf,fs) ; disp_hhs(im,[],fs)
    
13. %       disp_hhs(im,t,[],fs) ; disp_hhs(im,t,inf,fs)
    
14. %
    
15. %
    
16. % See also
    
17. %  emd, hhspectrum, toimage
    
18. %
    
19. % G. Rilling, last modification 3.2007
    
20. % gabriel.rilling@ens-lyon.fr
    
21. function disp_hhs(varargin)
    
22. error(nargchk(1,4,nargin));
    
23. fs = 0;
    
24. inf = -20;
    
25. im = varargin{1};   
    
26. t = 1:size(im,2);   
    
27. switch nargin
    
28.   case 1
    
29.     %raf
    
30.   case 2
    
31.     if isscalar(varargin{2})   
    
32.       inf = varargin{2};        
    
33.     else   
    
34.       t = varargin{2};         
    
35.     end
    
36.   case 3
    
37.     if isvector(varargin{2})   
    
38.       t = varargin{2};
    
39.       inf = varargin{3};
    
40.     else
    
41.       inf = varargin{2};
    
42.       fs = varargin{3};
    
43.     end
    
44.   case 4
    
45.     t = varargin{2};
    
46.     inf = varargin{3};
    
47.     fs = varargin{4};
    
48. end
    
49. if isempty(inf)
    
50.   inf = -20;
    
51. end
    
52. if inf > 0
    
53.   inf = -inf;
    
54. elseif inf == 0
    
55.   error('inf must be nonzero')
    
56. end
    
57. M=max(max(im));
    
58. warning off
    
59. im = 10*log10(im/M);
    
60. warning on
    
61. figure
    
62. if fs == 0
    
63.   imagesc(t,[0,0.5],im,[inf,0]);
    
64.   ylabel('normalized frequency')
    
65. else
    
66.   t=t/fs;
    
67.   imagesc(t,[0,0.5*fs],im,[inf,0]);
    
68.   ylabel('frequency')
    
69. end
    
70. set(gca,'YDir','normal')
    
71. xlabel('time')
    
72. title('Hilbert-Huang spectrum')
    
73. 
    

复制代码

是否应该是这样的？

[ 本帖最后由 bcyangbc 于 2010-4-20 16:02 编辑 ]

franciscowu

原来代码里error(nargchk(1,4,nargin))处，是 error(nargchk(1,3,nargin)).但楼主输入的是三参数，应该此处没有影响吧

bcyangbc

是的，可以用3个或4个参数，对结果没有影响。我想这里也就是3楼出现错误的原因

[ 本帖最后由 bcyangbc 于 2010-4-20 20:20 编辑 ]

qsq123

else
  t=t/fs;
  imagesc(t,[0,0.5*fs],im,[inf,0]);
  ylabel('frequency')

此处 t=t/fs; 将时间轴还原， 楼主说的问题通过这个应该就解决了。

xiangyu537

想对IMF1和IMF2做Hilbert谱分析
程序如下  感觉不对 求高人指点
imf = emd(ss);
emd_visu(ss,imf)
for j=1:2
[A,f,tt]=hhspectrum(imf(j,:));
end
axis([0 1024 0 0.1])
[im,tt]=toimage(A,f);
disp_hhs(im,[],fs);
axis([0 512 0 6000])
set(gca,'YDir','normal')
colormap(flipud(hot));

dzkt

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);
%axis([0 512 0 6000])
%set(gca,'YDir','normal')
%colormap(flipud(hot));
end