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

第八章 常微分方程数值解 由常微分方程理论可知，我们只能求一些特殊类型的常微分方程。而实际上许多常微分方程求解非常困难。 本章主要讨论一阶常微分方程的初值问题： （8－1） 从理论上讲只要方程中的连续且关于满足李普希兹（Lipschitz）条件，即存在常数L，使 则常微分方程存在唯一解。 微分方程数值解：就是求微分方程的解在一系列离散节点 处的近似值（i=1，2，…，n） 称为由到的步长，通常取为常数h。 求数值解，首先将微分方程离散化，常用方法有： 用差商代替微商 若用向前差商代替微商，即 （i=1，2，…，n） 则得 即 数值积分法 利用数值积分法左矩形公式 = 可得同样算法 用泰勒（Taylor）公式 得离散化计算公式 §1 欧拉（Euler）方法 1.1欧拉方法 对一阶微分方程（8—1），等分区间为份，， 则 由以上讨论可知，无论用一阶向前差商，还是用数值积分法左矩形公式，或者用泰勒公式取前两项都可得到同样的离散化计算公式 代入初值则得到数值算法： （i=1，2，…，n－1） （8－2） 称其为欧拉方法。 几何上欧拉方法就是用一条折线近似表示曲线。（如图）： P0 Pi Pi+1 O 1.2欧拉方法的误差估计 定义1 局部截断误差：假设为准确值，用某数值算法计算产生的误差 ，称为该数值算法的局部截断误差。 定义2 整体截断误差：准确解与数值解的误差，。 设有二阶导数，由泰勒公式有： = = 所以 = ， （8－3） 当h充分小时，欧拉方法的局部截断误差与h2是同阶无穷小，称其为一阶方法。 定义3 如果一数值解法的局部截断误差为，则称该算法为阶算法。 1.3 改进的欧拉方法 由微分方程数值解的三种基本构造方法知，若取不同的差商（如向后），不同的数值积分公式（如梯形公式），以及泰勒公式取前三项、四项等可得不同的算法。 如果用梯形公式计算积分： （8－4） 且 = （8－5） 由于此方程为的隐式方程，不易求解。一般将其与欧拉方法联合使用。 可得算法 （8－6） （k=0，1，2，…；i=1，2，…，n－1） 实际计算中，当比较小时，常取一次迭代后的近似值为，于是有改进的欧拉方法 （i=0,1,2,…,n-1） 例1 用欧拉方法和改进的欧拉方法求微分方程 的数值解（取h=0.1）。 解：由欧拉方法（8－2），得数值计算公式 =＋0.1× 计算结果如表8-1 由改进的欧拉方法（8－6），得数值计算公式 计算结果如表8-2 表8-1 xi 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 yi 1.0000 1.0069 1.0208 1.0391 1.0628 1.0923 1.1269 1.1643 误差 0.0000 0.0037 0.0077 0.0993 0.0120 0.0151 0.0189 0.0222 表8-2 xi 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 yi 1.0000 1.0033 1.0132 1.0292 1.0506 1.0773 1.1079 1.1422 误差 0.0000 0.0000 0.0000 0.0002 0.0020 0.0101 0.0122 0.0053