声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2152|回复: 2

[分形与混沌] 关于Lyapunov指数计算程序的问题

[复制链接]
发表于 2011-10-28 17:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
小弟初学混沌,有些问题请教各位。最近看到一段小数据量法计算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;
[tt1,y1]=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;

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2011-11-6 16:01 | 显示全部楼层
delt_t1 = 0.001是步长,通常取得很小。
发表于 2012-7-1 09:48 | 显示全部楼层
我的理解是在用龙格库塔计算混沌方程时的步长,离散系统应该取1,这在论坛里也有前辈说过
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-10 23:47 , Processed in 0.076719 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表