计算方法第一章误差资料.ppt

(5) 要选用数值稳定的计算公式 对于一个数学问题的求解往往有多种数值方法在选择数值方法时，要注意所用的数值方法不应将计算过程中难以避免的误差放大的较快，造成计算结果完全失真。 例13 计算积分 并估计误差 解 容易得到递推公式 即 为 则准确的理论递推式 实际运算的递推式 两式相减有 这就是说,若 与 的误差为 = - ,即 ，则误差的递推规律为 于是 计算 时的误差被扩大了 倍,显然算法是数值不稳定的。 如果将递推公式 变换一种形式 准确的理论递推式 实际运算的递推式 从而有 即 于是有 则这个算法的误差传递规律为 即每计算一步的误差的绝对值是上一步的十分之一，误差的传播逐步缩小，得到很好的控制，这个算法是数值稳定的 本章小结 误差在数值计算中是不可避免的，误差的传播和积累直接影响到计算结果的精度。在研究算法的同时，必须注重误差分析，使建立起来的算法科学有效。 按照误差产生的来源可分为模型误差、观测误差，截断误差、和舍入误差等。 误差的表示法有绝对误差和相对误差两种。 在表示一个近似数时,要用到有效数字的概念,这在数值计算中非常有用,有效数字是由绝对误差决定的 通常用函数的泰勒展开对误差进行估计 *  * test * * test * * * * * * * 1.3 .3 有效数字 例5. 3.142作为π的近似值时有几位有效数字 解： 3.141592…= 0.3141592…× 3.142 = 0.3142× p = 1 |π-3.142 |=|0.3141592…× -0.3142×? | ＜ 0.000041× ＜ 0.0005= × p –n =1–n =-3 所以 n =4，具有4位有效数字 例6. 当取3.141作为?的近似值时 ??-3.141?=?0.3141592…?101 -0.3141?101? ≤0.0000592 ?101 0.005=1/2 ? 10-2 p-n=1-n=-2 所以n=3具有3位有效数字 推论 如果近似数x*误差限是某一位的半个单位, 由该位到x*的第一位非零数字一共有n位 x*就有n位有效数字,也就是说准确到该位 再如3.1416作为?的近似值时 ? ?-3.1416 ?=? 0.3141592…?101-0.31416?101 ? ≤0?101 ≤ 000005 0.5 ? 10-4 p-n=1-n=-4 所以 n=5 x*= 3.1416有5位有效数字 例7 若 是 的具有六位 有效数字的近似值，求其误差限。 解 由于 具有六位有效数字，即 近似代替 时，准确到末位，于是误差限为 若 是 的具有五位有 效数字的近似值，求其误差限。 解 由于 具有五位有效数字， 即 近似代替 时，准确到末位，于是误差限为 关于有效数字说明 ① 用四舍五入取准确值的前n位x*作为近似值,则 x*必有n位有效数字。如3.142作为 ?的近似值 有4位有效数字，而3.141为3位有效数字 ② 有效数字相同的两个近似数，绝对误差不一定 相同。例如，设x1*=12345,设x2*=12.345,两者 均有5位有效数字但绝对误差不一样 ?x- x1* ?=?x- 12345? ≤ 0.5= 1/2 ? 100 ?x- x2* ?=?x- 12.345?≤0.0005=1/2?10-3 ③ 把任何数乘以10p(p=0,?1,…)不影响有效位数 ④ 准确值具有无穷多位有效数字,如三角形面积 S=1/2ah=0.5ah 因为0.5是真值,没有误差 ?*=0,因此n??,准确值具有无穷位有效数字 1.3.4 有效数字与相对误差 定理1.1 若近似数 具有 n 位有效数字，则其相对误差限 证明 由定义知 是具有n位有效数字的近似值，因此 相对误差限 而