常微分方程线性多步法研究.pptVIP

求解常微分方程初值问题线性多步法研究 陈婧 2009/6/15 信息与计算科学 南京工业大学浦江学院 目录 第一章 引言与综述 第二章 常微分方程初值问题的数值解法 第三章 基于数值积分的构造方法 第四章 基于Taylor展开的待定系数方法  第五章 线性多步法的算例 第六章 总结 第一章 引言与综述 现在，常微分方程作为微分方程的基本类型之一，在很多学科领域内有着重要的应用，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等。很多问题的数学表述都可以归结为常微分方程的定解问题，实际生活中很多问题的数学模型都是微分方程。这些问题都可以化为求常微分方程的解，或者化为研究解的性质的问题。 但在许多情况下，首先找到问题的解析解，然后再进行相关的计算往往非常困难，有时甚至是行不通的，基于此理由，我们可以避免求解析解而直接求相应的精度较高的数值解。在常微分方程的级数解法中，逐步逼近法就是近似解法。这些方法可以给出解的近似表达式，通常称为近似解析方法。还有一类近似方法称为数值方法，它可以给出解在一些离散点上的近似值。在解常微分方程时，我们经常使用的就是数值解法。 关于数值解法主要有如下几种方法： 第二章 常微分方程初值问题的数值解法 通常常微分方程初值问题的数值方法可以分为两类: （1）单步法——计算 在 处的值仅取决于 处的应变量及其导数值. （2）多步法——计算 在 处的值需要应变量及其导数在 之前的多个节点的值. 单步法： 1.欧拉法 若将函数 在点 处的导数 用两点式代 替，即： ，再用 近似值 代替,则 初值问题变为： 上式就是著名的欧拉(Euler)公式.以上方法称为欧拉 法或欧拉折线法. 2.龙格—库塔公式 由差商 和微分中值定理 并利用微分方程 得： 其中 称作区间 上的平均斜 率，记作 。 3. 线性多步法 求解初值问题的单步法在计算时只用到前面一步的结果，所以要提高精度时，需要增加中间函数值的计算，这就加大了计算量。如果在计算 时，不仅用到上 的近似值 ，还用到前面若干节点 上的近似值 称这种方法为多步法。 线性多步法的基本思想：在逐步推进的求解过程中， 计算 之前事实上已经求出了一系列的近似值 如果充分利用前面多步的信息来预测，则 可期望获得 较高的精度，这就是构造多步法的基本思想。 构造线性多步法有两种途径。一种是基于数值积分 的方法，对此我们着力介绍两种方法，一种是亚当斯显 式（四阶显式公式），还有一种是亚当斯隐式公式（四 阶隐式公式）。第二种是基于Taylor展 开的方法。 第四章 基于Taylor展开的待定系数方法 设想要构造如下形式的线性 步公式： 其中 和 为待定常数。将其局部截断误差在点 处作Taylor展开，得： 由上述公式我们可以推出： 2步隐式公式Simpson公式： 4步显式公式Milne公式： 3步隐式公式Hamming公式： 第五章 线性多步法的算例 在本论文中我们研究了两个算例，对于这 两个算例我们分别采用了单步法中的龙格——库塔 四阶方法和线性多步法中的Adams四阶显、隐式公式、 十阶外插公式、Milne公式来进行计算。 这两个算例分别为： 两个实验所得到的误差结果图 我们对上述两个算例采用了不同的方法求出解 与其精确解进行比较。 通过对数据进行比较我们可以发现以下结论： 1、对于同样的方法来说取不同的步长所产生的效果 是不一样的。对于同样的方法，所取步长越小的得 到的值就越精确，迭代的值的误差就越稳定，产生 的效果也就越好。但是计算量会增加，运行速度也 会降低。 2、对于相同的步长来说，选取多步法的阶数的不 同，产生的效果也是不同的。对于相同的步长来说 ，阶数越高的多步法所得出的值就越精确，效果越 好。 3、而对于相同步长，相同阶数来说，隐式公式的效 果要比显式公式