我正在研究Lyapunov指数,并且Rosenstein算法可用于计算最大指数。我在python中使用nolds包。我也知道单独的Lyapunov指数不能用来确定系统是混乱的还是周期性的。以下是我使用的以下示例代码:
import numpy as np
import nolds
rwalk = np.cumsum(np.random.random(1000))
nolds.lyap_r(rwalk)
这给出了以下输出:
array([5.94316837e+00, 9.23215127e-06])
为什么有两个不同的数字?我已阅读文档,但我不明白为什么输出是这样的。
答案 0 :(得分:-3)
我认为您正在估计Lyapunov指数的整个范围。要估计最大的Lyapunov指数,可以使用nolds,lyap_e()