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