矩阵分解-Cholesky分解

  • 矩阵分解
    • 程序

矩阵分解

矩阵分解是将矩阵拆解为数个矩阵的乘积,可分为三角分解、满秩分解、QR分解、Jordan分解和SVD(奇异值)分解等。

Cholesky分解法是求解对称正定线性方程组最常用的方法之一。

程序

可采用如下命令:
R=chol(A):产生一个上三角阵R,使R’R=A。若A为非对称正定,则输出一个出错信息。

[R,p]=chol(A):这个命令格式将不输出出错信息。当A为对称正定的,则p=0,R与上述格式得到的结果相同;否则p为一个正整数。如果A为满秩矩阵,则A为一个阶数为q=p-1的上三角阵,且满足R’R=X(1:q,1:q)。

Cholesky分解后,线性方程组Ax=b变成R‘Rx=b,所以x=R(R’b)。

也可采用如下程序:
function [X]=Decomchol(A,b)
[N,N]=size(A);
X=zeros(N,1);
Y=zeros(N,1);
for i=1:N
A(i,i)=sqrt(A(i,i)-A(i,1:i-1)*A(i,1:i-1)’);
if A(i,i)==0
‘A is singular. no unique solution’
break;
end;
for j=i+1:N
A(j,i)=(A(j,i)-A(j,1:i-1)*A(i,1:i-1)’)/A(i,i);
end

A
b
前代法
for j=1:N
Y(j)=(b(j)-A(j,1:j-1)*Y(1:j-1))/A(j,j);
end
Y;
A=A’
for k=N: -1 :1
X(k)=(Y(k)-A(k,k+1:N)*X(k+1:N))/A(k,k);
end

矩阵分解-Cholesky分解相关推荐

  1. 矩阵的Cholesky分解

    原文链接:矩阵的Cholesky分解 首先来复习线性代数中几个重要的概念. 1)如果一个复矩阵A = A*(共轭转置),则A称为Hermitian矩阵.(注意,矩阵A转置后仍为其本身,显然A一定是方阵 ...

  2. 『矩阵论笔记』详细介绍矩阵的三角分解(LR分解)+平方根分解(Cholesky分解)

    详细介绍矩阵的三角分解(LR分解)+平方根分解(Cholesky分解)! 文章目录 一. 三角分解(LR分解) 1.1. 方阵的两个重要分解 1.2. 上(下)三角阵的性质 1.3. 三角分解的概念 ...

  3. 矩阵的Cholesky分解的Matlab实现

    版权声明:本文为博主原创文章,未经博主允许不得转载.https://blog.csdn.net/weixin_38451800/article/details/88933683 1.Cholesky分 ...

  4. Matlab与线性代数--矩阵的Cholesky分解

    本图文介绍了Matlab对正交矩阵的Cholesky分解操作.

  5. 三十分钟理解:矩阵Cholesky分解,及其在求解线性方程组、矩阵逆的应用

    写一篇关于Cholesky分解的文章,作为学习笔记,尽量一文看懂矩阵Cholesky分解,以及用Cholesky分解来求解对称正定线性方程组,以及求对称正定矩阵的逆的应用. 文章目录 直接Choles ...

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

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

  7. 线性代数笔记: Cholesky分解

    1 介绍 当一个实矩阵A是对称正定矩阵的时候,它可以分解成一个下三角矩阵L以及它的转置的乘积,即: 1.1 矩阵半正定的情况 如果矩阵是正定的话,那么L唯一确定:如果矩阵是半正定的话,那么也可以分解, ...

  8. C++实现Cholesky分解

    题目: 编制程序求解矩阵 A 的 Cholesky 分解,并用程序求解方程组 Ax=b,其中 代码实现: #include <iostream> #include <math.h&g ...

  9. 视觉SLAM中的数学——解方程AX=b与矩阵分解:奇异值分解(SVD分解) 特征值分解 QR分解 三角分解 LLT分解

    前言 本博客主要介绍在SLAM问题中常常出现的一些线性代数相关的知识,重点是如何采用矩阵分解的方法,求解线性方程组AX=B.主要参考了<计算机视觉--算法与应用>附录A以及Eigen库的方 ...

  10. 矩阵分解——三角分解(Cholesky 分解)

    (1)一个对角元素都是1的下三角矩阵,称为单位下三角矩阵. (2)上(下)三角矩阵的乘积仍是上(下)三角矩阵: (3)一般来说,矩阵的三角分解不唯一. (4)实对称正定矩阵 AA,Δk>0\De ...

最新文章

  1. npm 引用子项目模块_Java / Web项目中的NPM模块Browser-Sync
  2. Python获取2019-nCoV疫情实时追踪数据
  3. $dbms=mysql_Oracle dbms
  4. [原]用三行代码实现对音量的控制,实现增大,减小,静音
  5. SQL 之后,GQL 成为 ISO/IEC 国际标准数据库语言项目
  6. 软件工程概论第十三周学习进度
  7. 分享一个在线301跳转代码生成工具
  8. oracle 弱口令扫描器,超级弱口令检查工具
  9. 普及一下行业尖端知识——腾讯自研分布式数据库TBase
  10. 记录一次解决后端接口设置cookie设置不上去经过,一级域名可以设置上去cookie,二级域名设置不上cookie
  11. 四咪唑基四苯乙烯,四苯基吡嗪(tetraphenylpyrazine,TPP)
  12. python最终目标_Python晋级之路-工欲善其事必先利其器
  13. 21天战拖记——Day8:猴子法则and继续复习!(2014-05-11)
  14. 用python实现erp出入库_ERP采购入库成品入库流程(精)
  15. 计算机二级第十九套题电子表格,计算机二级MS-OFFICE套题知识点总结
  16. Python爬虫--爬取厦大电费
  17. python 存储bmp格式图片
  18. VM安装mac及dmg文件转换iso
  19. linux桌面版远程控制软件,ubuntu64位ARM版本飞腾可用
  20. 阿里云总裁胡晓明:保护客户数据隐私是第一原则

热门文章

  1. 超强PHP集成开发工具:phpDesigner v7.2.1 最新版下载(附注册机)
  2. 2012中国软件业务收入百强企业发布 (zz.IS2120)
  3. 执行SQL语句时出现问题操作必须使用一个可更新的查询错误的解决方法
  4. 山寨qq java_MTK手机QQjava版本山寨qq2009通用版
  5. 安全管家安卓_联想电脑的安全软件为什么更推荐用联想电脑管家?
  6. 研华工控机u盘启动安装linux系统,研华工控机610L如何安装win7系统
  7. m序列的原理以及verilog实现
  8. java中移动光标怎么写_java 程序当中无法做到光标跟随
  9. OSX 安装配置Maven
  10. adb shell中的am pm命令