特征矩阵广播操作

时间:2016-04-25 19:19:20

标签: c++ eigen

我有一个MatrixXf变量和一个VectorXf变量。我想使用我的矩阵上的Vector执行行划分。有可能做这样的事吗?

#include <iostream>
#include "Eigen/Dense"
using namespace std;
using namespace Eigen;

int main() {
    MatrixXf mat(3, 2);
    mat << 1, 2,
           3, 4,
           5, 6;
    VectorXf vec(2);
    vec << 2, 3;
    mat = mat.rowwise() / vec;
    cout << mat << endl;
    return 0;
}

我期望得到一个价值矩阵[0.5,0.667; 1.5,1.333; 2.5,2]。

非常感谢!

1 个答案:

答案 0 :(得分:6)

您需要将矩阵和向量用作数组(而不是线性代数对象,请参阅docs)。为此,您需要将相关行重写为:

mat = mat.array().rowwise() / vec.transpose().array();
cout << mat << endl; // Note that in the original this was vec

transpose是必需的,因为根据定义VectorXf是一个列向量,你想要一个行向量。