HHT问题请教

chunhuajia

HHT高手们：请帮我解释一下下面程序出错的原因，我找了半天还是找不出来，谢谢了！
[x,fs]=wavread('c1.wav')；
x=x(130000:1:149600);
t=1:length(x);
imf=emd(x,t);%模态分解
[A,fa,tt]=hhspectrum(imf);
[im,tt1]=toimage(A,fa,tt,length(tt));
运行出错：??? Out of memory. Type HELP MEMORY for your options.

Error in ==> toimage at 65
im(splx,sply) = 0;
可以帮我说说toimage.m中下面这段的具体意思吗？为什么im的初始位置要这样设置？
for i=1:size(f,1)
  for j = 1:lt
    ff=floor(f(i,j)*2*(sply-1))+1;
    if ff <= sply % in case f(i,j) > 0.5
      im(floor(j*(splx-1)/lt)+1,ff)=im(floor(j*(splx-1)/lt)+1,ff)+A(i,j);
    end
  end  
end

cboboc

嗯。。顶一下，我也想知道