2012数值分析 -lec1--误差x.pptVIP

避免除数的绝对值远小于被除数的绝对值 两种算法与真实值的比较 0 0.1823 0.1823 0.182322 1 0.08850 0.08839 0.0883922 2 0.05750 0.05804 0.0580389 3 0.04583 0.04314 0.0431387 4 0.02085 0.03431 0.0343069 5 0.09575 0.02847 0.0284684 6 -0.3121 0.02433 0.0243249 7 1.703 0.02123 0.0212326 8 -8.392 0.01881 0.0188369 说明 在上表中， 是算法1计算的值， 是算法2计算的值，而 是真实值的一个近似。从上表我们不难直观的得出结论：随着n的增大，算法一的出来的值是越来越偏离真实值，我们可以说，算法1是不稳定的。 定义：对于某个算法，若输入数据的误差在计算过程中迅速增长而得不到控制，则称该算法是数值不稳定的，否则是数值稳定的。 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 4. 设计数值算法时的几个原则 数值运算总是在一个预先设计好的算法中进行的，所谓算法就是一个有限的基本运算序列。这个序列预定了怎样从输入数据去计算出问题的解。由于运算是在计算机上进行的，而计算机的字长有限，因而产生舍入误差。为减小舍入误差的影响，设计算法时应遵循以下一些原则： 要避免除数绝对值远远小于被除数的绝对值的除法 要避免两相近数相减 要防止大数“吃掉”小数 注意简化计算步骤，减少运算次数 在这种情况下，计算结果对分母的扰动很敏感，而分母 通常是近似值，所以计算结果不可靠； 很小的数作除数，有时会造成计算机的溢出而停机。 (会耗失许多有效数字,可以用数学公式化简后再做). 例：各有五位有效数字的23.034与22.993相减. 23.034-22.993=0.041 0.041只有两位有效数字,有效数字的耗失,说明准确度减小, 因此, 在计算时需要加工计算公式,以免这种情况发生. 例 当 x 较大时,计算 避免相近的两数相减 防止”大数”吃”小数” 当两个绝对值相差很大的数进行加法或减法运算时,绝对值小的数有可能被绝对值大的数吃掉从而引起计算结果很不可靠. 例:求 两者结果不同,因为计算机计算时做加减法要 “对阶”,“对阶”的结果使大数吃掉了小数.产生了误差.为了避免由于上述原因引起的计算结果严重失真,可以根据一些具体情况,存在需要把某些算式改写成另一种等价的形式. 注意简化计算步骤，减少运算次数 例 已知 a0, a1, a2 ,…, an, x, 计算多项式： 直接计算：运算量（乘） 秦九韶算法（1247年）： 运算量 作业 13: 习题：1，3，5，7，8. 阅读《Numerical Recipes in C》： 21. 关于误差和稳定性的内容。 * * * 计算量：一个算法所需的乘除运算总次数，单位是flop. 计算量是衡量一个算法好坏的重要标准。 石家庄铁道学院 数理系 数值分析 朱立永 北京航空航天大学 数学与系统科学学院 “诸位在校，有两个问题应该自己问问， 第一，到浙大来做什么？ 第二，将来毕业后做什么样的人？” 竺可桢 老校长的两句话刻在浙大紫金港校区的一块大石上 这一讲的主要内容 数值分析是做什么的？ 数值分析这门课程的主要内容 这门课程的特点及学习方法 数值分析中的基本概念: 误差 数值计算（分析）是做什么的？ 1、天体力学中的Kepler方程 x是行星运动的轨道，它是时间t 的函数。 非线性方程的数值解法! 2、已经测得在某处海洋不同深度处的水温如下： 深度（M） 466 741 950 1422 1634 水温（oC） 7.04 4.28 3.40 2.54 2.13 根据这些数据，希望合理地估计出其它深度 （如500米，600米，1000米…）处的水温 插值法! 3、铝制波纹瓦的长度问题 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而