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库。