误差分析的基本概念

1、第一章 误差分析的基本概念§1 误差的来源1 误差概念 ：精确值与近似值之差称为误差，也叫绝对误差。2 产生误差的主要原因 模型误差：在解决实际问题时，在一定条件下抓住主要因素将现实系统理想化的数学描述称为实际问题的数学模型，这种数学描述常常是近似的，数学模型与实际系统之间存在误差，这种误差称为模型误差。 观测误差：数学模型中往往含有一些由观测得到的物理量（如温度、电阻、长度）或由物理量估算出的模型参数，这些观测物理量或模型参数常常与实际数据存在误差。这种由观察产生的误差称为观测误差。 截断误差：数值计算中用有限运算近似代替无穷过程产生的误差。例如计算一个无穷次可微函数的函数值时，理

2、论上只要能算出这个函数的泰勒级数值即可，但是实际工程上仅用泰勒级数中前面有限项来近似计算函数值，而舍去高阶无穷小量。这个被舍的高阶无穷小量正是截断误差。 舍入误差：计算中按四舍五入进行舍入而引起的误差或因计算机字长有限，数据在内存中存放时进行了舍入而引起的误差。 3.举例说明例1 设一根铝棒在温度t 时的实际长度为Lt ,在 t=0时的实际长度为L0,用来表示铝棒在温度为t时的长度计算值，并建立一个数学模型：,其中是由实验观察得到的常数 （0.0000238±0.0000001）1/，称为模型误差，0.0000001/是的观测误差。这个问题中模型误差产生的原因是：实际上与t2有微弱关

3、系，也就是说模型未能完全反映物理过程。例2 已知在 x=0 处展开的泰勒级数为：为了计算近似值,可取前面有限项计算.如取前面五项计算,计算过程中与计算结果都取五位小数得e1+1+1/2+1/6+1/242.7083，e取五位小数时的准确值为=2.71828，于是截断误差为:这表明：只要在计算中采用了有限步运算近似代替无限步运算的方法，截断误差就一定存在。 例3.=3.1415926；=1.41421356，在计算机上运算时只能用有限位小数，如果我们取小数点后四位小数则：=-3.1416 =-0.0000074；=-1.4142=0.000013就是舍入误差。另外值得一提的是十进制数转化为二进制

4、数时有时也引起循环小数，因计算机上浮点数存储位数限制而舍弃尾部部分小数，如 存储时会引起舍入误差。这个数制转化问题表明：只要计算机内部采用二进制运算，无论计算机发展的多完善，这个舍入误差理论问题永远存在。总的来说，误差一般有：模型误差；观测误差；截断误差；舍入误差。在计算方法这门课程中，截断误差和舍入误差是误差的主要研究对象，讨论它们在计算过程中的传播和对计算结果的影响，并找出误差的上下界，对分析和改进算法都有重大的实际意义。§2 绝对误差 相对误差 有效数字 定义1：设x为准确数，为x的近似值，记e=x-x 称e为x与x 的误差，也叫x与x的绝对误差。显然，x= x+ e 即近似值

5、加误差就是准确值，因此把e也叫做近似值x的修正值，或者说近似值加上修正值就是准确值。误差可正可负，且有量纲单位，当误差为负时，近似值偏大，叫做“强近似”，当误差为正时，近似值偏小，叫做“弱近似”。例1 x=3.14159265 按四舍五入的原则保留不同位数的小数，计算其误差。用一位数字近似表示 用三位数字近似表示 用五位数字近似表示 用六位数字近似表示 定义2：如果 就叫做近似值x的“误差限”，也叫绝对误差限。误差限一定是一个正数。我们常用来表示近似值的精确度或准确值所在的范围（）。现在引入有效数字的概念。如果近似值的误差限是某一位上的半个单位，该位到的第一位非零数字共有n位，我们就说有“n位

6、有效数字”，或者说准确到该位。用四舍五入法取准确值的前n位作为近似值，则有n位有效数字。 以下观察有效数字的位数n与误差限之间的关系 3位有效数字 5位有效数字 6位有效数字 定义3：若用表示x的近似值，并将表示成=±, （及p 为整数, ;, ）若其误差限为就称近似值具有n位有效数字.利用定义3，由有效数字位数n和近似值可以确定误差限： 。注意，首先需要特别指出的是，在有效数字的记法中，有效数字0.123×10-3 和0.1230×10-3是有区别的，前者只有三位有效数字，后者却有四位有效数字；其次，如果只知道x* =300000的绝对误差限不超过500=，则应

