计算方法的稳定性

在实际数值计算过程中,由于不可避免地存在和不断产生各种误差,因此计算结果不是绝对精确的。如果误差使得计算结果和实际情况有较大差别或者出现错误的结果,则数值 计算便失去了价值和意义。因此,分析数值计算过程中误差的来源和传递规律,设法控制和减小误差。

1.误差的来源

来源:固有误差(模型误差、观测误差)和计算误差(截断误差、舍入误差)

舍入误差

设s是r进制数,p是r进制正负整数或零,则形如
微小偏差可以不进行的工作 微小误差舍去原则_进制
并满足
微小偏差可以不进行的工作 微小误差舍去原则_舍入误差_02
的数x,称为r进制浮点数,且s和p分别称为浮点数的尾数和阶数。

任何一种计算机只能用有限的位数来表示浮点数的尾数和阶数。设
微小偏差可以不进行的工作 微小误差舍去原则_舍入误差_03
其中微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_04为正整数,它们主要由计算机用多少位数来表示阶数而决定。如果尾数的小数尾数为t(t一般比m,M的位数大若干倍),则计算机的数系由一切阶数满足微小偏差可以不进行的工作 微小误差舍去原则_进制_05的t位r进制浮点数的集合F组成。可见,在计算机数系中,数的个数有限,数系中的每一个数都是有理数,且阶数相等的数以相等的距离分布在数轴的某一段上。

由于计算机数系是有限集,不仅无理数微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_06等不属于计算机数系,而且一些有理数(如微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_07)也属于计算机数系,因此常常用计算机数系中和它们接近的数来表示它们。同时,在利用计算机进行计算时,由于字长限制,参与计算的数的长度也是有限的,而由此产生的误差,称为舍入误差

例1

浮点数F的集合可以用以下4个参数来描述:
微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_08
其中,微小偏差可以不进行的工作 微小误差舍去原则_进制_09为基数,t是精度参数,整数L与U是阶码E(范围)的下限和上限微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_10

这样,F中的每一个浮点数x的值可表示为:
微小偏差可以不进行的工作 微小误差舍去原则_浮点数_11
式中的整数微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_12满足微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_13,同时又微小偏差可以不进行的工作 微小误差舍去原则_浮点数_14(E是整数)。

如果对F中每个非零的x,有微小偏差可以不进行的工作 微小误差舍去原则_浮点数_15,则称浮点数系F为规格化浮点数系。括号中的部分微小偏差可以不进行的工作 微小误差舍去原则_进制_16称为尾数。我们知道,一个实数常用【整数+小数点+尾数】的形式表示,它们在计算机中对应的浮点数微小偏差可以不进行的工作 微小误差舍去原则_浮点数_17则常用某种整数表示方式(例如以原码、反码或补码的形式)存储。

例2:求十进制数系与计算机采用的二进制数系之间的差别

:以在许多算法中常被选作步长的十进制0.1为例。在微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_18或者为2的幂的浮点数系中,10个0.1的步长并不刚好等于一个1.0的步长。事实上,当把微小偏差可以不进行的工作 微小误差舍去原则_浮点数_19转换成为以微小偏差可以不进行的工作 微小误差舍去原则_舍入误差_20为底的幂的有限项展开式时,有:
微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_21
用下标表示数基,如果微小偏差可以不进行的工作 微小误差舍去原则_微小偏差可以不进行的工作_18,则有:
微小偏差可以不进行的工作 微小误差舍去原则_进制_23
如果微小偏差可以不进行的工作 微小误差舍去原则_舍入误差_24,则
微小偏差可以不进行的工作 微小误差舍去原则_舍入误差_25
由于字长限制,等号右边的值只能取7位。很明显,当把10个这样的数相加时,其结果并不正好是1.0。这就是舍入误差造成的。所以,在计算机上进行的浮点运算(四则运算)只能是近似计算。

观测误差和数据的舍入误差虽然来源不同,但对计算结果的影响完全一样。在数值计算中涉及的误差一般指舍入误差(包含初始数据误差)和截断误差。

2 计算方法的稳定性

计算方法的稳定性是指数值计算中是否稳定的问题。在数值计算过程中,数值解是逐步计算出来的。由于计算机的字长有限,每一步计算都有误差存在,且前一步的舍入误差必然要影响下一步的近似解。如果运算序列的舍入误差不增长。误差的积累或传递对计算结果的影响是可控的,则该算法是数值稳定的,否则是数值不稳定的。

3. 数值计算的基本原则

评价一个数值计算方法优劣的标准:

  • 计算时间复杂度(运算次数或计算时间),包括收敛性问题
  • 计算空间复杂度(占用计算机存储空间)
  • 计算结果精确度(包括稳定性问题)

构造和选择一个好的计算方法:

  1. 避免两个相近的数相减

在数值计算过程中,两个相近的数相减,会严重损失有效数字,从而使相对误差变大。

如果两个相近的数相减,常采用变换的公式进行计算。如果计算公式不能改变,则采用增加有效数字位数的方法。

  1. 避免使用绝对值很小的数作分母
  2. 两个相差很大的数进行运算时,防止大数“吃掉”小数
  3. 简化计算步骤,减少运算次数