码迷,mamicode.com
首页 > 其他好文 > 详细

Lecture 4:LU分解

时间:2021-03-02 12:15:23      阅读:0      评论:0      收藏:0      [点我收藏+]

标签:math   没有   amp   进入   情况   end   app   矩阵分解   需要   

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

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

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

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

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

\[\begin{pmatrix} 2 & 1 \8 & 7 \end{pmatrix} = \begin{pmatrix} 1 & 0 \4 & 1 \end{pmatrix} \begin{pmatrix} 2 & 0 \0 & 3 \end{pmatrix} \begin{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\)次操作

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

Lecture 4:LU分解

标签:math   没有   amp   进入   情况   end   app   矩阵分解   需要   

原文地址:https://www.cnblogs.com/miraclepbc/p/14465551.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!