本文共 1123 字,大约阅读时间需要 3 分钟。
是一个方便使用的c++矩阵运算库,只要包含Eigen的源码头文件,就能使用。 下面的例子是计算矩阵的特征值特征向量,并把特征向量矩阵实部和虚部分开。
#include#include using namespace std;using namespace Eigen;int main(){ Matrix3d A; //定义3*3 double 型矩阵 A << 1, 2, 3, 5,-5,3,4,6,9; cout << "Here is a matrix, A:" << endl << A << endl << endl; EigenSolver es(A); cout << "The eigenvalues of A are:" << endl << es.eigenvalues() << endl; MatrixXcd vect = es.eigenvectors();//特征向量矩阵 cout << "复数形式输出:" << endl << vect << endl; MatrixXd realvect(vect.rows(), vect.cols()); MatrixXd imagvect(vect.rows(), vect.cols()); cout << "虚实拆分:" << endl; /*两种拆分方式等价 for (size_t i = 0; i < vect.rows(); i++) { for (size_t j = 0; j < vect.cols(); j++) { realvect(i, j) = real(vect(i, j)); imagvect(i, j) = imag(vect(i, j)); } }*/ for (size_t i = 0; i < vect.size(); i++) { *(realvect.data() + i) = real(*(vect.data() + i)); *(imagvect.data() + i) = imag(*(vect.data() + i)); } cout << "realvect:" << endl << realvect << endl; cout << "imagvect:" << endl << imagvect << endl;}
转载地址:http://yrqii.baihongyu.com/