Eigen库求取最大特征值和特征向量

最新推荐文章于 2024-09-26 19:02:06 发布

倔强不倒翁

最新推荐文章于 2024-09-26 19:02:06 发布

阅读量6.7k

点赞数 1

分类专栏： Eigen 文章标签： Eigen

Eigen 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

原文链接：http://blog.csdn.net/wcsgzc/article/details/53946345

Eigen库中有求取矩阵特征值和特征向量的函数EigenSolver，用起来很方便。
但是官网说明文档里，求取特征向量后仅仅是输出来表示，如何使用求得的这些参数，我花了一点时间琢磨出来。
举例如下：

    Matrix<double, 4, 4>K= MatrixXd::Random(4,4);
    EigenSolver<Matrix<double, 4, 4>> es(K);
    MatrixXcd evecs = es.eigenvectors();//获取矩阵特征向量4*4，这里定义的MatrixXcd必须有c，表示获得的是complex复数矩阵
    MatrixXcd evals = es.eigenvalues();//获取矩阵特征值 4*1
    MatrixXd evalsReal;//注意这里定义的MatrixXd里没有c
    evalsReal=evals.real();//获取特征值实数部分
    MatrixXf::Index evalsMax;
    evalsReal.rowwise().sum().maxCoeff(&evalsMax);//得到最大特征值的位置
    Vector4d q;
    q << evecs.real()(0, evalsMax), evecs.real()(1, evalsMax), evecs.real()(2, evalsMax), evecs.real()(3, evalsMax);//得到对应特征向量