第八章常微分方程数值解.docVIP

第八章 常微分方程数值解 由常微分方程理论可知，我们只能求一些特殊类型的常微分方程。而实际上许多常微 SKIPIF 1 0 分方程求解非常困难。 本章主要讨论一阶常微分方程的初值问题： SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 （8－1） 从理论上讲只要方程中的 SKIPIF 1 0 连续且关于 SKIPIF 1 0 满足李普希兹（Lipschitz）条件，即存在常数L，使 SKIPIF 1 0 则常微分方程存在唯一解 SKIPIF 1 0 。 微分方程数值解：就是求微分方程的解 SKIPIF 1 0 在一系列离散节点 SKIPIF 1 0 处的近似值 SKIPIF 1 0 （i=1，2，…，n） SKIPIF 1 0 称为由 SKIPIF 1 0 到 SKIPIF 1 0 的步长，通常取为常数h。 求数值解，首先将微分方程离散化，常用方法有： SKIPIF 1 0 用差商代替微商 SKIPIF 1 0 若用向前差商代替微商，即 SKIPIF 1 0 （i=1，2，…，n） 则得 SKIPIF 1 0 SKIPIF 1 0 即 SKIPIF 1 0 数值积分法 SKIPIF 1 0 利用数值积分法左矩形公式 SKIPIF 1 0 = SKIPIF 1 0 可得同样算法 SKIPIF 1 0 用泰勒（Taylor）公式 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 得离散化计算公式 SKIPIF 1 0 §1 欧拉（Euler）方法 1.1欧拉方法 对一阶微分方程（8—1），等分区间 SKIPIF 1 0 为 SKIPIF 1 0 份， SKIPIF 1 0 ， 则 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 由以上讨论可知，无论用一阶向前差商，还是用数值积分法左矩形公式，或者用泰勒 SKIPIF 1 0 公式取前两项都可得到同样的离散化计算公式 SKIPIF 1 0 代入初值则得到数值算法： SKIPIF 1 0 （i=1，2，…，n－1） （8－2） 称其为欧拉方法。 几何上欧拉方法就是用一条折线近似表示曲线 SKIPIF 1 0 。（如图）： P0 Pi Pi+1 O SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 1.2欧拉方法的误差估计 定义1 局部截断误差：假设 SKIPIF 1 0 为准确值，用某数值算法计算 SKIPIF 1 0 产生的误差 SKIPIF 1 0 ，称为该数值算法的局部截断误差。 定义2 整体截断误差：准确解 SKIPIF 1 0 与数值解 SKIPIF 1 0 的误差， SKIPIF 1 0 。 设 SKIPIF 1 0 有二阶导数，由泰勒公式有： SKIPIF 1 0 = SKIPIF 1 0 = SKIPIF 1 0 所以 SKIPIF 1 0 = SKIPIF 1 0 ， SKIPIF 1 0 （8－3） 当h充分小时，欧拉方法的局部截断误差与h2是同阶无穷小，称其为一阶方法。 SKIPIF 1 0 定义3 如果一数值解法的局部截断误差为 SKIPIF 1 0 ，则称该算法为 SKIPIF 1 0 阶算法。 1.3 改进的欧拉方法