李雅普诺夫指数

Posted studyer_domi

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了李雅普诺夫指数相关的知识,希望对你有一定的参考价值。

1. 李雅普诺夫指数

Lyaponuv(李雅普诺夫)指数表示相空间相邻轨迹的平均指数发散率的数值特征。又称李雅普诺夫特征指数,是用于识别混沌运动若干数值的特征之一。

李雅普诺夫指数常常被用来判定一个系统的混沌性,通过图像可以直观地看出某个系统或者映射是否是混沌系统或映射。

2. Lyaponuv指数性质

我们假定Lyaponuv指数用 来表示,那么对于的不同取值有以下不同的含义:

当lambda >0时,系统运动会进入混沌状态,对应的映射叫做混沌映射;

当 lambda<0时,系统的运动状态会趋于稳定且此时对系统的初始状态不敏感,也就是说此时的映射对初始值不敏感;

当 lambda=0时,系统则处于稳定状态。

3. Lyaponuv指数的计算


上式中lambda就是所要计算的LE氏指数。

Python实现

 


import numpy as np
import matplotlib.pyplot as plt
# show plots in notebook
% matplotlib inline

result = []
lambdas = []
maps = []

# define range of r 
rvalues = arange(0, 2, 0.01)

# loop through r
for r in rvalues:
    x = 0.1
    result = []
    # iterate system 100 times
    for t in range(100):
        x = x + r - x**2
        # calculate log of the absolute of the derivative
        result.append(log(abs(1 - 2*x)))
    # take average
    lambdas.append(mean(result))
    # for the map ignore first 100 iterations as transient time and iterate anew
    for t in range(20):
        x = x + r - x**2
        maps.append(x)    
    
fig = plt.figure(figsize=(10,7))
ax1 = fig.add_subplot(1,1,1)

xticks = np.linspace(0, 2, 4000)
# zero line
zero = [0]*4000
ax1.plot(xticks, zero, 'g-')
# plot map
ax1.plot(xticks, maps, 'r.',alpha = 0.3, label = 'Map')
ax1.set_xlabel('r')
# plot lyapunov
ax1.plot(rvalues, lambdas, 'b-', linewidth = 3, label = 'Lyapunov exponent')
ax1.grid('on')
ax1.set_xlabel('r')
ax1.legend(loc='best')
ax1.set_title('Map of x(t+1) = x(t) + r - x(t)^2 versus Lyapunov exponent')

 

 

以上是关于李雅普诺夫指数的主要内容,如果未能解决你的问题,请参考以下文章

Matlab求解混沌系统最大李雅普诺夫指数

matlab洛伦兹混沌系统时间序列李雅普指数计算

matlab洛伦兹混沌系统时间序列李雅普指数计算

回顾:系统的能控性、能观性和稳定性及李雅普诺夫方法

分数阶混沌系统李雅普指数和分岔图

数值计算 --求解连续微分系统和混沌系统