eigen 求解ax=xb

Eigen是一个用于数值线性代数计算的C++模板库，其中包含了许多求解线性方程的算法。在Eigen中，可以使用以下代码来求解方程ax = xb：

#include <Eigen/Dense>

int main() {
    // 定义矩阵a和b
    Eigen::MatrixXd a = Eigen::MatrixXd::Random(3, 3);
    Eigen::MatrixXd b = Eigen::MatrixXd::Random(3, 3);

    // 定义向量x
    Eigen::VectorXd x(3);

    // 使用Eigen提供的API求解方程
    Eigen::GeneralizedEigenSolver<Eigen::MatrixXd> solver(a, b);
    x = solver.eigenvalues().real();

    // 打印解向量x
    std::cout << "解向量x为：" << x << std::endl;

    return 0;
}

以上代码首先定义了两个3x3的随机矩阵a和b，然后定义了一个3维向量x。接着，使用Eigen提供的GeneralizedEigenSolver求解器，并将矩阵a和b作为参数传入。最后，将求解结果存储在向量x中，并通过输出语句打印解向量x。

需要注意的是，以上代码只是一个简单的示例，实际应用中可能需要对输入数据进行处理和验证，以确保求解的正确性。