zoukankan      html  css  js  c++  java
  • 线性代数复习(4)------方阵的LU分解

    线性代数复习(4)------方阵的LU分解

    矩阵乘法的逆

    首先假设A,B为可逆矩阵 A − 1 A^{-1} A1 B − 1 B^{-1} B1为A与B的逆矩阵, ( A B ) − 1 = B − 1 A − 1 (AB)^{-1}=B^{-1}A^{-1} (AB)1=B1A1
    即: ( A B ) ( B − 1 A − 1 ) = A ( B B − 1 ) A − 1 = A E A − 1 = E (AB)(B^{-1}A^{-1})=A(BB^{-1})A^{-1}=AEA^{-1}=E (AB)(B1A1)=A(BB1)A1=AEA1=E(E为单位矩阵)

    逆与转置

    转置矩阵:将原先的矩阵进行变化 a i j = a j i a_{ij}=a_{ji} aij=aji所得到的矩阵,记为 A T A^T AT
    矩阵乘法的转置有 ( A B ) T = B T A T (AB)^T=B^TA^T (AB)T=BTAT
    逆的转置 ( A A − 1 ) T = E T = ( A − 1 ) T A T (AA^{-1})^T=E^T=(A^{-1})^TA^T (AA1)T=ET=(A1)TAT
    逆的转置等于转置的逆 ( A − 1 ) T = ( A T ) − 1 (A^{-1})^T=(A^T)^{-1} (A1)T=(AT)1

    LU分解

    上三角矩阵U:主对角线以下都是零的方阵称为上三角矩阵。
    下三角方阵L:主对角线以上都是零的方阵称为下三角矩阵。
    LU分解:将原有的矩阵A写成LU相乘的形式。 A = L U A=LU A=LU
    举个例子:
    [ 2 1 8 7 ] egin{bmatrix}2&1 \8&7\end{bmatrix} [2817]在实际转换为上三角形式中我们实际的矩阵变化就是下三角矩阵的逆矩阵。
    [ 1 0 − 4 1 ] egin{bmatrix}1&0\-4&1\end{bmatrix} [1401] [ 2 1 8 7 ] = egin{bmatrix}2&1 \8&7\end{bmatrix}= [2817]= [ 2 1 0 3 ] egin{bmatrix}2&1 \0&3\end{bmatrix} [2013] ( L − 1 A = U ) (L^{-1}A=U) (L1A=U)

    转换一下有:
    [ 2 1 8 7 ] = egin{bmatrix}2&1 \8&7\end{bmatrix}= [2817]= [ 1 0 4 1 ] egin{bmatrix}1&0\4&1\end{bmatrix} [1401] [ 2 1 0 3 ] egin{bmatrix}2&1 \0&3\end{bmatrix} [2013] ( A = L U ) (A=LU) (A=LU)

    我们也可以从右往左来试试:
    [ 2 1 8 7 ] egin{bmatrix}2&1 \8&7\end{bmatrix} [2817] [ 1 − 1 2 0 1 ] = egin{bmatrix}1&-frac 12\0&1\end{bmatrix}= [10211]= [ 2 0 8 3 ] egin{bmatrix}2&0 \8&3\end{bmatrix} [2803] ( A U − 1 = L ) (AU^{-1}=L) (AU1=L)

    转换一下有:
    [ 2 1 8 7 ] = egin{bmatrix}2&1 \8&7\end{bmatrix}= [2817]= [ 2 0 8 3 ] egin{bmatrix}2&0 \8&3\end{bmatrix} [2803] [ 1 1 2 0 1 ] egin{bmatrix}1&frac 12\0&1\end{bmatrix} [10211] ( A = L U ) (A=LU) (A=LU)

    把两个合并一下有:
    [ 2 1 8 7 ] = egin{bmatrix}2&1 \8&7\end{bmatrix}= [2817]= [ 1 0 4 1 ] egin{bmatrix}1&0\4&1\end{bmatrix} [1401] [ 2 0 0 3 ] egin{bmatrix}2&0 \0&3\end{bmatrix} [2003] [ 1 1 2 0 1 ] egin{bmatrix}1&frac 12\0&1\end{bmatrix} [10211] ( A = L D U ) (A=LDU) (A=LDU)
    中间的D矩阵为对角矩阵。
    简单的总结一下( E i j E_{ij} Eij表示对i行(或列)与j行(或列)之间的变换)
    E 21 E 31 E 32 A = U E_{21}E_{31}E_{32}A=U E21E31E32A=U => A = E 32 − 1 E 31 − 1 E 21 − 1 U A=E_{32}^{-1}E_{31}^{-1}E_{21}^{-1}U A=E321E311E211U

    置换矩阵

    这里我们主要注意是置换矩阵在矩阵消元时所起到的行交换的作用,以保证所求的主元不为0。
    即我们的LU分解时往往是这样的 P A = L U PA=LU PA=LU,其中P为若干的置换矩阵进行的行处理。
    置换矩阵的个数为: n ! n! n!(对映相依的阶数n的阶乘)(由每行列中1的排列可得)
    置换矩阵的转置与逆相等 P T = P − 1 P^{T}=P^{-1} PT=P1

    转置矩阵

    矩阵交换行列位置之后产生的矩阵为转置矩阵。
    A i j = A j i T A_{ij}=A_{ji}^T Aij=AjiT
    当转置运算之后矩阵不变时,这样的矩阵称为对称矩阵, 这里的对称指的是对于主对角线的对称。( A T = A A^T=A AT=A
    矩阵于其转置矩阵的乘积为对称矩阵。
    矩阵乘法的转置有: ( A B ) T = B T A T (AB)^T=B^TA^T (AB)T=BTAT
    于是对于上面的证明就有: ( A A T ) T = A A T (AA^T)^T=AA^T (AAT)T=AAT(转置偶数次对应不变),所以有矩阵与其转置矩阵的乘积为对称矩阵。

  • 相关阅读:
    UVa 1605
    UVa 120
    UVa 10384
    UVa 11694
    UVa 11846
    常用小函数
    【DP】:CF #319 (Div. 2) B. Modulo Sum
    类的无参方法
    类和对象
    七言
  • 原文地址:https://www.cnblogs.com/yanzs/p/13788245.html
Copyright © 2011-2022 走看看