计算方法-第七章常微分数值解教学课件.pptVIP

Runge-Kutta法的基本思想（3） 二、二阶龙格－库塔方法 三、三阶龙格－库塔方法 四、四阶龙格－库塔方法 两点说明: 五、变步长的龙格—库塔方法 六、R-K方法的绝对稳定区域 考虑一阶常微分方程的初值问题 只要 f (x, y) 在[a, b] ? R1 上连续，且关于 y 满足 Lipschitz 条件，即存在与 x, y 无关的常数 L 使 对任意定义在 [a, b] 上的 y1(x) 和 y2(x) 都成立，则上述初值问题存在唯一解。 现要算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 ，即数值解。 节点间距 为步长，通常采用等距节点，即取 hi = h (常数)。 7.1 微分方程数值解法 第七章 常微分方程数值解 求解最基本的方法是单步法 单步法:从初值 开始,依次求出 ,后一步的值 只依靠前一步的 典型的单步法是Euler(欧拉)方法,其计算格式是: 例:求解常微分方程初值问题 由此可见,Euler公式的近似值接近方程的精确值. x0 x1 向前差商近似导数 记为 以Euler法为例说明构造初值问题数值方法的三种基本途径 1. 数值微分法, 用差商代替微商 2. Taylor展开法 7.2 欧拉公式 7.2.1 显式 欧拉公式 忽略高阶项,取近似值可得到Euler公式 3. 数值积分法区间 将 区间 积分 向后差商近似导数 x0 x1 )) ( , ( ) ( 1 1 0 1 x y x f h y x y + ? 由于未知数 yn+1 同时出现在等式的两边，不能直接得到，故称为隐式 欧拉公式。 一般先用显式欧拉公式计算一个初值，再迭代求解。 7.2.4 隐式欧拉公式 梯形公式 — 显、隐式两种算法的平均 中点欧拉公式 中心差商近似导数 x0 x2 x1 7.2.6 改进的欧拉公式 Step 1: 先用显式欧拉公式作预测，算出 ) , ( n n n y x f h y + = 1 + n y Step 2: 再将 代入隐式梯形公式的右边作校正，得到 1 + n y )] , ( ) , ( [ 2 1 1 + + + + = n n n n n x f y x f h y y 1 + n y 7.2.5 梯形公式 注：此法亦称为预测-校正。一方面它有较高精度，同时可以看到它是个单步递推格式，比隐式公式的迭代求解过程简单。后面将看到，它的稳定性高于显式欧拉法。 局部截断误差:设 是初值问题(1)的解,用单步法计算到第n步没有误差,即 ,则 定义 1. 整体截断误差和局部截断误差 整体截断误差:数值解 和精确解 之差 整体截断误差除与 步计算有关外,还与 的计算 有关 分析计算中的某一步,显式单步法的一般形式可写为: 其中 称为增量函数。如对于Euler公式其增量函数 7.2.3 局部截断误差、7.4 单步法收敛性、7.5稳定性 称为单步法在点 处的局部截断误差。 定义 　　　若某算法的局部截断误差为 ，则称该算法有p 阶精度。 欧拉法的局部截断误差,由Taylor展开: 欧拉法具有 1 阶精度。 类似可以证明改进的Euler方法具有2阶精度（P134） 2. 收敛性和整体截断误差 　　　 若某算法对于任意固定的 x = x0 + n h，当 h?0 ( 同时 n ? ?) 时有 yn? y( xn )，则称该算法是收敛的。 定义 例：就初值问题 考察欧拉显式格式的收敛性。 解：该问题的精确解为 欧拉公式为 对任意固定的 x = xn = nh ，有 ? ? 关于整体截断误差与局部截断误差的关系,有如下定理 定理:对初值问题的单步法 , 若局部截断误差为 ,且函数 对y 满足