数值计算的误差PPT课件

在数值分析中，除了研究数学问题的算法外，还要研究计算结果的误差是否满足精度要求，这就是误差估计问题。在数值计算方法中，主要讨论的是截断误差和舍入误差。1.2.2 误差与有效数字定义 1.1 设 是某实数的精确值， 是它的一个近似值，则称 为近似值 的绝对误差，或简称误差。 称为 的相对误差。 AxAxx AxxxxA/ )（Axx用计算机做数值计算时，一般也不能获得数值计算公式的准确解，需要对原始数据、中间结果和最终结果取有限位数字。我们将计算过程中取有限位数字进行运算而引起的误差称为舍入误差。 例如， 如果我们取小数点后四位数字，则 就是舍入误差。.33333.03/1.000033. 03333. 03/ 1第1页/共11页当 时，相对误差没有意义。在实际计算中，精确值 往往是不知道的，所以通常把 作为 的相对误差。0 xxAAxxx/ )( Ax 定义1.2 设 是某实值的精确值， 是它的一个近似值，并可对 的绝对误差作估计 ，则称 是 的绝对误差界，或者称误差界。称 是 的相对误差界。AxAxAAxx AAxAAx/ Axx 例1.1 我们知道 若取近似值 ，则 ，可以估计绝对误差界为0.002,相对误差界为0.0006。 1415926. 3 14. 3 A 0015926. 0 A 解 因为实际问题中所截取的近似数，其绝对误差界一般不超过最小刻度的半个单位，所以当 时，有 ，其相对误差界为cmx954 cmA5 . 0 %.5300005241095450 xA例1.2 1.2 测量一木板长是954cm954cm，问测量的相对误差界是是多大？第2页/共11页 定义1.3 设 是 的一个近似值，将 写成AxxAxikAaaax21.010 (1.2.1)它可以是有限或无限小数的形式，其中 是0，1，9中的一个数字， 为整数。如果), 2 , 1( iaika, 01 ,105 . 0nkAxx 则称 为 的具有 位有效数字的近似值。Axxn可见，若近似值 的误差界是某一位的半个单位，该位到 的第一位非零数字共有 位有效数字的近似值。AxAxn第3页/共11页通常在 的准确值已知的情况下，若要取有限位数的数字作为近似值，就采用四舍五入得到的近似值，其绝对误差界可以取被保留的最后数位上的半个单位。x显然，近似值的有效数字位数越多，相对误差越小，反之也对。下面，我们给出相对误差界与有效数字的关系。 定理1.1 设 的近似值 有（1.2.1）的表达式。xAx（1）如果 有 位有效数字，则Axn;102111nAAaxxx （1.2.2）。32105.0142.3,105.014.3 按定义，3.14和3.142分别是具有三位和四位有效数字的近似值。例如第4页/共11页证 由（1.2.1）可得到。）（111110110 kAkaxa（1.2.4）所以，当 有 位有效数字时，Axn，nknkAAaaxxx 1111102110105 . 0即（1.2.3）得证。则 至少具 位有效数字。（2）如果，）（nAAaxxx 1110121（1.2.3）Axn由（1.2.3）和（1.2.4）有 ，）（）（nknkAaaxx 105 . 0101211011111即说明 有 位有效数字，（2）得证。Axn第5页/共11页 例1.3 已知近似数 的相对误差界为0.3%，问 至少有几位有效数字？AxAx1.2.3 函数求值的误差估计 对一元函数 , 自变量x的一个近似值为 ，以 近似 ，其误差界记作 。若 具有2阶连续导数，与 的比值不太大，则可忽略 的二次项，由Taylor展开式得到 的一个近似误差界：）（xf)( Axf）（xfAx）（Axf）（Axf ）（xf）Axf ( Axx ）（Axf)()()(AAAxxfxf 解 设 有 位有效数字，由于 的第一个有效数 没有具体给定，而我们知道 一定是1，2， ，9中的一个，由nAx1a1a1210)19(21102110003 AAxxx故由（1.2.3）式知 =2，即 至少有2位有效数字。nAxAx第6页/共11页其中 可以得到函数值的一个近似误差界：）。，（）（nAAAkAkxxxfxxf21 ，），（），（)()(12121kAkAnkknAAAnxxxfxxxfxxxf 对n元函数 ，自变量 的近似值分别为 ，则有），（nxxxf21nxxx， 21,21nAAAxxx ）。（）（），（kAnkAknAAAxxfxxxf121特别地，对 有2121xxxxf ），（）。（）（）（AAAAxxxx2121 同样，可以得到），），（）（）（AAAAAAxxxxxx122121 。，）（）（）（0222122121 AAAAAAAAxxxxxxxx 第7页/共11页解 这里 并且有 ，）（，）（2 . 02 . 0 AAdl 。，210800mdlSldSdlSAAA 于是有误差界2422 . 0902 . 0120mSA ）（ 相对误差界。）（）（%39. 01080042AAAArdlSSl 例1.4 设有长为 ,宽为 的某场地。现测得 的近似值 M,d 的近似值 =90M，并已知它们的差界为 试估计该场地面积 的误差界和相对误差界。ld120 AlAd.2 . 0,2 . 0mddmllAA lds l第8页/共11页例1.5 设有三个近似数，24. 293. 131. 2 cbabcap p它们都有三位有效数字。试计算 的误差界，并问 的计算结果能有几位有效数字？ 。6332. 624. 293. 131. 2 p解 于是有误差界 ）（）（）（bcap ）（）（）（bccba 02585.024.293.1005.0005.0 ）（相对误差界。）（）（%39. 06332. 602586. 0 pppr 因为 所以 能有两位有效数字。,05. 002585. 0 ）（p 6332. 6 p1.2.4 计算机中数的表示和舍入误差 任意一个非零实数用（1.2.1）表示，是规格化的十进制科学记数方法。在计算机中通常采用二进制的数系（或其变形的十六进制等），并且表示成与十进制类似的规格化形式，即浮点形式第9页/共11页十进制输入计算机时转换成二进制，并对 位后面的数做舍入处理，使得尾数为 位，因此一般都有舍入误差。两个二进制数作算术运算时，对计算结果也要作类似 的舍入处理，使得尾数为 位，从而也有舍入误差。ttt在实现计算时，计算的最后结果与计算的精确解之间的误差，从根本上说是由机器的舍入误差造成的，包括输入数据和算术运算的舍入误差。因此有必要对计算机数的浮点表示方法和舍入误差有一个初步的了解。有时为了分析某一个计算方法可能出现的误差现象，为了适应人们的习惯，我们会采用十进制实数系统进行误差分析。，tm 21. 02 这里整数m称为阶码，用二进制表示为 或1 , S是阶的位数。小数 称为尾数，其中 或 t是尾数部位的位数。S和t与具体的机器有关。021 jsaaaam，），（sj21 t 21. 0，11 ,3,21tj ，0j由于计算机的字长总是有限位的，所以计算机所能表示的数系是一个特殊的离散集合，此集合的数称为机器数。第10页/共11页感谢您的观看。第11页/共11页