1-1误差分析c2.pdfVIP

第1章 数值分析与科学计算引论 1.1 数值分析的对象、作用与特点 1.2 数值计算的误差 1.3 误差定性分析与避免误差危害 1.4 数值计算中算法设计的技术 1.5 数学软件（略） 1 郑勋烨 1.3 误差定性分析与避免误差危害 一个工程或科学计算问题往往要运算千万次，由于每 步运算都有误差，如果每步都做误差分析是不可能的，也 不科学. 因为误差积累有正有负，绝对值有大有小，都按最坏 情况估计误差限得到的结果比实际误差大得多，这种保守 的误差估计不反映实际误差积累. 2 郑勋烨 考虑到误差分布的随机性，有人用概率统计方法，将 数据和运算中的舍入误差视为适合某种分布的随机变量， 然后确定计算结果的误差分布，这样得到的误差估计更接 近实际，这种方法称为概率分析法. 20世纪60年代以后对舍入误差分析提出了一些新方法， 较重要的有威尔金森(Wilkinson )的向后误差分析法和 穆尔(Moore)的区间分析法两种. 3 郑勋烨 1.3.1 算法的数值稳定性 用一个算法进行计算，如果初始数据误差在计算中传播 使计算结果的误差增长很快，这个算法就是数值不稳定的. 1 例5 计算 I n e 1  x n e x dx (n 0,1,) 并估计误差. 0 (依 “反对幂指三”顺序分部!) 由分部积分可得计算 的递推公式 I n I n  1  n I n 1 ( n  1,2, ), （3.1） 1 I e 1 e x dx 1e 1 . 0  0 若计算出 ， I0 代入（3.2），可逐次求出I 1 , I 2 ,  的值. 4 郑勋烨 1 而要算出 就要先计算 . I e 0 若用泰勒多项式展开部分和 1 (1)2 (1)k e 1(1)   , 2! k ! 并取k  7，用4位小数计算，则得 e 1  0 .3 6 7 9 ，