LU分解线代讨论课.pptVIP

关于矩阵的LU分解 刘强 （2010001040005） 指导老师：黄廷祝 邵晋梁 题外话 题：将矩阵 分解为对角元为1的下 三角矩阵L和上三角矩阵U的乘积，即A=LU。 HOW TO REALIZE IT ? 定义: 将矩阵A分解为对角元为1的下三角矩阵L和上三 角矩阵U的乘积叫做矩阵A的LU分解。 其中 i=1,2,……，n) 现在的问题在于如何确定 ( )和 ( )的值 方法一 高斯消元 构造矩阵 ，其中 且 。我们考虑到对矩阵A作行初等变换化为行阶梯 矩阵的过程中，若令 则 以此类推 令 则可得 我们令 则分解完成 L矩阵的元为 U矩阵的元为 但有些情况下，某些矩阵是不可用此方法分解的，比如 A= 。要其能够分解，由上面分解过程可以看出， 要求 （ ），其实质是 （ 其中 ）。 即在所有A（k） 可逆的情况下，矩阵A能进行LU分解，但这只是一个充分条件。 方法二 解方程 即根据A=LU，先求U矩阵第1行元，再求L矩阵第1列元，再求U矩阵第2行元，再求L矩阵第2列元……依此类推。比如： 类推可得 求解时只需分别令 ，2 ，3 …… 即可。并规定 。 下面讨论A矩阵可分解的条件： 1．先前已经说明， （ ）时，即在所有 可逆的情况下，其可以通过方法一不断作行变换，进而分解。如果能分解，由方法二可知，其分解必唯一（L矩阵和U矩阵的每一个元的解是唯一的）。 2．若存在 不可逆，而 可逆，则矩阵A必不能分解。 证明：用定义 的方式类似地定义 ， 。假设矩阵 能分解为 ，易知 。因为 不可逆， 可逆，易知 必不可逆。故 中某对角元 。则矩阵 也不可逆，推出 不可逆，矛盾。故原命题得证。 3。最后一种情况比较复杂：即 在 时可逆且在 时皆不可逆( )，或者所有 均不可逆。此时矩阵A可能可逆，也可能不可逆。鉴于时间关系，没能作出更深入的思考。 对于非方阵情况，可将其人为地添加行数或列数使之成为方阵后分解，再作拆开即可，不多作赘述。 谢谢观赏！