emd分解完画频谱图怎么老出错啊！！！帮看看

bgpz2007

源程序：
clear;clc;
fs=2000;
N=1024;
t=0:1/fs:(N-1)/fs;
x=(1+0.2*sin(2*pi*6*t)).*cos(2*pi*20*t+0.5.*sin(2*pi*10*t))+0.8.*sin(2*pi*80*t);
figure;
plot(t,x);
set(gcf,'color',[1 1 1]);
imf=emd(x);
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
fs=300;
disp_hhs(im,tt,fs);
colormap(flipud(gray));
set(gcf,'color',[1 1 1]);
ylim([0,110])
L_IMF=length(imf(:,1));
figure(3);
for i=1:L_IMF
    subplot(L_IMF,1,i);plot(t,imf(i,:),'k');
end
set(gcf,'color',[1 1 1]);


运算的这一步怎么老错：
disp_hhs(im,tt,fs);
??? Error using ==> disp_hhs at 35
inf doit etre < 0

帮解答一下，怎么能算通 啊！！！