常微分方程的数值解法实验报告.docVIP

常微分方程的数值解法 专业班级：信息软件 姓名：吴中原 学号：120108010002 一、实验目的 2、明确各种算法的精度与所选步长有密切关系；通过计算更加了解各种 算法的优越性。 二、实验题目 2、试分别取不同步长，考察某节点处数值解的误差变化情况； 3、试用不同算法求解某初值问题，结果有何异常； 4、分析各个算法的优缺点。 三、实验原理与理论基础 (6-1) (6-2) 用数值方法求解时，我们总是认为(6-1)、(6-2)的解存在且唯一。 欧拉法是解初值问题的最简单的数值方法。从(6-2)式由于y (x0) = y0已给定，因而可以算出 设x1 = h充分小，则近似地有： (6-3) 记 从而我们可以取 作为的近似值。利用及f (x1, y1)又可以算出的近似值： 一般地，在任意点处的近似值由下式给出 (6-4) 这就是欧拉法的计算公式，h称为步长。 ，它每一步计算的值一次，截断误差为。 改进的欧拉公式可以改写为： ，它每一步要计算的值两次，截断误差为。 改进的欧拉方法之所以比欧拉方法具有更高的精度，是因为在每一步它都比欧拉方法多计算了一次的值。因此，要进一步提高精度，可以考虑在每一步增加计算的次数。 如果考虑在每一步计算的值四次，则可以推得如下公式： 此公式称为标准四