• Lecture 4:LU分解


    对矩阵(A),有(E_{21}A = U)

    [egin{pmatrix} 1 & 0 \ -4 & 1 end{pmatrix} egin{pmatrix} 2 & 1 \ 8 & 7 end{pmatrix} = egin{pmatrix} 2 & 1 \ 0 & 3 end{pmatrix} ]

    现在要转换成(A = LU)的形式,这里的(L)即为(E_{21}^{-1})

    [egin{pmatrix} 2 & 1 \ 8 & 7 end{pmatrix} = egin{pmatrix} 1 & 0 \ 4 & 1 end{pmatrix} egin{pmatrix} 2 & 1 \ 0 & 3 end{pmatrix} ]

    这种分解形式相当于将一个矩阵分解为一个对角线全为(1)的下三角矩阵与一个对角线为主元的上三角矩阵的乘积。
    另外一种分解方式是(A = LDU),其中(D)是一个对角线为主元的对角矩阵(diagonal),如:

    [egin{pmatrix} 2 & 1 \ 8 & 7 end{pmatrix} = egin{pmatrix} 1 & 0 \ 4 & 1 end{pmatrix} egin{pmatrix} 2 & 0 \ 0 & 3 end{pmatrix} egin{pmatrix} 1 & frac{1}{2} \ 0 & 1 end{pmatrix} ]

    一般情况下,在没有行交换的情况下,所有的乘数都直接进入到(L)当中了。

    (n)阶矩阵(A)通过初等行变换,变为(U),总共大概需要(n^2 +(n - 1)^2 + dots ,1^2 approx frac{1}{3}n^3)次操作

    (LU)分解一个比较重要的意义在于加速(多个)线性方程组的求解。
    首先将(A)分解为(LU)的形式,然后再计算(LUx = b)
    计算(LUx = b),可以分为两步,即:(Ly = b, Ux = y)。这两个线性方程组的求解都是(O(n^2))的。

    置换矩阵(P)(P^{-1} = p^T)(n)阶置换总共有(n!)种可能性。

  • 相关阅读:
    JAVA 面向对象的扩展 内部类
    对于win10 更换JDK后 查询JDK路径还是原路径的解决办法
    懂得的懂
    稀疏数组转化二维数组
    Flume的安装配置
    CentOS7配置ip和ssh免密登录和hadoop环境
    AOP的使用和事务
    spring的个人理解
    单车月结算2-修改和删除功能
    单车月结算1
  • 原文地址:https://www.cnblogs.com/miraclepbc/p/14465551.html
Copyright © 2020-2023  润新知