非线性动力学之一瞥——Lorenz系统.doc

《非线性动力学之一瞥——Lorenz系统.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、非线性动力学非线性系统之一瞥——Lorenz系统2013-01-30-12-0前言0.1非线性系统动力学线性系统是状态变量和输出变量对于所有可能的输入变量和初始状态都满足叠加原理的系统；非线性系统就是这些量不满足叠加原理的系统。非线性系统在日常生活和自然界中不胜枚举，也远远多于线性系统。非线性动力学是研究非线性系统的各种运动状态的定性和定量变化规律，尤其是系统的长时期行为。研究的对象主要有分叉、混沌和孤立子等。0.2洛伦兹方程洛伦兹方程是美国气象学家洛伦兹在模拟天气这一非周期性现象时确定，这个方程的三个变量分别模拟温度

2、、湿度和压力。可以得出结论，初期微小的差别随着时间推移差别会越来越大，洛伦兹基于此提出长期的天气预报是不可能的。这也被视为研究非线性混沌理论的开始，所以洛伦兹系统在研究非线性系统中具有举足轻重的地位。本文借助洛伦兹系统对非线性进行简单的介绍。洛伦兹方程如下。x=σ-x+yy=-xz+μx-yz=xy-βz方程中，σ、μ和b都为实参数。实参不同，系统的奇点及数目也是不同的。1奇点和稳定性1.1奇点洛伦兹系统含有三个实参数，当参数变化，奇点的数目可能不同。首先，0，0，0一定是系统的奇点。β>0时，当μ≤1时，系统仅有0，

3、0，0一个奇-12-点；当μ>1时，系统还有另外两个奇点±βμ-1，±βμ-1，μ-1。下面仅解μ>1时的两个非原点奇点。令σ-x+y=0-xz+μx-y=0xy-βz=0方程第一式得x=y，第三式可得z=xyβ=x2β，将两式代入第二式得-x3β+μx-x=0即x=y=±βμ-1，z=x2β=μ-1。1.2奇点稳定性判别下面根据Liapunov稳定性判别方法，找出系统在原点处大范围渐进稳定的条件，取Liapunov函数V=x2+σy2+σz2/2。考虑σ>0，β>0的情况。则有dVdt=xx+σyy+σzz将洛伦兹方

4、程x=σ-x+yy=-xz+μx-yz=xy-βz代入上式，可得dVdt=-σx2-σy2-βσz2+σ+σμxy变换为二次型，系数矩阵为-σ12σ+σμ012σ+σμ-σ000-βσ已知σ>0，β>0，则系数矩阵负定的条件是μ<1。所以该系统是大范围渐进稳定的条件是μ<1，前提是σ>0，β>0。Liapunov函数V总是存在的，只要构造出合适的Liapunov函数，就可以通过Liapunov-12-稳定性定理直接判断奇点的稳定性，而不需要求解非线性方程组。有的Liapunov函数不易构造，则可以通过奇点处导算子的特征

5、值来判断：若所有的特征值实部都小于0，则方程组在该奇点是局部渐进稳定的；若特征值实部至少有一个为正，该奇点是不稳定的。仍以洛伦兹系统为例，求出导算子的特征值。Df0,0,0=-σσ0μ-1000-β特征矩阵的行列式（特征方程）为λ+σ-σ0-μλ+1000λ+β=λ+bλ2+σ+1λ+σ-σμ特征值λ1=-bλ2=-σ-1-σ2+4μ-2σ+12λ3=-σ-1+σ2+4μ-2σ+12显然，当σ>0，β>0时，λ1<0，λ2<0，要使方程在原点处渐进稳定，λ3必须小于0，因此σ+1>σ2+4μ-2σ+1两边同时平方可得1

6、-μσ>0因此μ<1显然μ<1使得系统在奇点0，0，0渐进稳定。1.3中心流形定理导算子的特征根的实部都不为0，它的Liapunov稳定性可由特征值实部来判定；若导算子的特征根实部有0存在，显然不能通过Liapunov定理判断，可以借助中心流形定理判断。特征根中，零实部特征根对应的特征向量构成的子空间比较特殊，从这个解子空间出发的轨线是周期轨。-12-前面计算得到，特征值中λ1<0，λ2<0，可以使λ3=-σ-1+σ2+4μ-2σ+12=0，此时μ=1，即三个特征值分别为λ1=-βλ2=-σ-1λ3=0特征向量分别为e

7、1=0，0，1Te2=σ，-1，0Te3=1，1，0T洛伦兹系统方程组不是标准形式，先将其化为标准形式，令xyz=0σ10-11100uvw=σv+ww-vu即uvw=0011σ+1-1σ+101σ+1σσ+10xyz=zxσ+1-yσ+1xσ+1+σyσ+1uvw=0011σ+1-1σ+101σ+1σσ+10xyz由洛伦兹方程中xyz=σ-x+y-xz+x-yxy-βz因此uvw=0011σ+1-1σ+101σ+1σσ+10σ-x+y-xz+x-yxy-βz因为标准形为xyz=σv+ww-vu因此-12-xyz=σ-

8、x+y-xz+x-yxy-βz=-1+σσv-σv+wu+σ+1vσv+ww-v-βu代入前一式可得uvw=σv+ww-v-βu1σ+1σv+wu-1+σv-σuσ+1σv+w若只判断奇点的稳定性，只需计算其中心流形Wlocc0，0，0Wlocc0，0，0=u，v，w

9、v=φw，u=ψw，φ0=φ0=Ψ0=Ψ0=0v=φw，u=ψ