Lorenz微分方程.docVIP

Lorenz微分方程 洛仑兹 E.N.Lorenz,1917- 混沌吸引子已成为混沌理论的徽标，代表着复杂性新科学，好 比行星轨道图代表着哥白尼、开普勒理论一样。英国的《新科学家》杂志曾办了一个专栏， 每一期的刊头都有洛仑兹吸引子标记。在国内清华大学曾国屏 1953- 撰写的《自组织的自 然观》一书封面一共选了三幅图：道家的阴阳鱼、行星轨道图和洛仑兹混沌图，三张图分别 位于一个魔方 也可视为骰子 外显的三个面上。这种设计是很有趣的，作者试图用三幅图分 别代表三种世界图景，第一种是各国早期朴素的辩证自然观，第二种是近代科学形成的还原 论的自然观，第三种是以自组织理论、复杂性理论为标志的新型自然观。 以前科学家们不自觉地认为微分方程的解只有那么几类：1 发散轨道；2 不动点；3 极限环 ；4 极限环面 torus，复数形式为tori 。除此以外，大概没有新的运动类型了，这是人们 的一种主观猜测，谁也没有给出证明。事实上这种想法是非常错误的。1963年美国麻省理工 学院气象科学家洛仑兹给出一个具体模型，清楚地展示了一种新型运动体制：混沌运动，轨 道既不收敛到极限环上也不跑掉。 洛仑兹是学数学出身的，1948年起在MIT作动力气象学博士后工作，1963年他在《大气科学 杂志》上发表的论文《确定性非周期流》是混沌研究史上光辉的著作，无论怎样评价都不算 过分。[有关此文的发表经过可参见洛仑兹的专著《混沌的本质》 The Essense of Chaos 第四章，中译本由刘式达等翻译，气象出版社1997年出版。]而略带腼腆的洛仑 兹在回顾自己的成就时总是十分谦虚。大家应当注意的是，他在那个时代就使用了计算机， 如果不用计算机他肯定发现不了确定性非周期流，不可能就“系统长期行为不可预测”作出 判断。现在科学家使用计算机已成为一种时髦，而那时却是一种冒险。 说了半天，早该给出洛仑兹方程的具体形式了，学数理的朋友们最喜欢公式 但学文的正相 反 。洛仑兹方程极其美妙，多少年来，人们并未逃出洛仑兹方程的框架，找来找去能出现 混沌的系统总是与洛仑兹系统差不多。洛仑兹方程形式很简单，只有三个方程： 其中是正的参数，这是一个自治的三阶方程组。此方程所 表示的向量场的散度 为divF 偏f_1/偏x+偏f_2/偏y+偏f_3/偏 所以洛仑兹系统是耗散的，相体积不断收缩。此方程虽然能一定程度上描述天气的复杂变化 过程，但它的真正意义并不在气象预报上!他首先是数学家，他用数学来思考问题。因而他 的模型及其从中所揭示出的新的运动体制的意义，就远远不止于气象学。 现在大家都已清楚，在二维连续系统中不可能出现混沌，三维是出现混沌所要求的最低维数 。洛仑兹模型恰好只有三维。 对于洛仑兹方程，一般是固定参数σ和b，单独考察r变化时，系统 行为的变化。当0＜r＜1时，有一个稳定不动点O 0,0,0 ；当1＜r ＜r^*≡1.34561…时 对于σ 10,b 8/3 ，又出现两个新的稳定不动 点A和B，这时共有3个不动点，不动点O已变为不稳定不动点。 A和B的性质总是相同的，因为方程在变换 x,y,z → -x,-y,z 下是不变的。当r继续增大到r_t σ σ+b+3 / σ-b-1 24.7368…时 对于σ 10, b 8/3 ，方程的3个不动点都变得不稳定了，r_t是系统行为变化的临界点， 这时就出现了洛仑兹发现的“确定性非周期流”。1963年洛仑兹研究时3个参数的取值为： σ 10,b 8/3,r 28。 这组参数值通常称标准情形 canonical case 。当年洛仑兹就是在这一组参数值下，采用计算机数值计算，发现了奇怪吸引子。当时还没 有这一概念，直到1971年吕埃尔 D.Ruelle 和塔肯斯 F.Takens 才提出“strange attractor”一词。 洛仑兹的伟大贡献是多方面的，可以轻易举出几条： 1 为非线性动力学研究贡献了一个绝好的数学模型，值得推崇的不是此模型在多大程度上 直接反映了某个具体的物理现象，而是它抓住了复杂性的本质，为进一步深入探索各种复杂 事物奠定了坚实的基础。 2 洛仑兹有一套清晰的研究方法，《确定性非周期流》一文运用的一系列合理的“工作程序 ”，已被70年代末以来的混沌学家广泛采用。 3 定性的数学分析与计算机定量模拟相结合，开创了复杂性研究的新方法。在现在看来使 用计算机是件平常的事，但在60年代初，这是不平凡的。特别是，他识破了复杂系统中“对 初始条件的敏感依赖性”，能够正确地理解数值计算结果，显示了一位杰出科学家的洞察力 。 苏联著名学者安德罗诺夫 A.A.Andronov，1901-1952 一生“工作在相平面上”，对非线性 振动有深刻理解，但最终没有发现混沌，没有完成质的飞跃。关键性的一步似乎注定要让