lorenz_lyapunov.m

来自「用小数据量方法求lyapunov指数」· M 代码 · 共 22 行

clear
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[tt1,y1]=ode45(@LorenzDifEqn1,[0,75],[0.1,0.1,0.1]);
xx1 = y1(:,1)';
delt_t1 = 0.001;
t1 = 0:delt_t1:75;
x1 = spline(tt1, xx1, t1);
Data = xx1(1000:1:1400);
m=10;
tau=14;
N=401;
P=ave_T(Data,N);
[Y_num,Y]=small_data_sets_Lyapunov(Data,N,m,tau,P,delt_t1);
for i=1:Y_num
    YY(i)=Y(i);
end
XX=1:Y_num;
pp=polyfit(XX,Y(XX),1);
lambda_1=pp(1)
yp=polyval(pp,XX);
plot(XX,Y(XX),'-o',XX,yp,'--');