数值分析第9章常微分方程数值解法

《数值分析第9章常微分方程数值解法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第十章常微分方程数值解法（NumericalMethodsforOrdinaryDifferentialEquations）问题驱动：蝴蝶效应洛伦兹吸引子(Lorenzattractor)是由MIT大学的气象学家EdwardLorenz在1963年给出的，他给出第一个混沌现象——蝴蝶效应。图10.1.1蝴蝶效应示意图洛伦兹方程是大气流体动力学模型的一个简化的常微分方程组：该方程组来源于模拟大气对流，该模型除了在天气预报中有显著的应用之外，还可以用于研究空气污染和全球侯变化。洛伦兹借助于这个模型，将大气流体运动的强度x与水平和垂直方向的温

2、度变化y和z联系了起来。参数称为普兰特数，是规范化的瑞利数，和几何形状相关。洛伦兹方程是非线性方程组，无法求出解析解，必须使用数值方法求解上述微分方程组。洛伦兹用数值解绘制结果图10.1.1，并发现了混沌现象。§1引言微分方程数值解一般可分为：常微分方程数值解和偏微分方程数值解。自然界与工程技术中的许多现象，其数学表达式可归结为常微分方程（组）的定解问题。一些偏微分方程问题也可以转化为常微分方程问题来（近似）求解。Newton最早采用数学方法研究二体问题，其中需要求解的运动方程就是常微分方程。许多著名的数学家，如Bernoulli（家族

3、），Euler、Gauss、Lagrange和Laplace等，都遵循历史传统，研究重要的力学问题的数学模型，在这些问题中，许多是常微分方程的求解。作为科学史上的一段佳话，海王星的发现就是通过对常微分方程的近似计算得到的。本章主要介绍常微分方程数值解的若干方法。一、初值问题的数值解法1、一阶常微分方程初值问题的一般形式常微分方程的数值解法分为（1）初值问题的数值解法（2）边值问题的数值解法(2)一般构造方法：离散点函数值集合+线性组合结构→近似公式2.迭代格式的构造(1)构造思想：将连续的微分方程及初值条件离散为线性方程组加以求解。由于

4、离散化的出发点不同，产生出各种不同的数值方法。基本方法有：有限差分法（数值微分）、有限体积法（数值积分）、有限元法（函数插值）等等。(3)如何保证迭代公式的稳定性与收敛性?3.微分方程的数值解法需要解决的主要问题(1)如何将微分方程离散化，并建立求其数值解的迭代公式？(2)如何估计迭代公式的局部截断误差与整体误差？称在区域D上对满足Lipschitz条件是指:记4、相关定义二、初值问题解的存在唯一性考虑一阶常微分方程的初值问题/*Initial-ValueProblem*/:则上述IVP存在唯一解。只要在上连续,且关于y满足Lipsc

5、hitz条件，即存在与无关的常数L使对任意定义在上的都成立，求函数y(x)在一系列节点a=x0

6、x1亦称为欧拉折线法/*Euler’spolygonalarcmethod*/在假设yi=y(xi)，即第i步计算是精确的前提下，考虑的截断误差Ri=y(xi+1)yi+1称为局部截断误差/*localtruncationerror*/。定义若某算法的局部截断误差为O(hp+1)，则称该算法有p阶精度。定义欧拉法的局部截断误差：Ri的主项/*leadingterm*/欧拉法具有1阶精度。例1:用欧拉公式求解初值问题取步长。解:应用Euler公式于题给初值问题的具体形式为：其中。计算结果列于下表：可用来检验近似解的准确程度。进行计算，

7、数值解已达到了一定的精度。这个初值问题的准确解为，从上表最后一列，我们看到取步长欧拉公式的改进：隐式欧拉法/*implicitEulermethod*/向后差商近似导数x0x1))(,()(1101xyxfhyxy+由于未知数yi+1同时出现在等式的两边，不能直接得到，故称为隐式/*implicit*/欧拉公式，而前者称为显式/*explicit*/欧拉公式。一般先用显式计算一个初值，再迭代求解。隐式欧拉法的局部截断误差：即隐式欧拉公式具有1阶精度。梯形公式/*trapezoidformula*/—显、隐式两种算法的平均注：的

8、确有局部截断误差，即梯形公式具有2阶精度，比欧拉方法有了进步。但注意到该公式是隐式公式，计算时不得不用到迭代法，其迭代收敛性与欧拉公式相似。中点欧拉公式/*midpointformula*/中心差商近似导