常微分方程初值问题数值解法的讨论与比较.docVIP

《微分方程数值解法》论文 常微分方程初值问题数值解法的讨论与比较 一、一阶常微分方程的初值问题 科学技术中常常要求解常微分方程的定解问题，这类问题最简单的形式是一阶方程的初值问题 我们知道，只要函数适当光滑，譬如关于满足Lipschitz条件 （1.3） 理论上就可以保证初值问题(1.1),(1.2)的解存在并且唯一。 虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，实际问题中归结出来的微分方程主要靠数值解法。 二、欧拉法 我们知道，在平面上，微分方程(1.1)的解称作它的积分曲线。积分曲线上一点的切线斜率等于函数的值，如果按函数在平面上建立一个方向场，那么，积分曲线上每一点的切线方向均与方向场在该点的方向相一致，基于上述几何解释，从初始点出发，先依方向场在该点的方向推进到上一点，然后再从依方向场的方向推进到上一点，循此前进推出一条折线，一般地，设已做出该折线的顶点，过依方向场的方向再推进到，显然两个顶点，的坐标有关系 即 (2.1) 这就是著名的欧拉(Euler)公式。若初值已知，则依公式(2.1)可逐步算出 例1 求解初值问题 解 欧拉公式的具体形式为 取步长，计算结果如下表： 表1 计算结果对比 0.1 0.2 0.3 0.4 0.5 1.1000 1.1918 1.2774 1.3582 1.4351 1.0954 1.18