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

原创 2015年07月09日 19:09:19

1、矩阵的四则运算

A、通过”+”、“-”、”*”、”/”、”%”分别进行矩阵的加、减、乘、点除、点乘运算。
B、示例

#include <iostream>
#include <armadillo>
using namespace arma;
int main()
{
    //1、构造矩阵A、B
    mat A(2,2),B(2,2);
    for (int i = 0; i < 4; i++)
    {
        A(i)=i+1;//i:以列为排列顺序,矩阵的第i个元素
        B(i)=i+5;
    }
    A.print("矩阵A为:");
    B.print("矩阵B为:");
    //2、矩阵A与矩阵B相加
    mat AAddB = A + B;
    AAddB.print("A + B =");
    //3、矩阵A与矩阵B相减
    mat AMinusB = A - B;
    AMinusB.print("A - B =");
    //4、矩阵A与矩阵B相乘
    mat AMulB = A * B;
    AMulB.print("A * B =");
    //5、矩阵A与矩阵B点除
    mat ADiviB = A / B;
    ADiviB.print("A / B =");
    //6、矩阵A与矩阵B点乘
    mat ADotMulB = A % B;
    ADotMulB.print("A % B =");
}

C、结果
这里写图片描述

2、矩阵的逆与伪逆

A、通过”inv(A)”、“pinv(A)”分别求矩阵A的逆与伪逆矩阵。
B、示例

#include <iostream>
#include <armadillo>
using namespace arma;
int main()
{
    //1、产生随机矩阵A,大小为5x5,每个数的范围为:(0,10)
    mat A = randu<mat>(5, 5)*10;
    mat A_inv = inv(A);//当A为方阵时,用inv()
    mat A_pinv = pinv(A);//当A不为方阵时,用pinv()
    mat AMulA_inv = A * A_inv;//验证结果
    A.print("原矩阵A:");
    A_inv.print("矩阵A的逆矩阵为:");
    AMulA_inv.print("A * A^-1为:");
    A_pinv.print("矩阵A的伪逆矩阵为:");
}

C、结果
这里写图片描述

3、解线性方程组

A、通过”solve(A,B)”来解线性方程组:AX=B。
B、示例

#include <iostream>
#include <armadillo>
using namespace arma;
int main()
{
    //1、产生随机矩阵A,大小为5x5,每个数的范围为:(0,10)
    mat A = randu<mat>(5, 5)*10;
    //2、获取矩阵的第2列
    colvec B = A.col(1);
    //3、解线性方程组AX=B
    mat X = solve(A,B);
    X.print("X = ");
}

C、结果
这里写图片描述

Armadillo矩阵库的使用(二)之API接口

一、获取矩阵的某个元素通过nameMat(i,j)获取矩阵nameMat的第(i,j)个元素。 #include #include using namespace arma; int main()...

矩阵伪逆介绍及C++/OpenCV/Eigen的三种实现

矩阵伪逆介绍及C++/OpenCV/Eigen的三种实现

Cpp类似Matlab的矩阵运算库Armadillo

在C++中数值计算需要对矩阵进行操作,今天就介绍一款类似Matlab的矩阵运算库Armadillo。真的是太美了。下面介绍Ubuntu14.04下的安装及其使用。1 安装按照源码中的README.tx...

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
  • 3250

C++中armadillo矩阵库使用说明

在http://blog.csdn.net/piaoxuezhong/article/details/58055709博文中介绍了eigen矩阵库的使用,这里介绍另一种矩阵库:armadillo~ ...

C++线性代数库Armadillo

做机器学习的同学一定对Matlab矩阵运算的简洁高效映像深刻,并且已经习以为常这种便捷。但实际工作中,我们经常需要将matlab代码移植成C/C++。       C++下的线性代数库诸如lapack...

Armadillo (C++ library)

Latest Files armadillo-4.419.alpha1.tar.gz (日期: 2014-09-03, 大小: 5.2 MB)armadillo-4.400.1.tar.gz (...

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
  • 1713

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

1、矩阵的SVD分解A、通过调用”svd(X)”函数进行矩阵X的奇异值分解(svd)。 B、示例 //1、产生随机矩阵A,大小为5x5,每个数的范围为:(0,10) mat A = r...
  • jnulzl
  • jnulzl
  • 2015年07月09日 22:33
  • 1676

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
  • 1982
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Armadillo C++ linear algebra library 学习笔记(4)——矩阵的运算
举报原因:
原因补充:

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