用定义求最大Lyapunov指数遇到的问题?

yibanlixue

在论坛上找到如下程序计算最大le指数，对这个程序有几点不明，请高手指点
1，d0可以设置成别的吗?其中z1=1+d0可以变成中x1=1+d0吗？在ode45的括号中的1表示积分时间是1秒吗？为什么设置成1呢?
本人初学最大Lyapunov指数，问题或很简单，请高手指点迷津！多谢了
Z=[];
d0=1e-8;
for a=linspace(32,40,80)
lsum=0;
x=1;y=1;z=1;
x1=1;y1=1;z1=1+d0;
for i=1:10000
   [T1,Y1]=ode45('Chen',1,[x;y;z;a;3;28]);
   [T2,Y2]=ode45('Chen',1,[x1;y1;z1;a;3;28]);
   n1=length(Y1);n2=length(Y2);
   x=Y1(n1,1);y=Y1(n1,2);z=Y1(n1,3);
   x1=Y2(n2,1);y1=Y2(n2,2);z1=Y2(n2,3);
   d1=sqrt((x-x1)^2+(y-y1)^2+(z-z1)^2);
   x1=x+(d0/d1)*(x1-x);
   y1=y+(d0/d1)*(y1-y);
   z1=z+(d0/d1)*(z1-z);
   if i>5000
       lsum=lsum+log(d1/d0);
   end
end
Z=[Z lsum/(i-5000)];
end
a=linspace(32,40,80);
plot(a,Z,'-');
title('Chen 系统最大lyapunov指数')
xlabel('parameter a'),ylabel('lyapunov exponents')

[ 本帖最后由 yibanlixue 于 2008-5-21 20:52 编辑 ]

yibanlixue

没人回复啊，在下真的很着急，请oct，liliangbiao等高手指点啊，非常感谢啊

yizhe

有人明白么，我也想请教