论述微分方程数值解法.pptVIP

　　　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri = y(xi+1) ? yi+1 称为局部截断误差 /* local truncation error */。 定义2.2 　　　 若某算法的局部截断误差为O(hp+1)，则称该 算法有p 阶精度。 定义2.1 2、欧拉法的局部截断误差 .精品课件. * ? 欧拉法的局部截断误差： Ri 的主项 /* leading term */ 欧拉法具有 1 阶精度。 .精品课件. * 例1: 用欧拉公式求解初值问题 取步长 。 解: 应用Euler公式于题给初值问题的具体形式为： 其中 。 计算结果列于下表： .精品课件. * .精品课件. * 可用来检验近似解的准确程度。 进行计算，数值解已达到了一定的精度。 这个初值问题的准确解为 ， 从上表最后一列，我们看到取步长 .精品课件. * 3、 欧拉公式的改进： ? 隐式欧拉法 /* implicit Euler method */ 向后差商近似导数 x0 x1 )) ( , ( ) ( 1 1 0 1 x y x f h y x y + ? .精品课件. * 由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。 .精品课件. * 一般先用显式计算一个初值，再迭代求解。 隐式欧拉法的局部截断误差： 即隐式欧拉公式具有 1 阶精度。 .精品课件. * ? 梯形公式 /*trapezoid formula */ — 显、隐式两种算法的平均 注：梯形公式的局部截断误差 ， 即梯形公式具有2 阶精度，比欧拉方法有了进步。 但注意到该公式是隐式公式，计算时不得不用到 迭代法，其迭代收敛性与欧拉公式相似。 .精品课件. * ?中点欧拉公式 /* midpoint formula */ 中心差商近似导数 x0 x2 x1 假设 , 则可以导出 即中点公式具有 2 阶精度。 .精品课件. * 方 法 ? ? 显式欧拉 隐式欧拉 梯形公式 中点公式 简单 精度低 稳定性最好 精度低, 计算量大 精度提高 计算量大 精度提高, 显式 多一个初值, 可能影响精度 .精品课件. * ? 改进欧拉法 /* modified Euler’s method */ Step 1: 先用显式欧拉公式作预测，算出 Step 2: 再将 代入隐式梯形公式的右边作校正，得到 1 + n y .精品课件. * .精品课件. * 注:此法亦称为预测-校正法 /* predictor-corrector method */ 可以证明该算法具有 2 阶精度，同时可以看到它 是个单步递推格式，比隐式公式的迭代求解过程 简单。后面将看到，它的稳定性高于显式欧拉法。 改进的欧拉法 .精品课件. * 在实际计算时，可将欧拉法与梯形法则相结合， 计算公式为 应用改进欧拉法,如果序列 收敛, 它的极限便满足方程 .精品课件. * 改进欧拉法的截断误差 因此，改进欧拉法公式具有 2 阶精度 .精品课件. * 例2: 用改进Euler公式求解例1中的初值问题， 取步长 。 解：对此初值问题采用改进Euler公式， 其具体形式为 计算结果列于下表： 例1: 用欧拉公式求解初值问题 .精品课件. * 改进的Euler法 Euler法 .精品课件. * 通过计算结果的比较可以看出，改进的Euler方法 的计算精度比Euler方法要高。 .精品课件. * 第6章 常微分方程数值解法 §6.1 引 言 §6.2 欧拉方法 §6.3 龙格—库塔方法 .精品课件. * §6.1 引 言 微分方程数值解一般可分为：常微分方程数值解和偏微分 方程数值解。自然界与工程技术中的许多现象，其数学表达式 可归结为常微分方程（组）的定解问题。一些偏微分方程问题 也可以转化为常微分方程问题来（近似）求解。Newton最早采 用数学方法研究二体问题，其中需要求解的运动方程就是常微 分方程。许多著名的数学家，如 Bernoulli（家族），Euler、 Gauss、Lagrange和Laplace等，都遵循历史传统，研究重要 的力学问题的数