以下是在本论坛找到的计算李雅普诺夫指数的程序,我对matlab不懂,请问:1、果将这个程序用于计算一个新的时间序列,需要修改程序中的哪几个参数?
2、其中linear=input('请输入线形部分的长度')是指什么的长度,如何确定?
谢谢
clear
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[tt1,y1]=ode45(@LorenzDifEqn1,[0,200],[0.1,0.1,0.1]);
xx1 = y1(:,1)';
delt_t1 = 0.01;
t1 = 0:delt_t1:75;
x1 = spline(tt1, xx1, t1);
Data = x1(1000:1:6000);
m=8;
tau=14;
N=5001;
P=ave_T(Data,N);
[Y_num,Y]=small_data_sets_Lyapunov1(Data,N,m,tau,P,delt_t1);
XX=1:Y_num;
plot(XX,Y(XX),'.','markersize',1);hold on;
linear=input('请输入线形部分的长度');
XX1=1:linear;
pp=polyfit(XX1,Y(XX1),1);
lambda_1=pp(1)
yp=polyval(pp,XX1);
plot(XX1,yp,'r--');