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

原文链接：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);//得到对应特征向量