一点点数学!开干!

参考书籍:《矩阵分析与计算》李继根 张新发编著

矩阵的LU分解:

LU分解定理:如果n阶方阵A的各阶顺序主子式≠0(K=1、2、3,…,n),即A的各阶顺序主子式矩阵都可逆,则存在唯一的单位下三角矩阵L与唯一的非奇异上三角矩阵U,使得A=LU

上面这个LU分解定理就是判断一个矩阵是否能LU分解的判断条件

下面用一个实例说明如何对一个矩阵进行LU分解:

一点讲解:也许你会有疑问,为什么

其实如果你线性代数学的扎实是很容易理解的,因为并起来然后进行行变换,右边这个单位矩阵I实际就相当于记录了A是经过怎样的行变换变成U的,也就是说

那么问题又来了,如果一个矩阵A不满足所有顺序主子式都不等于零,那么它就一定不能LU分解了吗?
一个矩阵A如下

显而易见,这个矩阵A的二阶顺序主子式为0,所以A不能进行LU分解,那么遇见这种情况我们是不是可以交换矩阵A的第二行和第三行,然后再验算一下,是不是就发现A的所有顺序主子式都不为零,可以进行LU分解了!

这就是带行变换矩阵P的LU分解,定义如下

带P的LU分解说白了就是将可逆方阵A重新排列,使其满足LU分解的条件,进行LU分解。

同时我们对带P的LU分解提出一个定理 (因为可以通过矩阵P重新排列行,所以LU分解的条件放宽成了矩阵A可逆)

简单来说只要矩阵A可逆,那么就存在带P的LU分解。那么我们现在对A进行带P的LU分解吧:

如果细心的小朋友估计会对这个带P的LU分解提出疑问,那就是为啥要交换第二行和第三行,交换第一行和第三行不也可以嘛!确实可以,而且基于这个问题我们还提出了一个方法,那就是列选主元法!

什么是列选主元法?

列选主元法:在每一步消元前,在主元所在的一列的对角线以及对角线下方的元素选取绝对值最大的元素作为主元素(就是把第一列绝对值最大的那一行放在第一行)。利用置换矩阵处理行变换。列选主元法可保证消去时的乘子都不超过1,可以抑制数据误差的传播。

所以对于上面这个带P的LU分解如果考虑列选主元法我们应该交换第一行和第三行,而不是交换第二行和第三行!

当然LU分解还有一些推广形式,比如说LDU分解

LDU分解就是通过对角矩阵D将U进一步简化为单位上三角矩阵。这里就不展开了,很简单。

矩阵的LU分解有何用?

LU 分解的意义在于,将矩阵的分解与方程的求解分离。

在解线性方程组的时候Ax=b 与求解LUx=b(L和U都是三角矩阵很容易求解)

两者求解的时间复杂度是不一样的

前者的时间复杂度是0的3次,后者的时间复杂度是0的2次,不是一个数量级的

