C++中Eigen线性代数库的使用

参考：【技术与比赛】C++中Eigen线性代数库的使用——421施公队Clang出品_哔哩哔哩_bilibili

Eigen库安装

安装Eige

https://www.cnblogs.com/lihanwen/p/9928063.html

安装命令：

sudo apt-get install libeigen3-dev

重设路径：

sudo cp -r /usr/include/eigen3/Eigen /usr/include

#include <eigen3/Eigen/Dense>
using namespace Eigen;

Eigen库的使用

#include<eigen3/Eigen/Dense>
using namespace Eigen

template <Typename T> using Mat2 = Matrix<T,2,2>;

int main()
{
    Mat2<int> mat1;
    mat1 << 1,1,
            1,1;

    Vector2d vec1;
    vec1<< 1,2;

    MatrixXd mat2(3,3);
    mat2 <<1,2,3,
           4,5,6,
           7,8,9;

    MatrixXd mat3 = MatrixXd::Zero(3,3);        //全零矩阵
    MatrixXd mat4 = MatrixXd::Ones(3,3);        //全1矩阵
    MatrixXd mat5 = MatrixXd::Identity(3,3);    //单位矩阵
    MatrixXd mat6 = MatrixXd::Random(3,3);      //随机矩阵

    mat3.rows();            //mat3的行数    区别于row(N)
    mat3.cols();            //mat3的列数    区别于col(N)

    mat3*mat4;              //矩阵乘法
    MatrixXd mat7 = mat3.transpose();       //得到mat3的转置矩阵
    MatrixXd mat7 = mat3.conjugate();       //得到mat3的共轭矩阵
    MatrixXd mat7 = mat3.inverse();         //得到mat3的逆矩阵
    MatrixXd mat7 = mat3.adjoint();         //得到mat3的伴随矩阵

    mat3.transposeInPlace();                //让mat3变为自身的转置矩阵
    mat3 = mat3.inverse();                  //让mat3变为自身的逆矩阵

    //向量块操作 block operation
    ArrayXd vec(9);
    vec << 1,2,3,4,5,6,7,8,9;
    vec.head(3);            //取vec的前3个元素
    vec.tail(3);            //取vec的后3个元素
    vec.segment(3,3);       //取vec第4个元素开始的三个元素，即4,5,6

    //矩阵块操作 block operation
    MatrixXd mat8(3,3);
    mat8<<1,2,3,
          4,5,6,
          7,8,9;
    mat8.col(2);            //取得mat8的第三列
    mat8.row(1);            //取得mat8的第二行
    mat8.block(1,1,2,2);    //取得mat8第二行第二列开始的后一行后一列，即(5,6;8,9)
    mat8.block<2,2>(1,1);   //同上
}