2013数值分析lec1误差.ppt

避免除数的绝对值远小于被除数的绝对值 说明 在上表中， 是算法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、铝制波纹瓦的长度问题 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的. 假若要求波纹瓦长4英尺,每个波纹的高度(从中心线)为1英寸,且每个波纹以近似2π英寸为一个周期. 求制做一块波纹瓦所需铝板的长度L. 这个问题就是要求由函数 f (x)=sin x 给定的曲线从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 上述积分称为第二类椭圆积分,它不能用初等函数表示出来，只能通过无穷级数表示. 数值积分! 数学 模型 求解复杂问题或运算 如 数值 计算 方法 计算机 近似解 数值分析是做什么用的？ 研究(构造）使用计算机求解各种科学 与工程计算问题的数值方法 对