4.矩阵的LU分解
转置:将矩阵中各个元素的行列互换(即将矩阵的行变为列,列变为行)得到新矩阵,新矩阵被称作原矩阵的转置矩阵,记作 A T A^{T} AT。当矩阵为方阵时,即将整个矩阵关于主对角线对称。
-
A A A的逆的转置等于 A A A的转置的逆,即 A A − 1 = I A A^{-1} = I AA−1=I,则 ( A − 1 ) T A T = I (A^{-1})^T A^T = I (A−1)TAT=I(其中 A A A为广义上的可逆矩阵)
证明: i i , j = r o w i o f A → ⋅ c o l u m n j o f A − 1 → = c o l u m n i o f A T → ⋅ r o w j o f ( A − 1 ) T → i_{i,j} = \overrightarrow{row\ i\ of\ A} \cdot \overrightarrow{column\ j\ of\ A^{-1}} = \overrightarrow{column\ i\ of\ A^T} \cdot \overrightarrow{row\ j\ of\ (A^{-1})^T} ii,j=row i of A⋅column j of A−1=column i of AT⋅row j of (A−1)T
因而 I = ( A − 1 ) T A T = A A − 1 I = (A^{-1})^T A^T = A A^{-1} I=(A−1)TAT=AA−1
-
假设 A A A在消元过程中不需要换行,已知 E A = U , A = L U EA = U , A = LU EA=U,A=LU,则 L = E − 1 L = E^{-1} L=E−1
例: [ 1 0 − 4 1 ] [ 2 1 8 7 ] = [ 2 1 0 3 ] E A U \begin{matrix} \begin{bmatrix} 1 & 0 \\ -4 & 1 \end{bmatrix} & \begin{bmatrix} 2 & 1 \\ 8 & 7 \end{bmatrix} & = & \begin{bmatrix} 2 & 1 \\ 0 & 3 \end{bmatrix} \\ E & A & & U \end{matrix} [1−401]E[2817]A=[2013]U
[ 2 1 8 7 ] = [ 1 0 4 1 ] [ 2 1 0 3 ] = [ 1 0 4 1 ] [ 2 0 0 3 ] [ 1 1 2 0 1 ] A L U L D U \begin{matrix} \begin{bmatrix} 2 & 1 \\ 8 & 7 \end{bmatrix} & = & \begin{bmatrix} 1 & 0 \\ 4 & 1 \end{bmatrix} & \begin{bmatrix} 2 & 1 \\ 0 & 3 \end{bmatrix} & = & \begin{bmatrix} 1 & 0 \\ 4 & 1 \end{bmatrix} & \begin{bmatrix} 2 & 0 \\ 0 & 3 \end{bmatrix} & \begin{bmatrix} 1 & \frac{1}{2} \\ 0 & 1 \end{bmatrix} \\ A & & L & U & & L & D & U \end{matrix} [2817]A=[1401]L[2013]U=[1401]L[2003]D[10211]U
在某些时候 L L L比 E E E能更加直观地展现信息
并不是所有可逆方阵都可以化为一个下三角阵和一个上三角阵的积的形式,只有在消元过程中不需要换行的才可以,否则都要在此基础上乘上一个置换矩阵
证明 L L L为下三角阵且主对角线元素为 1 1 1:
每次消元操作都是用某一行减去另一行的若干倍,被减行对应系数一直为 1 1 1,所以 E E E的主对角线元素一定为 1 1 1, 并且每次消元都是用某一行减去在它之前的另一行的若干倍,所以 E E E是一个下三角阵,又 L E = I LE = I LE=I,积的第 i i i行 第 i i i列为 1 1 1且该行和该列的其他元素为 0 0 0,所以用 L L L的行对 E E E的行进行线性组合时要确保 E E E中末尾连续 0 0 0个数较小 的行(即第 i i i行以后的行)对应的变量为 0 0 0且第 i i i行对应的变量为 1 1 1,这样 L L L就为下三角阵且主对角线元素为 1 1 1
即要使 L [ 1 0 ⋯ 0 ⋮ 1 ⋯ 0 ⋮ ⋮ ⋱ ⋮ ⋯ ⋯ ⋯ 1 ] = [ 1 0 ⋯ 0 0 1 ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ 1 ] L \begin{bmatrix} 1 & 0 & \cdots & 0 \\ \vdots & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ \cdots & \cdots & \cdots & 1 \end{bmatrix} = \begin{bmatrix} 1 & 0 & \cdots & 0 \\ 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & 1 \end{bmatrix} L 1⋮⋮⋯01⋮⋯⋯⋯⋱⋯00⋮1 = 10⋮001⋮0⋯⋯⋱⋯00⋮1 , L L L一定要满足需证明的结论,可以想象一下不满足会 怎么样
打赏
制作不易,若有帮助,欢迎打赏!