常微分方程数值解法.doc

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

1、第八章常微分方程数值解法教学目的1.掌握解常微分方程的单步法：Euler方法、Taylor方法和Runge-Kutta方法；2.掌握解常微分方程的多步法：Adams步法、Simpson方法和Milne方法等；3.了解单步法的收敛性、相容性与稳定性；多步法的稳定性。教学重点及难点重点是解常微分方程的单步法：Euler方法、Taylor方法和Runge-Kutta方法和解常微分方程的多步法：Adams步法、Simpson方法和Milne方法等；难点是理解单步法的收敛性、相容性与稳定性及多步法的稳定性。教学时数20学时教学过程§1基本概念1．1常微分方程初值问题的一般提

2、法常微分方程初值问题的一般提法是求函数，满足其中是已知函数，是已知值。假设在区域上满足条件：（1）在上连续；（2）在上关于变量满足Lipschitz条件：，（1.3）其中常数称为Lipschitz常数。我们简称条件（1）、（2）的基本条件。由常微分方程的基本理论，我们有：定理1当在上满足基本条件时，一阶常微分方程初值问题（1.1）、（1.2）对任意给定存在唯一解在上连续可微。定义1方程（1.1）、（1.2）的解称为适定的，若存在常数和，对任意满足条件及的和，常微分方程初值问题（1.4）存在唯一解，且适定问题的解连续依赖于（1.1）右端的和初值。由常微分方程的基本理

3、论，还有：定理2当在上满足基本条件时，微分方程（1.1）、（1.2）的解是适定的。我们在本章中假设在上满足基本条件，从而（1.1）、（1.2）的解存在且适定。一般的一阶常微分方程组初值问题是求解（1.5）（1．5）的向量形式是（1.5）其中记。类似于定理1和定理2，我们有：定理3若映射满足条件(1)在上是从到上的连续映射;(2)在上关于满足Lipschits条件;任意。则常微分方程组初值问题（1.5）存在的唯一的连续可微解而且解是适定的。高阶常微分方程初值问题一般为（1.6）其中是给定多元函数，为给定值。引进新的变量函数（1.7.）则初值问题（1.6）化成了一阶常

4、微分方程组初值问题通过求解（1.8）得到(1.6)的解。1.2初值问题数值解基本概念初值问题的数值解法，是通过微分方程离散化而给出解在某些节点上的近似值。在上引入节点称为步长。在多数情况下，采用等步长，即。记（1.1）,(1.2)的为准确解为,记的近似值为,记为.。求值问题数值解的方法是步进法，即在计算出后计算。数值的方法有单步与单步法之分。单步法在计算时只利用而多步法在计算时不仅要利用还要利用前面已算出的若干个。我们称要用到的多步法为步方法。单步法可以看作多步法，但两者有很大差别。步方法只能用于的计算，要用其它的方法计算；而且在稳定性上单性法比的多步法容易分析；

5、此外单步法容易改变步长。单步法和多步法又都有显式方法和稳式方法之分。单步显式法的计算公式可写成（1.9）隐式单步法的计算公式可写成（1.10）在(1.10)中右端项显含。从而（1.10）是的方程式，要通过解方程求出。显式多步法计算公式为（1.11）而隐式多步法计算公式为(1.12)右端项含。多步法中一类常用方法是线性多步法（1.13）其中是独立于和的常数。时（1.13）是显式的，时是隐式的。数值解法及到方法构造、误差分析、稳定性分析等内容。一些概念和定义在后面的论述中逐步引入。§2Euler的方法Euler方法是常微分方程初值问题数值方法中最简单的方法。Euler

6、方法精度低，较少有直接使用。但我们通过Euler方法介绍离散化途径、数值解法中的基本概念、术语和加速方法等。2.1显式Euler方法设节点为。初值问题（1.1）、（1.2）的显式Euler方法为（2.1）其中。显式Euler方法可以用多种途径导出。导出方法1：Taylor展开法。将在点进行Taylor展开,得（2.2）忽略这一阶项,分别用近似，和，得。结合初值条件即得（2.1）。导出方法2：向前差分近似微分法。用向前差分近似微分，得（2.3）将近似号改作等号，用近似，，并结合初值条件即得（2.1）。方法3：左矩数值积分法。将（1.1）两边从到积分得（2.4）用，近

7、似、，数值积分采用左矩公式得，从而亦得（2.1）Euler方法有几何意义，如图8-1，式（1.1）,(1.2)的解曲线过点，且具斜率。从出发以为斜率作直线段，交于，显然。式（1.1）过的解曲线具有斜率从出发以为斜率作直线要交于，余类推。这样我们得到了一条折线，它在点的右侧具有斜率，与（1.1）过的解曲线相切。我们取折线，作为（1.1）、（1.2）解曲线的近似曲线，所以Euler方法又称折线法。2.2隐式Euler方法和梯形方法若将在展开、忽略项，用和分别近似，及，可以得另一计算公式（2.5）（2．5）式称为隐式Enler方法。隐式Euler方法也可以利用向后差分近

8、似微分或用