求助：如何求三维分数阶系统的最大李指数呢？
最近在看关于分数阶系统的李指数方面的资料，也参考了论坛上的一些程序，多是针对整数阶系统。如下面的就是求整数阶Chen 系统的最大李指数，能不能将这个程序修改，以求分数阶系统的最大李指数呢？或者有没有其他的方法求分数阶系统的李指数呢？希望高人指点迷津！:@) 谢谢！
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') 你的这个很专业的，还是靠自己比等别人帮忙更实际一些，毕竟这玩意儿懂的人比较少。
回复 沙发 sogooda 的帖子
对！ 自己努力才最实际！但刚接触这块，还是请有这方面经验的高手，路过请给些建议或一些参考文献！ 从定义上好像是可行的，但是就不能再用ode45,看过的一些文献使用euler方法解出来的。。 在整数阶和其他经典方法相比较而言，这个程序我算的效果不是很好。
页:
[1]