第十章 线性方程组的数值解法.docVIP

数值计算中应注意的几个原则 关于数值稳定性的算法。 一个程序往往要进行大量的四则运算才能得出结果,每一步的运算均会产生舍入误差。在运算过程中,舍入误差能控制在某个范围内的算法称之为数值稳定的算法,否则就称之为不稳定的算法。 例: , 用分部积分公式得递推式:。 用四位有效数字计算: , , , ,, ,, ,. 可以估计出 故 ,。 于是与精确值已经面目全非,一位有效数字也没有。这是由于如果有误差,不计中间再产生的舍入误差,该误差随着计算过程分别乘以,到时已经变成了,误差扩大了4万倍。因而该算法不是稳定的。 如果递推式改为,由,逐步计算直到。计算结果有四位有效数字,如果有误差,其传播到所引起的误差仅为。故该算法是稳定的。 注意避免两个相近数的相减。 两个相近的数相减,有效数字会大大损失。前面一个例子已说明问题,这里再举一例: ,如用四位有效数字计算: ,结果只有一位有效数字;如改为: ,有四位有效数字。新算 法避免了两个相近数的相减。 简化计算步骤,减少运算次数。 例: 计算多项式的值: ,每取有次乘法运算,因此计算共需次乘法和加法运算。如将写成,用递推算法:,最终,共需次乘法和次加法运算。 一般地要注意,能在循环外计算,就不要放在循环内计算。 4.避免除数的绝对值远小于被除数的绝对值。 ,当时,