解线性方程组地直接法和迭代法.docVIP

数值分析方法中方程求解的直接法和迭代法 第3章 解线性方程组的直接法 消元法 高斯消元法(加减消元)：首先将A化为上三角阵，再回代求解。 步骤如下： 第一步： 第二步： 类似的做下去，我们有： 第k步：。 n－1步以后，我们可以得到变换后的矩阵为： 注意到，计算过程中处在被除的位置，因此整个计算过程要保证它不为0。 所以，Gauss消元法的可行条件为：。 就是要求A的所有顺序主子式均不为0，即 因此，有些有解的问题，不能用Gauss消元求解。 另外，如果某个很小的话，会引入大的误差。 列主元消元法 在Gauss消元第k步之前，做如下的事情： 若交换k行和j行 其实就按行变换，不改变方程组的解，同时又有效地克服了Gauss消元地缺陷。 Gauss-Jordan消元法 将在Gauss消元第k步，变为（遇上=0的则需要换一行）： 将该行上、下三角地部分都变为0，即对角阵，同时主元系数也要消，这样就不用回代了。 直接分解法 Gauss消元法的第k步： 从矩阵理论来看，相当于左乘矩阵： 。 因此，整个Gauss消元法相当于左乘了一个单位下三角阵。 所以有，L为单位下三角阵，U为上三角阵。 因此 我们可以通过2次反代过程求解方程组，分解的理论由Gauss消元得出，因此分解能够进