lorenz_lyapunov.m

小数据量法求Lyapunov的优化程序,把网上的源程序优化过了

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--');