非线性动力学.docxVIP

非线性动力学 非线性系统之一瞥Lorenz系统 0前言0.1非线性系统动力学 线性系统是状态变量和输出变量对于所有可能的输入变量和初始状态都满足叠加原理的系统；非线性系统就是这些量不满足叠加原理的系统。非线性系统在日常生活和自然界中不胜枚举，也远远多于线性系统。 非线性动力学是研究非线性系统的各种运动状态的定性和定量变化规律，尤其是系统的长时期行为。研究的对象主要有分叉、混沌和孤立子等。 0.2洛伦兹方程 洛伦兹方程是美国气象学家洛伦兹在模拟天气这一非周期性现象时确定，这个方程的三个变量分别模拟温度、湿度和压力。可以得出结论，初期微小的差别随着时间推移差别会越来越大，洛伦兹基于此提出长期的天气预报是不可能的。这也被视为研究非线性混沌理论的开始，所以洛伦兹系统在研究非线性系统中具有举足轻重的地位。本文借助洛伦兹系统对非线性进行简单的介绍。洛伦兹方程如下。 仲=o（-x+y）\y=-xz+|ix-y[z=xy-Pz 方程中，。、□和b都为实参数。实参不同，系统的奇点及数目也是不同的。 1奇点和稳定性1.1奇点 洛伦兹系统含有三个实参数，当参数变化，奇点的数目可能不同。首先，（0，0，0）一定是系统的奇点。P0时，当^1时，系统仅有（0，0，0）一个奇点；当口1时，系统还有另外两个奇点(±血|1-1),±血|1-1),(|1-1))。 下面仅解口1时的两个非原点奇点。令o(-x+y)=0-xz+px-y=xy-阻=0 方程第一式得厂〃第三式可得z弋弋，将两式代入第二式得x3 x3 px-x=0 即x=y=±j6(p-1)，z=#=p-1。 1.2奇点稳定性判别下面根据Liapunov稳定性判别方法，找出系统在原点处大范围渐进稳定的条件，取Liapunov函数V=(x2+oy2+oz2)/2。考虑o0，60的情况。则dV■^t=xx+oyy+ozz 将洛伦兹方程仲=o(-x+y) \y=-xz+px-y (z=xy-四 代入上式，可得oy2 oy2—oz2+(o+op)xy ---=ox2dt 变换为二次型，系数矩阵为-o1 -o 1 -(o+op)0 1-o2(o+op)一阮)(0 -o 一阮) 已知o0，60，则系数矩阵负定的条件是p1。所以该系统是大范围渐进稳定的条件是p1，前提是o0，60。 Liapunov函数V总是存在的，只要构造出合适的Liapunov函数，就可以通过Liapunov稳定性定理直接判断奇点的稳定性，而不需要求解非线性方程组。有的Liapunov函数不易构造，则可以通过奇点处导算子的特征值来判断：若所有的特征值实部都小于0，则方程组在该奇点是局部渐进稳定的；若特征值实部至少有一个为正，该奇点是不稳定的。仍以洛伦兹系统为例，求出导算子的特征值。 TOC\o1-5\h\z—oa0Df（0,0,0）=（|1—10）00—6 特征矩阵的行列式（特征方程）为入+a—a0—p入+10=（入+b）［入2+（a+1）入+a—ap_］00入+6 特征值入］=—b—a—1—Ja2+（4口—2）a+1入=22—a—1+Ja2+（4口—2）a+1｛入3=2 显然，当a0，60时，A10，入20，要使方程在原点处渐进稳定，入3必须小于0，因此a+1Ja2+（4口—2）a+1 两边同时平方可得（1—|i）a0 因此|i1 显然p1使得系统在奇点（0，0，0）渐进稳定。 1.3中心流形定理 导算子的特征根的实部都不为0，它的Liapunov稳定性可由特征值实部来判定；若导算子的特征根实部有0存在，显然不能通过Liapunov定理判断，可以借助中心流形定理判断。特征根中，零实部特征根对应的特征向量构成的子空间比较特殊，从这个解子空间出发的轨线是周期轨。 前面计算得到，特征值中入10，入20，可以使=-°T+Jb2;（叩一2）°+1=0，此时|1=1，即三个特征值分别为气=-6■入=―。一1 、入3=0 特征向量分别为e1=（0，0，1）Te2=（。，—1，0）T ■e3=（1，1，0）T 洛伦兹系统方程组不是标准形式，先将其化为标准形式，令u(v)=wx0(y)=(0za—101u1)( u (v)=w x0(y)=(0z a—10 1u 1)(v)=(0 av+ww—v)u 0 1 a+1 1 (a+1 0—1 a+1a a+1 0) x (y)=z zxy—a+1a+1 x+ay(a+1a+1) u(v)w 0 1 a+1 1 (a+1 0—1 a+1a a+1 0) x (y)z 由洛伦兹方程中Xa(—x+y) (y)=(—xz+x—y) Zxy—6z 因此 因为标准形为 0 1 O+1(寿0-1 G+1a a+1 。(￡dy) 0)xy-6z/ /x\/av+w\ (y)=(w-v)