数值计算中的误差教学课件

数值计算中的误差课件数值计算中的误差课件1 1、舟遥遥以轻飏，风飘飘而吹衣。、舟遥遥以轻飏，风飘飘而吹衣。、舟遥遥以轻飏，风飘飘而吹衣。、舟遥遥以轻飏，风飘飘而吹衣。2 2、秋菊有佳色，裛露掇其英。、秋菊有佳色，裛露掇其英。、秋菊有佳色，裛露掇其英。、秋菊有佳色，裛露掇其英。3 3、日月掷人去，有志不获骋。、日月掷人去，有志不获骋。、日月掷人去，有志不获骋。、日月掷人去，有志不获骋。4 4、未言心相醉，不再接杯酒。、未言心相醉，不再接杯酒。、未言心相醉，不再接杯酒。、未言心相醉，不再接杯酒。5 5、黄发垂髫，并怡然自乐。、黄发垂髫，并怡然自乐。、黄发垂髫，并怡然自乐。、黄发垂髫，并怡然自乐。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作1、例子、例子问题：问题：计算算式计算算式 的近似值：的近似值：可用下列四个式子进行计算：可用下列四个式子进行计算：分别采用近似值：分别采用近似值：和和内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作错误的结果错误的结果内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作 因此，在研究算法的同时，还必须因此，在研究算法的同时，还必须正确掌握误差的基本概念，以及误差在正确掌握误差的基本概念，以及误差在近似值运算中的传播规律，误差分析、近似值运算中的传播规律，误差分析、估计的基本方法和算法的数值稳定性概估计的基本方法和算法的数值稳定性概念。否则，一个合理的算法也可能会得念。否则，一个合理的算法也可能会得出一个错误的结果来。出一个错误的结果来。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作三、三、绝对误差和相对误差绝对误差和相对误差（一）、绝对误差和相对误差限1、绝对误差、绝对误差设某一个量的准确值（称之为真值）为设某一个量的准确值（称之为真值）为 ，其近似值，其近似值为为 ，则，则 与与 的差的差称为近似值称为近似值 的的绝对误差绝对误差，简称误差。当，简称误差。当 时，时，称为称为亏亏近似值或近似值或弱弱近似值，反之则称为近似值，反之则称为盈盈近似值或近似值或强强近似值。近似值。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作2、绝对误差限，或精度、绝对误差限，或精度此 称为近似值 的绝对误差限绝对误差限，或，或精度精度。由于真值往往是未知或无法知道的，因此 的准确值（真值）也就是无法求出。但一般可估计出此绝对误差 的上限，也即可以求出一个正数 ，使内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作（二）、相对误差和相对误差限（二）、相对误差和相对误差限1、为什么要讨论相对误差为什么要讨论相对误差2、相对误差、相对误差定义：绝对误差与真值之比定义：绝对误差与真值之比4、绝对误差与相对误差的关系：3、相对误差限绝对误差与相对误差比较还有一个差别：量纲之差。绝对误差与相对误差比较还有一个差别：量纲之差。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作5、相对误差的其它定义相对误差的其它定义 因为一个量的真值往往是不可能求出的，所以在求相对误差时，常用绝对误差与近似值的比来描述。于是有：百分误差：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作（三）、有效数字及其与误差的关系（三）、有效数字及其与误差的关系1、有效数字、有效数字引子：末位的半个单位引子：末位的半个单位有效数字的通俗理解有效数字的通俗理解分析：当近似值分析：当近似值 的误差限是其某一位上的误差限是其某一位上的半个单位时，就称其的半个单位时，就称其“准确准确”到这一位，且到这一位，且从该位起直到前面从该位起直到前面第一位第一位非零数字为止的所非零数字为止的所有数字都称为有数字都称为有效数字有效数字。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作例如例如 的五、六位有的五、六位有效数字分别为：效数字分别为：数字的规格化形式数字的规格化形式一般说，设有一个数一般说，设有一个数 ，其近似值，其近似值 的规格化形式的规格化形式式中：式中：都是都是0,1,2,3,4,5,6,7,8,9中的一个数中的一个数字，字，；n是正整数；是正整数；m是整数。是整数。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作有效数字有效数字则称则称 为具有为具有n位有效数字的有效数，或称为它位有效数字的有效数，或称为它精度到精度到 。其中每一位数字。其中每一位数字 都都 是的有效数字。是的有效数字。若若 的误差限为的误差限为内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作有效数尾部的零的作用存疑数字：具有n位有效数字的有效数与真值x精确到第n位的近似值在同一位可能相同或相差可能为1。有效位数的长短受到计算机字长的限制。几点注意几点注意203（3），0.0203（3），0.0203（3），0.020300（5）准确值准确值准确值准确值，近似值，近似值，近似值，近似值存凝数字存凝数字内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作2、有效数字与误差的关系、有效数字与误差的关系由上式，可以从有效数字算出近似值得绝对误差限；有效数字的位数越多，其绝对误差限也就越小。绝对误差限相对误差限参考：易大义，参考：易大义，计算方法计算方法，浙江大学，浙江大学内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作绝对误差和相对误差的计算以及有效数字？绝对误差和相对误差的计算以及有效数字？例1 当用 来表示 的近似值时，它的相对误差是多少？解：具有五位有效数字，由(7)有计算题内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作例2 要使积分 的近似值 的相对误差不超过0.1，问至少取几位有效数字？解：可以知道，这样 ，可得可得n=3，即，即 只要取三位有效数字只要取三位有效数字 ，就能保证就能保证 的相对误差不大于的相对误差不大于0.1%。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作四、误差的传播与估计四、误差的传播与估计1、误差估计的一般公式、误差估计的一般公式在实际的数值计算中，参与运算的数据往往都在实际的数值计算中，参与运算的数据往往都是些近似值，带有误差。这些数据误差在多次是些近似值，带有误差。这些数据误差在多次运算过程中会进行传播，使计算结果产生误差。运算过程中会进行传播，使计算结果产生误差。而确定计算结果所能达到的精度，显然是十分而确定计算结果所能达到的精度，显然是十分重要的，但这往往也是件很困难的事。但做一重要的，但这往往也是件很困难的事。但做一些有用的估计还是可以做到的。这里介绍一种些有用的估计还是可以做到的。这里介绍一种常用的误差估计的一般公式，它是利用函数的常用的误差估计的一般公式，它是利用函数的泰勒泰勒(Taylor）展开得到的。）展开得到的。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作以二元函数为例绝对误差：绝对误差：分别是 和 对 的绝对误差增长因子，它们分别表示绝对误差 经过传播后增大或缩小的倍数。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作相对误差：相对误差：分别是分别是 和和 对对 的绝对误差增长因子，它的绝对误差增长因子，它们分别表示绝对误差们分别表示绝对误差 经过传播后增大或缩小的倍数。经过传播后增大或缩小的倍数。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作2、误差在算术运算中的传播、误差在算术运算中的传播在具体应用时，应注意分析加、减、乘、除、乘在具体应用时，应注意分析加、减、乘、除、乘方和开方等算术运算对数据误差的传播规律。方和开方等算术运算对数据误差的传播规律。1、加、减运算近似值之和的绝对误差等于各近似值的绝对误差近似值之和的绝对误差等于各近似值的绝对误差的代数和。的代数和。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作即：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作例如，当要求计算例如，当要求计算 ，结果精确到第，结果精确到第五位数字时，至少取到（八位）五位数字时，至少取到（八位）才能达到具有五位有效数字的要求。如果变换算式（五位）：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作3、例题分析、例题分析问题：利用问题：利用 计算代数式计算代数式 的值，的值，并分析对计算结果的影响。并分析对计算结果的影响。内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作计算结果计算结果内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作分析算法对数值计算结果有重要的影响。分析算法对数值计算结果有重要的影响。由于两个近似数相减，使计算结果的有效由于两个近似数相减，使计算结果的有效数字位数显著减少，以第二种算法尤为严数字位数显著减少，以第二种算法尤为严重。而后两种算法中，则有效数字的损失重。而后两种算法中，则有效数字的损失较少。又由于近似值的较少。又由于近似值的p次乘方的相对误差次乘方的相对误差是该近似值本身的相对误差的是该近似值本身的相对误差的p倍，因此，倍，因此，在后两种算法中以最后一种为最佳。在后两种算法中以最后一种为最佳。目的：目的：分析：分析：内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作l应选用数值稳定的计算算法，避开不稳定的算式；l注意简化计算步骤，减少运算次数；l大数“淹没”小数的现象发生；l应避免两相近数相减（变换）；l绝对值太小的数不宜作为除数；l注意计算过程中误差的传播与积累。五、防止误差传播的若干方法五、防止误差传播的若干方法内江师范学院数学与信息科学学院内江师范学院数学与信息科学学院 吴开腾吴开腾 制作制作介绍了误差理论的基本概念，误差在近似值运算中的介绍了误差理论的基本概念，误差在近似值运算中的传播规律以及估算方法，以及数值稳定性的概念；传播规律以及估算方法，以及数值稳定性的概念；误差的表示方法：绝对误差和相对误差；误差的表示方法：绝对误差和相对误差；误差产生的原因：过失误差和非过失误差；误差产生的原因：过失误差和非过失误差；计算机计算受计算机的字长的限制：因而有效数字的计算机计算受计算机的字长的限制：因而有效数字的概念很重要；概念很重要；常用的误差估计方法：泰勒展开方法；常用的误差估计方法：泰勒展开方法；防止误差传播的几个常用方法。防止误差传播的几个常用方法。本章小结本章小结作业：作业：P11 1,2,7,8,9,12,1341、学问是异常珍贵的东西，从任何源泉吸收都不可耻。阿卜日法拉兹42、只有在人群中间，才能认识自己。德国43、重复别人所说的话，只需要教育；而要挑战别人所说的话，则需要头脑。玛丽佩蒂博恩普尔44、卓越的人一大优点是：在不利与艰难的遭遇里百折不饶。贝多芬45、自己的饭量自己知道。苏联