C++矩阵处理库--Eigen初步使用 (转载自 CSDN cyxcw1博客).docxVIP

C++矩阵处理库--Eigen初步使用 (转载自: CSDN cyxcw1的博客)??2013-07-16 21:20:43|??分类：?计算机视觉|举报|字号?订阅项目要进行比较多的矩阵操作，特别是二维矩阵。刚开始做实验时，使用了动态二维数组，于是写了一堆Matrix函数，作矩阵的乘除加减求逆求行列式。实验做完了，开始做代码优化，发现Matrix.h文件里适用性太低，而且动态二维数组的空间分配与释放也影响效率，于是寻找其他解决方案。首先考虑的是与Matlab混合编程，折腾了半天把Matlab环境与VS2010环境之后，发现Matlab编译出来的函数使用起来也比较麻烦，要把数组转化成该函数适用的类型后才能使用这些函数。我的二维数组也不是上千万维的，估计这个转化的功夫就牺牲了一部分效率了。（如果谁有混合编程的心得，求帮忙，囧。。。）接着想到使用一维数组的方法，或者把一维数组封装在一个类里边。想着又要写一堆矩阵操作函数头就大，索性谷歌了一下矩阵处理库，除了自己之前知道的OpenCV库（之前由于转化cvarr麻烦，于是放弃），还有Eigen, Armadillo。/houston11235/article/details/8501135该博客对这三个库的效率做了一个简单的评测，OpenCV库的矩阵操作效率是最低的，还好我没使用。Eigen速度最快，与自己定义数组的操作效率相当（- -，才相当吗？我本来还想找个更快的呢）。于是选择使用Eigen。进入正题。安装：/index.php?title=Main_Page这里是官网，直接把包下载下来，不大，也就几M，我是直接放在自己项目文件夹（考虑项目封装时，这样比较方便），放在VS2010 INCLUDE文件夹。简单使用：看了一下官方文档，Eigen库除了能实现各种矩阵操作外，貌似还提供《数学分析》中的各种矩阵操作（包括L矩阵U矩阵）。目前我使用到的还是简单的矩阵操作，如加减乘除，求行列式，转置，逆，这些基本操作只要:[cpp]?view plaincopyprint?#include?Eigen/Eigen??using?namespace?Eigen;??就能实现，别忘了名空间Eigen。包含的类型：MatricesArraysMatrixfloat,Dynamic,Dynamic =?MatrixXfMatrixdouble,Dynamic,1 =?VectorXdMatrixint,1,Dynamic =?RowVectorXiMatrixfloat,3,3 =?Matrix3fMatrixfloat,4,1 =?Vector4fArrayfloat,Dynamic,Dynamic = ArrayXXfArraydouble,Dynamic,1 = ArrayXdArrayint,1,Dynamic = RowArrayXiArrayfloat,3,3 = Array33fArrayfloat,4,1 = Array4f如上表，主要包括两种类型，Matrices与Arryays，接着是这两种类型的派生类型。现在我用到的是Matrices（我不明白这两种类型在效率间有什么差距，囧。。。），其中Matrix代表二维矩阵，Vector代表列向量RowVector代表行向量。如果后面跟着X，则代表是动态的数组，运行时可以根据需求改变，如果是数字，则代表是静态的（根据实验，最多能建立4维的静态矩阵或者数组，- -，为嘛不是6维，实验正好需要）。i代表int类型，f代表float类型，d代表double。对应关系：Matrix二维矩阵Vector列向量RowVector行向量X动态固定数字n静态，4=n=1iintffloatddoubleArrays类型的话也跟Matrices差不多。基本操作，定义，初始化，矩阵操作：[cpp]?view plaincopyprint?#include?iostream??#include?Eigen/Eigen??using?namespace?std;??using?namespace?Eigen;????void?foo(MatrixXf?m)??{??????Matrix3f?m2=Matrix3f::Zero(3,3);??????m2(0,0)=1;??????m=m2;??}??int?main()??{??????/*?定义，定义时默认没有初始化，必须自己初始化?*/??????MatrixXf?m1(3,4);???//动态矩阵，建立3行4列。??????MatrixXf?m2(4,3);???//4行3列，依此类推。??????MatrixXf?m3(3,3);??????Vector3f?v1;????????//若是静态数组，则不用指定