Lorenz系统最大Lyapunov指数程序程序通用其他系统

Lorenz系统最大Lyapunov指数程序（程序通用其他系统）

%%计算lymax
% function jeffcott（）
global a;
Z=[];
d0=1e-8;
ti=0.1;
X1=[];
for a=linspace（85012）
lsum=0;
x=0;y=1;z=0;
x1=0;y1=1;z1=d0;
for i=1:1000
%    [T1Y1]=ode45（@Lorenz1[x;y;z;a;3;28]）; 
%    [T2Y2]=ode45（@Lorenz1[x1;y1;z1;a;3;28]）;
x0=[x;y;z];
x01=[x1;y1;z1];
Y1=solveLorenz（x0ti）;
Y2=solveLorenz（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（85012）;
figure;plot（aZ）;
title（‘Lorenz 系统最大lyapunov指数‘）
xlabel（‘parameter a‘）ylabel（‘lyapunov exponents‘）
% figure;plot（aX1）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-12-17 09:19  Lyapunov_max_Lorenz\
     文件         834  2015-12-17 09:18  Lyapunov_max_Lorenz\Lorenz_lyapunov.m
     文件         373  2015-12-17 09:16  Lyapunov_max_Lorenz\solveLorenz.m