armadillo库之对称矩阵的特征值与特征向量
2015-05-03 17:09
926 查看
# 上面的测试用例能正常运行
#include <iostream>
#include <armadillo>
#include <vector>
// using namespace std;
// using namespace arma;
int main(int argc, char** argv)
{
// mat A = randu<mat>(4,5);
// mat B = randu<mat>(4,5);
std::vector<int> ivec;
ivec.push_back(0);
ivec.push_back(1);
ivec.push_back(2);
arma::vec x(3);
x[0] = 0;
x[1] = 1;
x[2] = 2;
//x.print();
//x = -1*x;
arma::vec y(3);
y[0] = 1;
y[1] = 2;
y[2] = 3;
//x = x + y;
arma::rowvec z(3);
z[0] = 0;
z[1] = 1;
z[2] = 2;
//(x*z).print();
arma::mat a = x*z;
//a.print();
arma::vec eigval;
arma::mat eigvec;
eig_sym(eigval, eigvec, a);
eigval.print();
eigvec.print();
//std::cout<<a.n_rows<<std::endl;
return 0;
}
#include <iostream>
#include <armadillo>
#include <vector>
// using namespace std;
// using namespace arma;
int main(int argc, char** argv)
{
// mat A = randu<mat>(4,5);
// mat B = randu<mat>(4,5);
std::vector<int> ivec;
ivec.push_back(0);
ivec.push_back(1);
ivec.push_back(2);
arma::vec x(3);
x[0] = 0;
x[1] = 1;
x[2] = 2;
//x.print();
//x = -1*x;
arma::vec y(3);
y[0] = 1;
y[1] = 2;
y[2] = 3;
//x = x + y;
arma::rowvec z(3);
z[0] = 0;
z[1] = 1;
z[2] = 2;
//(x*z).print();
arma::mat a = x*z;
//a.print();
arma::vec eigval;
arma::mat eigvec;
eig_sym(eigval, eigvec, a);
eigval.print();
eigvec.print();
//std::cout<<a.n_rows<<std::endl;
return 0;
}
相关文章推荐
- 数组分析 jacobi法求实对称矩阵全部特征值与特征向量 MATLAB实现
- Jacobi法确定对称矩阵特征值与特征向量
- 实对称阵属于不同特征值的的特征向量是正交的
- 特征值与特征向量
- 再谈特征值与特征向量的求取
- 特征值与特征向量概念
- 特征值与特征向量的求法
- 第一次上课:特征值特征向量的几何直观
- 特征向量和特征值的几何意义
- Jacobian矩阵、Hessian矩阵、特征值、特征向量
- 说文解字----矩阵分析(二)特征值特征向量 奇异值分解(SVD)
- PCA、协方差矩阵、特征向量、特征值之间的关系
- 特征向量与特征值
- 特征值与特征向量理解
- 线性代数(三十) :特征值与特征向量的性质
- 特征值、特征向量的物理、几何意义
- 特征值与特征向量
- 浅浅地聊一下矩阵与线性映射及矩阵的特征值与特征向量
- 特征值与特征向量的几何意义
- 证明:实对称矩阵中,属于不同特征值的特征向量相互正交