Vector4d B(2,4,4,9); Vector4d x = A.colPivHouseholderQr().solve(B); Vector4d x2 = A.llt().solve(B); Vector4d x3 = A.ldlt().solve(B); std::cout << "The solution is:\n" << x <<"\n\n"<
}
除了colPivHouseholderQr、LLT、LDLT，还有以下的函数可以求解线性方程组，请注意精度和速度： 解小矩阵（4*4）基本没有速度差别
最小二乘求解
最小二乘求解有两种方式，jacobiSvd或者colPivHouseholderQr，4*4以下的小矩阵速度没有区别，jacobiSvd可能更快，大矩阵最好用colPivHouseholderQr
#include
#include
using namespace std;
using namespace Eigen;
int main()
{
MatrixXf A1 = MatrixXf::Random(3, 2);
std::cout << “Here is the matrix A:\n” << A1 << std::endl;
VectorXf b1 = VectorXf::Random(3);
std::cout << “Here is the right hand side b:\n” << b1 << std::endl;
//jacobiSvd 方式:Slow (but fast for small matrices)
std::cout << “The least-squares solution is:\n”
<< A1.jacobiSvd(ComputeThinU | ComputeThinV).solve(b1) << std::endl;
//colPivHouseholderQr方法:fast
std::cout << “The least-squares solution is:\n”
<< A1.colPivHouseholderQr().solve(b1) << std::endl;
}
稀疏矩阵
稀疏矩阵的头文件包括：
#include
typedef Eigen::Triplet T;
std::vector tripletList;
triplets.reserve(estimation_of_entries); //estimation_of_entries是预估的条目
for(…)
{
tripletList.push_back(T(i,j,v_ij)); //第 i,j个有值的位置的值
}
SparseMatrixType mat(rows,cols);
mat.setFromTriplets(tripletList.begin(), tripletList.end());
// mat is ready to go!
2.直接将已知的非0值插入
SparseMatrix mat(rows,cols);
mat.reserve(VectorXi::Constant(cols,6));
for(…)
{
// i,j 个非零值 v_ij != 0
mat.insert(i,j) = v_ij;
}
mat.makeCompressed(); // optional
稀疏矩阵支持大部分一元和二元运算:
sm1.real() sm1.imag() -sm1 0.5*sm1
sm1+sm2 sm1-sm2 sm1.cwiseProduct(sm2)
二元运算中，稀疏矩阵和普通矩阵可以混合使用
//dm表示普通矩阵
dm2 = sm1 + dm1;
也支持计算转置矩阵和伴随矩阵
参考以下链接
点击这里跳转查看更多
第三部分：
其他相关博客：
1、单独下载与安装：https://blog.csdn.net/augusdi/article/details/12907341
2、一篇较详细的教程：https://blog.csdn.net/wzaltzap/article/details/79501856
3、计算特征值特征向量：https://blog.csdn.net/wokaowokaowokao12345/article/details/47375427
Eigen介绍及简单使用（PCL库实现） 转载自：https://blog.csdn.net/fengbingchun/article/details/47378515，

Eigen是可以用来进行线性代数、矩阵、向量操作等运算的C++库，它里面包含了很多算法。它的License是MPL2。它支持多平台。 Eigen采用源码的方式提供给用户使用，在使用时只需要包含Eigen的头文件即可进行使用。之所以采用这种方式，是因为Eigen采用模板方式实现，由于模板函数不支持分离编译，所以只能提供源码而不是动态库的方式供用户使用。 我是使用PCL做的，对头文件 略有改动。

一、 矩阵的定义
Eigen中关于矩阵类的模板函数中，共有六个模板参数，常用的只有前三个。其前三个参数分别表示矩阵元素的类型、行数和列数。

矩阵定义时可以使用Dynamic来表示矩阵的行列数为未知。Eigen中无论是矩阵还是数组、向量，无论是静态矩阵还是动态矩阵都提供默认构造函数，也就是定义这些数据结构时都可以不用提供任何参数，其大小均由运行时来确定。矩阵的构造函数中只提供行列数、元素类型的构造参数，而不提供元素值的构造，对于比较小的、固定长度的向量提供初始化元素的定义。

二、 矩阵类型
Eigen中的矩阵类型一般都是用类似MatrixXXX来表示，可以根据该名字来判断其数据类型，比如”d”表示double类型，”f”表示float类型，”i”表示整数，”c”表示复数；Matrix2f，表示的是一个2*2维的，其每个元素都是float类型。
三、 数据存储
Matrix创建的矩阵默认是按列存储，Eigen在处理按列存储的矩阵时会更加高效。如果想修改可以在创建矩阵的时候加入参数，如:
Matrix Acolmajor;
Matrix 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)、矩阵相乘、矩阵向量相乘
使用操作符*，共有和=两种操作符；
下面这一段参考博客：https://blog.csdn.net/yang_q_x/article/details/52383289
需要注意的是，在做矩阵计算是，用的是循环语句，所以是串行计算的方式。

