Rufter2007 发表于 2009-12-23 21:28

关于Hilbert边际谱的问题

clear
load sj;
fs=8000;
N=8001;
a1=s(1:N,1);
imf=emd(a1);
=hhspectrum(imf);
=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
=toimage(A,fa,tt,length(tt));”
以此程序对0.5s的话音进行处理却没问题,跟数据序列的大小有关么?还是哪里出了问题,不知道怎么改。哪位高人点拨下,感激。

[ 本帖最后由 Rufter2007 于 2009-12-23 21:41 编辑 ]
页: [1]
查看完整版本: 关于Hilbert边际谱的问题