非线性动力学之一瞥_Lorenz系统

非线性动力学
非线性系统之一瞥——Lorenz系统
2013-01-30
0 前言
0.1非线性系统动力学
线性系统是状态变量和输出变量对于所有可能的输入变量和初始状态都满足叠加原理的系统；非线性系统就是这些量不满足叠加原理的系统。

非线性系统在日常生活和自然界中不胜枚举，也远远多于线性系统。

非线性动力学是研究非线性系统的各种运动状态的定性和定量变化规律，尤其是系统的长时期行为。

研究的对象主要有分叉、混沌和孤立子等。

0.2洛伦兹方程
洛伦兹方程是美国气象学家洛伦兹在模拟天气这一非周期性现象时确定，这个方程的三个变量分别模拟温度、湿度和压力。

可以得出结论，初期微小的差别随着时间推移差别会越来越大，洛伦兹基于此提出长期的天气预报是不可能的。

这也被视为研究非线性混沌理论的开始，所以洛伦兹系统在研究非线性系统中具有举足轻重的地位。

本文借助洛伦兹系统对非线性进行简单的介绍。

洛伦兹方程如下。

方程中，、和都为实参数。

实参不同，系统的奇点及数目也是不同的。

1 奇点和稳定性
1.1 奇点
洛伦兹系统含有三个实参数，当参数变化，奇点的数目可能不同。

首先，一定是系统的奇点。

时，当时，系统仅有一个奇点；当时，系统还有另外两个奇点。

下面仅解时的两个非原点奇点。

令
方程第一式得，第三式可得，将两式代入第二式得
即，。

1.2 奇点稳定性判别
下面根据Liapunov稳定性判别方法，找出系统在原点处大围渐进稳定的条件，取Liapunov函数。

考虑，的情况。

则有
将洛伦兹方程
代入上式，可得
变换为二次型，系数矩阵为
已知，，则系数矩阵负定的条件是。

所以该系统是大围渐进稳定的条件是，前提是，。

Liapunov函数V总是存在的，只要构造出合适的Liapunov函数，就可以通过Liapunov稳定性定理直接判断奇点的稳定性，而不需要求解非线性方程组。

有的Liapunov函数不易构造，则可以通过奇点处导算子的特征值来判断：若所有的特征值实部都小于0，则方程组在该奇点是局部渐进稳定的；若特征值实部至少有一个为正，该奇点是不稳定的。

仍以洛伦兹系统为例，求出导算子的特征值。

特征矩阵的行列式（特征方程）为
特征值
显然，当，时，，，要使方程在原点处渐进稳定，必须小于0，因此
两边同时平方可得
因此
显然使得系统在奇点渐进稳定。

1.3中心流形定理
导算子的特征根的实部都不为0，它的Liapunov稳定性可由特征值实部来判定；若导算子的特征根实部有0存在，显然不能通过Liapunov定理判断，可以借助中心流形定理判断。

特征根中，零实部特征根对应的特征向量构成的子空间比较特殊，从这个解子空间出发的轨线是周期轨。

前面计算得到，特征值中，，可以使
，此时，即三个特征值分别为
特征向量分别为
洛伦兹系统方程组不是标准形式，先将其化为标准形式，令
即
由洛伦兹方程中
因此
因为标准形为
因此
代入前一式可得
若只判断奇点的稳定性，只需计算其中心流形
代入变换后的方程组
下面用近似的方法求中心流形，即分别将和表示成幂级数的形式
由条件可以使幂级数化为
代入第一式可得
展开可得
比较和的系数可得
因此
且
得
因此
因此中心流上的解满足
因此
1）：
时，当时；时，当时。

因此时，奇点是稳定的。

2）：
时，当时，远离原点；时，当
时，远离原点。

因此时，奇点是不稳定的。

图 1.1 Lorenz系统中心流形
2 分叉
2.1 概念
分叉是指在含有参数的系统中，参数值连续地变化到某一值时系统的定性或定量行为发生了质的变化。

这种变化包括奇点的数目和稳定性变化、闭轨数目的变化以及同宿轨、异宿轨的出现或消失。

2.2 洛伦兹系统的奇点分叉
通过前面的讨论可以发现，洛伦兹系统随参数的变化奇点的数目和奇点稳定性将发生改变，这就是奇点分叉的实例。

（1）奇点数目的改变：
前面计算奇点数目时发现，奇点的数目与无关，而与和有关。

若，时有三个奇点和；时就有一个奇点。

若，恰与之相反。

若，也只能有一个奇点。

（2）奇点稳定性的变化
因为有三个参数，奇点不唯一且变化，所以讨论起来比较麻烦。

下面以奇点为例分析。

特征值为
显然，从负变为正时，从正变为负，奇点一定从不稳定变为稳定。

但是时，出现了零特征根的情形，在这一点是否稳定需要通过中心流形来判断，方法同前面讨论时的一样。

（3）hopf分叉
Hopf分叉是指导算子的特征值沿复平面的上方或下方（即不是通过实轴）穿过虚轴。

下面试讨论洛伦兹系统在处的Hopf分叉。

对于，显然是实数，变化经过实轴穿过虚轴，显然不是Hopf分叉。

只有和有可能不经过实轴。

不经过实轴时，，实部的变化与有关。

令，得，时实部大于0。

图 2.1 的正负
分情况讨论（如图 2.1）：
时，要求
时，要求
时，，因此不存在Hopf分叉。

图 2.2 分叉类型3 混沌吸引子
用数值方法计算出的lorenz吸引子
图 3.1 ，，
图 3.2 和（红色），，
图 3.2 ，（红色）和（黄色），
，
C语言代码
#include "stdafx.h"
#include "stdafx.h"
#define N 4001
void main()
{
double sigma,mu,beta,delta;
double x[N],y[N],z[N];
int i;
sigma=10;
mu=28;
beta=8/3;
delta=0.01;
x[0]=1;
y[0]=0;
z[0]=0;
for(i=1;i<=N-1;i++)
{
x[i]=x[i-1]+delta*sigma*(-x[i-1]+y[i-1]);
y[i]=y[i-1]+delta*(-x[i-1]*z[i-1]-y[i-1]+mu*x[i-1]);
z[i]=z[i-1]+delta*(x[i-1]*y[i-1]-beta*z[i-1]);
}
FILE*out=fopen("lorenz.csv","w");
for(i=0;i<=N-1;i++)
{
fprintf(out,"%f,",x[i]);
}
fprintf(out,"\n");
for(i=0;i<=N-1;i++)
{
fprintf(out,"%f,",y[i]);
}
fprintf(out,"\n");
for(i=0;i<=N-1;i++)
{
fprintf(out,"%f,",z[i]);
}
fprintf(out,"\n");
fclose(out);
}。