Eigen还提供了一些常用的变换函数：transpose()，转置操作：conjugata()，共轭操作；adjoint()，共轭转置操作；需要注意的是这些是非原位操作，例如a = [1, 2; , 3, 4] 那么a = a.transpose() 得到的a为[1, 2, 2, 4]，这样就产生了所粗。原位转置操作可以使用transposeInPlace() 函数。

矩阵乘法是一个特殊的运算，为了防止上述的错误，矩阵与矩阵相乘都会使用临时变量，例如m = m* m; 在实际编译的时候是tem = m* m; m = tem; 如果在乘法中你可以确定不会有错误的产生，那么可以使用noaliasd()函数来避免产生临时变量，例如：c.noalias() += a * b;

Eigen还提供了dot()点乘函数和cross()叉乘函数。其中dot可以针对任意两个长度相等的向量，而叉乘只能是两个三维向量，例如Vector3d v(1, 2, 3); Vector3d w(0, 1, 2); 那么v.dot(w) 得到的结果是8（等价于v.adjoint() * w），v.corss(w)得到的结果是(1; -2; 1)。

(4)、矩阵的块操作
有两种使用方法：

matrix.block(i,j, p, q) : 表示返回从矩阵(i, j)开始，每行取p个元素，每列取q个元素所组成的临时新矩阵对象，原矩阵的元素不变；matrix.block(i, j) :可理解为一个p行q列的子矩阵，该定义表示从原矩阵中第(i, j)开始，获取一个p行q列的子矩阵，返回该子矩阵组成的临时矩阵对象，原矩阵的元素不变；

(5)、向量的块操作
获取向量的前n个元素：vector.head(n);

获取向量尾部的n个元素：vector.tail(n); 获取从向量的第i个元素开始的n个元素：vector.segment(i,n); Map类：在已经存在的矩阵或向量中，不必拷贝对象，而是直接在该对象的内存上进行运算操作。

1. 从http://eigen.tuxfamily.org/index.php?title=Main_Page下载最新稳定版本3.2.5，解压缩；

   
   

2. 新建一个vs2013 TestEigen控制台工程，将Eigen文件所在目录加入到工程属性的C/C++附加包含目录中，这样就可以使用Eigen中的函数了；

   
   

3. TestEigen.cpp文件中的内容为：

   
   

//#include “stdafx.h”//原博客中用的是这个，但是我的电脑无法找到，于是用下面这三行替换。
#include
#include “stdlib.h”
#include
#include
#include
//十、eigen的使用
using namespace std;
using namespace Eigen;
template
static void matrix_mul_matrix(T* p1, int iRow1, int iCol1, T* p2, int iRow2, int iCol2, T* p3)
{
if (iRow1 != iRow2) return;

//列优先 //Eigen::Map< Eigen::Matrix > map1(p1, iRow1, iCol1); //Eigen::Map< Eigen::Matrix > map2(p2, iRow2, iCol2); //Eigen::Map< Eigen::Matrix > map3(p3, iCol1, iCol2); //行优先 Eigen::Map< Eigen::Matrix > map1(p1, iRow1, iCol1); Eigen::Map< Eigen::Matrix > map2(p2, iRow2, iCol2); Eigen::Map< Eigen::Matrix > map3(p3, iCol1, iCol2); map3 = map1 * map2;

}
int main(int argc, char* argv[])
{
//1. 矩阵的定义
cout << “1. 矩阵的定义” << endl;
Eigen::MatrixXd m(2, 2);
Eigen::Vector3d vec3d;
Eigen::Vector4d vec4d(1.0, 2.0, 3.0, 4.0);

//2. 动态矩阵、静态矩阵 cout << "2. 动态矩阵、静态矩阵" << endl; Eigen::MatrixXd matrixXd; Eigen::Matrix3d matrix3d; //3. 矩阵元素的访问 cout << "3. 矩阵元素的访问" << endl; m(0, 0) = 1; m(0, 1) = 2; m(1, 0) = m(0, 0) + 3; m(1, 1) = m(0, 0) * m(0, 1); std::cout << m << std::endl << std::endl; //4. 设置矩阵的元素 cout << "4. 设置矩阵的元素" << endl; m << -1.5, 2.4, 6.7, 2.0; std::cout << m << std::endl << std::endl; int row = 4; int col = 5; Eigen::MatrixXf matrixXf(row, col); matrixXf << 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20; std::cout << matrixXf << std::endl << std::endl; matrixXf << Eigen::MatrixXf::Identity(row, col); //Identity是单位矩阵的意思：Identity matrix std::cout << matrixXf << std::endl << std::endl; //5. 重置矩阵大小 cout << "5. 重置矩阵大小" << endl; Eigen::MatrixXd matrixXd1(3, 3); m = matrixXd1; std::cout << m.rows() << "" << m.cols() << std::endl << std::endl; //6. 矩阵运算 cout << "6. 矩阵运算" << endl; m << 1, 2, 7, 3, 4, 8, 5, 6, 9; std::cout << m << std::endl; matrixXd1 = Eigen::Matrix3d::Random(); m += matrixXd1; std::cout << m << std::endl << std::endl; m *= 2; std::cout << m << std::endl << std::endl; std::cout << -m << std::endl << std::endl; std::cout << m << std::endl << std::endl; //7. 求矩阵的转置、共轭矩阵、伴随矩阵 cout << "7. 求矩阵的转置、共轭矩阵、伴随矩阵" << endl; std::cout << m.transpose() << std::endl << std::endl; std::cout << m.conjugate() << std::endl << std::endl; std::cout << m.adjoint() << std::endl << std::endl; std::cout << m << std::endl << std::endl; m.transposeInPlace(); std::cout << m << std::endl << std::endl; //8. 矩阵相乘、矩阵向量相乘 cout << "8. 矩阵相乘、矩阵向量相乘" << endl; std::cout << m*m << std::endl << std::endl; vec3d = Eigen::Vector3d(1, 2, 3); std::cout << m * vec3d << std::endl << std::endl; std::cout << vec3d.transpose()*m << std::endl << std::endl; //9. 矩阵的块操作 cout << "9. 矩阵的块操作" << endl; std::cout << m << std::endl << std::endl; std::cout << m.block(1, 1, 2, 2) << std::endl << std::endl; std::cout << m.block<1, 2>(0, 0) << std::endl << std::endl; std::cout << m.col(1) << std::endl << std::endl; std::cout << m.row(0) << std::endl << std::endl; //10. 向量的块操作 cout << "10. 向量的块操作" << endl; Eigen::ArrayXf arrayXf(10); arrayXf << 1, 2, 3, 4, 5, 6, 7, 8, 9, 10; std::cout << vec3d << std::endl << std::endl; std::cout << arrayXf << std::endl << std::endl; std::cout << arrayXf.head(5) << std::endl << std::endl; std::cout << arrayXf.tail(4) * 2 << std::endl << std::endl; //11. 求解矩阵的特征值和特征向量 cout << "11. 求解矩阵的特征值和特征向量" << endl; Eigen::Matrix2f matrix2f; matrix2f << 1, 2, 3, 4; Eigen::SelfAdjointEigenSolver eigenSolver(matrix2f); //形式就是这个形式，遇到高维的就直接用高维的替换掉matrix2f。 if (eigenSolver.info() == Eigen::Success) { std::cout << eigenSolver.eigenvalues() << std::endl << std::endl; std::cout << eigenSolver.eigenvectors() << std::endl << std::endl; }//12. 类Map及动态矩阵的使用 cout << "12. 类Map及动态矩阵的使用" << endl; int array1[4] = { 1, 2, 3, 4 }; int array2[4] = { 5, 6, 7, 8 }; int array3[4] = { 0, 0, 0, 0 }; matrix_mul_matrix(array1, 2, 2, array2, 2, 2, array3); for (int i = 0; i < 4; i++) std::cout << array3[i] << std::endl; return 0;

}
参考文献：

