关于Hilbert边际谱的问题

Rufter2007

clear
load sj;
fs=8000;
N=8001;
a1=s(1:N,1);
imf=emd(a1);
[A,fa,tt]=hhspectrum(imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
title('Hilbert边际谱');
以上程序欲求解数据序列sj.mat的边际谱，其中，sj.mat是一段随机的1s话音序列。但是运行结果提示“??? Out of memory. Type HELP MEMORY for your options.
Error in ==> fliplr at 17
y = x(:,end:-1:1);
Error in ==> toimage at 79
im=fliplr(im)';
Error in ==> test7 at 8
[E,tt1]=toimage(A,fa,tt,length(tt));”
以此程序对0.5s的话音进行处理却没问题，跟数据序列的大小有关么？还是哪里出了问题，不知道怎么改。哪位高人点拨下，感激。

[ 本帖最后由 Rufter2007 于 2009-12-23 21:41 编辑 ]