数值分析-第5章-线性方程组的直接方法.pptVIP

例5.13 线性方程组 的系数矩阵带误差，成为如下方程组 求方程组系数矩阵的条件数, 并说明方程组的性态 解 因为 所以 因此方程组是良态的 5.7.2 精度分析 求得方程组Ax=b的一个近似解以后,希望判断其精度，检验精度的一个简单办法是将近似解再回代到原方程组去求出余量r. r = b-A 如果r很小，就认为解是相当精确的。 定理5.14 设 是方程组Ax=b的一个近似解,其精确解记为 ,r为 的余量。则有 证明略 例5.14 设A为正交矩阵，证明：cond2(A)=1 分析：由正交矩阵和条件数的定义便可推得 解：因为A是正交矩阵， 故ATA= AAT=I， A-1= AT，从而 例9 求Hilbert矩阵H3的条件数. Hilbert 阵 cond (H2)? = 27 cond (H3)? ? 748 cond (H6)? = 2.9 ? 106 cond (Hn)?? ? as n ? ? 注：一般判断矩阵是否病态，并不计算A?1，而由经验出。 ? 行列式值很大或很小（如某些行、列近似相关）； ? 元素间相差大数量级，且无规则； ? 主元消去过程中出现小主元； ? 特征值相差大数量级。 小结 本章介绍了解线性方程组的直接法。直接法是一种计算量小而精度高的方法。直接法