计算方法第1章预备知识与误差分析

计算方法北京工业大学应用数理学院杨中华第一章预备知识与误差分析1.2误差分析

所谓误差就是真实值或精确值与近似值之差。用计算机解决数学问题的过程中会产生多种误差，而引起误差的原因是多方面的，下面我们通过考察解决实际问题的过程来具体说明误差的来源。

1.误差的来源及误差类型

一般使用计算机解决实际问题须经过如下几个过程：实际问题数学模型数值方法程序设计计算结果

(1)模型误差、根据实际问题建立数学模型的过程中通常会忽略某些次要因素而对问题进行简化，由此产生的误差称为模型误差。例如自由落体运动规律的公式

是忽略了空气阻力的影响，因此产生的误差就是模型误差。

误差分析(1.6)

(2)参数误差、很多数学模型都含有一些参数，而有些参数往往又是观测得到的近似值，如此取得的近似参数与真实参数值之间的误差称为参数误差或观测误差。

在自由落体运动规律的公式中重力加速度通常取g=9.8，实际上参数g的值与落体所在位置的地球纬度和海拔高度有关，此近似参数将导致出现参数误差。

(3)截断误差、通常在选定数学模型及参数后还要选择数值算法，很多数值算法往往是将一个无限过程截断为有限过程，此类误差称为截断误差或方法误差。例如前述的差商近似导数的计算方法实质上是在如下展开式中“截断”上式中

的二阶以上无穷小项得到的计算公式，由此产生的误差就是截断误差。误差分析

(4)舍入误差、因为计算机表示浮点数的字长有限(通常4字节或8字节)，除极少量能够准确表示的数据外绝大部分是在超过界定的某位4舍5入所得的近似值，比如圆周率等存入计算机内存时均需要在某位进行4舍5入，此种误差称为舍入误差。

虽然舍入误差在单步运算时或许微不足道，但是在一个较为复杂且连贯的数值算法中，舍入误差可能会积累、传播，所以我们必须加以关注。误差分析

2.绝对误差

定义1.1

设x是准确值，x*是它的一个近似值，称

为x的绝对误差，简称误差，误差的绝对值的上限ε*

称为x*的绝对误差限，简称误差限。

例如在使用圆周率时通常取为π*=3.14，其绝对误差及误差限为

在工程技术或商品规格中通常用如下方式表示产品的误差限：误差分析(1.8)(1.7)

3.相对误差

度量误差的另一种方式是相对误差，例如从10±0.1和1000±1中我们肯定不会认为绝对误差限小的比绝对误差限大的更精确。

定义1.2设x为准确值,x*是它的一个近似值，称比值

为近似值x*的相对误差，相对误差的绝对值的上限记为

称为x*的相对误差限，即有

误差分析(1.9)(1.10)

对于x=10±0.1和1000±1，我们可以计算出它们的相对误差分别为：

从绝对误差看，后者误差限更大些，但从相对误差的角度看，后者更为精确。

在(1.9)式中，在可以确定准确值x的前提下也可以取为x作为分母。误差分析

4.有效数字

定义1.3设x为准确值，其近似值为

这里

为整数，且

则称近似值x*有p位有效数字。

例1.3圆周率

，如果取3.1415和3.14159作为其近似值，问分别有几位有效数字？

解：(1)取

则3.1415有4为有效数字，最后一位数字5是无效数字。误差分析(1.11) (2)取

则取3.14159有6为有效数字。

我们在近似计算中四舍五入保留小数点以后n位数字时有效数字的个数，如同样保留小数点后4位数字 1.414213… 近似值1.4142有5位有效数字 0.012353… 近似值0.0124有3位有效数字 28.08918… 近似值28.0892有6位有效数字

总之：对精确值进行四舍五入所得结果后所有非零数字均是有效数字。

误差分析

定理1.1设近似值x*具有(1.11)式的形式，如果有n位有效数字或是由其准确值4舍5入后得到的近似值，则x*的相对误差限为：

证明：如果

有n位有效数字，则其绝对误差限为

相对误差限为

如果x*是由其准确值4舍5入后得到的近似值，则此近似值必然有n位有效数字，从而结论也必然成立。误差分析

5.近似计算中减少误差的几个策略 1)避免两个相近的数相减

我们知道一个近似值其有效数字愈多愈为精确，如果两个相近的数相减势必减少其有效数字的个数，使得相对误差变大。

例1.4已知

x=1.5846,y=1.5839，求x-y的值。

解：x–y=1.5846–1.5839=0.0007，显然x，y均有5位有效数字，而

x–y仅有1位有效数字，使得其有效数字大大的减少。误差分析

例1.5已知

x=18.496，y=18.493，取4位有效数字，计算x-y的近似值，并估计其相对误差。

解：先取4为有效数字，x*=18.50，y*=18.49，再计算其差

x*–y*=18.50–18.49=0.01

实际上按准确值计算 x–y=18.496–18.493=0.003

相对误差为：误差分析

在编程序时，可采取以下措施减少误差尽量使用双精度定义变量当两个接近的量相减时，最好做等价变换避免损失有效数字，如当x的值较大时令：当x与y非常接近时令：误差分析

