数值第一章.ppt

* 所以 x* = 3.14 作为π的近似值,有3位有效数字； 而取x*=3.1416 时, 所以 x* = 3.1416 作为 π 的近似值,有5位有效数字。 例 π = 3.1415926535?, 取 x* = 3.14 时， |x*-x|0.002≤0.005 |x*-x| ≤ 0.000008≤0.00005 * 下面给出有效数字的另一等价定义 用 表示x 的近似值，并将 表示成 若其误差限 ,则称 具有 n 位有效 数字, 这里 m 是整数, a1, a2 ,?, ak 为 0~r-1 中的一个数字, 且a1 ? 0. 定义： 自左向右看, 第一位非零数字 x*=±0.a1a2…an×rm n m x x - * × ￡ - r 2 1 * 例 = 3.1415926535 ?, 取 = 3.14时， 即 m- n = - 2, m=1, n = 3, 所以 = 3.14 作为 近似值时, 就有3 位有效数字。 2 10 2 1 005 . 0 002 . 0 - * × = ￡ ￡ - x x 例 按四舍五入原则，写出下列和数具有5位有效数字的近似数 187.9325 0 8.000033 2.7182818 * 习题: = 3.1415926535 ?, 取 = 3.1415时，5不准应该进位 所以有效数字不准 有几位有效数字,证明你的结论. * 四、 误差限与有效数字的关系 证明 , 故: 此定理说明，相对误差是由有效数字决定的。 定理1 有n 位有效数字， 。则其相对误差限为 设近似值 1 1 10 2 1 ) ( + - * × ￡ n r a x e 0 1 1 a m n a a a x 10 . 0 2 1 × ± = * L m n a a a x 10 . 0 2 1 × ± = * L 1 1 1 1 10 ) 1 ( | | 10 - * - × + ￡ ￡ × m m a x a 1 1 1 1 10 2 1 10 10 5 . 0 ) ( + - - - * * * × = × × ￡ - = n m n m r a a |x | x x x e * 已知 定理2 设近似值 的相对误差 则它至少有n 位有效数字。 不大于 故 至少有n 位有效数字。 证明 m n a a a x 10 . 0 2 1 × ± = * L , 10 ) 1 ( 2 1 1 1 + - × + n a 1 10 -n+1 ￡ 2(a1+1) × n m - × = 10 2 1 1 10 -n+1 × 2(a1+1) × × 10 m-1 ≤(a1+1) * * 例 解 由于 所以 由定理1有 即 得 要使 的近似值的相对误差限小于0.1% , 要取几位有效数字。 故只要对 的近似数取4 位有效数字,其相对误差 因此,可取 就可小于0.1%, * 绝对误差，相对误差，有效位数是度量近似数精度的常用三种量。 实际计算时最终结果以有效数给出, 也就隐含了绝对误差限和相对误差限。 * * 1.4 误差的传播 误差分析的重要性 例 计算 ，n=0,1,2,…,并估计误差 解：因为 可得递推关系 取I0=0.095310,e(I0) ≤0.2×10-6, 计算结果见表1-1. 计算结果不正确! 在数值计算方法中，除了研究数学问题的算法外，还要研究计算结果的误差是否满足精度要求，这就是误差估计问题。 *  一、算术运算的误差 1. 由于 的误差 可看作是x 的微分, 即 则： 可见, 和、差的误差是误差之和、差, 但是因为 * 所以和或差的误差限是误差限之和。 以上的结论适用于任意多个近似数的和或差。 相对误差有： 需讨论两种情况：（1）x和y同号 （2）x和y异号 * 同理可得：乘、除运算的误差,以两数为例有： * 例2 求解二次方程x2-26x+1=0，并估计误差 解：利用二次方程的求解公式得 取 有 于是有 * 相对误差 尽管 的绝对误差比较小，但相对误差却很大，原因是计算 时有效数字的丢失比较多，所以，如果把计算 的公式改为 则有 * 这是 比较好的近似，因为 此时，x2*与x1*的相对误差基本一样。 * 二、函数运算误差 设 f (x)在(a,b)内连续可微, x 的近似值为 , f (x)的近似值为 , 其误差为 , 误差限为 取绝对值得 其中