求最大李雅普诺夫指数的matlab程序

程序一function dx=Lorenz(t,x);dx(1,1)=10*(x(2)-x(1));dx(2,1)=x(1)*(30-x(3))-x(2);dx(3,1)=x(1)*x(2)-8

程序一 function dx=Lorenz(t,x); dx(1,1)=10*(x(2)-x(1)); dx(2,1)=x(1)*(30-x(3))-x(2); dx(3,1)=x(1)*x(2)-8/3*x(3); dx(4,1)=0; dx(5,1)=0; dx(6,1)=0; function lambda_1=lyapunov_wolf1(data,N,m,tau,P) %Lyapunov --Wolf 该函数用来计算时间序列的最大指数方法 %m: 嵌入维数 %tau: 时间延迟 %data: 时间序列 %N: 时间序列长度 %P:,I 时间序列的平均周期选择演化相点距当前点的位置差，即若当前相点为，则演化相点 |IJ|>P 只能在－的相点中搜寻 %lambda_1:lyapunov 返回最大指数值 %************************************************************************** %ode 计算整数阶系统的时间序列 %****************************************************************** delt_t1 =0.001; t1 =0:delt_t1:60; [tt1,y1]=ode45(@lorenz,t1,[-1,0,1]); xx1 =y1(:,1)'; x1 =spline(tt1, xx1, t1); data= x1(20000:10:60000);% 采样 N=length(data); m=3; tau=11; %*****************************************************