Eigen库矩阵的基本类型

1.Eigen: 是C++开源矩阵计算工具.

2.Eigen用源码的方式提供给用户使用,在使用时只需要包含Eigen的头文件即可进行使用。之所以采用这种方式,是因为Eigen采用模板方式实现,由于模板函数不支持分离编译,所以只能提供源码而不是动态库的方式供用户使用,不过这也更方面用户使用和研究。

3.Eigen中矩阵类的模板函数中,共有6个模板参数,但是目前常用的只有前三个,如下所示:

template<typename _Scalar,

                 int _Rows,

                 int _Cols,

                 int _Options,

                 int _MaxRows,

                 int _MaxCols>  

struct traits<Matrix<_Scalar,  _Rows,  _Cols,  _Options,  _MaxRows,  _MaxCols> >  

其中,前三个参数分别表示矩阵元素的类型,行数和列数。

矩阵定义时可以使用Dynamic来表示矩阵的行列数为未知,例如:

typedef Matrix<double,Dynamic, Dynamic> MatrixXd;

在Eigen中也提供了很多常见的简化定义形式,例如:

typedef Matrix< double , 3 , 1> Vector3d  //3维列向量

4.Eigen中无论是矩阵还是数组、向量,无论是静态矩阵还是动态矩阵都提供默认构造函数,也就是你定义这些数据结构时都可以不用提供任何参数,其大小均由运行时来确定。

5.矩阵的构造函数中只提供行列数、元素类型的构造参数,而不提供元素值的构造,对于比较小的、固定长度的向量提供初始化元素的定义,例如:

Vector2d a(5.0, 6.0);  

Vector3d b(5.0, 6.0, 7.0);  

Vector4d c(5.0, 6.0, 7.0, 8.0);  



6.动态矩阵是指其大小在运行时确定,静态矩阵是指其大小在编译时确定,在Eigen中并未这样称呼矩阵。

//代码1

#include <iostream>  

#include <Eigen/Dense>  

using namespace Eigen;  

using namespace std;  

int main()  

{  

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

m = (m + MatrixXd::Constant(3,3,1.2)) * 50;  

cout <<"m =" << endl << m << endl;  

VectorXd v(3);  

v << 1, 2, 3;  

cout <<"m * v =" << endl << m * v << endl;  


//代码2

#include <iostream>  

#include <Eigen/Dense>  

using namespace Eigen;  

using namespace std;  

int main()  

{  

Matrix3d m = Matrix3d::Random();  

m = (m + Matrix3d::Constant(1.2)) * 50;  

cout <<"m =" << endl << m << endl;  

Vector3d v(1,2,3);  

cout <<"m * v =" << endl << m * v << endl;  

}  

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

2) 代码段2中Matrix3d表示元素类型为double大小为3*3的矩阵变量,其大小在编译时就知道;

3)上例中向量的定义也是类似,不过这里的向量时列优先,在Eigen中行优先的矩阵会在其名字中包含有row,否则就是列优先。

4)向量只是一个特殊的矩阵,其一个维度为1而已,如:typedef Matrix< double , 3 , 1> Vector3d


7.访问矩阵元素: 在矩阵的访问中,行索引总是作为第一个参数,需注意Eigen中遵循大家的习惯让矩阵、数组、向量的下标都是从0开始。矩阵元素的访问可以通过()操作符完成,例如m(2,3)即是获取矩阵m的第2行第3列元素(注意行列数从0开始)。可参看如下代码:

#include <iostream>  

#include <Eigen/Dense>  

using namespace Eigen;  

int main()  

{  

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);  

std::cout <<"Here is the matrix m:\n" << m << std::endl;  

VectorXd v(2);  

v(0) = 4;  

v(1) = v(0) - 1;  

std::cout <<"Here is the vector v:\n" << v << std::endl;  

}  


8.矩阵元素赋值: 在Eigen中重载了"<<"操作符,通过该操作符即可以一个一个元素的进行赋值,也可以一块一块的赋值。另外也可以使用下标进行复制,例如下面代码:

9.设置矩阵的大小: 当前矩阵的行数、列数、大小可以通过rows(),cols()和size()来获取,对于动态矩阵可以通过resize()函数来动态修改矩阵的大小.

10.固定大小的矩阵是不能使用resize()来修改矩阵的大小;

11.resize()函数会析构掉原来的数据,因此调用resize()函数之后将不能保证元素的值不改变。

12.使用“=”操作符操作动态矩阵时,如果左右边的矩阵大小不等,则左边的动态矩阵的大小会被修改为右边的大小。例如下面的代码段:


13.如何选择动态矩阵和静态矩阵?

Eigen对于这问题的答案是:

对于小矩阵(一般大小小于16)的使用固定大小的静态矩阵,它可以带来比较高的效率,

对于大矩阵(一般大小大于32)建议使用动态矩阵。

还需特别注意的是:如果特别大的矩阵使用了固定大小的静态矩阵则可能造成栈溢出的问题.


参考:

1.https://blog.csdn.net/liufengl138/article/details/78405652

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 230,501评论 6 544
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 99,673评论 3 429
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 178,610评论 0 383
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 63,939评论 1 318
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 72,668评论 6 412
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 56,004评论 1 329
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 44,001评论 3 449
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 43,173评论 0 290
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 49,705评论 1 336
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 41,426评论 3 359
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 43,656评论 1 374
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 39,139评论 5 364
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 44,833评论 3 350
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 35,247评论 0 28
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 36,580评论 1 295
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 52,371评论 3 400
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 48,621评论 2 380

推荐阅读更多精彩内容

  • 容器 在实际的开发过程中, 数据结构本身的重要性不会逊于操作于数据结构的算法的重要性, 当程序中存在着对时间要求很...
    编程小兔崽阅读 1,103评论 0 1
  • Eigen3的安装方式一般用两种,一是简单的命令安装,二是用源码编译安装; 1.简单命令安装 sudo apt-g...
    AI秘籍阅读 39,394评论 0 2
  • C++基础 模板及标准模板库 模板的作用模板使程序员能够快速的建立具有类型安全得库集合和函数集合,它的实现,方便了...
    I踏雪寻梅阅读 521评论 0 4
  • 技术交流QQ群:1027579432,欢迎你的加入! 1.Cpp中的数组 C++支持数组数据结构,它可以存储一个固...
    CurryCoder阅读 1,606评论 0 1
  • 今天用来晓岛看书。就是高晓松开的书店。没有想象中大。需要提前一周预约,限200人,所以不会太拥挤。 到了朝阳大悦城...
    Mia面阅读 499评论 3 1