数值分析教程.ppt

考察第n步的误差 我们有责任改变。 造成这种情况的是不稳定的算法 /* unstable algorithm */ 迅速积累，误差呈递增。 可见初始的小扰动 ? 公式 注意此公式与公式一 在理论上等价。 方法：先估计一个IN ,再反推要求的In ( n N )。 可取 取 We just got lucky? 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 /* stable algorithm */ 在我们今后的讨论中，误差将不可回避，算法的 稳定性会是一个非常重要的话题。 定义：一个算法如果输入数据有扰动（即有误差），而计算过程中舍入误差不增长，则称此算法是数值稳定的，否则称此算法为不稳定的。 1. 要避免除数绝对值远远小于被除数绝对值的除法; 2。避免两个相近的数相减; 3. 要防止大数“吃掉”小数; 2。应选用数值稳定性的计算方法; 2。简化计算步骤和公式，设法减少运算次数。 避免误差危害的若干原则 3.3 避免误差危害的若干原则 避免两个相近的数相减 当遇到两个相近的数相减时，参与运算的数应当多保留几位有效数字或者变换原来公式以避免这种情况的发生 。 由前面公式可知 可以看到，如果两个相近的数相减，则 而相对误差限就会比较大，故有效数字位会大大减少。 较小， 例 给定 若使用计算机计算有 ，应如何变换公式使有效数字位增加？ ，若使用计算器取四位有效数字计算 解 使用计算器计算取四位有效数字得 从而得到 但由于 而使用计算器取四位有效数字得 所以有 这说明变换公式后能使有效数字位由1位增加到3位。 ? 几种经验性避免方法： 当 | x | 1 时： 取右端的有限项近似代替左端。 要防止小数被大数“吃掉”而使有效数字位损失 例 求一元二次方程 在数值运算中，如果两个参与运算的数相差太大，则小数有可能被大数“吃掉”而使有效数字位损失，从而影响计算结果的可靠性。 的根。 远远大于 解 求一元二次方程的根可以使用公式 有可能 可能损失有效数字位，使计算结果出现错误。 按新的求根公式计算得到方程两个准确根为 例如，在只有7位有效数字的计算机系统上使用求根公式解方程 得到的两个根为 要避免这种错误的发生，可以修改求根公式为 ， 要注意减少运算的次数 对于一个计算问题，如果能减少运算次数的话，我们不仅能减少计算时间，提高运行的速度，而且还可以减少误差的积累。 如果把原式子改写为 解 按公式直接计算每一项后，再把每一项求和，就要进行 则计算n次多项式的算法可以是 按秦九韶算法计算n次多项式的值，只需要n次乘法和n次加法。 的值。 例 计算n次多项式 次乘法和n次加法。 例 计算 ln2 的近似值，要求误差小于10??. 解： ①计算量太大②各项的舍入误差会损失和的有效数字 (b) 用级数 计算 用前 9 项（即取 m = 8）计算就能达到精度要求: (a) 用级数 计算 分母接近零的数会产生溢出错误,因而产生大的 误差, 此时可以用数学公式化简后再做. 避免做除数绝对值远远小于被除数绝对值的除法 利用等价变换使下列表达式计算比较精确. 例 答案 误差的种类 模型误差： 观测误差 截断误差 舍入误差 绝对误差 相对误差 误差的表示法 内容回顾 算法设计遵循的条件： （5）绝对值太小的数不宜作为除数 。 （1）应选用数值稳定性的计算方法 ； （2）简化计算步骤和公式，设法减少运算次数 ； （3）合理安排运算顺序，防止大数淹没小数 ； （4）避免两相近数相减 ； * 建立数学模型后，为什么必须用计算机去求解呢？ * 模型误差：在物理中，比如说光滑的水平面，忽略摩擦力，自由落体运动，忽略空气阻力等； 截断误差：近似解与精确解之间的差别，如泰勒级数展开，忽略一个高阶项的时候，就会带入截断误差； 舍入误差：计算机字长有限引入的误差，比如pi的值，现在有效位数已经能够算出多少亿位，但计算机只能取近似值，这时就引入舍入误差； 测量误差：