文档介绍：中北大学
《数值分析》
常微分方程初值问题的数值解法
专业:
班级:
学号:
姓名:
日期:
常微分方程初值问题的数值解法
摘要
微分方程的数值解法在科学技术及生产实践等多方面应用广泛. 文章分析了构造常微分方程初值问题数值解法的三种常用基本方法,差商代替导数法,数值积分法及待定系数法,推导出了Euler系列公式及三阶龙格-库塔公式,指出了各公式的优劣性及适用条件,并对Euler公式的收敛性、稳定性进行了分析。
Abstract
The numerical solution of differential equations is widely used in science, technology, production practices and many other fields. This paper analyzed three kinds of basic methods for constructing numerical solutions for initial value problem of ordinary differential equations:difference quotient instead of derivative method, numerical integral method and undetermined coefficients method. At the same time, the paper deduces the Euler series formula and the classical third order Runge-Kutta formula. In addition, the paper pointed out the advantages and disadvantages of each formula and application condition, it also analyzed the convergence and stability of the Euler formula.

科学技术及实际生产实践中的许多问题都可归结为微分方程的求解问题,使用较多的是常微分方程初值问题的求解。对于一阶常微分方程的初值问题,其中f为已知函数,是初始值。如果函数f关于变量y满足Lipschitz条件,则初值问题有唯一解。只有当f是一些特殊类型的函数时,才能求出问题的解析解,但一般情况下都满足不了生产实践与科学技术发展的需要,因此通常求其数值解法。

数值解法是一种离散化的方法,可以求出函数的精确解在自变量一系列离散点处的近似值。基本思想是离散化,首先要将连续区间离散化,对连续区域进行剖分,为步长;其次将其函离散
化,基本方法有差商代替导数法、积分插值方法、待定系数法;最后研究其稳定性、收敛性质。

分别用一阶向前、向后、中心差商近似代替在处的导数,可以求出Euler系列公式.

用一阶向前差商代替导数,即
则
因为,公式可以简记为
即为显式Euler公式。Euler公式是最简单的一种数值解法,一阶的,精度较差,可直接求解,有明显的几何意义,也称为折线法,但其方法对于更复杂的情况有着较为普遍的意义。

用一阶向后差商代替导数,即
则
简记为
即为隐式Euler公式或后退Euler公式,计算比显式麻烦,但稳定性好。

二阶中心差商代替导数
则
即为两步Euler公式,稳定性较单步好。

在区间上求积分,,对右端积分采用积分插值公式,即得,为梯形公式,也可看成显式Euler公式与隐式Euler公式的算术平均,可以证明梯形方法是二阶方法。
实际计算时,初始近似,可由Euler法求解,则梯形法每步完整的计算公式为
又称为改进的Euler公式。由Euler公式给出预测值,再由梯形进行校正,收敛速度较快,一般只需求出两次迭代即可满足要求。

同样,对在区间上求积分,
并对右端采用Gauss公式,得
也称为两步Euler公式。
对上式右端积分采用Simpson公式
即
称为隐式两步法。

R-K基本思想即用位于上的若干个点处的线性组合来近似它,将三阶R-K公式一般形式表示为:
是未知常数,我们的主要目的要目的是估算和F3,,可以从下列正交逼近中得出
此正交公式可以使多项式的精度达到2阶,待定系