Eigen的下载安装不说了,这个很简单,不懂的查一下就好了。
在Eigen中,所有矩阵和向量都是矩阵模板类的对象。向量只是矩阵的一种特殊情况,有1行或1列。
Matrix类共有六个模板参数,这里只介绍前三个参数,剩下的三个参数有其默认值,我们不需要更改。矩阵的三个强制性模板参数是:
Matrix<typename Scalar, int RowsAtCompileTime, int ColsAtCompileTime>
RowsAtCompileTime
和ColsAtCompleTime
是在编译时已知的矩阵的行数和列数Eigen提供了许多方便的typedef来涵盖常见的情况。例如,Matrix4f是浮点的4x4矩阵。以下是Eigen对它的定义:
typedef Matrix<float, 4, 4> Matrix4f;
在Eigen中,向量只是矩阵的一种特殊情况,具有1行或1列。它们有1列的情况是最常见的;这样的矢量称为列矢量,通常简称为矢量。在它们有1行的另一种情况下,它们被称为行向量。
// 包含3个浮点的(列)向量
typedef Matrix<float, 3, 1> Vector3f;
// 包含2个的整数(行)向量
typedef Matrix<int, 1, 2> RowVector2i;
Eigen并不局限于在编译时维数已知的矩阵。RowsAtCompileTime
和ColsAtCompleTime
模板参数可以采用特殊值Dynamic
,该值表示在编译时大小未知,因此必须作为运行时变量处理。在特征术语中,这种大小被称为动态大小(dynamic size);而在编译时已知的大小被称为固定大小(fixed size)。
// 创建一个双精度的动态矩阵
typedef Matrix<double, Dynamic, Dynamic> MatrixXd;
// 创建一个整型列向量
typedef Matrix<int, Dynamic, 1> VectorXi;
// 创建拥有固定数量的行和动态数量的列的矩阵
Matrix<float, 3, Dynamic>
Eigen提供了默认的构造函数,他不会提供动态内存的分配,也不会初始化任何矩阵的值。Eigen类型可以这样使用:
// 一个3 * 3 的矩阵,矩阵的元素都没有被初始化
Matrix3f a;
// 一个动态矩阵,大小是0 * 0,也就是说还没有为该矩阵分配内存
MatrixXf b;
构造函数提供指定矩阵大小的重载。对于矩阵,总是先传递行数。对于向量,只需传递向量大小。它会分配矩阵或向量所需内存的大小,但是不会初始化他们的值。
// 初始化一个10 * 15的矩阵,内存进行分配,但是没有初始化
MatrixXf a(10,15);
// 初始化一个大小为30的向量,但是没有初始化
VectorXf b(30);
Eigen使用圆括号()
访问矩阵或者向量的元素,序号从0开始。
Eigen不支持使用方括号[]
访问矩阵的元素(向量除外)。
// 矩阵不可以使用方括号进行初始化
// 例如m[0, 0] = 3;
// 这是不合法的
MatrixXd m(2, 2);
m(0, 0) = 3;
// 使用圆括号或方括号进行初始化向量都是合法的
VectorXd b(3);
b(0) = 1;
b[1] = 2;
矩阵和向量也可以从系数列表中初始化。
在 C++11 之前,此功能仅限于小型固定大小的列或大小不超过 4 的向量:
Vector2d a(5.0, 6.0);
Vector3d b(5.0, 6.0, 7.0);
Vector4d c(5.0, 6.0, 7.0, 8.0);
如果启用C++11,则可以通过传递任意数量的系数来初始化任意大小的固定大小的列或行向量:
Vector2i a(1, 2);
Matrix<int, 5, 1> b {1, 2, 3, 4, 5};
Matrix<int, 1, 5> c = {1, 2, 3, 4, 5};
在具有固定或运行时大小的矩阵和向量的一般情况下,系数必须按行分组,并作为初始化器列表的初始化器列表传递
MatrixXi a { // 构造一个 2x2 矩阵
{1, 2}, // 第一行
{3, 4} // 第二行
};
Matrix<double, 2, 3> b {
{2, 3, 4},
{5, 6, 7},
};
对于列或行向量,允许隐式换位。这意味着可以从一行初始化列向量:
VectorXd a {
{1.5, 2.5, 3.5}}; // 具有 3 个系数的列向量
RowVectorXd b {
{1.0, 2.0, 3.0, 4.0}}; // 具有 4 个系数的行向量
逗号表达式初始化:
Matrix3f m;
m << 1, 2, 3,
4, 5, 6,
7, 8, 9;
std::cout << m;
矩阵的当前大小可以通过rows()
、cols()
和size()
来检索。这些方法分别返回行数、列数和矩阵元素数。
通过resize()
方法重置矩阵的大小,如果矩阵大小没有变化,那么resize()
操作没有任何影响。如果矩阵的大小改变了,那么resize()
操作没有任何影响。如果矩阵的大小改变了,那么矩阵的值可能会改变。如果你想在重设大小过程中不改变矩阵的值使用conservativeResize()
使用赋值符号“=”,Eigen自动调整左侧矩阵的大小,使其与右侧矩阵的大小相匹配。
当然,如果左侧是固定大小的,则不允许调整其大小。
**加减法:**重载C++中“+”
、“-”
、“+=”
、“-=”
操作符,要求左右操作数的相同。不允许一个向量加上或者减去一个数。只要维度一样,向量和矩阵是可I加的,这也说明了向量是特殊的矩阵。
如果需要实现矩阵中的每一个元素都加/减同一个元素的操作,在Eigen中是不可以用矩阵直接加/减一个数的。我们需要用array()
将矩阵转换为数组然后在加/减数字。【详见第二部分,数组类章节】
MatrixXd m(3, 5);
m.setZero();
m = m.array() + 5;
**数乘与数除:**重载C++中“\*”
、“/”
、“*=”
、“/=”
操作符,支持矩阵和向量乘以或者除以一个数。
如果我们并不希望两矩阵按照矩阵乘法的形式相乘,而是希望对应位置的元素相乘,则我们需要用array()
将矩阵转换为数组然后再进行乘法操作。【详见第二部分,数组类章节】
Vector3i m;
Vector3i n;
m << 1, 2, 3;
n = m;
m = m.array() * n.array();
cout << m << endl;
**转置、共轭、共轭转置、逆矩阵:**转置、共轭、共轭转置和求逆分别通过transpose()
、conjugate()
、adjoint()
和inverse()
实现。对于实数而言,共轭没有任何影响,共轭转置等价于转置。
MatrixXcf a = MatrixXcf::Random(2, 2);
cout << "矩阵 a\n " << a << endl;
cout << "转置 a^T\n " << a.transpose() << endl;
cout << "共轭 a\n " << a.conjugate() << endl;
cout << "共轭转置 a^*\n" << a.adjoint() << endl;
cout << "逆矩阵 a^-1\n " << a.inverse() << endl;
// 如果矩阵不可逆,会返回无穷值,已经不能再进行计算了
使用a=a.transpose()
可能会出现错误,这是因为Eigen在进行转置或者共轭操作时,会同时写左操作数,从而得到意想不到的结果。要实现这种功能可以使用a.transposeInPlace()
。类似的,也支持adjointInPlace()
。
Matrix2i a; a << 1, 2, 3, 4;
cout << "矩阵 a:\n" << a << endl;
a = a.transpose(); // !!! 不可以这么做 !!!
cout << "混叠效应的结果:\n" << a << endl;
MatrixXf a(2,3); a << 1, 2, 3, 4, 5, 6;
cout << "矩阵 a:\n" << a << endl;
a.transposeInPlace();
cout << "转置后:\n" << a << endl;
**点乘和叉乘:**点乘又可称为内积,Eigen分别使用dot()
和cross()
来实现内积和向量积。叉乘只适用于三维向量。
Vector3d v(1, 2, 3);
Vector3d w(0, 1, 2);
std::cout << "点乘: " << v.dot(w) << std::endl;
std::cout << "叉乘:\n" << v.cross(w) << std::endl;
基础的代数运算:
#include <iostream>
#include <Eigen/Dense>
using namespace std;
int main()
{
Eigen::Matrix2d mat;
mat << 1, 2,
3, 4;
cout << "求和: " << mat.sum() << endl;
cout << "连乘积: " << mat.prod() << endl;
cout << "平均值: " << mat.mean() << endl;
cout << "最小值: " << mat.minCoeff() << endl;
cout << "最大值: " << mat.maxCoeff() << endl;
cout << "迹: " << mat.trace() << endl;
}
计算最大值和最小值的函数支持返回最大值和最小值的位置。
Matrix3f m = Matrix3f::Random();
std::ptrdiff_t i, j;
float minOfM = m.minCoeff(&i, &j);
cout << "矩阵 m:\n" << m << endl;
cout << "最小元素 (" << minOfM
<< ") 在下标处 (" << i << "," << j << ")\n\n";
RowVector4i v = RowVector4i::Random();
int maxOfV = v.maxCoeff(&i);
cout << "向量 v: " << v << endl;
cout << "最大元素 (" << maxOfV
<< ") 在下标处 " << i << endl;
Array类提供通用数组,而Matrix类用于线性代数。此外,Array类提供了一种执行按系数运算的简单方法,这种运算不具有线性代数意义,例如向Array中的每个系数添加一个常数或按系数乘以两个数组。
Array是一个类模板,采用与Matrix相同的模板参数。与Matrix一样,前三个模板参数是强制性的:
Array<typename Scalar, int RowsAtCompileTime, int ColsAtCompileTime>
参数的含义和上面Matrix相同,所以不多解释。
Eigen 还为一些常见情况提供了 typedef,其方式类似于矩阵typedef,但有一些细微差别,因为“数组”一词用于一维和二维数组。 我们采用 ArrayNt 形式的typedef代表一维数组的约定,其中 N 和 t 是大小和标量类型, 对于二维数组,使用 ArrayNNt 形式的类型定义。 下表显示了一些示例:
Type | Typedef |
---|---|
Array<float,Dynamic,1> | ArrayXf |
Array<float,3,1> | Array3f |
Array<double,Dynamic,Dynamic> | ArrayXXd |
Array<double,3,3> | Array33d |
访问 Array 中的值的方式和 Matrix 是一样的。下面介绍和 Matrix 不一样的操作:
**加减乘除:**将对应位置元素相加(减、乘、除)。
**平方根、绝对值:**调用 sqrt()
、abs()
求两个矩阵的最小值:a.min(b)
,要求 a 和 b 的维度相同,求对应位置的两个 数的最小值。
那么,如何选择 Matrix 和 Array 呢? 不能对数组应用矩阵运算,也不能对矩阵应用数组运算。因此,如果你需要做线性代数运算,比如矩阵乘法,那么你应该使用矩阵;如果需要进行系数运算,那么应该使用数组。如果又要支持线性代数操作,又要支持数运算操作,那么就可以使用.array()
和.matrix()
实现类型的转换,这种转换是没有计算代价的。这些操作也不会改变调用该函数的矩阵或者数组本身,而是返回一个副本。
#include <Eigen/Dense>
#include <iostream>
using Eigen::MatrixXf;
int main()
{
MatrixXf m(2,2);
MatrixXf n(2,2);
MatrixXf result(2,2);
m << 1,2,
3,4;
n << 5,6,
7,8;
result = m * n;
std::cout << "-- Matrix m*n: --\n" << result << "\n\n";
result = m.array() * n.array();
std::cout << "-- Array m*n: --\n" << result << "\n\n";
result = m.cwiseProduct(n);
std::cout << "-- With cwiseProduct: --\n" << result << "\n\n";
result = m.array() + 4;
std::cout << "-- Array m + 4: --\n" << result << "\n\n";
}
/*
-- Matrix m*n: --
19 22
43 50
-- Array m*n: --
5 12
21 32
-- With cwiseProduct: --
5 12
21 32
-- Array m + 4: --
5 6
7 8
*/
Eigen 支持的数组元操作见下图所示。
块是矩阵或数组的矩形部分。块表达式既可以用作右值,也可以用作左值。与Eigen表达式一样,只要让编译器进行优化,这种抽象的运行时成本为零。
最常用的是 block()操作,它有两个版本:
// 动态大小的块
matrix.block(i,j,p,q);
// 固定大小的块
matrix.block<p,q>(i,j);
i、j 表示块的起始位置(块的左上角元素在矩阵中的角标),p、q 表示块的大小。 固定大小的块在运行时效率要高。
Eigen::MatrixXf m(4, 4);
m << 1, 2, 3, 4,
5, 6, 7, 8,
9, 10, 11, 12,
13, 14, 15, 16;
cout << "中间的块:" << endl;
cout << m.block<2, 2>(1, 1) << endl << endl;
for (int i = 1; i <= 3; ++i)
{
cout << "块的大小: " << i << "x" << i << endl;
cout << m.block(0, 0, i, i) << endl << endl;
}
/*
中间的块:
6 7
10 11
块的大小: 1x1
1
块的大小: 2x2
1 2
5 6
块的大小: 3x3
1 2 3
5 6 7
9 10 11
*/
block()
操作支持任意的块操作,对于一些特殊的访问 Eigen 也提供了 API。
ith row * matrix.row(i);
jth column * matrix.col(j);
col()
和row()
的参数是要访问的列或行的索引,col()
和row()
操作既可以是左值又可以是右值。
矩阵支持的块操作:
向量支持的块操作:
文章浏览阅读2.1k次。 需求如下:该搜索框是对整个页面的input检索 ,但与弹出层中的input冲突 博主几经辗转 简单处理 解决问题,思路如下:排除掉特定class的input。代码如下:$('input:not(.pop)', this.footer()).on('keyup change', function () { if (that.search() !== th..._input排他选择器
文章浏览阅读5.6k次,点赞6次,收藏20次。看到别人有个1024的勋章,特意留了一篇在今年的10.24日,看看会不会获得。在日常开发中可能涉及接口之间的相互调用,虽然在现在微服务的理念推广下,很多公司都采用轻量级的JSON格式做为序列化的格式,但是不乏有些公司还是有一些XML格式的报文,最近就在对接某个合作方的时候遇到了XML报文。在JSON报文爽快的转换下很难试用一个一个的拿报文参数,还是希望能直接将报文转换成Bean。接下来就了解到..._jaxb 泛型
文章浏览阅读1.2k次。numpy的主要数据对象是多维数组,其中包含相同类型的元素,通常是数字类型,每个元素都有一个索引。使用numpy前通常要导入包。import numpy as np目录类型维度创建运算索引和切片类型numpy的数组被称为ndarray。numpy.array只处理一维数组,而ndarray对象才提供更多功能。a = np.array([[1, 2, 3], [4, 5, 6]])type(a) # <class 'numpy.ndarray'>dtype属性可以获得元素的数_ndarray的位置
文章浏览阅读1.6w次。还在苦于网上找到的一些指令已经不适用了吗?还在苦于有些地方的指令有误吗?还在苦于有些地方整理的指令不够全面吗?那么你来对地方了!小编为大家整理了《我的世界》原版游戏常用的指令,这些基本足以满足各位的基本需求了!大家来一起看看吧!注:表示的是必须输入的部分,[方括号]表示的是可选择性输入的部分基本命令列表命令描述/?/help的替代命令,提供命令使用帮助。/ban + 玩家名字将玩家加入封禁列表。/..._gamemode指令java
文章浏览阅读1.5w次,点赞3次,收藏3次。Spring Boot 结合shiro做第三方登录验证1、首先,说一下我的具体实现思路。在做spring boot拦截器的过程中,开始我准备用spring security来实现,但是研究了一段时间之后发现spring security的集成度太高,需要修改的东西比较多,而且对它本身的使用方法不是很了解,后来转而使用Apache shiro。由于是第三方登录,是不需要我来验证密码的。最开始,我陷入了_shiro 第三方token登录
文章浏览阅读1.9k次,点赞4次,收藏4次。[INFO ] __init__:get_config:71 - Loading config file from:C:\Users\xxx\.labelmercTraceback (most recent call last): File .... line 191, in <module> main() File ...., line 145, in main config = get_config(config_file_or_yaml, config_fro_file "c:\rgzn\labelme-main\setup.py", line 91, in main if sys.argv[1] == "re
文章浏览阅读8.7k次。代码错误的原因和调试方法_代码报错
文章浏览阅读5.2k次,点赞9次,收藏40次。---恢复内容开始---1.认识游戏 1.1什么是游戏 1.1.1游戏的定义 任何人类正常生理需求之外的活动均可称为游戏 1.1.2游戏的分类 RPG角色扮演游戏、ACT动作游戏、AVG冒险游戏、FPS第一人称视角射击游戏、TPS第三人称视角射击游戏、FTG格斗游戏、SPT体育游戏、RAC竞速游戏、RTS即时战略游戏、STG..._深度解析java游戏服务器开发
文章浏览阅读4k次。CSRF是什么我就不解释了,百度一搜全是,比波姐的片源还要多,千篇一律都他么是复制粘贴。那为什么这个令牌(token)操作可以防范CSRF呢?下面我就随便说说说错了大家不要介意。首先我们要知道令牌是存储在session里面的,这个很重要 php代码如下<?php namespace app\index\controller; //我直接允许跨域,因为伪装..._tp5 开启csrf令牌
文章浏览阅读1.7k次,点赞2次,收藏6次。市盈率PE市盈率 = 市值/净利润概念解析:买入一家公司,几年回本,年化收益率:净利润/市值(市盈率的倒数)举例:砖头10万买个砖头,每年拍人带来1万利润,需要10年回本市盈率:10/1 = 10年化收益率:1/10 = 10%市净率PB市净率 = 市值/净资产净资产 = 总资产 - 负债举例:张三便利店,净资产:120万市值:1..._净资产收益率和股息率
文章浏览阅读737次。教育部昨举行「102年国立馆所文创商品设计比赛」颁奖典礼,台北科技大学创新设计研究所硕士生谢镇宇,为TW艺术教育馆设计「墨器」杯垫,取「默契」谐音,用5片压克力板,展现水墨画层层渲染效果,增加立体视觉感受,并在杯架后方加入LED光源,获评审肯定夺特优奖和奖金10万元。台南应用科技大学商品设计系学生高郁翔,为国立自然科学博物馆设计「恐龙化石钉书机」,他认为小朋友把钉书机钉下去的那一刻,会觉得像暴龙準_杯垫文创设计说明
文章浏览阅读404次。XML与对象,集合的相互转化 今天小伙伴在群里问了一下关于XML与对象之间的相互转换,作为菜鸟的我正好趁着闲着的时间学习了一波,直接上代码了,有疑问或者有错误的地方还请大家指正,谢谢。。。。 1 using System; 2 using System.Collections.Generic; 3 using System.IO; 4 using System...._c# xml转集合