对数值计算中误差分析.docVIP

对数值计算中误差分析 （一）问题背景： 随着科学技术的突飞猛进，无论是工农业生产还是国防尖端技术，例如机电产品的设计、建筑工程项目的设计、气象预报和新型尖端武器的研制、火箭的发射等，都有大量复杂的数值计算问题亟待解决。他们的复杂程度已达到远非人工手算所能解决的地步。数字电子计算机的出现和飞速发展大大推动了数值计算方法的进展，许多复杂的数值计算问题现在都可以通过电算得到妥善解决。 利用计算机、电子计算机等计算工具来求出数学问题的数值解的全过程，称为数值计算。 关于数值计算中误差的产生与传播以及如何分析与控制各种误差的方法与过程。数据近似值与精确值之差是衡量数据可靠性和精确度的重要方面。应用数值方法在计算机上求解实际问题时，由于模型、测量手段和计算工具等方面的限制，以及计算方法的差异，所得结果往往不是所考虑对象的准确值，而是近似值。 误差按其来源可分为模型误差、观测误差、截断误差和舍人误差等。 1　模型误差 用数值计算方法解决实际问题时，首先必须建立数学模型。由于实际问题的复杂性，在对实际问题进行抽象与简化时，往往为了抓住主要因素而忽略了一些次要因素，这样就会使得建立起来的数学模型只是复杂客观现象的一种近似描述，它与实际问题之间总会存在一定的误差。 2　测量误差 在数学模型中往往包含一些由观测或实验得来的物理量，如电阻、电压、温度、长度等，由于测量工具精度和测量手段的限制，它们与实际量大小之间必然存在误差，这种误差称为测量误差。上面近似公式中地球半径是要经过测量得到，然而无论使用什么工具，其误差是无法避免的。 3　截断误差 由实际问题建立起来的数学模型，在很多情况下要得到准确解是困难的，通常要用数值方法求出它的近似解。例如常用有限过程逼近无限过程，用能计算的问题代替不能计算的问题。这种数学模型的精确解与由数值方法求出的近似解之间的误差称为截断误差，由于截断误差是数值计算方法固有的，故又称为方法误差。 4　舍入误差 无论用计算机、计算器计算还是笔算，都只能用有限位小数来代替无穷小数或用位数较少的小数来代替位数较多的有限小数。在上面的近似公式中的，因为是一个无理数，在计算机中无法精确表示，只能取有限位，一般取3.14159，而将后面无穷多位舍弃。不仅无理数，即便是十分简单的有理数如1/3，也只能用有限位的计算机数近似地表示为0.333333（保留6位）。因此在用计算机进行数值计算时，由于计算机的位数有限，在数值计算时只能近似地表示这些数字，由此而产生的误差称为舍入误差。 舍入地方法比较多，有收尾法（只入不舍）、去尾法（只舍不入）和四舍五入法等，一般常用人们所熟知的四舍五入法。 当然在计算过程中，这类误差往往是有舍有入的，而且单从一次的舍入误差来看也许是微不足道的，但应当注意的是，在数值计算中，当计算机上完成了千百万次运算之后，舍入误差的积累却可能是十分惊人的，这些误差一经迭加或传递，对精度可能有较大的影响。所以，在做数值计算时，对舍入误差应予以足够的重视。 综上所述，数值计算中除了可以完全避免的过失误差外，还存在难以回避的模型误差、观测误差、截断误差和舍入误差。显然，四类误差都会影响计算结果的准确性，而在这四种误差来源的分析中，前两种误差是客观存在的，称为固有误差，而固有误差往往是计算工作者不能独立解决的，是需要与各有关学科的科学工作者共同研究的问题；后两种误差是由计算方法所引起的，称为计算误差，计算误差将是数值计算方法的主要研究对象。 （二）方法综述 ： 误差估计的一般公式，本文介绍一种常用的误差估计的一般公式，它是利用函数的泰勒展开得到的。 设二元可微函数中的自变量*1，*2相互独立，又1，2是自变量*1，*2的近似值，则的近似值。 将函数在点(1，2)处作泰勒展开，并略去其中的高阶无穷小项，即可得到y*的近似值y的绝对误差的估计式为： ＝ ≈ 　　　　　　　　 ＝ 其中两个偏导数应该是在点(1，2)处的值。 近似值y的相对误差的估计式为： 　　　　　　　　　＝ 以上式子中的及分别为各个对的绝对误差和相对误差的增长因子，分别表示绝对误差和相对误差经过传播后增大或缩小的倍数。 由以上两个公式，很容易导出两个近似值和与差的绝对误差和相对误差的估计式： 对于绝对误差有：。即和与差的绝对误差不大于各绝对误差之和。 对于相对误差，考虑最坏的情况是所有相对误差同号，当时，可得： 即和的相对误差不超过各单项中的最大相对误差。 　　同理可得两个近似值之积、商绝对误差和相对误差的估计式： 　　　 在算数的数值稳定性中本文介绍一种方法利用根与系数的关系式，在利用根与系数关系式求第二根时，必须先算出绝对值较大的一个根，然后再求另一个根，才能得到精度较高的结果。 （三）程序设计： 1．编程解一元二