emd分解完画频谱图怎么老出错啊！！！帮看看
源程序：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',);
imf=emd(x);
=hhspectrum(imf);
=toimage(A,f);
fs=300;
disp_hhs(im,tt,fs);
colormap(flipud(gray));
set(gcf,'color',);
ylim()
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',);
运算的这一步怎么老错：
disp_hhs(im,tt,fs);
??? Error using ==> disp_hhs at 35
inf doit etre < 0
帮解答一下，怎么能算通 啊！！！
页:
[1]