二元函数的误差传播公式-Read.PPTVIP

* 引论 误差 §1 误差的来源与分类 从实际问题中抽象出数学模型 —— 模型误差 通过测量得到模型中参数的值 —— 观测误差 将无穷过程截断 ——截断误差 机器字长有限 —— 舍入误差 如估算 如估算 两个天体的距离100万光年 ±1光年. 注：e* 理论上讲是唯一确定的，可能取正，也可能取负。 e* 0 不唯一，当然 e* 越小越具有参考价值。 §2 误差的基本概念 ? 绝对误差 其中x为精确值，x*为x的近似值。？ ，例如： 工程上常记为 ，称为绝对误差限 ， 的上限记为 零件直径是 20cm?1cm. ? 相对误差定义 x 的相对误差上限 定义 但精确值往往不知道，公式中x常替换为x* 注：严格的说法是， 与 ， 与 是否反映了同一数量级的误差？ 例: 解: 可见,经四舍五入取近似值,其绝对误差限将不超过其末位数字的半个单位 解法之一：将 作Taylor展开后再积分 S4 R4 取 则 称为截断误差 | 舍入误差 | 由截去部分 引起 0.743 由留下部分 引起 如果近似值x*的误差限是某一数位的半个单位,就说准确到该位,从这一位起向左到前面第一个非零的数字为止,所有的数字为x*的有效数字。 有4位有效数字 有6位有效数字 有8位有效数字 只有4位有效数字 ？？ §3 有效数字与误差 ?有效数字文字描述 ?有效数字形式化描述 （其中 ）。若 （即截取按四舍五入规则），则称 为有n 位有效数字，精确到 。 例 问： 有几位有效数字？请证明你的结论。 证明： 有 位有效数字，精确到小数点后第 位。 4 3 用科学计数法，记 （其中 ）。若 （即截取按四舍五入规则），则称 为有n 位有效数字，精确到 。 ?有效数字与相对误差的关系 ? 有效数字 ? 相对误差限 已知 x* 有 n 位有效数字，则其相对误差限为 ? 相对误差限 ? 有效数字 已知 x* 的相对误差限可写为 则 可见 x* 至少有 n 位有效数字。 例：为使 的相对误差小于0.001%,至少应取几位有效数字？ 解：假设 ?* 取到 n 位有效数字，则其相对误差上限为 要保证其相对误差小于0.001%，只要保证其上限满足 已知 a1 = 3，则从以上不等式可解得 n 6 ? log6，即 n ? 6，应取 ?* = 3.14159。 §4 误差的传播 例：蝴蝶效应 —— 纽约的一只蝴蝶翅膀一拍，风和日丽的北京就刮起台风来了？！ NY BJ 以上是一个病态问题 关于本身是病态的问题，我们还是留给数学家去头痛吧！ ※ ※ 例：计算 ? 公式一： 注意此公式精确成立 记为 初始误差 ? ?? ? ! ! ! What happened?! §1 Introduction: Spread Accumulation 考察第n步的误差 我们有责任改变。 造成这种情况的是不稳定的算法 迅速积累，误差呈递增走势。 可见初始的小扰动 ? 公式二： 注意此公式与公式一 在理论上等价。 方法：先估计一个IN ,再反推要求的In ( n N )。 可取 §1 Introduction: Spread Accumulation 取 We just got lucky? 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 在我们今后的讨论中，误差将不可回避， 算法的稳定性是一个非常重要的话题。 ?一元函数的误差传播公式 问题：对于 y = f (x)，若用 x* 取代 x，将对y 产生什么影响？ 分析：e*(y) = f (x*) ? f (x) = f ’(? )(x* ? x) x* 与 x 非常接近时，可认为 f ’(? ) ? f ’(x*) ，则有： |e*(y)| ? | f ’(x*)|·|e*(x)| 即：x*产生的误差经过 f 作用后被放大/缩小了| f ’(x*)|倍。故称| f ’(x*)|为放大因子 /* amplification factor */ 或 绝对条件数 /* absolute condition number */. 相对误差条件数 /* relative