Ch5解线性方程组的直接方法导论.ppt

下面再分析H3的情况 (2) 考虑方程组 (1) 计算H3的条件数cond(H3)∞ ||H3||∞=11/6, ||H3-1||∞=408, 所以cond(H3)∞=748. 同样可计算cond(H6)∞=2.9×107, cond(H7)∞=9.85×108. 由此看出当n越大时, Hn矩阵病态越严重. H3 x=(11/6, 13/12, 47/60)T=b, 设H3及b都有微小误差(取3位有效数字)有 简记为(H3+?H3)(x+?x)=b+?b. 方程组H3x=b与(5.8)的精确解分别为x=(1, 1, 1)T, x+?x=(1.089512538, 0 1.491002798)T. 于是 ?x=(0.0895, -0.5120, 0.4910)T, 这就是说H3与b相对误差不超过0.3%, 而引起解的相对误差超过50%. 由上面的讨论, 要判别一个矩阵是否病态需要计算条件数cond(A) =||A-1||||A||,而计算A-1是比较费劲的, 那么在实际计算中如何发现病态情况呢? (1) 如果在A的三角约化时(尤其是用主元素消去法解(5.3)时)出现小主元, 对大多数矩阵来说, A是病态矩阵, 例如用选主元的直接三角分解法