Ch1-1 -绪论.ppt

2 误差与有效数字 /* Error and Significant Digits */ ? 绝对误差 /* absolute error */ 其中x为精确值，x*为x的近似值。 Hey isn’t it simple? Oh yeah? Then tell me the absolute error of Oops! ，例如： 工程上常记为 ，称为绝对误差限 /* accuracy */， 的上限记为 注：e* 理论上讲是唯一确定的，可能取正，也可能取负。 e* 0 不唯一，当然 e* 越小越具有参考价值。 I can tell that this part’s diameter is 20cm?1cm. I can tell that distance between two planets is 1 million light year ±1 light year. Of course mine is more accurate ! The accuracy relates to not only the absolute error, but also to the size of the exact value. 相对误差是与量纲无关的量 相对误差限： 注意：相对误差衡量的是近似程度的好坏。 如：1百万+1 和 20+1 （绝对误差一样） 1百万光年+1 光年 和 20cm+1cm ? 相对误差 /* relative error */ ?有效数字 /* significant digits */ 用科学计数法，记 （其中 ）。若 （即 的截取按四舍五入规则），则称 为有n 位有效数字，精确到 。 例： 问： 有几位有效数字？请证明你的结论。 证明： 有 位有效数字，精确到小数点后第 位。 4 3 注：0.2300有4位有效数字，而00023只有2位有效。12300如果写成0.123?105，则表示只有3位有效数字。 数字末尾的0不可随意省去！ ?有效数字与相对误差的关系 ? 有效数字 ? 相对误差限 已知 x* 有 n 位有效数字，则其相对误差限为 ? 相对误差限 ? 有效数字 已知 x* 的相对误差限可写为 则 可见 x* 至少有 n 位有效数字。 例：为使 的相对误差小于0.001%,至少应取几位有效数字？ 解：假设 ?* 取到 n 位有效数字，则其相对误差上限为 要保证其相对误差小于0.001%，只要保证其上限满足 已知 a1 = 3，则从以上不等式可解得 n 6 ? log6，即 n ? 6，应取 ?* = 3.14159。 §3 函数的误差估计 /*Error Estimation for Functions*/ 问题：对于 y = f (x)，若用 x* 取代 x，将对y 产生什么影响？ 分析：e*(y) = f (x*) ? f (x) e*(x) = x* ? 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 condition number*/ f 的条件数在某一点是小\大，则称 f 在该点是好条件的 /* well-conditioned */ \坏条件的 /* ill-conditioned */。 注：关于多元函数 的讨论，请参阅教材第8页。 几点注意事项 /* Remarks */ 2. 避免相近二数相减 (详细分析请参阅教材p.10) 例：a1 = 0.12345，a2 = 0.12346，各有5位有效数字。 而 a2 ? a1 = 0.00001，只剩下1位有效数字。 ? 几种经验性避免方法： 当 | x | 1 时： 1. 避免小分母 : 分母小会造成浮点溢出 /* over flow */ 3. 避免大数吃小数 例：用单精度计算