Lyapunov指数谱计算程序

基于Lyapunov指数谱计算工具箱
工具箱下载网站为 http://www.math.rsu.ru/mexmat/kvm/matds/。这里不在赘述。下载工具箱后，只要根据简化Lorenz系统编写下列函数。

function f=SimLorenz_ly(t,X)
global c
x=X(1); y=X(2); z=X(3);
Y= [X(4), X(7), X(10);
    X(5), X(8), X(11);
    X(6), X(9), X(12)];
f=zeros(12,1);
f(1)=10*(y-x);
f(2)=-x*z+(24-4*c)*x+c*y;
f(3)=x*y-8/3*z;
 Jac=[-10,10,0;
     24-4*c-z,  c,  -x;
      y,    x,   -8/3];
  f(4:12)=Jac*Y;

计算系统的Lyapunov指数谱，新建一个m文件，
编写代码如下:

clc
clear
global c
c=2;
[T,Res]=lyapunov(3,@SimLorenz_ly,@ode45,0,0.5,500,[ 0.1 0 0.1],5);
Ly=Res(end,:)