矩阵的LU分解,LU分解的推广,LU分解有什么意义,为什么要用LU分解。相关推荐

  1. 分解得到的时频域特征_【推荐文章】基于变分模态分解和广义Warblet变换的齿轮故障诊断...

    <机械传动>2018年  第42卷   第7期 文章编号:1004-2539(2018)07-0157-05 DOI:10.16578/j.issn.1004.2539.2018.07.0 ...

  2. 机器学习(十一)——机器学习中的矩阵方法(1)LU分解、QR分解

    http://antkillerfarm.github.io/ 因子分析的EM估计(续) 去掉和各参数无关的部分后,可得: ∑i=1mE[logp(x(i)|z(i);μ,Λ,Ψ)]=∑i=1mE[1 ...

  3. crout分解计算例题_矩阵与数值计算(2)——矩阵三角分解LU、PALU、Cholesky三角分解、QR分解...

    前言 矩阵分解是设计算法的主要技巧,通过分解可以将复杂问题转化为几个简单问题求解,通常完成这一转化任务的主要技巧就是矩阵分解.例如,我们知道上三角矩阵和下三角矩阵是容易求解的,或者对角矩阵是最理想的求 ...

  4. lu分解法matlab_MIT 18.065—机器学习中的矩阵方法02 矩阵乘法与矩阵分解

    数据分析.信号处理和机器学习中的矩阵方法 第02讲 矩阵乘法与矩阵分解 新MIT 线性代数|机器学习(中英机翻字幕)18.065 by Gilbert Strang_哔哩哔哩 (゜-゜)つロ 干杯~- ...

  5. [Math]常见矩阵分解及复杂度 Cholesky QR LU

    本篇尝试介绍几个常见的矩阵分解及其数值方法,探究其算法复杂度及不稳定性来源,能够让我们更好的理解工具,使用工具.常见的矩阵分解可以在维基中找到 文章目录 Cholesky 定义 数值方法 数值方法的数 ...

  6. matlab lud矩阵分解,MIT线性代数总结笔记——LU分解

    MIT线性代数总结笔记--LU分解 矩阵分解 矩阵分解(Matrix Factorizations)就是将一个矩阵用两个以上的矩阵相乘的等式来表达.而矩阵乘法涉及到数据的合成(即将两个或多个线性变换的 ...

  7. 几种矩阵分解算法: LU分解,Cholesky分解,QR分解,SVD分解,Jordan分解

    目录 1.LU分解 2. LDLT分解法 3. Cholesky分解的形式 4. QR分解 5.SVD分解 5.1 SVD与广义逆矩阵 6. Jordan 分解 参考文章: ---------我只是搬 ...

  8. 矩阵的三角分解(LU分解)

    矩阵的三角分解将矩阵分解成一个下三角矩阵L和一个上三角矩阵U的乘积. 定义:如果n阶矩阵A能够分解成一个下三角矩阵L和一个上三角矩阵U的乘积,则称这种分解为三角分解或LU分解,如果n阶矩阵A能够分解为 ...

  9. Cholesky和LU矩阵分解

    1.Cholesky分解 在线性代数中,矩阵分解是将矩阵分解为矩阵的乘积.有许多不同的矩阵分解.其中之一就是Cholesky分解. Cholesky分解是将Hermitian正定矩阵分解为下三角矩阵及 ...

最新文章

  1. 服务器php将视频转为m3u8,php应用ffmpeg对mp4文件转换并生成m3u8视频流文件
  2. Windows Server 2012关闭Server Manager开机自启动
  3. 我的第一个appium+Python自动化实例
  4. python—操作MySQL增加数据
  5. php 获取权限,php 获取权限列表函数_PHP教程
  6. 腾讯重磅发布全栈机器学习平台Angel 3.0
  7. oracle新增表单,Oracle Retail(Retek)中Form生成列表方法总结
  8. 想Get热搜同款?GitHub开源神器让父亲重返18岁!
  9. python写入指定行_python文件操作如何写在指定的行
  10. exception java doc_在javadoc中,標簽@throws和@exception之間有什么區別?
  11. day26-2 基于TCP协议的套接字编程
  12. freebsd mysql tmp_FreeBSD下安装MySQL与MySQLdb笔记
  13. java oracle11g jar_oracle11g驱动jar包下载
  14. 切比雪夫多项式MATLAB绘图
  15. 用 MeGUI 压制 DVDrip 入门
  16. uva11045(最大二分图匹配)
  17. 电脑误删的文件怎么恢复?分享90%的人都会的这2招
  18. 让我们看看往年的1024传智播客怎么安排的
  19. 推广网店的12个秘诀
  20. 关于字符串的输入几种方法

热门文章

  1. iOS苹果相册“地点”相簿效果实现浅谈
  2. 第07篇 部署区块链浏览器(2)
  3. houdini环境变量服务器文件读不了,Windows下在普通命令行窗口里初始化Houdini环境...
  4. Spring之ApplicationContext介绍
  5. 固态硬盘使用的优化——适用于任何版本的固态
  6. 私钥,公钥,钱包地址,助记词,keyStore的区别
  7. java、C语言数组反转
  8. 水下动作捕捉的两种实现方式
  9. 数学通道的应用(十五)-电阻的测量
  10. freeswitch + dbl 呼入设置