Eigen是C++中可以用来调用并进行矩阵计算的一个库,里面封装了一些类,需要的头文件和功能如下:
1.Eigen下载
Eigen的官网下载地址:http://eigen.tuxfamily.org/index.php?title=Main_Page#Download
另外,CSDN资源里也可以下载,http://download.csdn.net/detail/piaoxuezhong/9767251
2.Eigen的简介
矩阵的定义:Eigen中关于矩阵类的模板函数中,共有六个模板参数,常用的只有前三个。其前三个参数分别表示矩阵元素的类型、行数和列数。
矩阵定义时可以使用Dynamic来表示矩阵的行列数为未知。
Eigen中无论是矩阵还是数组、向量,无论是静态矩阵还是动态矩阵都提供默认构造函数,也就是定义这些数据结构时都可以不用提供任何参数,其大小均由运行时来确定。矩阵的构造函数中只提供行列数、元素类型的构造参数,而不提供元素值的构造,对于比较小的、固定长度的向量提供初始化元素的定义。
矩阵类型:Eigen中的矩阵类型一般都是用类似MatrixXXX来表示,可以根据该名字来判断其数据类型,比如”d”表示double类型,”f”表示float类型,”i”表示整数,”c”表示复数;Matrix2f,表示的是一个2*2维的,其每个元素都是float类型。
数据存储:Matrix创建的矩阵默认是按列存储,Eigen在处理按列存储的矩阵时会更加高效。如果想修改可以在创建矩阵的时候加入参数,如:
Matrix<int,3, 4, ColMajor> Acolmajor;
Matrix<int,3, 4, RowMajor> Arowmajor;
动态矩阵和静态矩阵:动态矩阵是指其大小在运行时确定,静态矩阵是指其大小在编译时确定。
MatrixXd:表示任意大小的元素类型为double的矩阵变量,其大小只有在运行时被赋值之后才能知道。
Matrix3d:表示元素类型为double大小为3*3的矩阵变量,其大小在编译时就知道。
在Eigen中行优先的矩阵会在其名字中包含有row,否则就是列优先。Eigen中的向量是一个特殊的矩阵,其维度为1。
矩阵元素的访问:在矩阵的访问中,行索引总是作为第一个参数,Eigen中矩阵、数组、向量的下标都是从0开始。矩阵元素的访问可以通过”()”操作符完成。例如m(2, 3)既是获取矩阵m的第2行第3列元素。
针对向量还提供”[]”操作符,注意矩阵则不可如此使用。
设置矩阵的元素:在Eigen中重载了”<<”操作符,通过该操作符即可以一个一个元素的进行赋值,也可以一块一块的赋值。另外也可以使用下标进行赋值。
重置矩阵大小:当前矩阵的行数、列数、大小可以通过rows()、cols()和size()来获取,对于动态矩阵可以通过resize()函数来动态修改矩阵的大小。注意:
(1)、固定大小的矩阵是不能使用resize()来修改矩阵的大小;
(2)、resize()函数会析构掉原来的数据,因此调用resize()函数之后将不能保证元素的值不改变;
(3)、使用”=”操作符操作动态矩阵时,如果左右两边的矩阵大小不等,则左边的动态矩阵的大小会被修改为右边的大小。
如何选择动态矩阵和静态矩阵:
对于小矩阵(一般大小小于16)使用固定大小的静态矩阵,它可以带来比较高的效率;对于大矩阵(一般大小大于32)建议使用动态矩阵。注意:如果特别大的矩阵使用了固定大小的静态矩阵则可能会造成栈溢出的问题。
矩阵和向量的算术运算:在Eigen中算术运算重载了C++的+、-、*
(1)、矩阵的运算:提供+、-、一元操作符”-”、+=、-=;二元操作符+/-,表示两矩阵相加(矩阵中对应元素相加/减,返回一个临时矩阵);一元操作符-表示对矩阵取负(矩阵中对应元素取负,返回一个临时矩阵);组合操作法+=或者-=表示(对应每个元素都做相应操作);矩阵还提供与标量(单一数字)的乘除操作,表示每个元素都与该标量进行乘除操作;
(2)、求矩阵的转置、共轭矩阵、伴随矩阵:可以通过成员函数transpose()、conjugate()、adjoint()来完成。注意:这些函数返回操作后的结果,而不会对原矩阵的元素进行直接操作,如果要让原矩阵进行转换,则需要使用响应的InPlace函数,如transpoceInPlace()等;
(3)、矩阵相乘、矩阵向量相乘:使用操作符*,共有*和*=两种操作符;
(4)、矩阵的块操作:有两种使用方法:
matrix.block(i,j, p, q) : 表示返回从矩阵(i, j)开始,每行取p个元素,每列取q个元素所组成的临时新矩阵对象,原矩阵的元素不变;
matrix.block<p,q>(i, j) :<p, q>可理解为一个p行q列的子矩阵,该定义表示从原矩阵中第(i, j)开始,获取一个p行q列的子矩阵,返回该子矩阵组成的临时矩阵对象,原矩阵的元素不变;
(5)、向量的块操作:
获取向量的前n个元素:vector.head(n);
获取向量尾部的n个元素:vector.tail(n);
获取从向量的第i个元素开始的n个元素:vector.segment(i,n);
3.配置
直接在属性页的包含目录中加入eigen的目录即可
4.验证
实例1:(动态矩阵设置)
#include <iostream>
#include <Eigen/Dense>
using Eigen::MatrixXd;
using namespace std;
int main(int argc, char* argv[])
{
MatrixXd m(2, 2);
m(0, 0) = 3;
m(1, 0) = 2.5;
m(0, 1) = -1;
m(1, 1) = m(1, 0) + m(0, 1);
cout << m << endl;
return 0;
}
运行结果:
实例2:(静态矩阵设置,及简单矩阵计算)
#include <iostream>
#include <Eigen/Dense>
using namespace std;
using namespace Eigen;
int main()
{
Matrix2d mat;
mat << 1, 2,
3, 4;
Vector2d u(-1, 1), v(2, 0);
cout << "the mat is: \n" << mat << endl;
cout << "the u is:\n" << u << endl;
cout << "the v is:\n" << v << endl;
cout << "Here is mat*mat:\n" << mat*mat << endl;
cout << "Here is mat*u:\n" << mat*u << endl;
cout << "here is u^T\n" << u.transpose() << endl;
cout << "Here is u^T*v:\n" << u.transpose()*v << endl;
cout << "Here is u*v^T:\n" << u*v.transpose() << endl;
}
运行结果:
实例3:(Map矩阵与数组转换)
#include <iostream>
#include "Eigen\Dense"
using namespace Eigen;
using namespace std;
int main()
{
double *Mat1 = new double[6];
for (int i = 0; i < 6; i++)
{
Mat1[i] = i;
}
//静态矩阵,map使用
Eigen:Map<Matrix<double, 2, 3>>staMat(Mat1);
cout << "显示静态矩阵:" << endl;
for (int i = 0; i < staMat.size(); i++)
cout << *(staMat.data() + i) << " ";
cout << endl;
//动态矩阵,map使用
Map <MatrixXd> dyMat(Mat1, 2, 3);
cout << "显示动态矩阵:" << endl;
for (int i = 0; i < dyMat.size(); i++)
cout << *(dyMat.data() + i) << " ";
cout << endl;
}
实例4:(转置矩阵,逆矩阵,伴随矩阵,特征值等求解)
#include <iostream>
#include "Eigen\Dense"
using namespace Eigen;
using namespace std;
int main()
{
Matrix3d Mat1;
Mat1 << 1, 2, 3,
4, 6, 8,
7, 9, 9;
cout << "Mat1=\n" << Mat1 << endl;
cout << "Mat1转置矩阵:\n" << Mat1.transpose() << endl;
cout << "Mat1伴随矩阵:\n" << Mat1.adjoint() << endl;
cout << "Mat1逆矩阵:\n" << Mat1.inverse() << endl;
cout << "Mat1行列式:\n" << Mat1.determinant() << endl;
SelfAdjointEigenSolver<Matrix3d>eigensover(Mat1);
if (eigensover.info() != Success) abort();
cout << "特征值:\n" << eigensover.eigenvalues() << endl;
cout << "特征向量:\n" << eigensover.eigenvectors() << endl;
}
需要注意,在运行时需要设置:
在工程上右键-》属性-》c/c++-》代码生成-》运行库
改成(release为MT,debug为MTD)
否则,可能会出现如下报错:
error LNK2038: 检测到“RuntimeLibrary”的不匹配项: 值“MT_StaticRelease”不匹配值“MD_DynamicRelease”
参考:
http://blog.csdn.net/chenbang110/article/details/12304123
http://blog.csdn.net/abcjennifer/article/details/7781936
http://blog.csdn.net/houjixin/article/details/8477522
http://blog.csdn.net/houjixin/article/details/8490941
http://blog.csdn.net/fengbingchun/article/details/47378515
http://blog.csdn.net/augusdi/article/details/12907341 (eigen较为详细的介绍)
http://blog.csdn.net/wpc320/article/details/8496957