6958192460 发表于 2011-10-28 17:20

关于Lyapunov指数计算程序的问题

小弟初学混沌,有些问题请教各位。最近看到一段小数据量法计算Lyapunov指数的程序,请问其中的delt_t1指的是什么,怎样确定?文章中从来没提到过,请大家赐教
下面是程序:
function lambda_1=lyapunov_wolf1(data,N,m,tau,P)
%该函数用来计算时间序列的最大Lyapunov 指数--Wolf 方法
%m: 嵌入维数
%tau:时间延迟
%data:时间序列
%N:时间序列长度
%P:时间序列的平均周期,选择演化相点距当前点的位置差,即若当前相点为I,则演化相点只能在|I-J|>P的相点中搜寻
%lambda_1:返回最大lyapunov指数值
%**************************************************************************
%ode计算整数阶系统的时间序列
%******************************************************************
delt_t1 = 0.001;
t1 = 0:delt_t1:60;
=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;

qzm020811 发表于 2011-11-6 16:01

delt_t1 = 0.001是步长,通常取得很小。

zoulei87 发表于 2012-7-1 09:48

我的理解是在用龙格库塔计算混沌方程时的步长,离散系统应该取1,这在论坛里也有前辈说过
页: [1]
查看完整版本: 关于Lyapunov指数计算程序的问题