No 3-6 参考阅读：北大学生作业 Lorenz.docVIP

数值代数上机作业第六次 关于Lorenz问题与混沌的简单讨论 谭昌汇北京大学数学科学学院科学与工程计算系，100871 完成于2007年6月6日 摘 要 (Abstract) 常微分方程的数值算法有很多中，诸如Euler算法，以及Runge-Kutta算 法等等。用这些方法可以求解一些常微分方程的数值解。 本文使用四级四阶古典Runge-Kutta公式对Lorenz方程 ( dx dt = σ(y ? x) ( ( ( ( dy dt = ρx ? y ? xz ( ( ( ( dz dt = xy ? βz 进行数值求解，并分析解曲线的性质，从而对“混沌”现象有一个初步的 认识和分析。同时也是对常微分方程的数值方法的具体应用。 关键词：Runge-Kutta方法，Lorenz方程，混沌，常微分方程数值方法，数 值分析 1 介绍 (Introduction) 1.1 常微分方程的数值方法 (Numerical method for di?erential equations) 常微分方程的求解是数学上的一大重要问题，多年来，数学家都在为寻求 求解常微分方程的普遍的方法而努力。然而，众所周知的是，一些常微分方程 是没有初等形式的积分解的，这对求解常微分方程来说是一个重大的打击。 在工程应用中，我们并不像纯数学家那样追求结论的精确，往往一个足够 令人满意的近似解就可以起到很好的作用。为此，计算数学学家提出了很多常 微分方程的数值方法，诸如Guass方法、Runge-Kutta方法，线性多步法，以及 辛几何的算法等等。这些方法都很好的对常微分方程的积分解做了较为成功的 数值估计，也在工程上有着广泛的应用。 1.2 Lorenz 沌 (Lorenz problem and chaos) 问题与混 Lorenz方程是一个很简单的三维常微分方程。其微分形式如下： ( ( ( ( ( dx dt = σ(y ? x) ( dy dt = ρx ? y ? xz ( ( ( ( dz ( dt = xy ? βz 这个看似十分普通的方程其实非常的麻烦，其大致的三维图像如下图（具 体求作方法见下面的讨论，这里只是一个概述） 1FIGURE 3D graphic of the equation (N=5000, y(0)=(1,1,1),h=0.01) 50 40 30 20 10 0 40 20 20 0 0 ?20 ?40 ?20 三维图形可能显得挤了一点，我们来看一下的是其在xOy平面上的投影， 其大致形状如下图。 FIGURE Projection on xOy of the equation (N=15000, y(0)=(1,1,1),h=0.01) 30 20 10 0 ?10 ?20 ?30 ?20 ?15 ?10 ?5 0 5 10 15 20 可以发现，这个投影由两个“耳朵”组成。函数不停地在两个“耳朵”的 地方缠绕，作近似周期性的运动。 通过后面的分析，可以知道，当参数特定选取时，初值稍有变化，函数的 2绕法就发生变化。也就是说，微分系统对初值的敏感度相当的大，这就形成了 所谓的“混沌”。 混沌现象在实际生活中是相当常见和有用的。比如大气情况和天气预报就 与混沌相关。因此，对混沌的研究是相当必要的。这里只对混沌现象做出一个 简略的引入。 2 Runge-Kutta 解Lorenz 方法求 问题 (Solving Lorenz problem using Runge-Kutta method) 2.1 Runge-Kutta 法 (Runge-Kutta method) 方 Runge-Kutta方法是解决常微分方程数值解问题的一个重要的方法。其想法 其实很简单，就是通过待定系数法确定一些权因子，使得最终的迭代可以由一 下的方法得出（具体论证可见参考文献[1]） yn +1 = yn + h(c1K1 + · · · K1 = f(xn, yn), + cmkm), K2 = f(xn + a2h, yn + hb21K1), . . . . . . . . . . . . Xm?1 Km = f(xn + amh, yn + h bm1K1). i=1 通过计算可以选取适当的参数使得估计的阶数尽量高，从而尽可能提升算法的 精度。 2.2 Lorenz 解 (Numerical result of Lorenz problem) 问题的数值 首先由上面的待定系数法可以求出一个四级四阶的古典Runge-Kutta公式、 ( yn+1 = yn + h6 (K1 + 2K2 + 2K3 + K4) ( ( ( ( ( ( K1 = f(xn, yn) ( ( ( K2