方程组的误差估计.ppt

第四章方程组的直接解法,4.3.2 方程组的误差估计,4.3.1 矩阵的条件数,4.3 方程组的性态和误差估计,解Ly=b得y=(6,1,-1),T,解L,T,x=D,-1,y得x=(2,1,-1),T,4.3.1 矩阵的条件数,定义4.1,如果方程组Ax=b中，矩阵A和b右端的微小变化，引起解向量x的很大变化，则称A为关于解发才组和矩阵求逆的,病态矩阵,，称相应的方程组为,病态方程组,否则，称A为,良态矩阵,，称相应的方程组为,良态方程组,的准确解是（1，1）,T,若A及b作微小的变化，考扰动后的方程组,其准确解为（-2，10）,T,先看一个例子，说明方程组的解对或的扰动的敏感性问题例4.9 方程组,在上例中，A和b的微小变化引起x很大的变化，x对A和b的扰动是敏感的这种现象的出现完全是有方程组的性态决定的方程组右端所引起的解向量的相对误差就可能越大可见，量,是相对误差,的倍增因子，该量越大，,又由于，,x=,A,-1,b,，,即得,我们需要一种能刻画矩阵和方程组病态程度的标准暂且不考虑矩阵A的扰动，仅须考虑b 的扰动对方程组的影响，设方程组Ax=b的扰的方程组为A（x+x）=b+b，则,为矩阵A的条件数。

如果矩阵范数取2范数，则记,定义4.2,设AR,nn,为可逆矩阵，按算子范数，称,cond（A）=,同样可以定义cond,（A）和cond,1,（A）按（4.3.1），,（1）其中cond（A）1，cond（A）=cond（A,-1,），cond（A）=cond（A），其中 R，0,（2）若U为正交矩阵，即U,T,U=I则,cond,2,（U）=1，cond,2,（A）=cond,2,（AU）=cond,2,（UA）3）设 与 为A按绝对值最的和最小的特征值，则,若A对称，则cond,2,（A）=,cond（A）,设A,-1,存在，条件数有如下一些性质：,例4.10 下列Hilbert矩阵是一族著名的病态矩阵：,它是一个nn的对称矩阵，可以证明是正定的计算条件数有cond,2,（H,4,）=1.5514 10,4,，cond,2,（H,6,）=1.4951 10,7,，cond,2,（H,8,）=1.525 10,10,由此可见，随着n的增加，H,n,的病态可能越严重H,n,常常在数据拟合和函数逼近中出现对于实际问题，条件数一般是很难计算的下列现象可能表示方程组Ax=b是病态的1）如果矩阵A的按绝对值最大特征值和最小特征值之比很大，则A是病态的。

2）如果系数矩阵A的元素间数量级很大，并且无一定规则，则A可能病态3）如果系数矩阵A的莫些行或列是近似相关的，或系数矩阵的行列式值相对说很小，则A可能病态4）如果在A的三角化过程中出现小指元或采用选用选主远技术，主元素数量级相差悬除时，则A可能病态对于病态方程组，数值求解必须小心进行，否则达不到所要求的准确度有时可以用高精度（如双精度或扩充精度）的运算，以改善或减轻方程组的病态程度，有时也可以对圆方程组作预处理，以降低,系数矩阵的条件数，即选择非奇异矩阵P和Q，一般选着为对角阵或三角矩阵，使,cond（PAQ）cond(A),然后，求解等价方程组PAQ y=P y，y=Q,-1,x例如，对矩阵,有cond,10,5,若进行预处理,则cond,（B）=4，条件数的改善4.3.2 方程组的误差估计,定理4.9,设Ax=b，A为非奇异矩阵，b为非零向量，A和b分别有扰动，A和 b,若 1,则有误差估计式,由于舍入误差，我们解方程组往往得到的是近似解下面利用条件数给出近似解的事前误差估计，即计算之前和计算之后的误差估计将上式两端取范数，则有,证.将代入扰动方程组,，整理后有,经整理后得,由于,，即得所证。

再利用,，则有,若,时，则由（）有,定理4.10,设Ax=b，b0，则对方程组的近似解,有误差估计式,该定理说明，当cond(A)很大时，即使方程组余量r 的相对误差已经很小，但近似解的相对误差仍然可能很大如果用直接解法得到的近似解 误差很大，我们可以用迭代改善的办法对近似解进行修正设r=b-A ，x为修正量，为新的近似解这样，我们可以通过求解 A x=r得到 ，显然，在准确运算下有,（4.3.3）,=A（x-）,证 由Ax=b有 r=Ax-A,又由x=A,-1,b，有,定理得证其中r=b-A 为剩余向量方程组直接解法的稳定性是应当注重的问题如果通过直接计算每一步设入误差对解的影响来获得近似解的误差界，那将是非常困难的等人提出了“向后误差分析法”，其基本思想是把计算过程中设入误差对解的影响归结为原始数据对解的影响下面给出一个定理来说明这方面的结果然而，再实际计算时，方程组（4.3.3）不大可能求解，所以解（4.3.3）只能提供有限的修正因此，需要反复求解为（4.3.3）的方程组，不断对所得的近似解进行改进这种近似值逐进接近真解的过程称为迭代解法为了节省计算量，可事先对矩阵A进行LU分解，把反复解形为（4.3.3）的方程组改为反复解形为Ly=r，Ux=y 的方程组。

为了保证计算精度，计算剩余向量r可采用高精度计算定理4.11,设AR,nn,，A为非奇异矩阵，用列主元法或全主元法解方程组,Ax=b，其计算解,满足,记计算机尾数字长,是消去过程中A,(k),中的元,为t，且n2-t 0.01素，则有,该定理说明，矩阵A的阶数越高条件数越大矩阵元素的增长因子 越大和计算机字长越短，则舍入误差对解的影响越严重因此计算精度取决于矩阵的规模、方程组的性态、所选取的算法和计算机字长1）若A的LU分解计算结果为,则有,（2）,（3）计算解有精度估计：,。