请教相空间重构matlab程序

jyxhappy

在网上看了一篇论文：孙海云，曹庆杰.相空间重构延迟时间的确定.山东工业大学学报.2000.30（2），有个例子方程是这样的

dx/dt=A-(B+1)*x+x^2*y+a*cos(x*t);

dy/dt=B*x-x^2*y;

然后用matlab用最基本的方法和相空间重构的方法分布进行了相空间绘制，我把基本的相图作法程序贴出，但重构的那种方法我一直做不出来请高手指教

Brusselator.m文件为

function dy=Brusselator(t,y);

dy=zeros(2,1);

dy(1)=0.4-2.2*y(1)+y(1)^2*y(2)+0.08*cos(0.86*t);

dy(2)=1.2*y(1)-y(1)^2*y(2);

运行下面代码

% [T,Y]=ode45('Brusselator',[0 1000],[1 0]);

% plot(Y(:,1),Y(:,2))

请问重构的程序怎么做，我想画出如文章中的图

[ 本帖最后由 jyxhappy 于 2008-7-12 14:47 编辑 ]

无水1324

不明白，你这个是确定的系统，怎么还需要重构呢？

jyxhappy

这篇文章的用意就是用重构的方法做出的相图和确定的系统做出的相图形状基本相同，以证明相空间重构方法的正确。文章只是这么写，但程序也没带，具体怎么重构的方法也没提。
我从网上看到一个很流行的代码，我把他也贴上来了。压缩文件中有两个文件，LorenzData.mat和PhaSpaRecon.m
运行代码
load LorenzData.mat%将LorenzData.mat中的数据读入
[xn,dn]=PhaSpaRecon(x,6,3)%进行重构
plot(xn(1,:),dn)%绘制相图
即可绘出相图。我想是不是上面我说得那个问题是他通过确定的方程求解出x，y的值，然后也存成一组数据，用我上传的代码的方法来做呢。

[ 本帖最后由 jyxhappy 于 2008-7-13 16:20 编辑 ]

douermi2004

看到你的相空间重构的程序,我有问题请教:为什么len-1-(m-1)*tau,而不是len-(m-1)*tau,我觉得后者是对的,不知道你的想法是怎样的,能加我QQ:745724181,共同讨论吗?

s06092062

地板说的len-（m-1）*tau，我认为是对的，楼主你的程序改改看:victory:

李锦鹏2008268

求你了，给我发一份，急用，谢谢
421069750@163.com

yuling

1. function Data=reconstitution(data,m,tau)
   
2. %该函数用来重构相空间
   
3. % m:嵌入空间维数
   
4. % tau:时间延迟
   
5. % data:输入时间序列
   
6. % Data:输出,是m*n维矩阵
   
7. 
   
8. N=length(data); % N为时间序列长度
   
9. M=N-(m-1)*tau; %相空间中点的个数
   
10. Data=zeros(m,M);
    
11. for j=1:M
    
12.   for i=1:m           %相空间重构
    
13.     Data(i,j)=data((i-1)*tau+j);
    
14.   end
    
15. end

复制代码

lvtaomei

重构图像怎么做呢？

yuling

一般系统的嵌入维数为十几甚至几十，这样通过重构获得的相空间维数较高，而我们只能作出二维和三维的相图，这就要求我们在保留原始数据主要信息的前提下，将原有的高维数据投影到较低维的数据空间，K-L变换正是进行这项操作的理想工具

mathtang

回复 yuling 的帖子

请问下：知道延迟时间，给出一组数据，如何画二维吸引子图，用相空间重构的方法，。谢谢了 ！！能否给我个源程序，多谢了！我邮箱mathtang@163.com。。。急用！！多谢了 啊！！

mathtang

回复 mathtang 的帖子

哦，谢谢了 帖子发了，还没有回

rosserneri

回复 3 # jyxhappy 的帖子

牛叉！！！！！！！！！

欧罗巴领主

遇到了相似的问题。。。T,T

mphzzia

哈哈哈哈哈哈哈

mphzzia

能不能下载连接