最大lyapunov指数matlab程序

求解Chen系统的最大lyapunov指数。求解方法为定义法。两条相轨线的步长为初设距离d0的基础上加上相对分量。求解时直接运行chen_lyapunov.m即可。可移植性强，比如换求其他系统的最大lyapunov指数，只需要改变变量数目，方程即可。

%%计算lymax
% function jeffcott（）
global a;
Z=[];
d0=1e-8;
ti=0.1;
X1=[];
for a=linspace（255026）
lsum=0;
x=0;y=1;z=0;
x1=0;y1=1;z1=d0;
for i=1:1000
%    [T1Y1]=ode45（@Chen1[x;y;z;a;3;28]）; 
%    [T2Y2]=ode45（@Chen1[x1;y1;z1;a;3;28]）;
x0=[x;y;z];
x01=[x1;y1;z1];
Y1=solvechen（x0ti）;
Y2=solvechen（x01ti）;
   n1=length（Y1）;n2=length（Y2）;
   x=Y1（n11）;y=Y1（n12）;z=Y1（n13）;%    X=x;
   x1=Y2（n21）;y1=Y2（n22）;z1=Y2（n23）;
   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>500
       lsum=lsum+log（d1/d0）;
   end
end
% X1=[X1 X];
Z=[Z lsum/（i-500）];
end
a=linspace（255026）;
figure;plot（aZ）;
title（‘Chen 系统最大lyapunov指数‘）
xlabel（‘parameter a‘）ylabel（‘lyapunov exponents‘）
% figure;plot（aX1）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        826  2014-02-24 12:31  chen_lyapunov.m

     文件        371  2014-02-23 22:24  solvechen.m

----------- ---------  ---------- -----  ----

                 1197                    2