数值分析01-误差.ppt

1、阜师院数科院第一章 误差,1-1,主编 余孝华 王谦,数值分析,阜师院数科院第一章 误差,1-2,第一章,误 差,阜师院数科院第一章 误差,1-3,第一章目录,1 误差来源 1.1 模型误差 1.2 观察误差 1.3 舍入误差 1.4 截断误差 2 绝对误差、相对误差和有效数字 2.1 绝对误差与相对误差 2.2 有效数字 3 基本运算中的误差估计,阜师院数科院第一章 误差,1-4,第一章 误 差,半个世纪以来计算机还给我们这个世界的诸多烦恼中，误差问题最为突出。小到银行利率的错算，大到导弹的错误发射，除了操作人员的疏忽、机器的故障引起的过失误差外，计算机在处理数据过程中还存在计算误差。这是计

2、算机机器数系所引起的，这一数系的特点是有限、离散、支离破碎；这和数学上常用的实数系无限、稠密、连续的特点完全不同。机器数的表示方法通常采用浮点数形式，即：,数值计算方法就是“研究用于求得数学问题近似解的方法和过程”，由于算法的实现必须在计算机上进行，虽然计算机是非常准确且快捷的计算工具，但计算机并不是象一般人想象哪样可以解决一切问题而不出差错。,阜师院数科院第一章 误差,1-5,误 差 （ 续1）,其中 ，且 都是整数09中的任一个数。 称为尾数，尾数的位数n是有限正整数； 中的m称为阶数，阶数也是有界的数。所以，机器数中有最大的数，也有最小的数。用机器数表示实数时，很多情况下都带有误差 。,

3、在2400多年前，古希腊人提出了被称为几何三大问题的古典难题。这说明在历史上，人类就常被误差所困扰。下面问题就是三大难题之一。,阜师院数科院第一章 误差,1-6,例 题,解 不妨设已知立方体体积为1。要作的立方体体积为2，则所求方立体高度应该为 ，用计算机计算出 ，（15位数）。尽管精确度相 当高，但仍是近似值。下面的表1-1列出了对h取前有限位数时，计算所得体积的误差。,例1 立方倍积问题。作一个立方体，使其体积 为已知立方体的二倍 。,阜师院数科院第一章 误差,1-7,例 1（续）,表1-1 立方倍积问题的计算,由上表可知，计算机机器数的有限位特点使这一问题只能在满足一定的精度条件下解决，

4、误差是无法消除的。,阜师院数科院第一章 误差,1-8,1 误差来源,（2）在给出的数学模型中往往涉及一些根据观测得到的物理量，如电压、电流、温度、长度等，而观测难免不带误差，这种误差称为观测误差。,一个物理量的真实值和我们算出的值往往不相等，其差称为误差。引起误差的原因是多方面的。,（1）从实际问题转化为数学问题，即建立数学模型时，对被描述的实际问题进行了抽象和简化，忽略了一些次要因素，这样建立的数学模型虽然具有“精确”、“完美”的外衣，其实只是客观现象的一种近似。这种数学模型与实际问题之间出现的误差称为模型误差。,阜师院数科院第一章 误差,1-9,方法误差与舍入误差,（4）在计算中遇到的数据

5、可能位数很多，也可能是无 穷小数，如， ， 等，由于计算机数系是间断的且有界，即计算时只能对有限位数进行运算，因此必须进行四舍五入，这样产生的误差称为舍入误差 。,（3）在计算中常常遇到只有通过无限过程才能得到的结果，但实际计算时，只能用有限过程来计算。如无穷级数求和，只能取前面有限项求和来近似代替，于是产生了有限过程代替无限过程的误差，称为截断误差，这是计算方法本身出现的误差，所以也称方法误差，这种误差是需要特别重视的。,阜师院数科院第一章 误差,1-10,有时，带有误差的数据也被人们频繁使用。例如，在某次人口普查，经统计我国某省的人口数为7123万，这就是一个近似数，其舍入误差不超过0.5

6、万。 用3.1415926来代替圆周率，其舍入误差为,舍入误差,在对收敛的无穷级数计算中，常取有限项代替无穷项。如对于正弦函数:,取 ，作近似计算，则 为其截断误差。,阜师院数科院第一章 误差,1-11,条 件 问 题,计算方法中有一类问题称为条件问题，条件问题是一个算法 （公式）由于初始数据或者中间某些数据微小摄动对计算结果产生影响的敏感性的问题。舍入误差、观测误差都属初始数据的摄动。研究坏条件问题的计算方法是十分重要的课题，有的时候，一些问题的条件并不坏，但由于算法不恰当，初始数据的微小摄动或舍入误差在计算过程中不断被放大，而可能导致计算结果的精度大大降低，甚至使计算失去意义。,阜师院数科

