数值分析第五版计算实习题_数值计算笔记-误差分析

1.1 概述

1. 定义数值计算目标: 寻找一个能迅速完成的（迭代算法）算法，同时估计计算结果的准确度。

1.2 误差分析基础

1. 误差来源：截断误差、舍入误差、数学建模时的近似、测量误差（数据误差）

2. 误差的分类：

绝对误差

；误差限

相对误差

或者

；相对误差限

3. 定义有效数字：从左到右第一位非零数字开始的所有数字

定理：设x与其近似值

的第一位有效数字相同，均为

，若

有p位正确的有效数字，则其相对误差满足：

定理：设对x保留p位有效数字后得到近似值

，则相对误差满足：

定理：设x的第一位有效数字为

，若近似值

的相对误差满足

则

具有p位正确的有效数字，或者在保留p位有效数字后

定理：若x的近似值在

相对误差满足

，则

至少有p位正确的有效数字，或者在保留p位有效数字后

应用：可以不严谨的说如果相对误差不超过

怎有p位正确的有效数字

4. 区分：精度（precision）：有效数字的位数有关

准确度（accuracy）：与准确的有效数字的位数有关

5. 数据传递误差与计算误差：考虑

计算误差：计算过程中的近似引起的误差，例

数据传递误差：单纯由输入数据误差引起的计算结果的误差，例

6. 截断误差：实际结果，算法准确计算得到的结果的差

舍入误差：算法精确计算得到的结果，算法经有限精度计算的得到结果的差

7. 敏感性：输入数据的扰动对问题得影响程度得大小

条件数：

计算公式：

推导可得：

反问题（字面理解）：原问题是y = f(x)，反问题是给定y确定x。反问题的条件数就是探讨反函数

的敏感性。可以化简得出反问题得条件数和原问题得条件数互为倒数，即

8. 算法稳定性

理解：算法很多是迭代得，每一步算法都会不断放大进而传导误差，如果传导误差不会特别大则算法稳定。

分析算法稳定性工具为误差限的四则运算，

的误差限为

，误差限的四则运算如下：

简单的例题：

计算黄金比例的20次方

，初始的近似值x = 0.618034，对于误差限有一下迭代公式：

可以看到系数在以斐波那契数列增长，因此算法不稳定。

1.3 减少误差的方法

避免“大数吃小数”：在编程中很常见的

这种问题；

避免大小相近的数相减：例如1923.05 - 1921.37 = 1.68，这种情况下前两个大数末尾几位已经不精确，在计算过程中又“丢去”高位，这样结果更加不精确。

总结上面的所有对误差的讨论，可以得到对一个问题误差分析如下图：