Armadillo C++ linear algebra library 学习笔记(5)——矩阵的分解

原创 2015年07月09日 22:33:03

1、矩阵的SVD分解

A、通过调用”svd(U,S,V,X)”函数进行矩阵X的奇异值分解(svd)。X = USV’
B、示例

#include <iostream>
#include <armadillo>
using namespace arma;
int main()
{
    //1、产生随机矩阵A,大小为5x5,每个数的范围为:(0,10)
    mat A = randu<mat>(5, 5)*10;
    //2、对矩阵A进行svd分解,USV' = A
    mat U,V;//U,V均为正交矩阵
    colvec S;//S为奇异值构成的列向量
    svd(U,S,V,A);//进行奇异值分解
    A.print("原矩阵A:");
    U.print("正交矩阵U=\n");
    S.print("奇异值列向量S=\n");
    V.print("正交矩阵V=\n");
}

结果
这里写图片描述

2、矩阵的Cholesky分解

A、通过调用”chol(X)”函数进行矩阵X的Cholesky分解。X = R.t()*R
B、示例

#include <iostream>
#include <armadillo>
using namespace arma;
int main()
{
    //1、产生随机矩阵A,大小为5x5,每个数的范围为:(0,10)
    mat B = randu<mat>(5, 5)*10;
    mat A = B.t()*B;//A为对称矩阵
    mat R = chol(A);//对矩阵A进行Cholesky分解,R.t()*R = A,矩阵A必须是对称矩阵;
    A.print("原对称矩阵A:");
    R.print("矩阵ACholesky分解后的R为:\n");
}

结果
这里写图片描述

3、矩阵的QR分解

A、通过调用”qr(Q, R, X)”函数进行矩阵X的QR分解。X = QR
B、示例

#include <iostream>
#include <armadillo>
using namespace arma;
int main()
{
    //1、产生随机矩阵A,大小为5x5,每个数的范围为:(0,10)
    mat A = randu<mat>(5, 5)*10;
    mat Q,R;
    qr(Q,R,A);//QR=A
    A.print("原矩阵A:");
    Q.print("QR分解后的矩阵Q为:\n");
    R.print("QR分解后的矩阵R为:\n");
}

结果
这里写图片描述

4、矩阵的LU分解

A、通过调用”lu(L, U, P, X)”函数进行矩阵X的LU分解。PX = LU
B、示例

#include <iostream>
#include <armadillo>
using namespace arma;
int main()
{
    //1、产生随机矩阵A,大小为5x5,每个数的范围为:(0,10)
    mat A = randu<mat>(5, 5)*10;
    mat P,L,U;
    lu(L,U,P,A);//PA = LU
    A.print("原矩阵A:");
    L.print("LU分解后的矩阵L为:\n");
    U.print("LU分解后的矩阵U为:\n");
    P.print("LU分解后的矩阵P为:\n");
}

结果
这里写图片描述

矩阵奇异值分解简介及C++/OpenCV/Eigen的三种实现

矩阵奇异值分解简介及C++/OpenCV/Eigen的三种实现
  • fengbingchun
  • fengbingchun
  • 2017年06月03日 19:00
  • 2162

奇异值分解 SVD C++代码

/** 函数原型: bool svd(vector > A, int K, std::vector > &U, std::vector &S, std::vector > &V); 输入 待分解矩阵A...
  • qxessay
  • qxessay
  • 2014年10月30日 20:17
  • 1815

矩阵的奇异值分解

定义设A∈Cm×nA\in C^{m\times n},则矩阵AHAA^{H}A的nn个特征值λi\lambda _i的算术平方根δi=λi−−√\delta _{i}=\sqrt {\lambda ...
  • T_27080901
  • T_27080901
  • 2016年03月24日 20:44
  • 639

矩阵操作比较:Armadillo,Eigen,OpenCV

有的时候写出来的程序慢也许并不是算法有问题,而是用的库比较慢;也可能并不是库本身慢,而只是你的写法不够高效。在经历了无数次令人蛋疼的等待后,我决定比较一下这几个所谓的高效的线性代数库(OpenCV虽然...
  • zyh821351004
  • zyh821351004
  • 2015年05月30日 11:07
  • 836

Armadillo C++ linear algebra library 学习笔记(3)——矩阵的基本操作(2)

1、重置矩阵为空矩阵A、通过调用”.reset()”来重置原始矩阵为空矩阵。 B、示例#include #include using namespace arma; int main() { ...
  • jnulzl
  • jnulzl
  • 2015年07月09日 18:33
  • 1879

Armadillo C++ linear algebra library 学习笔记(2)——矩阵的基本操作(1)

1、获取矩阵的某个元素A、通过nameMat(i,j)获取矩阵nameMat的第(i,j)个元素。 B、示例#include #include using namespace arma; int...
  • jnulzl
  • jnulzl
  • 2015年07月08日 21:20
  • 3544

Armadillo C++ linear algebra library 学习笔记(4)——矩阵的运算

1、矩阵的四则运算A、通过”+”、“-”、”*”、”/”、”%”分别进行矩阵的加、减、乘、点除、点乘运算。 B、示例 //1、构造矩阵A、B mat A(2,2),B(2,2); ...
  • jnulzl
  • jnulzl
  • 2015年07月09日 19:09
  • 2635

Armadillo C++ linear algebra library 学习笔记(6)——生成矩阵

1、主对角线为1的矩阵A、通过调用”eye(rows, cols)”函数生成主对角线为1的矩阵。 B、示例#include #include using namespace arma; int...
  • jnulzl
  • jnulzl
  • 2015年07月17日 16:40
  • 2105

Armadillo C++ linear algebra library 学习笔记(7)——统计相关函数

1、求矩阵行、列的最大值A、通过调用”max(A, dim)”函数求矩阵A的行(dim=1)、列(dim=0)最大值。 B、示例#include #include using namespace...
  • jnulzl
  • jnulzl
  • 2015年07月17日 17:12
  • 1837

Armadillo C++ linear algebra library 学习笔记(1)——简介与安装

1、Armadillo C++ 简介 Armadillo C++ Library是一种C++的线性代数库(矩阵数学)以取得良好的平衡速度与易用性。整数,浮点,而复杂的数字支持,以及一个子集,三角和统...
  • jnulzl
  • jnulzl
  • 2015年07月08日 20:58
  • 3769
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Armadillo C++ linear algebra library 学习笔记(5)——矩阵的分解
举报原因:
原因补充:

(最多只允许输入30个字)