求助：disp_hhs.m

fgsunny

经过研究 ，我的程序终于没有错误了，但是边际谱出来了，可是HHT谱却没有显示。（只有运行出了，图谱版但是上面什么也没有：没有坐标轴，也没有坐标）
请求高手指点，能够简要的帮我解释一下disp_hhs.m是什么原理画谱的，
function disp_hhs(varargin)
error(nargchk(1,3,nargin));
fs = 0;
inf = -20;
im = varargin{1};
t = 1:size(im,2);
switch nargin
  case 1
    %raf
  case 2
    if isscalar(varargin{2})
      inf = varargin{2};
    else
      t = varargin{2};
    end
  case 3
    if isvector(varargin{2})
      t = varargin{2};
      inf = varargin{3};
    else
      inf = varargin{2};
      fs = varargin{3};
    end
  case 4
    t = varargin{2};
    inf = varargin{3};
    fs = varargin{4};
end
if isempty(inf)
  inf = -20;
end
if inf > 0
  inf = -inf;
elseif inf == 0
  error('inf must be nonzero')
end
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]);
  ylabel('normalized frequency')
else
  imagesc(t/fs,[0,0.5*fs],im,[inf,0]);
  ylabel('frequency')
end
set(gca,'YDir','normal')
xlabel('time')
title('Hilbert-Huang spectrum')
研究了一下，可是发现基本没看明白：:@(
小妹在此万分感谢！:@)

fgsunny

和实际品频率和采样点个数有关吗？
我分析的实验测得的数据，在读取数据仿真是设定其
         Ns=length(y2);  %采样点的长度 （y2为采样数据）
        fs=1000;                  % Sampling frequenc
        t=0:(1/fs):((Ns-1)/fs);   % time series
        freq=[1:Ns/2]*fs/Ns;      % frequency axis



请求论坛的朋友们给予帮助，谢谢了:@)