用定义求最大Lyapunov指数遇到的问题?
在论坛上找到如下程序计算最大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
=ode45('Chen',1,);
=ode45('Chen',1,);
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=;
end
a=linspace(32,40,80);
plot(a,Z,'-');
title('Chen 系统最大lyapunov指数')
xlabel('parameter a'),ylabel('lyapunov exponents')
[ 本帖最后由 yibanlixue 于 2008-5-21 20:52 编辑 ] 没人回复啊，在下真的很着急，请oct，liliangbiao等高手指点啊，非常感谢啊 有人明白么，我也想请教
页:
[1]