求助：如何求三维分数阶系统的最大李指数呢？

wh_Helen

最近在看关于分数阶系统的李指数方面的资料，也参考了论坛上的一些程序，多是针对整数阶系统。如下面的就是求整数阶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
   [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')

sogooda

你的这个很专业的，还是靠自己比等别人帮忙更实际一些，毕竟这玩意儿懂的人比较少。

wh_Helen

对！ 自己努力才最实际！
但刚接触这块，还是请有这方面经验的高手，路过请给些建议或一些参考文献！

lihaitao123

从定义上好像是可行的，但是就不能再用ode45,看过的一些文献使用euler方法解出来的。。

lihaitao123

在整数阶和其他经典方法相比较而言，这个程序我算的效果不是很好。