7、院第一章 误差,1-12,递 推 算 法,递推算法是解决实际问题中使用相当普遍的一种算法，它的数学描述是带初值的递推关系式。,例2 小猴吃桃问题。有一天小猴摘下了若干个桃子，当即吃掉了一半，还觉得不过瘾，又多吃了一个。第二天接着吃了剩下的一半，又多吃了一个。以后每天都是吃掉尚存的桃子的一半零一个。到第十天早上，小猴准备吃桃子时，看到只剩下1个桃子了。问小猴第一天共摘下了多少个桃子？,解 设第k天的桃子数为pk，则桃子数目变化规律为,阜师院数科院第一章 误差,1-13,递 推 算 法（续 1）,这是正向递推的关系式，解之，可得逆向递推关系式,由初值 ，根据上式设计算循环算法计算 出 即第一天的桃

8、子数为1534。,上例中仅涉及整数序列递推，根据初值条件来选择正向递推或逆向递推使实际问题得以解决。尽管正向递推和逆向递推公式在数学上完全等价，却导致两种完全不同的算法。对于实数序列的递推由于初始误差的存在，可以一种方向的递推会使误差扩大，而另一方向的递推会使得误差逐步减小。在设计（选用）算法时要用使初始误差不增长的算法。,阜师院数科院第一章 误差,1-14,解: 当n=0时,由此可得 出递推计 算公式:,于是可设计如下两种算法：,递 推 算 法（续 2）,例3,阜师院数科院第一章 误差,1-15,两种算法,算法1,算法2,，由（1-2）可得：,阜师院数科院第一章 误差,1-16,表1-1,由

9、表中结果可见，按算法1得到 ，这显然是错的。,阜师院数科院第一章 误差,1-17,说 明,因为对任意n0均有：,以及 且 时， 。而按算法2计 算，尽管 取值精度不高，其误差 但递推计算得到的 却有8位有效数字，为什么会出现这样的现象？下面的分析说明，这是舍入误差在计算过程中传播所引起的后果。,设 有舍入误差（可能由计算机自动舍入引起），假 定计算过程中不产生新的舍入误差，则由式（1-2）有：,阜师院数科院第一章 误差,1-18,说 明（续1）,阜师院数科院第一章 误差,1-19,说明（续2）,而对算法2，以 计算 应有 从而有: 因此从 出发计算到 时，其误差已缩小 倍。,上例说明，对于同一

10、问题，不同的算法对初始数据的误差（或计算过程中某一步的舍入误差）的传播是不同的，一个算法，经过指定次数计算后，若仍能将初始数据的误差的影响限于一定范围之内，这个算法的稳定性就好，反之稳定性差，上例中，算法2具有数值稳定性、而算法1则是数值不稳定的。显然，只有选用数值稳定性好的算法，才能求得较准确的结果。,阜师院数科院第一章 误差,1-20,2 绝对误差、相对误差和有效数字,2.1 绝对误差与相对误差 设 x *为准确值的近似值，记,一般情况下，准确值是不知道的，从而也不能算出绝对误差e的准确值，但往往可以根据测量工具或计算的情况估计出e 的取值范围，即估计出绝对误差的一个上界 :,这样的称为x

11、 *的绝对误差限或误差限。 显然，误差限不是唯一的。,阜师院数科院第一章 误差,1-21,误差限的意义,容易看出，经过四舍五入得到的数，其误差必定不超过被保留的最后数位上的半个单位，即最后数位上的半个单位为其误差限。例如若取的近似值为3.14，则：,有误差限及近似值，就 可以得到准确值x的范围：,即准确值必定在区间 x*, x*+内， 也常记作：x= x*,阜师院数科院第一章 误差,1-22,因此，要刻划近似值的精确程度，不仅要看绝对误差的大小，还必须考虑所测量值本身的大小，这就是相对误差 er。,误差限的大小不能完全反映近似值的准确程度。例如测量百米跑道长时，误差不超过10厘米，而测量黑板长

12、时得其长度为3米，误差不超过1厘米。就误差限而言，前者为后者的10倍，但由于前者误差只占所量长度的千分之一，而后者误差则占所量长度的三百分之一，显然测量百米跑通的结果更为精确。,相对误差,由于准确值x未知，故一般取相对误差为:,阜师院数科院第一章 误差,1-23,相对误差（续）,可以证明，当|er|很小时， 是er 的高阶无穷小， 可以忽略不计。所以，取绝对误差与近似值之比为相对误差是合理的。,同样，相对误差也只能估计其 上限，如果存在正数r ，使得：,则称r为x*的相对误差限，显然， 可作为x*的一个相对误差限，例如，由实验测得光速近似值为C*=2.997925105公里/秒，其误差限为0.

