如何通过索引设置Eigen矩阵/向量的值？

嗨，你问的这个问题其实Eigen本身就提供了很直观的方式来实现，而且和你想要的matrix[i][j]/vector[i]操作非常接近，我给你详细拆解下：

一、向量的元素赋值 #

Eigen的向量（不管是列向量还是行向量）都支持两种索引访问方式来赋值：

• 使用圆括号()：这是Eigen官方推荐的方式，语法为vec(index)
• 使用方括号[]：和C++数组的访问方式一致，语法为vec[index]

注意：Eigen的索引都是从0开始的，和C++原生容器保持一致。

举个简单的例子：

#include <Eigen/Core>
#include <iostream>

int main() {
    // 声明一个含5个元素的动态列向量
    Eigen::VectorXd col_vec(5);
    // 用()赋值
    col_vec(0) = 1.0;
    col_vec(1) = 2.0;
    // 用[]赋值
    col_vec[2] = 3.0;
    col_vec[3] = 4.0;
    col_vec[4] = 5.0;

    // 行向量的操作完全一致
    Eigen::RowVectorXd row_vec(5);
    row_vec(0) = 10.0;
    row_vec[1] = 20.0;

    std::cout << "列向量:\n" << col_vec << "\n\n行向量:\n" << row_vec << std::endl;
    return 0;
}

二、矩阵的元素赋值（含行主序场景） #

你提到了行主序（RowMajor），首先要明确：Eigen中矩阵的逻辑索引（行号、列号）和存储顺序无关，不管是默认的列主序还是你指定的行主序，都可以用以下两种方式赋值：

方式1：直接使用matrix(row, col)（官方推荐） #

这是最直接高效的方式，直接通过行号和列号定位元素，语法为mat(i, j)，代表第i行第j列的元素。

方式2：使用matrix[row][col]（类似你想要的操作） #

这种方式等价于mat.row(row)[col]，先获取第row行的表达式，再访问该行的第col个元素，也可以正常赋值。

下面是行主序矩阵的示例：

#include <Eigen/Core>
#include <iostream>

int main() {
    // 声明一个3x3的行主序动态矩阵
    Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> mat(3, 3);

    // 方式1：用(i,j)赋值
    mat(0, 0) = 1;
    mat(0, 1) = 2;
    mat(0, 2) = 3;

    // 方式2：用[i][j]赋值
    mat[1][0] = 4;
    mat[1][1] = 5;
    mat[1][2] = 6;

    // 混合两种方式
    mat(2, 0) = 7;
    mat[2][1] = 8;
    mat(2, 2) = 9;

    std::cout << "行主序矩阵:\n" << mat << std::endl;
    return 0;
}

几个需要注意的点 #

• 访问元素前，务必确保矩阵/向量已经分配了足够的空间（动态大小的容器可以用resize()或构造时指定大小，固定大小的容器则是编译时确定的）。
• 在debug模式下，如果索引超出容器的大小，Eigen会触发断言错误提示你；但在release模式下，越界访问会导致未定义行为，一定要避免。
• 虽然[]可以正常使用，但官方更推荐()来访问矩阵元素，因为matrix[i][j]会额外生成行表达式的临时对象，在极端性能敏感的场景下可能有微小的开销。

内容的提问来源于stack exchange，提问作者tangy