武汉大学 计算方法Chapter1_1.pptVIP

§1.2.2 误差与有效数字 (Error and Significant Digits) ? 绝对误差 /* absolute error */ 其中 x*为精确值，x为x*的近似值。 例如： 工程上常记为 的上限记为 , 称为绝对误差限 /* accuracy*/， ? 相对误差 ( relative error ) x 的相对误差上限 定义为 ? 有效数字 则说x*近似表示x准确到小数后第k位，并从这第k位起 直到最左边的非零数字之间的一切数字(共有n位)都称为 有效数字。 定义1.2：如果 （1.7） 由上述定义 有效数位为3位 有效数位为5位 有效数位为4位 把四舍五入原则抽象成数学语言如下: 定义: 设 是 的近似值,将 写成: (1) 其中ai是0,1,…,9中的一个数,且n为正数,m为整数,且x的绝对误差限满足不等式, 则称具有n位有效数字. 例1:e的近似值为2.718,它的有效数字是几位? 解: 即有效数字为4位,绝对误差限为 例2:按四舍五入的原则写出下列各数具有五位有效数字的近似值. 187.9352,08.000033,2.7182818,0.0002816651. ?有效数字与绝对误差和相对误差的关系： （1）有效数字越多, 绝对误差就越小. （2）有效数字越多, 相对误差就越小. 定理1:若 具有n位有效数字,(x用(1.4)式表示)则其相对误差限为 ,即 证明: 说明:此定理的意义还在于只要知道了近似值的有效数字的位数n和第一个非零数字,就可以估计出它的相对误差限. 定理2:若近似值的相对误差限为  则x至少有n位有效数字. 证明:由于 故x至少有n位有效数字. 1．要使用数值稳定的算法 例4：求 （n = 0, 1, 2, …, 8）的值。 解：由于 初值 §1.3 数值计算中应该注意的一些原则 递推公式 （1.8） 注意此公式精确成立 按 (1.8) 式就可以逐步算出 What happened?! 不稳定的算法 ！ 由递推公式（1.8）可看出，In-1的误差扩大了5倍后传给In， 因而初值I0的误差对以后各步计算结果的影响，随着n的增大 愈来愈严重。这就造成I4的计算结果严重失真。 这就是误差传播所引起的危害 ! 改变公式： 将公式 变为 不妨设I9 ? I10，于是由 可求得I9 ? 0.017，按公式（1.9）可逐次求得 （1.9） I8 ? 0.019 I7 ? 0.021 I6 ? 0.024 I8 ? 0.028 I4 ? 0.034 I3 ? 0.043 I2 ? 0.058 I1 ? 0.088 I0 ? 0.182 稳定的算法 ！ 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 2．要避免两个相似数相减 在数值计算中，两个相近的数作减法时有效数字会损失。 （1.10） 的值。当x = 1000，y 的准确值为0.01580 1、直接相减 2、将（1.10）改写为 则 y = 0.01581 例5: 求 类似地 3. 避免大数吃小数 例7：用单精度计算 的根。 精确解为 ? 算法1：利用求根公式 在计算机内，109存为0.1?1010，1存为0.1?101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1 的指数部分须变为1010，则：1 = 0.0000000001 ? 1010，取单精度时就成为： 109+1=01010+0?1010=0?1010 算法2：先解出 注：求和时从小到大相加，可使和的误差减小。 例8：按从小到大、以及从大到小的顺序分别计算 1 + 2 + 3 + … + 40 + 109 再利用 5. 绝对值太小的数不宜作除数 例6： 如分母变为0.0011，也即分母只有0.0001的变化时 4. 先化简再计算，减少步骤，避免误差积累。 Thank you ! * * 数 值 计 算 方 法 （Numerical Computational Methods） 主 讲 教 师：邹秀芬教授 数学与统计学院 2012年9月 教材 (Text Book) 数值计算方法 朱方生等 编著