疑问：关于CC 方法计算时间延迟和嵌入维数

cqupenghao · 发表于 2011-8-23 21:10

回复 59 # tangyz_hust 的帖子

图片可能小了点，是我的失误，上面显示 a=b=4，
初始条件为 x（0）=0.1
Lyapunov指数值为 LE1=ln（2）=0.693.

tangyz_hust · 发表于 2011-8-24 11:05

cqupenghao 发表于 2011-8-23 21:07

登录/注册后可看大图

回复 59 # tangyz_hust 的帖子

给你看一个截图，是用LET工具箱计算的logistic的Lyapunov指数值

   谢谢你的回复！

   我想问的问题是：为什么网上公布的程序计算结果，与理论上的结果很不一致？有什么经验技巧没有？是不是程序有问题？
   比如Logistic映射：x(i+1)=4*x(i)*(1-x(i))，理论上它的LE=ln(2), 理论时延tau = 1; 理论嵌入维m = 2！我想看看那个计算Lyapunov指数（wolf 方法）的程序是否有问题？就做了如下验证：
先迭代方法产生10000个Logistic映射数据，再使用如下入口参数，代入程序计算，结果却与理论结果ln(2)相差很远！

这是怎么回事？盼高人指点一下？

附：入口参数

tau = 1;             % 时延
m = 2;                % 嵌入维
taumax = 30;          % 最大离散步进时间
P = 100;             % 序列平均周期
fs = 1;                % 采样频率

计算Lyapunov指数（wolf 方法）程序就不贴了？网上都有！随便搜搜就可以得到！

cqupenghao · 发表于 2011-8-24 16:07

回复 62 # tangyz_hust 的帖子

你考虑到一个问题没，理论时延是tau=1，理论嵌入维是m=2，那么这是针对多少个数据来的呢？？而你是用的迭代方法产生10000个数据，那么这个序列的实际时延tau和实际嵌入维m以及序列平均周期又各是多少呢？我个人觉得你可以就把这10000个数据组成的时间序列当成一个带球Lyapunov指数的时间序列，然后你按部就班的求它的Lyapunov指数就可以了，不要总想着它是logistic映射数据，你觉得呢？

tangyz_hust · 发表于 2011-8-26 12:04

谢谢你！我要去看看书消化一下！

wangcan149 · 发表于 2012-8-12 11:13

我下载了一个陆振波的工具箱里面有CC 法但是看不怎么懂啊。比如说，他的原始数据在哪里导入，也没有。它就直接出来了一个结果，好像不能理解。
注明：我是初学者。

冰雪荷花 · 发表于 2013-4-5 21:28

winterdij 发表于 2007-9-4 16:31

登录/注册后可看大图

我用Logistic和Lorenz做了，都正确，程序也是我从Chaos Toolbox Ver.2.0工具箱里找的，我认为没有问题。

您好！看了您的帖子，请问您能否把验算程序给我发一份，不胜感激！！！！！

ABC2014 · 发表于 2014-4-6 14:13

求你的程序一用

ABC2014 · 发表于 2014-4-6 14:14

求你的程序一用

gghhjj · 发表于 2014-4-8 15:10

ABC2014 发表于 2014-4-6 14:14
求你的程序一用

http://forum.vibunion.com/forum. ... 071&fromuid=222

账号		自动登录	找回密码
密码			我要加入

[分形与混沌] 疑问：关于CC 方法计算时间延迟和嵌入维数