数值误差教案解析.ppt

科学和工程计算 第一章 绪论 教材及参考资料 清华大学出版社《科学和工程计算基础 》 施妙根 顾丽珍 编著 清华大学出版社《数值分析》李庆扬 王能超 易大义 编 1.1 数值分析研究对象与特点 数值分析的学科别名 1.1 数值分析研究对象与特点应用实例 湖水在夏天会出现分层现象，接近湖面温度较高，越往下温度变低，这种上热下冷的现象影响了水的对流和混合过程，是的下层水域缺氧，导致水生鱼类的死亡。如果把水温看成深度的函数T(x)，有某个湖的观测数据如下： 环境工程师希望：根据给定的数据能求出T(x)。 1.1 数值分析研究对象与特点 什么是数值分析 1.1 数值分析研究对象与特点 计算机只能进行加减乘除四则运算和简单的函数计算 1.数值代数:求解线性和非线性方程的解法,分直接方法和间接方法。 2.插值和数值逼近。 3.数值微分和数值积分。 4.常微分方程和偏微分方程数值解法。 1.2 数值计算的误差 误差的基本理论 误差来源与分类 误差来源与分类 误差来源与分类 机器字长有限—舍入误差 误差来源与分类 小结:模型误差.观测误差不是数值分析讨论的内容,计算方法主要研究截断误差和舍入误差在计算过程中的传播和对计算结果的影响,以提高计算的精度. 1.2 数值计算的误差 误差的基本理论 误差与有效数字 绝对误差/* Absolute error */ 有效数字 有效数字 §2 Error and Significant Digits 例：为使 的相对误差小于0.001%,至少应取几位有效数字？ 解：假设 ?* 取到 n 位有效数字，则其相对误差上限为 要保证其相对误差小于0.001%，只要保证其上限满足 已知 a1 = 3，则从以上不等式可解得 n 6 ? log6，即 n ? 6，应取 ?* = 3.14159。 例3. 解: 则有定理3,相对误差满足 即应取4位有效数字,近似值的误差不超过0.1%. 问题：对于 y = f (x)，若用 x* 取代 x，将对y 产生什么影响？ 分析：e*(y) = f (x*) ? f (x) e*(x) = x* ? x Mean Value Theorem = f ’(? )(x* ? x) x* 与 x 非常接近时，可认为 f ’(? ) ? f ’(x*) ，则有： |e*(y)| ? | f ’(x*)|·|e*(x)| 即：x*产生的误差经过 f 作用后被放大/缩小了| f ’(x*)|倍。故称| f ’(x*)|为放大因子 /* amplification factor */ 或 绝对条件数 /* absolute condition number */. 数值运算的误差估计 * * 计算方法 科学与工程计算 “数值分析”就是研究在计算机上解决数学问题的理论和数值方法。 数值算法的构造 算法的理论分析 计算机解决实际问题的步骤 建立数学模型 选择数值方法 编写程序 上机计算 数值 分析 输入复杂问题或运算 ? ? ? ? 计算机 近似解 1 用计算机进行实际问题的数值计算时,往往求得的是问题的近似解,都存在误差 在建立数学模型过程中,要将复杂的现象抽象归结为数学模型,往往要忽略一些次要因素的影响,而对问题作一些简化,因此和实际问题有一定的区别.—模型误差 在建模和具体运算过程中所用的数据往往是通过观察和测量得到的,由于精度的限制,这些数据一般是近似的,即有观测误差 如: 若将前若干项的部分和作为函数值的近似公式, 由于以后各项都舍弃了,自然产生了误差。 Taylor展开 由于计算机的字长有限,只能对有限位数进行运算,超过的位数按一定规则舍入,产生“舍入误差”. 2 误差是不可避免的,既要允许误差,又要 控制误差.要重视误差分析,分析误差的来源,误差的传播及对误差作出估计 据说，美军 1910 年的一次部队的命令传递是这样的: 营长对值班军官: 明晚大约 8点钟左右，哈雷彗星将可能在这个地区看到，这种彗星每隔 76年才能看见一次。命令所有士兵着野战服在操场上集合，我将向他们解释这一罕见的现象。如果下雨的话，就在礼堂集合，我为他们放一部有关彗星的影片。 值班军官对连长: 根据营长的命令，明晚8点哈雷彗星将在操场上空出现。如果下雨的话，就让士兵穿着野战服列队前往礼堂，这一罕见的现象将在那里出现。 连长对排长: 根据营长的命令，明晚8点，非凡的哈雷彗星将身穿野战服在礼堂中出现。如果操场上下雨，营长将下达另一个命令，这种命令每隔76年才会出现一次。 排长对班长: 明晚8点，营长将带着哈雷彗星在礼堂中出现，这是每隔 76年才有的事。如果下雨的话，营长将命令彗星穿上野战服到操场上去。 班长对士兵