Lyapunov指数计算[共20页].pdfVIP

Lyapunov 指数的计算与仿真 摘要：Lyapunov 指数在研究动力系统的分岔、混沌运动特性中起着重要的作用，是衡量系 统动力学特性的一个重要定量指标，它表征了系统在相空间中相邻轨道间收敛或发散的平均 指数率。对于系统是否存在动力学混沌，可以从最大Lyapunov 指数是否大于零非常直观的 判断出来：一个正的Lyapunov 指数意味着在相空间中，无论初始两条轨线的间距多么小， 其差别都会随时间的演化而成为指数率的增加以致达到无法预测，这就是混沌现象。本文介 绍了Lyapunov 指数的定义、性质及计算原理和数值计算的实现方法，应用Matlab 软件、编 写了计算Lyapunov 指数程序，计算了Logistic 映射系统、Henon 系统及时间序列的Lyapunov 指数。实例的计算机仿真表明Lyapunov 指数是研究分岔、混沌运动，解决工程实际问题的 有效方法之一。 关键词：非线性动力系统；混沌；Lyapunov 指数；时间序列 1. Lyapunov 指数的定义及性质 1.1 Lyapunov 指数的定义 对于一维离散映射系统xn+1 F (xn ) ，在每次迭代过程中，初始两点是相互 分离还是相互靠拢由下式决定： dF 1，映射使两点分开 dx dF 1，映射使两点靠拢 dx dF 但是在不断地迭代过程中， 的值要不断变化。为了表示从整体看两相邻初始 dx 状态分离的情况，必须对时间（或迭代次数）取平均。为此，设平均每次迭代所 引起的指数分离中的指数为λ，则初始两点 与 相距 ，经n 次迭代后，ε x x +ε 0 0 n n nλ 点变为 与 ，相距 ，即： F (x0 ) F (x0 +ε) εe εenλ(x0 ) = n ( ) - n ( ) （1） F x +ε F x 0 0 取极限ε →0 ，n →∞，由上式得到： F n (x ε) F n (x ) 1 + − 0 0 λ(x 0 ) lim lim ln n→∞ →ε 0 n ε （2 ） 1 dF n (x) lim ln