v2-第9章-常微分方程数值解-wb.docxVIP

PAGE PAGE # PAGE PAGE # 第九章常微分方程数值解法 在工程和科学技术的实际问题中，常需要求解常微分方程，许多实际问题的数学模型都是常微分方程 或常微分方程的定解问题，如物体运动、电路振荡、化学反映及生物群体的变化等。例如：对于单摆问题 (如图9.1)，选取摆长I的函数，来描述单摆运动。到下面二阶常微分方程1和质量m 1，我们希望得到摆角q的关于时间 (如图9.1)，选取摆长I 的函数，来描述单摆运动。 到下面二阶常微分方程 1和质量m 1，我们希望得到摆角q的关于时间t 根据力矩与角加速度的关系, 经一定的简化后可得 d2q dt2 gsinq 其中g为重力加速度。当单摆在摆动开始时刻 t t0时的初始摆角q t0 qo 和初始角速度凹 dt t to q (t0) q0都确定时, 单摆运动规律q(t)惟一确定, 这就可以写成一个初值问题: d2 d2q dF gsinq q to qo, q (to) q。 虽然常微分方程广泛存在于众多研究和应用领域，但是常微分方程中往往只有少数较简单和典型的微 分方程(例如线性常系数常微分方程等)的解能用初等函数、特殊函数或它们的级数与积分表达 ，对于变系 数常微分方程的解析求解就比较困难，而一般的非线性常微分方程的求解困难就更不用说了。大多数情况 下，常微分方程只能用近似方法求解。这种近似解法可分为两大类：一类是近似解析法，如级数解法、逐 次逼近法等；另一类是数值解法，它给出方程在一些离散点上的近似解。 在具体求解常微分方程时，需具备某种定解条件，常微分方程和定解条件合在一起组成定解问题。定 解条件有两种：一种是给出积分曲线在初始点的状态，称为初始条件，相应的定解问题称为初值问题；另 一类是给出积分曲线首尾两端的状态，称为边界条件，相应的定解问题称为边值问题。本章主要介绍常微 分方程初值问题的基本数值方法、理论和算法。 9.1基本概念 9.1.1常微分方程初值问题 常微分方程可分为线性、非线性、高阶方程与方程组等类，其中线性方程包含于非线性类中，高阶方 程可化为一阶方程组，若方程组中的所有未知量看作一个向量，则方程组可写成向量形式的单个方程。因 此研究一阶常微分方程的初值问题dyf(x, y) 此研究一阶常微分方程的初值问题 dy f(x, y)， dx y(a) yo 的数值解法具有典型性，其中方程的解为 y : R R。 a x b (9-1) 只有保证问题(9-1)的解存在惟一的前提下，其数值解法的研究才有意义。由常微分方程的基本理论, 我们有：定理9.1如果(9-1)中的f(x,y)满足条件f (x, y)在区域 D {(x, y) a x b， y }上连续；f (x, y)在D上关于y满足Lipschitz条件，即存在常数 只有保证问题(9-1)的解存在惟一的前提下，其数值解法的研究才有意义。由常微分方程的基本理论, 我们有： 定理9.1如果(9-1)中的f(x,y)满足条件 f (x, y)在区域 D {(x, y) a x b， y }上连续； f (x, y)在D上关于y满足Lipschitz条件，即存在常数 L，使得 f (x, y) f (x, y) Ly y| 则初值问题(9-1)在区间［a,b］上存在惟一的连续解 y y(x)。 在本章的讨论中，我们总假定方程满足以上两个条件，从而方程总存在惟一的连续解。 在实际问题的研究中，许多问题最后可以归结为一阶常微分方程组的初值问题 Yi fi (x, Y1 , Y2 丄，Ym) (i 1 ,L , m) Yi(a) Yio 若记 Y ( Y1, Y2,L , Ym)T , Yo (Y10, Y20,L , ymo)T, f (f1, f2,L , fm)T，则初值问题(9-2)可写成如下向量 (9-2) 形式的单个方程 如果向量函数f(x,y)在区域D : a x b, y 使得对 x ［a,b］, Y1, Y2 Rm，都有 I f (x,yj 那么问题(9-3)在［a,b］是存在惟一解y二y(x)。 y f (x,y) y(a) y。 Rm连续，且关于y满足Lipschitz条件，即存在L 0， f(x,y2) L||y1 y2， 问题(9-3)与问题(9-1)形式上完全相同，故对初值问题 (9-1)所建立的各种数值解法都可应用于求解问题 (9-3)，只需将y换成向量y，f (x,y)换成f (x, y)即可。 对于高价常微分方程的初值问题的处理， 下：设有m阶常微分方程初值问题 可以通过变量代换化为一阶常微分方程组初值问题，过程如 (9-3) 引入新变量y1 y,y2 y ,l ,ym (m) y y(a) y(m y1 y2 y2 y3 M (m 1)、 f (x, Y