关于HHT频谱的问题

lk677976

Length=1000;
x=zeros(1,Length);
n1=200;
x1=0.8;
n=Length+n1;
for i=1:n
x2=3.7*x1*(1-x1);
if(i>n1)
    j=i-n1;
    x(j)=x2;
end
x1=x2;
end
t=1:length(x);
[imf,ort,nbits]=emd(x,t);
ximf=imf(:,:);
savefilename='ximf';
save(savefilename,'ximf');
[A,f,tt] = hhspectrum(ximf);
[a,b]=xcorr(x,'unbiased');
plot(b,a)这个显示的是对整个IMF的HHT频谱吗，还是对单阶的IMF变换啊，不是很清楚啊，那个toimage函数又是用来干嘛的呢