常微分方程的数值法解.docVIP

第六章 常微分方程的数值解法 §6.0 引言 §6.1 算法构造的主要途径 §6.2 Runge-Kutta Method算法 §6.3 线性多步法 §6.4 线性多步法的一般形式 §6.5 算法的稳定性、收敛性 §6.0 引 言 1． 主要考虑如下的一阶常微分方程初值问题的求解: 微分方程的解就是求一个函数y=y(x)，使得该函数满足微分方程并且符合初值条件。 例如微分方程: xy-2y=4x ；初始条件: y(1)=-3。 于是可得一阶常微分方程的初始问题 。 显然函数y(x)=x2-4x满足以上条件,因而是该初始问题 的微分方程的解。 但是,只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解，有的甚至无法用解析表达式来表示。因此，只能依赖于数值方法去获得微分方程的数值解。 4． 微分方程的数值解：设微分方程问题的解y(x)的存在区间是[a,b]，初始点x0=a，将[a,b]进行划分得一系列节点x0 , x1 ,...,xn，其中a= x0 x1… xn =b。 y(x)的解析表达式不容易得到或根本无法得到，我们用数值方法求得y(x)在每个节点xk的近似值y(xk)，即y≈y(xk)，这样y0 , y1 ,...,yn称为微分方程的数值解。如图所示: §6.1 算法构造的主要途径 1 欧拉公式 1