计算方法第二讲12.ppt

第二讲 数值计算中的精确度分析 误差来源与误差估计问题 算法的数值稳定性 病态问题和条件数 1 误差及有关概念 1.1误差的来源 在科学计算的各个环节都会产生误差，按其来源不同，误差可做以下分类。 模型误差 观测误差 截断误差或方法误差 舍入误差 1.2 绝对误差与相对误差 设x为真正值， 为近似值，称： 为近似值的绝对误差（简称误差） 通常我们要求绝对误差不能超过某个值?，?称为绝对误差限或误差限。 设x为真正值， 为近似值，称： 为 的相对误差。 如果存在?r，使得 ，称之为 相对误差限。 在实际计算中，相对误差限很小时，也取： 1.3有效位数与有效数字 如果 的误差限为0.5×10-n，即 则称其准确到小数后第n位，并称 的第一个非零数字到第n位的全部数字为 的有效数字。 问若写成标准表示， m，n，k的关系是什么？ 1.4 数据误差的影响 对两个数x1和x2，简单计算可得： 可见，当x1和x2同号时， 反之，当x1和x2异号时，尤其 这表明，大小接近的异号数相加或大小接近的同号数相减，会严重损失有效数字！ 乘数绝对值很大，或除数接近零时，可能会严重扩大绝对误差，减少精度！ 开方会减少相对误差，提高精度。 一般地，设数学问题的解为 , 近似解为： 则绝对误差为： 相对误差为： 和 起对误差的放大和缩小作用，其绝对值分别称为所求解的数学问题的绝对误差下的条件数和相对误差下的条件数。条件数很大时称该问题为病态问题或坏条件问题，它是问题固有的属性，与算法无关。但由于这类问题数据的微小变化会引起解的剧烈变化，对于这类问题的计算，一般要采用高精度计算，或改变问题的提法，降低条件数。 1.5舍入误差的影响 在计算机中，用浮点法表示的数（称为浮点数）的尾数，位数是固定的，称为字长。设计算机字长为t，任意数x十进制是按舍入原则表为浮点数 绝对误差： 则相对误差的绝对值 记 称?为计算机的相对精度。我们有： 则对于多数相加 相对误差 类似地有 在 大体相同情况下，如 则 于是可得，多数相加时，一般先加绝对值较小的数，相对误差较小！ 通常称舍入误差对计算结果影响不大的算法为稳定的算法，反之为不稳定的算法。 计算数学的特点： 1、面向计算机（构造计算机能用的算法）； 2、要有可靠的理论分析（指误差、算法的收 敛性、稳定性等）； 3、要有好的计算复杂性（指算法省时、节省 内存）； 4、要有数值试验 (检验算法及程序的正确性) 避免误差危害的若干原则： 1、避免除数绝对值远远小于被除数绝对值的 除法。（危害：导致舍入误差增大） 2、避免相近两数相减。（危害：引起有效数 字的严重损失） 3、防止大数吃掉小数。 4、简化计算步骤，减少运算次数。（可以节 省计算机计算时间，更可以减少舍入误差） 例1：求方程 根，如z10系数 ?210略有误差，为?210.000000119，则根20变为20.847，19和18变为19.502?1.94i . 例2：求解微分方程 某些问题的计算中，由于数据的微小变化引起解的剧烈变化，称这类问题为病态问题和坏条件问题。对于这类问题的计算，一定要采用高精度计算。 但对于非病态的良态问题，如算法不当，由于计算机的近似性，有时也可能得到不可靠的结果。 例3：如在尾数为4位的计算机上计算 其真正值为0但计算结果为： 0.0560，但如果先进行有理化在计算，结果为：0.05574，显然，后一种计算精度高。 例4：如在尾数为4位的计算机上计算 精确值为34.5612，计算时如先加前两项，再加后一项，结果为34.57，如先加后两项，再加前一项，结果为34.56，显然，后一种算法更好。 例5：如在尾数为4位的计算机上计算 按两种不同递推计算，结果为： 由此可见，舍入误差对计算有影响，影响小的算法称为数值稳定的算法。 有些算法具有递推性，称之为迭代法或逐次逼近法。 再计算一些复杂的函数的值，有时我们用一些简单的函数（如多项式、有理函数等）来近似之，这称为函数逼近。有时要求逼近函数与被逼近函数之间在某些点函数值及若干阶导数值相等，这种逼近称为插值；有时要求在逼近区间上的最大误差取极小，这称为最佳一致逼近；或在某些点