《计算方法》课件1第1章.pptVIP

1.2.4数据误差在运算中的传播

设x*、y*分别是初始数据x、y的近似值，即



其中，e(x*)、e(y*)分别为x*、y*的绝对误差。下面考察用x*、y*代替x、y时，函数值z=f(x，y)会产生怎样的误差。

假设e(x*)、e(y*)的绝对值都很小,而且函数z=f(x，y)可微，记此时z的近似值为z*=f(x*，y*)，则有



上式可近似地表示为

(1.11)而且

(1.12)从(1.11)式容易得到，在进行数值运算时，初始数据误差与计算结果产生的误差之间有以下关系：

(1)若f(x,y)=x±y，则

(1.13)

(2)若f(x,y)=xy,则

(1.14)

(3)若，则

(1.15)从(1.12)式容易得到以下关系:

(1)若f(x,y)=x±y,则

(1.16)

(2)若f(x,y)=xy,则

(1.17)

(3)若，则

(1.18)1.3数值计算中应注意的问题

数值计算中的误差分析是一个既复杂而又不可避免的问题。为了保证数值计算结果的正确性，就必须把误差控制在比较小的范围内，也就是要求掌握误差产生、传播的规律。

可是实际上，目前尚无有效的方法对误差做出定量的估计，所以，在解决具体数值问题时，往往首先需要对该问题做出定性的分析。1.3.1算法的数值稳定性

如果在用某种算法进行数值计算过程中，舍入误差在一定条件下能够控制在一定范围内(即舍入误差的增长不会影响结果的可靠性)，则称该算法是数值稳定的,否则称该算法为数值不稳定的。

例1.3对n=0，1，2，…，计算积分。

解由于取I0=ln1.2≈0.182，用如下递推公式计算:

①(n=1,2,…,8)(1.19)

得到(精确到小数点后第3位)令εn=In-I*n，则有εn=-5εn-1。可见，若从I*n-1计算I*n，其误差将以每步5倍的速度增长，从而有

这表明计算公式①是数值不稳定的。

现在换一种方案，将(1.19)式倒过来算。取I*9=0.018，得到递推公式:

②(n=8,7,…,1)



计算可得且此时

I*0≈0.182

由此可见,利用公式②得到的结果与精确值比较相近，这是因为虽然e(I*9)较大，但利用公式②计算时，每步误差将缩小为上步误差的1/5，故算法②是数值稳定的。

此例说明，在实际数值计算中，数值方法不稳定的算法是不能使用的。一个数值算法，如果对任何允许的初始值都稳定，则称此算法为无条件稳定的；若对某些初始值稳定，而对另一些初始值不稳定，则这样的算法称为条件稳定的。1.3.2避免误差危害的若干原则

数值计算中，除了要分清算法是否稳定外，还应尽量避免误差危害，防止有效数字的损失，下面给出若干原则。

1.避免两相近数相减

由(1.16)式看到，两相近的数字相减时，会将计算结果的相对误差变得很大，称此时的精度损失为相减相消。为了防止这种现象产生，最好是改变数值计算方法。例1.4求二次方程ax2+bx+c=0(a≠0)的两个根。

解求根x1,x2的公式为

，

若b0，且b2-4ac0，计算x1时就会产生