博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
矩阵分析-线性系统-3 LU分解
阅读量:4449 次
发布时间:2019-06-07

本文共 763 字,大约阅读时间需要 2 分钟。

1. LU分解

LU分解在本质上是的一种表达形式。实质上是将A通过初等行变换变成一个上三角矩阵,其变换矩阵就是一个单位下三角矩阵。这正是所谓的杜尔里特算法(Doolittle algorithm):从下至上地对矩阵A做初等行变换,将对角线左下方的元素变成零,然后再证明这些行变换的效果等同于左乘一系列单位下三角矩阵,这一系列单位下三角矩阵的乘积的逆就是L矩阵,它也是一个单位下三角矩阵。

 

下面以具体例子来说明。若AX=b是一个非奇异系统,那么高斯消元法将A化简为一个上三角矩阵。若主轴上没有0值,则无需交互行,因此只需进行第3类初等行变换(把第 i 行加上第 jk 倍)即可完成此变换。例如

                                                                             (1)

第3类行变换可以通过左乘相应的初等矩阵实现,对上例来说进行的3个变换就是相应初等矩阵的乘积。注意最右边是一个下三角矩阵L

                                  =L            (2)

从而有,即。因此A=LU,为一个下三角与一个上三角矩阵的乘积,因此称为LU分解。

注意

1)U是高斯消元的结果,且对角线上是主元

2)L对角线上是1,对角线下面的元素恰恰是在式1中用于消去(i,j)位置上元素的乘子。

 

2. LDU分解

LU分解存在不对称,因为L矩阵的主对角线元素为1,而U矩阵主对角线不为1。为了对此进行弥补,可U矩阵进一步做如下分解

                     = DU

这就是LDU分解,其中L和U是单位三角阵(主对角线都为1),D是对角阵。

 

下图说明了的LDU分解 

 

                             

 

 

3. LDLT分解

当A是对称阵时,LDU分解为A=LDLT 。例子如下

对称阵   , LU分解为, LDLT 分解为。

 

进一步,由于对角阵D的元素为正,可有,从而得到,其中U是上三角阵。

转载于:https://www.cnblogs.com/pegasus/archive/2011/11/21/2257258.html

你可能感兴趣的文章
inux中bin与sbin目录的作用及区别介绍
查看>>
USACO 3.1 Contact
查看>>
Office之什么是高内聚低耦合
查看>>
一些奇怪的问题求回答
查看>>
这些年踩过的坑
查看>>
iOS开发拓展篇——如何把项目托管到GitHub
查看>>
性能优化之数据库优化
查看>>
类的继承、菱形继承、派生、多态
查看>>
mysql约束
查看>>
javascript鼠标及键盘事件总结及案例
查看>>
mysql表之间的关系及级联操作
查看>>
mac 搭建virtualenv的那些坑
查看>>
多路复用IO模型
查看>>
并发、串行、并行及多道技术原理
查看>>
hashlib、pickle、hmac、logging模块使用
查看>>
javascript常用知识点总结
查看>>
2019秋招复习笔记--数据库基本操作
查看>>
2019秋招复习笔试--手写代码
查看>>
2019秋招复习笔记--智力题
查看>>
MySQL学习笔记
查看>>