7 常微分方程2 计算方法课件及实验 教学课件.pptVIP

几种方法比较 ? 最常用为四级4阶经典龙格-库塔法 /* Classical Runge-Kutta Method */ ： 注： ? 龙格-库塔法的主要运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系： 7 5 3 可达到的最高精度 6 4 2 每步须算Ki 的个数 ? 由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h 取小。 算例： 分别用Euler公式，改进的Euler公式，经典4阶R-K 公式计算一阶常微分方程初值问题。 并与准确解 比较。 解：Euler公式, 改进的Euler公式取步长h=0.1, 经典4阶R-K公式取步长h=0.2。 4阶R-K公式： 1.7321 1.7321 1.7379 1.7848 1.0 1.6733 1.6782 1.7178 0.9 1.6125 1.6125 1.6153 1.6498 0.8 1.5492 1.5525 1.5803 0.7 1.4832 1.4833 1.4860 1.5090 0.6 1.4142 1.4164 1.4351 0.5 1.3416 1.3417 1.3434 1.3582 0.4 1.2649 1.2662 1.2774 0.3 1.1832 1.1832 1.1841 1.1918 0.2 1.0954 1.0959 1.1000 0.1 1.0000 1.0000 1.0000 1.0000 0.0 准确值 4阶R-K公式 改进Euler公式 Euler公式 X 计算结果见下表： * 计 算 方 法 南昌航空大学信息工程学院 第七章 常微分方程的数值解法 孙成立 §7.1 基本离散方法 §7.2 龙格-库塔方法 §7.3 线性多步法 §7.4 算法的收敛性与稳定性 §7 常微分方程数值解法 考虑一阶常微分方程的初值问题 : 例如： 其解析解为： §7.1 基本离散方法 但是, 只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解。 因此，只能依赖于数值方法去获得微分方程的数值解。 例如： 其解析解为： 很难得到其解析解因此，只能依赖于数值方法去获得微分方程的数值解。 ? Taylor级数复习 函数 在 点作Taylor级数展开： 这里 x , x0 都可以是任意一点。 ?Taylor级数复习 则： 7.1 Euler公式 利用等距分割，数值微分来代替导数项，建立差分格式。 1、向前欧拉公式（或显式欧拉公式） 所以，可以构造差分方程 称为局部截断误差。显然，这个误差在逐步计算过程中会传播，积累。因此还要估计这种误差 Euler公式 定义:　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri = y(xi+1) ? yi+1 称为局部截断误差。 定义: 　若某算法的局部截断误差为O(hp+1)，则称该算法有p 阶精度。 yk 局部截断误差 p阶精度 O Tk+1 yk+1 y(x) y x h Euler公式的精度为1阶 欧拉公式几何意义 2、向后欧拉公式(或隐式欧拉公式) 可以由向前差商公式求出 由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 欧拉公式。可以用迭代求解右端yn+1估计值 可由迭代法求得 3、中心欧拉公式 是两步欧拉格式，2阶误差，该格式不稳定 需要2个初值 y0和 y1来启动递推 过程，这样的算法称为双步法 /* double-step method */，而前面的三种算法都是单步法 /* single-step method */。 对微分方程 积分有: 类似，可以算出其误差估计式： 2阶的方法 所以，有 是个隐式的方法， 要用迭代法求解 局部截断误差 4、梯形公式 方 法 ? ? 显式欧拉 隐式欧拉 梯形公式 中点公式 简单 精度低 稳定性最好 精度低, 计算量大 精度提高 计算量大 精度提高, 显式 多一个初值, 可能影响精度 改进欧拉方法 Step 1: 先用显式欧拉公式作预测，算出 Step 2: 再将 代入隐式梯形公式的右边作校正，得到 又称为预估－校正Euler方法 Euler方法、隐式Euler方法、梯形方法与单步法计算公式的 显式单步法 对应关系 隐式单步法 显式 Euler方法 隐式Euler方法 梯形方法(隐式) 7.1.3 总结