python求微分方程 洛伦兹_Lorenz微分方程

Lorenz

微分方程

洛仑兹

(E.N.Lorenz,1917- )

混沌吸引子已成为混沌理论的徽标，代表着复杂性新科

学，好

比行星轨道图代表着哥白尼、开普勒理论一样。英国的《新科学家》杂志曾办了一

个专栏，

每一期的刊头都有洛仑兹吸引子标记。

在国内清华大学曾国屏

(1953-

)

撰写的

《自

组织的自

然观》一书封面一共选了三幅图：道家的阴阳鱼、行星轨道图和洛仑兹混沌图，

三张图分别

位于一个魔方

(

也可视为骰子

)

外显的三个面上。这种设计是很有趣的，作者试

图用三幅图分

别代表三种世界图景，第一种是各国早期朴素的辩证自然观，第二种是近代

科学形成的还原

论的自然观，第三种是以自组织理论、复杂性理论为标志的新型自然观。

以前科学家们不自觉地认为微分方程的解只有那么几类：

1)

发散轨道；

2)

不动点；

3)

极限环

；

4)

极限环面

(torus

，复数形式为

tori)

。除此以外，大概没有新的运动类型了，

这是人们

的一种主观猜测，谁也没有给出证明。事实上这种想法是非常错误的。

1963

年美

国麻省理工

学院气象科学家洛仑兹给出一个具体模型，清楚地展示了一种新型运动体制：

混沌运动，轨

道既不收敛到极限环上也不跑掉。

洛仑兹是学数学出身的，

1948

年起在

MIT

作动力气象学博士后工作，

1963

年他在

《大气

科学

杂志》上发表的论文《确定性非周期流》是混沌研究史上光辉的著作，无论怎样评价

都不算

过分。

[

有关此文的发表经过可参见洛仑兹的专著《混沌的本质》

(

The Essense of

Chaos

)

第四章，

中译本由刘式达等翻译，

气象出版社

1997

年出版。

]

而略带腼腆的洛仑

兹在

回顾自己的成就时总是十分谦虚。大家应当注意的是，他在那个时代就使用了计算机，

如

果不用计算机他肯定发现不了确定性非周期流，不可能就“系统长期行为不可预测”作出

判断。现在科学家使用计算机已成为一种时髦，而那时却是一种冒险。

说了半天，早该给出洛仑兹方程的具体形式了，学数理的朋友们最喜欢公式

(

但学文的

正相

反

)

。

洛仑兹方程极其美妙，多少年来，

人们并未逃出洛仑兹方程的框架，找来找去能

出现

混沌的系统总是与洛仑兹系统差不多。洛仑兹方程形式很简单，只有三个方程：

;

/

;

/

);

(

/

bz

xy

dt

dz

xz

y

rx

dt

dy

y

x

dt

dx



















其中

b

r

,

,



是正的参数，这是一个自治的三阶方程组。此方程所

表示的向量场的散度

为

div

F

=

偏

f

_1/

偏

x

+

偏

f

_2/

偏

y

+

偏

f

_3/

偏

0

)

1

(

1



















b

b

z





.