7、把它写成300×103或3.00×105,如果仍记为300000，则表示它的误差限不超过0.5，这是因为前者有三位有效数字，后者有六位有效数字；再次，还需要指出的是，一个准确数字的有效位数，应当说有无穷多位。例如对于1/4=0.25不能说只有两位有效数字。例2 若=3587.64是x的具有六位有效字的近似值，那么它的误差限为 定义4：称为近似值的相对误差，当比较小时，有时也把称为近似值的相对误差。相对误差无量刚。相对误差可正可负。我们把相对误差绝对值的上界叫做相对误差限，记作=|, 其中是的误差限(也叫绝对误差限)。推论1. 近似数（n 、及p为整数, 19; 09, 2

8、in）有n位有效数字，则其相对误差限为：证明： 由于有n位有效数字，故与x的绝对误差限应为由相对误差限的定义得:而 由此可以看出，有效数字位数越多，相对误差限就越小。推论2：若近似数（ n, 及p为整数，19； 09，2in）的相对误差限满足：则至少有n位有效数字。证明： (高位进1,舍去尾数,其值变大)由定义3知道：近似数有n位有效数字。证毕。例3. 用来表示e具有三位有效数字的近似值，相对误差限是多少？解：=0.272×101 , n=3 , p=1 , =2 . 由推论1得： =0.0025例4. 为了使的近似值的相对误差小于，问至少要取几位有效数字？解： 由推论2 故按题目要

9、求 令. 则有 即n至少要取为4取n=4查数学用表，其相对误差小于 0.1%§3. 和 差 积 商的误差1. 和 差 积 商的误差设是x的近似值，是y的近似值，用来表示的近似值，则它的误差为(x±y)-(x*±y*)=(x-x*) ±(y-y*) （1-3-1）于是有如下结论：结论1： 和的误差是误差之和，差的误差是误差之差。 |(x±y)-(x*±y*)|x-x*| +|y-y*| （1-3-2）结论2： 两个数和或差的绝对误差限不超过各数绝对误差限之和。结论3： 任意多个数和或差的绝对误差限不超过各个数的绝对误差限之和。结论4：

10、若令则相对误差是对数函数的微分 （1-3-3）设u=xy 则lnu=lnx+lny dlnu=dlnx+dlny 于是有如下结论：结论5 乘积的相对误差是各乘数的相对误差之和。设u=x/y 则lnu=lnx-lny dlnu=dlnx-dlny 于是有如下结论：结论6： 商的相对误差是被除数的相对误差减去除数的相对误差。结论7： 任意多次连乘，连除所得计算结果的相对误差限不超过各乘数和除数的相对误差限之和。证明： 设 w=(uv)/(xy) 则 lnw=lnu+lnv-lnx-lny ; dlnw=dlnu+dlnv-dlnx-dlny |dlnw|dlnu|+|dlnv|+|dlnx|+|d

11、lny| 证毕。例1设y=f(x) 则的相对误差是 例2设则,因此.的相对误差是x的相对误差的n倍。2一般数值运算的误差估计设的近似值依次是，把近似值代入函数y=f()运算得，显然是的近似值，的误差、相对误差如何估计？如果函数y=f()在（ ）附近有连续的二阶偏导数，函数值的误差可用多元函数在()处的泰勒展开式得到。y=f()=令 于是y的误差： （1-3-4）按相对误差定义，y的相对误差为： （1-3-5） 例3 测得某桌面的长a的近似值a*=120cm，宽b的近似值b*=60cm，若已知|a-a*|0.2cm，|b-b*|0.1cm，试求近似面积s*=a*b*的绝对误差限与相对误差限。解：