1. http://blog.csdn.net/augusdi/article/details/12907341

   
   

2. http://www.docin.com/p-863098431.html

   
   

3. http://www.360doc.com/content/15/0325/15/21172899_457946100.shtml

   
   

4. http://www.360doc.com/content/15/0413/22/21172899_463003614.shtml

   
   

【C++小知识——VS2013 MFC基于对话框编程】GitHub：https://github.com/fengbingchun/Eigen_Test

推荐阅读

• 

  刑法报复陷害犯罪规定的刑事量刑档次是怎样的
• 

  台铃米柚电池是什么电池 这电池有什么作用
• 

  实锤！2023雷克萨斯RX“即将登场”，增加2.4T引擎，时尚又大气
• 

  佳能正版清零软件怎么开
• 

  2020-11-27晨间日记
• 

  恰遇惊蛰，我是惊醒来的
• 

  原始取得与继受取得有区别吗？
• 

  亲戚的戚怎么写的 亲戚的戚怎么写
• 

  C++|C++ std::string使用总结
• 

  武魂2活动玩法攻略视频 武魂2活动玩法攻略
• 

  关于生命的意义的非全面思考
• 

  手汗症|手汗“多如雨下”，你有这种困扰吗？五种方法送给你
• 

  怎么去鹅蛋的腥味
• 

  尼康d3200入门教程 尼康d3200图解
• 

  幺鸡是什么意思
• 

  煮生蚝能杀死寄生虫吗
• 

  遇到松下空调一开机显示2怎么办？3个处理方法帮助你
• 

  磁铁的形状有几种 常见的磁铁形状有哪三种
• 

  拍身份证照片可以化妆吗
• 

  香溢融通控股集团股份有限公司

• mfc|MFC简介
• 笔记|C++（继承和派生）
• Windows|Windows --- 常用消息以及消息附加参数
• VC++|VC++使用画刷绘图 简单画刷 FillRect函数 位图画刷 透明画刷
• VC++|VC++ 利用MFC的CWindowDC类实现画线功能 在桌面窗口中画线 绘制彩色线条 CPen nPenStyle nWidth crColor
• VC++|VC++简单绘图 MFC消息映射机制 MessageBox函数
• VC++|VC++ 类向导