数值分析ch1概要.ppt

数值分析 研究内容 数值逼近 数值微分 数值积分 非线性方程的数值解 常微、偏微的数值解 数值线性代数 研究方法 理论分析 算法分析 误差分析 收敛性分析 收敛速度 误差分析 误差传播与积累 * 提问：数值分析是做什么用的？ 数值 分析 输入复杂问题或运算 ? ? ? ? 计算机 近似解 来源与分类 从实际问题中抽象出数学模型 —— 模型误差 /* Modeling Error */ 通过测量得到模型中参数的值 —— 观测误差 /* Measurement Error */ 求近似解 —— 方法误差 (截断误差 /* Truncation Error */ ) 机器字长有限 —— 舍入误差 /* Roundoff Error */ §1 Introduction: Source Classification 大家一起猜？ 1 1 / e 解法之一：将 作Taylor展开后再积分 S4 R4 /* Remainder */ 取 则 称为截断误差 /* Truncation Error */ | 舍入误差 /* Roundoff Error */ | = 0.747… … 由截去部分 /* excluded terms */ 引起 由留下部分 /* included terms */ 引起 例：计算 ? 公式一： 注意此公式精确成立 记为 则初始误差 ? ?? ? ! ! ! What happened?! 考察第n步的误差 我们有责任改变。 造成这种情况的是不稳定的算法 /* unstable algorithm */ 迅速积累，误差呈递增。 可见初始的小扰动 ? 公式 注意此公式与公式一 在理论上等价。 方法：先估计一个IN ,再反推要求的In ( n N )。 可取 取 We just got lucky? 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 /* stable algorithm */ 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 误差与有效数字 /* Error and Significant Digits */ ? 绝对误差 /* absolute error */ 其中x为精确值，x*为x的近似值。 ，例如： 工程上常记为 ，称为绝对误差限 /* accuracy */， 的上限记为 注：e* 理论上讲是唯一确定的，可能取正，也可能取负。 e* 0 不唯一，当然 e* 越小越具有参考价值。 I can tell that distance between two planets is 1 million light year ±1 light year. §2 Error and Significant Digits ? 相对误差 /* relative error */ x 的相对误差上限 /* relative accuracy */ 定义为 注：从 的定义可见， 实际上被偷换成了 ，而后才考察其上限。那么这样的偷换是否合法？ 严格的说法是， 与 是否反映了同一数量级的误差？ ?有效数字 /* significant digits */ 用科学计数法，记 （其中 ）。若 （即 的截取按四舍五入规则），则称 为有n 位有效数字，精确到 。 例： 问： 有几位有效数字？请证明你的结论。 证明： 有 位有效数字，精确到小数点后第 位。 4 3 注：0.2300有4位有效数字，而00023只有2位有效。12300如果写成0.123?105，则表示只有3位有效数字。 数字末尾的0不可随意省去！ ?有效数字与相对误差的关系 ? 有效数字 ? 相对误差限 已知 x* 有 n 位有效数字，则其相对误差限为 ? 相对误差限 ? 有效数字 已知 x* 的相对误差限可写为 则 可见 x* 至少有 n 位有效数字。