12、 因为 s=ab ，由（1-3-4）和（1-3-5）式|e*(s*)|600.2|+|1200.1|=24cm2 |e(s*)|=故s*的绝对误差限为24 cm2，相对误差限为0.33% .§4 近似计算中需要注意的几个问题1. 要避免两个相近的数相减在数值计算中，两个相近的数相减，则这两个数的前几位相同的有效数字会在在它们之差中消失，有效数字位数大大减少。例如计算 时，当x接近于零则应变换为来计算；再例如计算，当x充分大时应变换为来计算；当x=1000时，若取4位有效数字计算，两者相减结果为0.02，这个结果只有一位有效数字。但用计算，则得0.01581,它有四位有效数字。这说明应

13、当尽量避免出现这类运算，改变计算方法可以避免两个相近的数相减而引起有效数字损失。通常根据具体情况采用一些数学上的恒等变形如因式分解、分子分母有理化、三角函数恒等式、Taylor展开式等计算公式。2. 两个相差很大的数进行计算时，要防止大数“吃掉”小数例1 计算二次方程解：因式分解得二次方程两根为：，；按求根公式：;其中 =。若计算机上只能表达到小数后八位，则对阶运算时1=在计算中将不起作用，因此，。类似的分析将有， （对阶运算时4ac作零处理）故求得两个近似根 ,；类似方程还有许多如 ：这表明无论计算机发展的多完善，这个大数“吃”小数的问题永远存在。因此，设计算法或编制程序时，一般不要将大小相

14、差非常悬殊的两个数放在一起来运算。3. 要注意计算步骤的简化，减少运算的次数简化计算公式十分重要，它直接影响着计算的速度和误差的积累，有时可以使一项无法实现的计算能够实现，快速富氏变换就是典型例子。下面我们以计算多项式的值为例来说明简化计算公式的重要性。例2 计算多项式 （1-4-1）的值，若直接用上面公式来计算，计算k次项的值需要进行k次乘法，所以计算多项式共需n(n+1)/2次乘法和n次加法才能得到批p (x)的值，但如果我们将公式（1-4-1）改写成下面的形式:令 （1-4-2）对k=1,2, n 反复执行算式（1-4-2）的第2式，则共需n次乘法和n次加法即可得到一个多项式值。这就是著

15、名的秦九韶算法。从上面简单的例子可以看出化简公式不仅能减少运算次数，提高计算速度，而且还能简化逻辑结构，减少误差积累。4. 使用递推关系要注意递推方向的选择，以控制误差的扩大例3 计算 n=0，1，2，7利用定积分的分部积分法，容易得出递推关系式： ，在已知之后可算而得到表(1-4-1)中的第一列。当然也可以按恒等形式的递推关系式：In-1=(1In)n ，在已知之后，可算得而得到表(1-4-1)中的第二列，这八个积分的精确值为表(1-4-1)中的第三列。 表(1-4-1) 两种递推算法对比表 第一种计算法第二种计算法真值I0I1I2I3I4I5I6I70.63210.36800.26400.

16、20800.16800.16000.4000.72000.63200.36800.26430.20730.17080.14550.12690.11240.63210.36790.26420.20730.17090.14550.12680.1124由表中看出，在第一种算法中，随着递推次数的增大，计算结果偏离真值越来越远；而在第二种算法中，随着递推次数的增大，计算结果能稳定地接近真值。我们称第一种算法是不稳定的递推算法，第二种算法是稳定的递推算法。两种算法仅仅只是递推顺序不同，为什么却会出现不同的误差传播呢？如果精确值的近似值有误差=；精确值，近似值，与有误差=；，精确值的近似值有误差=；这就是说

17、若有误差，则的误差的绝对值就是误差的绝对值的n!倍；类似分析可以知道若近似值有误差，则的误差的绝对值就是误差的绝对值的1/（n!）倍。这表明第一种算法计算过程中误差不断扩大，而第二种算法计算过程中误差不断被缩小，这正是算法是否稳定的实质。 通过上述几个问题的简单讨论，我们可以看出，即使有了数学模型，进一步甚至数学上已经有了完善的结果，但仍然存在能不能在计算机上解算和如何实现解算的问题。所以我们必须研究数值计算方法，寻求数学问题在计算机上的有效算法。 习题一1下列各近似数的绝对误差限是最末位的的半个单位，试指出各近似数的绝对误差限及其有效数字位数。 2用秦九韶法计算P(x)=2x3+7x2-9 在x=2处的值。3若a=1.1062, b=0.947是经四舍五入后得到的近似值，问a+b, a×b 有几位有效