zzw1983107 发表于 2008-5-6 08:49

关于emd的问题,请帮忙!

请各位大侠帮忙:
我在运行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==> =toimage(A,fa,tt,length(tt));

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

zhangnan3509 发表于 2008-5-6 09:09

回复 楼主 的帖子

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

eight 发表于 2008-5-6 09:49

本帖最后由 wdhd 于 2016-9-10 14:42 编辑

原帖由 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 发表于 2008-5-6 12:17

回复 3楼 的帖子

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

jinnian 发表于 2008-5-6 13:02

本帖最后由 wdhd 于 2016-9-10 14:42 编辑

原帖由 zzw1983107 于 2008-5-6 12:17 发表
首先谢谢楼上的两位,我的数据有6000个数,我把虚拟内存改大了,也还是不行啊。
请问你们在求边际谱的时候最大处理多少个数?
6000个数不算多啊,我试了一下,速度还挺快的。把程序贴上来看看,还有最好把matlab升级到7.1以上。

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

zzw1983107 发表于 2008-5-6 13:44

回复 5楼 的帖子

= 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集
%%%%%%%%%%%%%求边际谱
=hhspectrum(y);
=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 发表于 2008-5-6 14:47

本帖最后由 wdhd 于 2016-9-10 14:42 编辑

原帖由 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包是什么时候的版本,我的新版的是不能这样调用的:=toimage(A,fa,tt,length(tt))。正确的调用例子如下: = toimage(A,f); = toimage(A,f,t); = toimage(A,f,sply); = toimage(A,f,splx,sply); = toimage(A,f,t,splx,sply);当然新版的也用不着E=flipud(E)这一句,但这一句肯定不会引起内存溢出。6000点的数据量真的不算不多。我算过至少100000点的数据。

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

zzw1983107 发表于 2008-5-6 15:15

回复 7楼 的帖子

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

jinnian 发表于 2008-5-6 15:28

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

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

gordylau 发表于 2008-6-5 20:41

用clear清楚WORKSPACE里的变量,用完一个清一个(指比数据量比较大的变量),可以缓解这种情况。6000的数据量好像真的不大,要是6000*6000可就够大的了
页: [1]
查看完整版本: 关于emd的问题,请帮忙!