《数值计算中的误差课件》由会员分享，可在线阅读，更多相关《数值计算中的误差课件（31页珍藏版）》请在金锄头文库上搜索。
1、内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作数值计算中的误差课件Stillwatersrundeep.流静水深流静水深,人静心深人静心深Wherethereislife,thereishope。有生命必有希望。有生命必有希望内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作 数值计算方法数值计算方法，是指将所欲求解的数学模型（数学问，是指将所欲求解的数学模型（数学问题）简化成一系列算术运算和逻辑运算，以便在计算机上题）简化成一系列算术运算和逻辑运算，以便在计算机上求出问题的数值解，并对算法的收敛性、稳定性和误差进求出问题
2、的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。行分析、计算。一、误差分析一、误差分析1 数值计算方法数值计算方法2 误差的含义及其理解误差的含义及其理解误差无处不在。一个合理的算法也可能得出错误的结果。误差无处不在。一个合理的算法也可能得出错误的结果。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作3 算法的算法的数值稳定性数值稳定性 算法算法选得不恰当，不仅影响到计算的速度和效率，还会选得不恰当，不仅影响到计算的速度和效率，还会由于计算机计算的近似性和误差的传播、积累直接影响到计由于计算机计算的近似性和误差的传播、积累直接影响到计算结果的精度，有
3、时甚至直接影响到计算的成败。不合适的算结果的精度，有时甚至直接影响到计算的成败。不合适的算法会导致计算误差达到不能容许的地步，而计算最终失败，算法会导致计算误差达到不能容许的地步，而计算最终失败，这就是算法的数值稳定性问题。这就是算法的数值稳定性问题。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作二、误差的种类及其来源二、误差的种类及其来源模型误差模型误差观测误差观测误差截断误差截断误差舍入误差舍入误差(凑整误差凑整误差)过失误差或疏忽误差过失误差或疏忽误差非过失误差非过失误差内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作
4、制作1、例子、例子问题：计算算式问题：计算算式 的近似值：的近似值：可用下列四个式子进行计算：可用下列四个式子进行计算：分别采用近似值：分别采用近似值：和和内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作错误的结果错误的结果内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作 因此，在研究算法的同时，还必须因此，在研究算法的同时，还必须正确掌握误差的基本概念，以及误差在正确掌握误差的基本概念，以及误差在近似值运算中的传播规律，误差分析、近似值运算中的传播规律，误差分析、估计的基本方法和算法的数值稳定性概估计的基本方法和算法的数
5、值稳定性概念。否则，一个合理的算法也可能会得念。否则，一个合理的算法也可能会得出一个错误的结果来。出一个错误的结果来。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作三、三、 绝对误差和相对误差绝对误差和相对误差（一）、 绝对误差和相对误差限1、绝对误差、绝对误差设某一个量的准确值（称之为真值）为设某一个量的准确值（称之为真值）为 ，其近似值，其近似值为为 ，则，则 与与 的差的差称为近似值称为近似值 的的绝对误差绝对误差，简称误差。当，简称误差。当 时，时，称为称为亏亏近似值或近似值或弱弱近似值，反之则称为近似值，反之则称为盈盈近似值或近似值或强强近似值。
6、近似值。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作2、绝对误差限，或精度、绝对误差限，或精度此 称为近似值 的绝对误差限绝对误差限，或，或精度精度。由于真值往往是未知或无法知道的，因此 的准确值（真值）也就是无法求出。但一般可估计出此绝对误差 的上限，也即可以求出一个正数 ，使内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作（二）、相对误差和相对误差限（二）、相对误差和相对误差限1、 为什么要讨论相对误差为什么要讨论相对误差2 、相对误差、相对误差定义：绝对误差与真值之比定义：绝对误差与真值之比4、 绝对误差与相对误
7、差的关系：3、相对误差限绝对误差与相对误差比较还有一个差别：量纲之差。绝对误差与相对误差比较还有一个差别：量纲之差。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作5、 相对误差的其它定义相对误差的其它定义 因为一个量的真值往往是不可能求出的，所以在求相对误差时，常用绝对误差与近似值的比来描述。于是有：百分误差：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作（三）、有效数字及其与误差的关系（三）、有效数字及其与误差的关系1、有效数字、有效数字引子：末位的半个单位引子：末位的半个单位有效数字的通俗理解有效数字的通俗理解分析
8、：当近似值分析：当近似值 的误差限是其某一位上的误差限是其某一位上的半个单位时，就称其的半个单位时，就称其“准确准确”到这一位，到这一位，且从该位起直到前面且从该位起直到前面第一位第一位非零数字为止的非零数字为止的所有数字都称为所有数字都称为有效数字有效数字。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作例如例如 的五、六位有的五、六位有效数字分别为：效数字分别为：数字的规格化形式数字的规格化形式一般说，设有一个数一般说，设有一个数 ，其近似值，其近似值 的规格化形式的规格化形式式中：式中： 都是都是0,1,2,3,4,5,6,7,8,9中的一个数中的一个
9、数字，字， ；n是正整数；是正整数；m是整数。是整数。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作有效数字有效数字则称则称 为具有为具有n位有效数字的有效数，或称为它位有效数字的有效数，或称为它精度到精度到 。其中每一位数字。其中每一位数字 都都 是的有效数字。是的有效数字。若若 的误差限为的误差限为内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作有效数尾部的零的作用存疑数字：具有n位有效数字的有效数与真值x精确到第n位的近似值在同一位可能相同或相差可能为1。有效位数的长短受到计算机字长的限制。几点注意几点注意203（
10、3），0.0203（3），0.0203（3），0.020300（5）准确值准确值准确值准确值，近似值，近似值，近似值，近似值存凝数字存凝数字内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作2、有效数字与误差的关系、有效数字与误差的关系由上式，可以从有效数字算出近似值得绝对误差限；有效数字的位数越多，其绝对误差限也就越小。绝对误差限相对误差限参考：易大义，参考：易大义，计算方法计算方法，浙江大学，浙江大学内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作绝对误差和相对误差的计算以及有效数字？绝对误差和相对误差的计算以及有效数字
11、？例1 当用 来表示 的近似值时，它的相对误差是多少？解： 具有五位有效数字， ，由(7)有计算题内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作例2 要使积分 的近似值 的相对误差不超过0.1，问至少取几位有效数字？解：可以知道， 这样 ，可得可得n=3，即，即 只要取三位有效数字只要取三位有效数字 ，就能保证就能保证 的相对误差不大于的相对误差不大于0.1%。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作四、误差的传播与估计四、误差的传播与估计1、误差估计的一般公式、误差估计的一般公式在实际的数值计算中，参与运算的数
12、据往往都在实际的数值计算中，参与运算的数据往往都是些近似值，带有误差。这些数据误差在多次是些近似值，带有误差。这些数据误差在多次运算过程中会进行传播，使计算结果产生误差。运算过程中会进行传播，使计算结果产生误差。而确定计算结果所能达到的精度，显然是十分而确定计算结果所能达到的精度，显然是十分重要的，但这往往也是件很困难的事。但做一重要的，但这往往也是件很困难的事。但做一些有用的估计还是可以做到的。这里介绍一种些有用的估计还是可以做到的。这里介绍一种常用的误差估计的一般公式，它是利用函数的常用的误差估计的一般公式，它是利用函数的泰勒泰勒(Taylor）展开得到的。）展开得到的。内江师范学院数学与
13、信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作以二元函数为例绝对误差：绝对误差：分别是 和 对 的绝对误差增长因子，它们分别表示绝对误差 经过传播后增大或缩小的倍数。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作相对误差：相对误差：分别是分别是 和和 对对 的绝对误差增长因子，它的绝对误差增长因子，它们分别表示绝对误差们分别表示绝对误差 经过传播后增大或缩小的倍数。经过传播后增大或缩小的倍数。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作2、误差在算术运算中的传播、误差在算术运算中的传播在具体应用时
14、，应注意分析加、减、乘、除、乘在具体应用时，应注意分析加、减、乘、除、乘方和开方等算术运算对数据误差的传播规律。方和开方等算术运算对数据误差的传播规律。1、加、减运算近似值之和的绝对误差等于各近似值的绝对误差近似值之和的绝对误差等于各近似值的绝对误差的代数和。的代数和。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作即：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作例如，当要求计算例如，当要求计算 ，结果精确到第，结果精确到第五位数字时，至少取到（八位）五位数字时，至少取到（八位）才能达到具有五位有效数字的要求。如果变换
15、算式（五位）：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作3、例题分析、例题分析问题：利用问题：利用 计算代数式计算代数式 的值，的值，并分析对计算结果的影响。并分析对计算结果的影响。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作计算结果计算结果内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作分析算法对数值计算结果有重要的影响。分析算法对数值计算结果有重要的影响。由于两个近似数相减，使计算结果的有效由于两个近似数相减，使计算结果的有效数字位数显著减少，以第二种算法尤为严数字位数显著
16、减少，以第二种算法尤为严重。而后两种算法中，则有效数字的损失重。而后两种算法中，则有效数字的损失较少。又由于近似值的较少。又由于近似值的p次乘方的相对误差次乘方的相对误差是该近似值本身的相对误差的是该近似值本身的相对误差的p倍，因此，倍，因此，在后两种算法中以最后一种为最佳。在后两种算法中以最后一种为最佳。目的：目的：分析：分析：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作l应选用数值稳定的计算算法，避开不稳定的算式；l注意简化计算步骤，减少运算次数；l大数“淹没”小数的现象发生；l应避免两相近数相减（变换）；l绝对值太小的数不宜作为除数；l注意计算过程
17、中误差的传播与积累。五、防止误差传播的若干方法五、防止误差传播的若干方法内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作介绍了误差理论的基本概念，误差在近似值运算中的介绍了误差理论的基本概念，误差在近似值运算中的传播规律以及估算方法，以及数值稳定性的概念；传播规律以及估算方法，以及数值稳定性的概念；误差的表示方法：绝对误差和相对误差；误差的表示方法：绝对误差和相对误差；误差产生的原因：过失误差和非过失误差；误差产生的原因：过失误差和非过失误差；计算机计算受计算机的字长的限制：因而有效数字的计算机计算受计算机的字长的限制：因而有效数字的概念很重要；概念很重要；常用的误差估计方法：泰勒展开方法；常用的误差估计方法：泰勒展开方法；防止误差传播的几个常用方法。防止误差传播的几个常用方法。本章小结本章小结作业：作业：P11 1,2,7,8,9,12,13
