数值分析期末复习要点总结.pptVIP

* 故， 简单迭代法 27 * Newton 法 基本思想 将非线性方程线性化 设 xk 是 f (x)=0 的近似根，将 f(x) 在 xk 处 Taylor 展开 令： 条件： f’(x) ? 0 * Newton 法 x y x* xk xk+1 * 因为 得牛顿迭代公式为 例7 用Newton法求方程 在(1,2)的根，取初值 [解] Newton法与弦截法 根据 Newton迭代 故在[1,2]内方程有唯一实根 取初值 得 22 * 代入初值得 Newton法与弦截法 故取 23 * k = 0, 1, 2, . . . 迭代函数 牛顿法至少二阶局部收敛 * 例：用 Newton 法求 f(x) = x2 – C=0 的正根, 并验证这一方法的二阶收敛性 解： 对任意 x00， 总有 |q|1， 即牛顿法收敛 * 二阶收敛 关于二次收敛性，事实上 * * * 常微分方程数值解法 (9-2) 公式(9-2)称为显式欧拉公式. Euler法的局部截断误差也可表示为 Euler方法是一阶方法. 第九章 常微分方程数值解法 30 * (9-9) 这就是求解初值问题式（9-1）的梯形公式. 梯形公式（9-9）的局部截断误差为 常微分方程数值解法 故梯形公式为二阶方法. 31 * 常微分方程数值解法 一般地,RK方法设近似公式为 (9-13) 其中 都是参数, 确定它们的原则是使近似 公式在 　处的Taylor展开式与 　在 　处的 Taylor展开式的前面的项尽可能多地重合,这样就使近 似公式有尽可能高的精度。 32 * 常微分方程数值解法 这就是改进Euler公式,它是一种二阶龙格-库塔公式。 33 * 常微分方程数值解法 (9-11) 式（9-11）称为由Euler公式和梯形公式得到的预测— 校正系统，也叫改进Euler法，它是显式单步法。 预测 校正 改进Euler法的局部截断误差为 故它是二阶方法。 34 * 例1用欧拉法求初值问题 当h = 0.02时在区间[0, 0.10]上的数值解。 解 把 代入欧拉法计算公式。就得 具体计算结果如下表： n xn yn y(xn) ?n = y(xn) - yn 0 0 1.0000 1.0000 0 1 0.02 0.9820 0.9825 0.0005 2 0.04 0.9650 0.9660 0.0005 3 0.06 0.9489 0.9503 0.0014 4 0.08 0.9336 0.9354 0.0018 5 0.10 0.9192 0.923 0.0021 * 常微分方程数值解法 例 用改进的Euler法求初值问题 的数值解, 取 . 解 改进的Euler公式为 所以 35 * 常微分方程数值解法 将 代入上式得 将 代入上式得 将 代入上式得 36 * 常微分方程数值解法 例12 用二阶方程初值问题 化为一阶方程组初值问题,并写出欧拉方法求解的 解 令 计算公式. 可得 计算公式为 41 * 定义 　　　 若某算法对于任意固定的 x = xi = x0 + i h，当 h?0 ( 同时 i ? ?) 时有 yi ? y( xi )，则称该算法是收敛的。 定义 　　　若某算法在计算过程中任一步产生的误差在以后的计算中都逐步衰减，则称该算法是绝对稳定的 一般分析时为简单起见，只考虑试验方程 常数，可以是复数 当步长取为 h 时，将某算法应用于上式，并假设只在初值产生误差 ，则若此误差以后逐步衰减，就称该算法相对于 绝对稳定， 的全体构成绝对稳定区域。 h l h = h * 例4.3 对初值问题 证明用梯形公式求得的近似解为 , 并证明当步长h?0时,yn 收敛于精确解 证明: 解初值问题的梯形公式为: ∵ ∴ 整理成显式 反复迭代,得到: * 整理成显式 反复迭代,得到: ∵ ∴ 例 对初值问题 证明用梯形公式求得的近似解为 , 并证明当步长h?0时,yn 收敛于精确解 * 由于 ，有: ∴ 证毕 * 例：考察显式欧拉法 由此可见，要保证初始误差?0 以后逐步衰减， 必须满足： 0 - 1 - 2 例：考察隐式欧拉法 可见绝对稳定区域为： 2 1 0 注：一般来说，隐式欧拉法的绝对稳定性比同阶的显式法的好。 * 6．证明求解初值问题 的如下单步法 是二阶方法。 谢 谢！ 放映结束 感谢各位批评指导！ 让我们共