Eigen快速入门

最新推荐文章于 2025-04-26 20:52:37 发布

超速前行

最新推荐文章于 2025-04-26 20:52:37 发布

阅读量1.3k

点赞数

分类专栏： C++ 文章标签： Eigen

本文链接：https://blog.csdn.net/sn_gis/article/details/79009006

版权

C++ 专栏收录该内容

13 篇文章

订阅专栏

Eigen快速入门

一个简单的例子

#include <iostream>
#include "Eigen/Dense"

using namespace std;
using Eigen::MatrixXd;


void fun0()
{
      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;
}

Eigen头文件定义了许多类型，对于简单的应用来说使用MatrixXd就足够
它代表任意大小的矩阵，X代表此意，其中的每一个元素都是double
参照快速参考指引查看其它类型
Eigen/Dense头文件中包含了MatrixXd中所有成员函数的定义和相关类型
定义在该头文件中的所有类和函数，还有其它头文件中的，都位于Eigen命名空间
第一行声明MatrixXd变量，并指明行数和列数，元素没有被初始化
设置左上角元素为3，你需要小括号来引用该元素
索引都是从0开始
紧接着三行设置其它三个值，最后输出该矩阵

Dense头文件中包含的内容：

#include "Core"
#include "LU"
#include "Cholesky"
#include "QR"
#include "SVD"
#include "Geometry"
#include "Eigenvalues"

矩阵个向量

运行时设置大小

void fun1()
{
    MatrixXd m=MatrixXd::Random(3,3);
    m=(m+MatrixXd::Constant(3,3,1.2))*50;
    cout<<"m="<<endl;
    cout<<m<<endl;
    Eigen::VectorXd v(3);
    v<<1,2,3;
    cout<<"v:"<<endl;
    cout<<v<<endl;
    cout<<"m*v="<<endl<<m*v<<endl;

}

这个例子声明了一个3*3的矩阵，并且使用Random()进行初始化
第二行使用现行变换将值限制在10到110之间
MatrixXd::Constant(3,3,1.2)生成一个3*3的矩阵，所有元素都是1.2
引入了一个新的类型，VectorXd，代表任意大小的向量，它声明包含三个元素，没有初始化
使用都好初始化器进行初始化
矩阵m和向量v相乘并输出结果

编译时设置大小


void fun2()
{
    Eigen::Matrix3d m=Eigen::Matrix3d::Random();
    m=(m+Eigen::Matrix3d::Constant(1.2))*50;
    cout<<"m="<<endl;
    cout<<m<<endl;
    Eigen::Vector3d v(1,2,3);
    cout<<"v:"<<endl;
    cout<<v<<endl;
    cout<<"m*v="<<endl<<m*v<<endl;
}