在C++中实现矩阵运算（Eigen基础使用）

在C++中实现矩阵运算（Eigen基础使用）

很多使用过matlab的人肯定对于其中方便的矩阵运算有很深的印象，那么如何在C++上同样实现像matlab一样的便捷计算呢？最近新学到了一个名为Eigen的库，这里记录一下希望加深一下印象也希望可以帮助到他人。

什么是Eigen

其实说白了Eigen就是一个库，在里面为矩阵以及向量进行了定义，并且定义了很多矩阵运算，让你像使用matlab一样在C++中实现矩阵运算（或者说就是很直白的运算），我在这里介绍的也只是其中矩阵的部分，其他部分功能如图所示：

Eigen在windows下的安装大家可以参考：http://blog.csdn.net/abcjennifer/article/details/7781936；

Linux下的安装可以参考：https://zhuanlan.zhihu.com/p/36706885

矩阵

矩阵模板

矩阵模板函数中一共包含六个模板参数，前三个是比较常用的，分别表示矩阵元素的类型、行数、列数。在矩阵定义的时候可以使用 Dynamic 来表示行或者列数未知。

template<typename _Scalar, int _Rows, int _Cols, int _Options, int _MaxRows, int _MaxCols>
class Eigen::Matrix< _Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols >

Eigen::Matrix<int, 3, 4> mat1;              //  3x4 的 int 类型的矩阵 mat1
Eigen::Matrix<double, 3, Dynamic> mat2;     //  3x? 的 double 类型的矩阵 mat2
Eigen::Matrix<float, Dynamic, 4> mat3;      //  ?x4 的 float 类型的矩阵 mat3
Eigen::Matrix<long, Dynamic, Dynamic> mat4; //  ?x? 的 long 类型的矩阵 mat4

在 Eigen 中 typedef 了很多矩阵的类型，通常命名为 Matrix 前缀加一个长度为 1∼4的字符串 S 的命名——MatrixS。其中 S 可以用来判断该矩阵类型，数字 n 表示 n ∗ n，n 的范围是2∼4，字母 d、f、i、c 表示 double、float、int、complex,另外 X 表示行或者列数未知的矩阵。

typedef Matrix<std::complex<double>, 2, 2> Eigen::Matrix2cd;            //  2x2 的 cd 类型的矩阵
typedef Matrix<double, 2, 2> Eigen::Matrix2d;                           //  2x2 的 d 类型的矩阵
typedef Matrix<std::complex<double>, 2, Dynamic> Eigen::Matrix2Xcd;     //  2x? 的 cd 类型的矩阵
typedef Matrix<std::complex<float>, Dynamic, 2> Eigen::MatrixX2cf;      //  ?x2 的 cf 类型的矩阵
typedef Matrix<std::complex<double>, Dynamic, Dynamic> Eigen::MatrixXcd;//  ?x? 的 cd 类型的矩阵
typedef Matrix<int, Dynamic, Dynamic> Eigen::MatrixXi;                  //  ?x? 的 i 类型的矩阵

行优先/列优先

在Eigen库中的矩阵默认存储是列优先的，这一点和matlab很相似，而我们知道在C++中的存储方式是行有限，所以如果我们为了自己不搞混，我们可以使用矩阵模板中的第四个参数（默认是 ColMajor）

Matrix<int, 5, 6, RowMajor> matRow; //  行优先的 5x6 的 int 类型矩阵

静态矩阵和动态矩阵

在 Eigen 中可以用 Dynamic 表示行或者列数未知，所以在定义一个矩阵时并不能确定矩阵的大小，只有在运行时才可以确定大小，然后进行动态分配，而静态矩阵则是在定义时便明确给定了行数以及列数，在编译时就可以分配好内存。

MatrixXd m = MatrixXd::Random(3,3);
Matrix3d m = Matrix3d::Random();

MatrixXd表示任意大小的元素类型为double的矩阵变量，其大小只有在运行时被赋值之后才能知道；MatrixXd::Random(3,3)表示产生一个元素类型为double的3*3的临时矩阵对象。

Matrix