13、1 公里/秒，于是,所以410-7 是C * 的一个 相对误差限。,阜师院数科院第一章 误差,1-24,2.2 有效数字,一个数能表示大小，如果这个数是一个近似值x*，当然希望能指明它的精确程度，如8与8.000大小一样，但若作为近似值，在引进有效数字概念后，可知其精确程度不一样。,通常要将某个位数很多的数表示成一定的位数，用四舍五入的方法，如=3.14159265.，可表为3.14，3.1416等，这种表示方法的特点是：近似数的误差限为其最末一位的半个单位。即：,x1*=3.14为所有三位数中与相差最小的数，不超过末位（第三位，百分位）的半个单位，即0.510-2;,是所有五位数中与相差最小

14、的数，不超过末位（第五位）的半个单位即0.510-4。,阜师院数科院第一章 误差,1-25,有效数字的定义,定义1,按定义x1*=3.14可称为准确到第三位或有三位有效数字， 而 x2*=3.1416称为准确到第五位，或有五位有效数字。,如果近似值x*的误差限是它的某一位的半个单位，就说 x* “准确”到这一位，并且从这一位直到前面第一个非零数字为止的所有数字均称为有效数字。,也可以给出如下定义：,阜师院数科院第一章 误差,1-26,同样，x*2=3.1416有五位有效数字，因为 x*2=3.1416=0.31416101 。而：,有效数字的定义（续）,x*1=3.14有三位有效数字，是因为x

15、*1=3.14=0.314101，而：,按上述定义，有效数字的概念实际上是说：以x*近似x，如果x*从x依四舍五入规则得到，那么x*的每一位都是有效数字。因此，实际应用时：,阜师院数科院第一章 误差,1-27,有效数字定义的进一步解释,1. 若x已知，可根据四舍五入的原则得x*；,若x未知，则需从近似值的误差界来判断x*的有效位数；,阜师院数科院第一章 误差,1-28,4.有效数字越多，其绝对误差也越小，相对误差同样也越小；,并且：若x*有n位有效数字，则其相对误差限为 ， 若x*的相对误差限为 ，则x*有n位有效数字；,5. 0.0023与0.002300不同，前者最多为二位有效数字，而0.

16、002300则可能具有四位有效数字。,3. 记近似值x*=0.a1a2an10m，若要保留五位有效数字（这是 以后常会用到的），即要求误差限0.510m-n,则n=5；,有效数字定义的进一步解释（续）,阜师院数科院第一章 误差,1-29,有效数字举例,例4,取=3.141592653的近似值分别为3.14，3.141，3.142，3.141592时，其有效数字位数分别为3、3、4、6、6，而作为数0.0509966的近似值，其值分别为0.051、0.0510、0.05100、0.0509、0.05099时，其有效数字位数分别为2、3、4、2、3。,例5,阜师院数科院第一章 误差,1-30,例3

17、（续）,阜师院数科院第一章 误差,1-31,3 基本运算中的误差估计,这里主要讨论四则运算和常用函数的计算中数据误差的传播情况。,设原始数据x1，x2，xn,，y 与xi有关，是由xi计算 所得的解。若x1，x2，xn,的近似值为x1*，x2*， xn,*，那么相应的解也有一定的误差，记为y*，此时解的 绝对误差为:,阜师院数科院第一章 误差,1-32,基本运算中的相对误差,相对误差为:,我们可以利用这两个公式来估计按函数 f 的计算误差。给定 f 的具体形式，就可得到加减乘除及开方这几种基本运算中数据误差与计算结果误差间的关系：,紧 接下屏,阜师院数科院第一章 误差,1-33,具体误差估计,

18、如：对加法：,类似地有：,对乘法：,对除法：,对开方：,因此，有更细的估计总结分析如下：,阜师院数科院第一章 误差,1-34,更细的误差估计分析1,1）对加法：,即：和的绝对误差（或相对误差）不超过相加各 项的绝对误差（或相对误差）之和 。,而x1+x20表示，大小相近的x1，x2异号相加，大小相近的x1，x2同号相减，此时|er(x1+x2)|很大，x1+x2的有效数字会减少。应该避免上述情况出现。,阜师院数科院第一章 误差,1-35,更细的误差估计分析2,2）对乘法：,当x1或x2的绝对值很大时|e(x1x2)|可能很大；,3）对除法 ：,除数x2接近于零时，|e(x1/x2)|可能很大；,4）对开方 ：,通常会缩小相对误差，提高精度；,5）对乘方：,这表明xn的相对误差是x的相对