常微分方程数值.docVIP

第五章 常微分方程数值解法 科学技术中的许多问题，在数学中往往归结为微分方程的求解问题。例如天文学中研究星体运动,空间技术中研究物体飞行等，都需要求解常微分方程初值问题。 最简单的常微分方程初值问题是 定理 如果在区域内连续，且关于满足Lipschitz条件，那么初值问题的解存在且惟一的近似值，所以总假定微分方程的解存在惟一，即初值问题中的满足定理的条件。 除特殊情形外，初值问题一般求不出解析解，即使有的能求出解析解，其函数表示式也比较复杂，计算量比较大，况且实际问题往往只要求在某一时刻解的函数值，因此，有必要研究初值问题的数值解法。 所谓初值问题的数值解法，就是寻求解在区间上的一系列点 上的近似值. 记表示相邻两个节点的间距，称为步长。 求微分方程数值解的主要问题： (1) 如何将微分方程离散化,并建立求其数值解的递推公式; (2) 递推公式的局部截断误差、数值数与精确解的误差估计； (3) 递推公式的稳定性与收敛性. §1 解常微分方程初值问题的Euler方法 Euler法是最早的解决一阶常微分方程初值问题的一种数值方法，虽然它不够精确，很少被采用，但是它在某种程度上反映了数值方法的基本思想和特征，所以我们首先介绍它。 一、Euler方法 考虑初值