C++使用eigen库做本征分解

要在C++中使用Eigen库进行本征分解，首先需要包含Eigen头文件，然后可以使用Eigen::SelfAdjointEigenSolver类来进行本征分解操作。下面是一个简单的示例代码：

#include <iostream>#include <Eigen/Dense>int main() {    Eigen::Matrix3f matrix;    matrix << 1, 2, 3,              2, 5, 6,              3, 6, 9;    Eigen::SelfAdjointEigenSolver<Eigen::Matrix3f> eigensolver(matrix);    if (eigensolver.info() != Eigen::Success) {        std::cout << "Eigen decomposition failed!" << std::endl;        return 1;    }    Eigen::Vector3f eigenvalues = eigensolver.eigenvalues();    Eigen::Matrix3f eigenvectors = eigensolver.eigenvectors();    std::cout << "Eigenvalues:" << std::endl;    std::cout << eigenvalues << std::endl;    std::cout << "Eigenvectors:" << std::endl;    std::cout << eigenvectors << std::endl;    return 0;}

此代码首先创建一个3x3的矩阵(matrix)，然后使用Eigen::SelfAdjointEigenSolver对其进行本征分解。如果分解成功，可以通过eigensolver.eigenvalues()和eigensolver.eigenvectors()方法获得本征值和本征向量。接下来，该代码打印出本征值和本征向量。

注意，要编译此代码，需要将Eigen库的头文件路径添加到编译器的包含目录中，并链接Eigen库。