cc.m

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

clear
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[tt1,y1]=ode45(@LorenzDifEqn1,[0,75],[0.1,0.1,0.1]);
xx1 = y1(:,1)';
delt_t1 = 0.01;
t1 = 0:delt_t1:75;
x = spline(tt1, xx1, t1);
Data = x(1000:1:3999);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sigma=std(Data);
R=[sigma/2 sigma sigma*3/2 2*sigma];
LengthD=length(Data);
% Data是长度为N的时间序列。
t_max=200;
[delt_ave_s,ave_s,S_cor]=DAS(LengthD,Data,t_max,R);
figure,plot(1:200,delt_ave_s);title('delt ave s');
figure,plot(1:200,ave_s);title('ave s');
figure,plot(1:200,S_cor);title('S cor');
%subplot(3,1,1);plot(1:200,delt_ave_s);title('delt_ave_s');
%subplot(3,1,2);plot(1:200,ave_s);title('ave_s');
%subplot(3,1,3);plot(1:200,S_cor);title('S_cor');