10 矩阵的三角分解.docVIP

第十讲 矩阵的三角分解 一、 Gauss消元法的矩阵形式 n元线性方程组 设,设A的k阶顺序主子式为，若， 可以令 并构造Frobenius矩阵 计算可得 初等变换不改变行列式，故，若，则，又可定义 ，并构造Frobenius矩阵 依此类推，进行到第(r-1)步，则可得到 则A的r阶顺序主子式，若，则可定义，并构造Frobenius矩阵 直到第（n－1）步，得到 则完成了消元的过程 而消元法能进行下去的条件是 二、 LU分解与LDU分解 容易求出 为下三角矩阵 令为上三角矩阵，则 （L: lower U: upper L: left R: right） 以上将A分解成一个单位下三角矩阵与上三角矩阵的乘积，就称为LU分解或LR分解。 LU分解不唯一，显然，令D为对角元素不为零的n阶对角阵，则 可以采用如下的方法将分解完全确定，即要求 L为单位下三角矩阵 或 U为单位上三角矩阵 或 将A分解为LDU，其中L，U分别为单位下三角，单位上三角矩阵，D为对角阵，而 。 n阶非奇异矩阵A有三角分解LU或LDU的充要条件是A的顺序主子式 n个顺序主子式全不为零的条件实际上是比较严格的，特别是在数值计算中，很小时可能会带来大的计算误差。因此，有必要采取选主