关于emd的问题，请帮忙！

zzw1983107

请各位大侠帮忙：
我在运行emd得到IMF分量后，求边际谱时出现了下面的问题：

??? Out of memory. Type HELP MEMORY for your options.

Error in ==> D:\MATLAB6p5\toolbox\matlab\elmat\fliplr.m
On line 16  ==> y = x(:,n:-1:1);

Error in ==> F:\calculate program\hilbert-huang\toimage.m
On line 77  ==> im=fliplr(im)';

Error in ==> F:\calculate program\hilbert-huang\BJP.m
On line 76  ==> [E,tt1]=toimage(A,fa,tt,length(tt));

请大家帮忙看看！多谢了！

zhangnan3509

请把把程序贴出来，这样有问题也好找:handshake

eight

原帖由 zzw1983107 于 2008-5-6 08:49 发表
请各位大侠帮忙：
我在运行emd得到IMF分量后，求边际谱时出现了下面的问题：

??? Out of memory. Type HELP MEMORY for your options.

Error in ==> D:\MATLAB6p5\toolbox\matlab\elmat\fliplr.m
On line 16 ...

估计你的数据量太大，内存太少了。这个错误估计搜索论坛就找到如何解决了

zzw1983107

首先谢谢楼上的两位，我的数据有6000个数，我把虚拟内存改大了，也还是不行啊。
请问你们在求边际谱的时候最大处理多少个数？

jinnian

原帖由 zzw1983107 于 2008-5-6 12:17 发表
首先谢谢楼上的两位，我的数据有6000个数，我把虚拟内存改大了，也还是不行啊。
请问你们在求边际谱的时候最大处理多少个数？

6000个数不算多啊，我试了一下，速度还挺快的。把程序贴上来看看，还有最好把matlab升级到7.1以上。

[ 本帖最后由 jinnian 于 2008-5-6 13:05 编辑 ]

zzw1983107

[imf,ort,nbits] = emd(x); %imf中包含了各个IMF和趋势R分量
n=size(imf,1)
ort=io(x,imf)  %计算各分量的正交性
figure(2)
subplot(n+1,1,1)
plot(x);
for i=1:n
        subplot(n+1,1,i+1);
    plot(imf(i,:));
end
figure(3)
N=length(imf(1,:));
y=imf; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
  for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;
end
ft=(0:N-3)/N*(fs/2);
plot(ft,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
end
请高手给帮忙看看吧，电脑是1G内存，x为6000个数，得到的E应该是6000*6000的矩阵。算边际谱的时候就出现out of memory的错误了.
还有一个有的程序在    bjp(k)=sum(E(k,:))*1/fs*1/T;  那怎么理解？   谢谢各位了！

jinnian

原帖由 zzw1983107 于 2008-5-6 13:44 发表
= emd(x); %imf中包含了各个IMF和趋势R分量
n=size(imf,1)
ort=io(x,imf)  %计算各分量的正交性
figure(2)
subplot(n+1,1,1)
plot(x);
for i=1:n
        subplot(n+1,1,i+1);
    plot(imf(i,:));
end
figure(3 ...

不知楼主用的emd包是什么时候的版本，我的新版的是不能这样调用的：[E,tt1]=toimage(A,fa,tt,length(tt))。正确的调用例子如下： [im,tt,ff] = toimage(A,f);[im,tt] = toimage(A,f,t);[im,tt,ff] = toimage(A,f,sply);[im,tt,ff] = toimage(A,f,splx,sply);[im,tt,ff] = toimage(A,f,t,splx,sply);当然新版的也用不着E=flipud(E)这一句，但这一句肯定不会引起内存溢出。6000点的数据量真的不算不多。我算过至少100000点的数据。

[ 本帖最后由 jinnian 于 2008-5-6 14:56 编辑 ]

zzw1983107

我的数据量太大了，因为这个没法算，挺苦恼的！

jinnian

好的，我发一个给你！等一下吧！……已发，请查收！matlab可能要用7.1以上版本。再说一句，你的数据量真的不算大。我的机器配置也不高的。

[ 本帖最后由 jinnian 于 2008-5-6 15:34 编辑 ]

gordylau

用clear清楚WORKSPACE里的变量，用完一个清一个（指比数据量比较大的变量），可以缓解这种情况。6000的数据量好像真的不大，要是6000*6000可就够大的了