2)避免绝对值太小的数作为分母

如

时，计算比值

会使绝对误差变大，通常也做等价变换，如

时：

3)防止大数“吃掉”小数

例1.6计算定积分

，其中N是一个非常大的正数。

解：根据Newton-Leibniz公式得到误差分析

但是如果按此式进行计算，由于N+1时1可能被N吃掉，所以实际计算时可能是N+1=N，因此计算结果可能是

但是实际上

结果严重失真，正确的计算方法应该是：误差分析

另一方面，根据积分中值定理有

我们对于

编程实际计算，结果如下：

因为

，断定后两个结果的精度优于第一个计算结果。误差分析1.3向量和矩阵的范数

如同误差分析是计算方法的预备知识一样，有关范数的概念也是学习计算方法所必须的基础知识。

通常用绝对值来度量实数的大小，用模来度量复数的大小，用范数来度量向量、矩阵的大小。

范数是论述、证明某些数值方法的收敛性和稳定性所必须的工具，本节将简单介绍范数的有关概念。 1.向量的范数

定义1.4如果定义在

上的一个实值函数

，对于任意的

都有 1)非负性：

且仅当

时等号成立 2)齐次性： 3)三角不等式：

则称该实值函数

为

上的一个向量范数。

满足定义1.4中三个条件的实值函数均是向量范数，但最常用的向量范数有： 1-范数 2-范数 ∞-范数向量和矩阵的范数(1.12)(1.13)(1.14)

对于2-范数与内积的关系，有著名的Cauchy-Schwarz(柯西—许瓦尔兹)不等式：

我们可以证明这三种向量范数均满足定义1.4中的三个条件，仅以2-范数为例加以证明，另两个留作习题。

证明：条件1)非负性是显然的，对于条件2)的齐次性

对于条件3)向量和矩阵的范数(1.16)

以上三种向量范数可以统一于p-范数定义之下，

p-范数定义如下：

当p=1、2时(1.17)式与(1.12)、(1.13)显然是一致的，对于p=∞情况，考虑下不等式向量和矩阵的范数

即有

对此式取p→∞时的极限，则有(1.17)例1.7求向量

前述的三种向量范数。

解：

如前所述向量范数可以度量向量的“大小”，而范数又不是唯一的，虽然对应同一个向量不同的范数值是不一样的，但是我们可以证明任意两种向量范数是等价的，首先给出向量范数等价性的定义。向量和矩阵的范数

定义1.5设有两种范数

，如果存在常数

使

成立，则称两种范数是等价的。

按此等价性的定义，常用的三种范数是等价的。事实上很容易证明以下三式成立

三种范数关系的综合不等式为：向量和矩阵的范数(1.18)

2.矩阵的范数

定义1.6如果定义在

上的一个实值函数

，对于任意的

都有 1)非负性：

且仅当

时等号成立 2)齐次性： 3)三角不等式：

4)相容性：

则称该实值函数

为

上的一个矩阵范数。

一个最常用而且便于计算的矩阵范数为Frobenius(弗罗贝纽斯)范数，定义如下：

该矩阵范数也简称为矩阵的F-范数，如果把矩阵按相邻行的首尾连接构成向量，则F-范数正是此向量的2-范数。向量和矩阵的范数(1.19)

以下讨论矩阵范数与向量范数之间的关系与相容性问题。

定义1.7对于给定

上的向量范数

和

上的矩阵范数

如果有

则称矩阵范数

与向量范数

是相容的。

我们也可以按相容性条件定义从属于向量范数的矩阵范数，即对于指定的向量范数

来定义矩阵的范数：

如此定义的矩阵范数

称为向量范数的从属范数，可以证明此从属范数满足矩阵范数定义的4个条件。向量和矩阵的范数(1.20)

从属于向量范数的矩阵1-范数、2-范数、

范数的表达式：

其中

是矩阵

的最大特征值。

对于以上3种从属范数我们均可以给出证明，此处我们仅就(1.21)式的1-范数给出证明向量和矩阵的范数(1.21)(1.22)(1.23)向量和矩阵的范数

证明：将矩阵A按列分块记为

假设第k列向量为范数最大的，即

从

可推得

特别的取

，则有

因此前述推导的上限，可以达到，则(1.21)式成立。例1.8

求矩阵

的各种矩阵范数。

解：先计算

其特征值为：

则可求出向量和矩阵的范数

3.矩阵的范数一些结论

以下不加证明的给出一些范数的有关定义和定理

定理1.2(范数等价定理)

设

与

是两种不同的向量范数和从属的矩阵范数，则存在常数

，使得对任何

，有

由于任何两种范数都是等价的，在本教材其后的某些章节讨论数值算法的收敛性、稳定性时我们可以以任一种范数对向量、矩阵的“大小”进行度量。向量和矩阵的范数

定义1.8设矩阵

，其特征值为

，则称

为矩阵A的谱半径。

矩阵2-范数与谱半径的关系：

对称矩阵：

正交矩阵：

定理1.3

设矩阵

，则对矩阵的任何范数

均有：

本定理说明，谱半径是所有范数的一个下界。

定理1.4设

的充要条件是