C++ Eigen庫計算矩陣特征值及特征向量
更新時間：2019年06月17日 11:33:57 作者：呂小豬不壞
這篇文章主要為大家詳細介紹了C++ Eigen庫計算矩陣特征值及特征向量，具有一定的參考價值，感興趣的小伙伴們可以參考一下
本文主要講解利用Eigen庫計算矩陣的特征值及特征向量并與Matlab計算結(jié)果進行比較。
C++Eigen庫代碼
#include <iostream>
#include <Eigen/Dense>
#include <Eigen/Eigenvalues>
using namespace Eigen;
using namespace std;
void Eig()
{
Matrix3d A;
A << 1, 2, 3, 4, 5, 6, 7, 8, 9;
cout << "Here is a 3x3 matrix, A:" << endl << A << endl << endl;
EigenSolver<Matrix3d> es(A);
Matrix3d D = es.pseudoEigenvalueMatrix();
Matrix3d V = es.pseudoEigenvectors();
cout << "The pseudo-eigenvalue matrix D is:" << endl << D << endl;
cout << "The pseudo-eigenvector matrix V is:" << endl << V << endl;
cout << "Finally, V * D * V^(-1) = " << endl << V * D * V.inverse() << endl;
}
int main()
{
Eig();
}
計算結(jié)果：
最大最小特征值及其索引位置
//maxCoeff //minCoeff int col_index, row_index; cout << D.maxCoeff(&row_index, &col_index) << endl; cout << row_index << " " << col_index << endl;
Matlab 代碼
clear all clc A = [1 2 3;4 5 6;7 8 9] [V,D] = eig(A)
Matlab計算結(jié)果
使用sort()函數(shù)對特征值排序
主成份分析以及許多應(yīng)用時候，需要對特征值大小排列。
A = magic(6); [V,D] = eig(A) [D_S,index] = sort(diag(D),'descend') V_S = V(:,index)
結(jié)果
V = 0.4082 -0.2887 0.4082 0.1507 0.4714 -0.4769 0.4082 0.5774 0.4082 0.4110 0.4714 -0.4937 0.4082 -0.2887 0.4082 -0.2602 -0.2357 0.0864 0.4082 0.2887 -0.4082 0.4279 -0.4714 0.1435 0.4082 -0.5774 -0.4082 -0.7465 -0.4714 0.0338 0.4082 0.2887 -0.4082 0.0171 0.2357 0.7068 D = 111.0000 0 0 0 0 0 0 27.0000 0 0 0 0 0 0 -27.0000 0 0 0 0 0 0 9.7980 0 0 0 0 0 0 -0.0000 0 0 0 0 0 0 -9.7980 D_S = 111.0000 27.0000 9.7980 -0.0000 -9.7980 -27.0000 V_S = 0.4082 -0.2887 0.1507 0.4714 -0.4769 0.4082 0.4082 0.5774 0.4110 0.4714 -0.4937 0.4082 0.4082 -0.2887 -0.2602 -0.2357 0.0864 0.4082 0.4082 0.2887 0.4279 -0.4714 0.1435 -0.4082 0.4082 -0.5774 -0.7465 -0.4714 0.0338 -0.4082 0.4082 0.2887 0.0171 0.2357 0.7068 -0.4082
結(jié)語
本人是在實驗中利用Eigen庫求取最小特征值對應(yīng)特征向量做PCA分析時使用，曾經(jīng)再不知道有Eigen庫的情況下自己寫過矩陣相關(guān)運算的模板類，現(xiàn)在接觸到Eigen庫，就把困擾過自己的問題今天做一個小小總結(jié)。
以上就是本文的全部內(nèi)容，希望對大家的學習有所幫助，也希望大家多多支持腳本之家。
相關(guān)文章
這篇文章主要為大家詳細介紹了C++實現(xiàn)機票預(yù)訂系統(tǒng)，具有一定的參考價值，感興趣的小伙伴們可以參考一下2019-01-01
這篇文章主要介紹了C/C++中多重繼承詳解及其作用介紹，本文給大家介紹的非常詳細，對大家的學習或工作具有一定的參考借鑒價值，需要的朋友可以參考下2021-09-09
C語言驅(qū)動開發(fā)內(nèi)核枚舉IoTimer定時器解析
這篇文章主要為大家介紹了C語言驅(qū)動開發(fā)內(nèi)核枚舉IoTimer定時器解析，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進步，早日升職加薪2022-10-10
大家好，本篇文章主要講的是C語言實現(xiàn)第一次防死版掃雷游戲，感興趣的同學趕快來看一看吧，對你有幫助的話記得收藏一下2022-01-01
C/C++實現(xiàn)動態(tài)數(shù)組的示例詳解
動態(tài)數(shù)組相比于靜態(tài)數(shù)組具有更大的靈活性，因為其大小可以在運行時根據(jù)程序的需要動態(tài)地進行分配和調(diào)整，本文為大家介紹了C++實現(xiàn)動態(tài)數(shù)組的方法，需要的可以參考下2023-08-08
C語言自定義類型詳解（結(jié)構(gòu)體、枚舉、聯(lián)合體和位段）
這篇文章主要給大家介紹了關(guān)于C語言中結(jié)構(gòu)體、枚舉、聯(lián)合體和位段自定義類型的相關(guān)資料,分別介紹了結(jié)構(gòu)體、枚舉、聯(lián)合體和位段等四種自定義類型,文中通過示例代碼介紹的非常詳細,需要的朋友可以參考下2021-08-08
這篇文章主要為大家介紹了C++面試八股文之STL標準模板庫使用詳解，<BR>有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進步，早日升職加薪2023-06-06
這篇文章主要為大家詳細介紹了C語言控制臺打磚塊小游戲，文中示例代碼介紹的非常詳細，具有一定的參考價值，感興趣的小伙伴們可以參考一下2022-05-05
以下是對exit與atexit的區(qū)別進行了詳細的分析介紹，需要的朋友可以過來參考下2013-09-09
關(guān)于C++智能指針shared_ptr和unique_ptr能否互轉(zhuǎn)問題
C++中的智能指針最常用的是shared_ptr和unique_ptr，C++新手最常問的問題是我從一個函數(shù)中拿到unique_ptr，但要轉(zhuǎn)成shared_ptr才能使用，要怎么轉(zhuǎn)換？同理是否能將shared_ptr轉(zhuǎn)換成unique_ptr，面對這些問題，跟隨小編一起看看吧2022-05-05