视觉SLAM十四讲第三章习题５

假设有一个大矩阵，想把它的左上角３×３的块取出来（这里没有将其赋值为单位矩阵）

#include <iostream>
#include <Eigen/Core>
#include <Eigen/Geometry>
#define MATRIX_SIZE 4
using namespace std;
/*提取大矩阵左上角３×３矩阵，有两种方式：
 * 1.直接从０－２循环遍历大矩阵的前三行和三列
 * ２．用矩阵变量.block(0,0,3,3) 从左上角００位置开始取３行３列*/
int main(int argc, char **argv) {
    //设置输出小数点后３位
    cout.precision(3);
    //用随机数来初始化矩阵matrix_NN
    Eigen::Matrix<double,MATRIX_SIZE,MATRIX_SIZE> matrix_NN =
            Eigen::MatrixXd::Random(MATRIX_SIZE,MATRIX_SIZE);
    //用单位矩阵来初始化矩阵matrix_3d
    Eigen::Matrix3d matrix_3d = Eigen::Matrix3d::Identity();
    //方法一：循环遍历矩阵的三行三列
    for (int i = 0; i < 3; ++i) {　　//遍历行
        for (int j = 0; j < 3; ++j) {　　//遍历列
            matrix_3d(i,j) = matrix_NN(i,j);　//将matrix_NN矩阵中前３行前３列的值依次赋值给matrix_3d
        }
    }
    cout << "大矩阵" <<endl;
    cout<<matrix_NN <<endl;   
    cout << "从大矩阵中取值后的矩阵为:"<<endl;
    cout<< matrix_3d <<endl;　　
    //方法二：使用.block()函数
    cout<<"大矩阵为"<<endl;
    cout<<matrix_NN<<endl;
    //使用.block(startRow,startCol,blockRows,blockCols)函数取出矩阵
    matrix_3d = matrix_NN.block(0,0,3,3);
    //matrix_3d = Eigen::Matrix3d::Identity();
    cout << "从大矩阵中提取３行３列后的矩阵为:"<< endl;
    cout << matrix_3d<<endl;
    return 0;
}

CMakeLists.txt

cmake_minimum_required(VERSION 3.17)
project(extractmatrix)

set(CMAKE_CXX_STANDARD 14)

include_directories("/usr/include/eigen3")
add_executable(extracting extracting.cpp)

参考自https://www.cnblogs.com/newneul/p/8306430.html