关于如何用代码表示矩阵及矩阵的一些算子（三维重建task1-1）

关于如何用代码表示矩阵及矩阵的一些算子（三维重建task1-1）

有一段日子没更新，参加了个数模美赛，弄了一些小项目的资料，过了个年，hh。
接下来的学习主要以代码解读为主，参考一下人家大佬写的代码，对自己缺少的知识进行补充，边学边写，有说的不对的地方欢迎大家指正。

代码呈上

#include <iostream>
#include <math/matrix_svd.h>
#include "math/matrix.h"
#include "math/vector.h"

int main(int argc, char *argv[])
{

    /*构建一个维度为4x5的矩阵，数据类型为double的矩阵*/
    math::Matrix<double, 4, 5> A;

    /*矩阵元素的设置和访问*/
    int id=0;
    for(int i=0; i< A.rows; i++){
        for(int j=0; j< A.cols; j++){
            A(i,j) = ++id;
            std::cout<<A(i, j)<<" ";
        }
        std::cout<<std::endl;
    }
    std::cout<<std::endl;

    /*取矩阵的列元素*/
    math::Vector<double, 4> col4 = A.col(4); // 取第5列元素
    std::cout<<"col4: "<<col4<<std::endl;

    /*取矩阵的行元素*/
    math::Vector<double, 5> row2 = A.row(2); // 取第3行元素
    std::cout<<"row2: "<<row2<<std::endl;

    // 向量的创建
    math::Vector<double, 5> v1;
    for(int i=0; i<v1.dim; i++){
        v1[i] = i;
    }

    std::cout<<"v1: ";
    for(int i=0; i<v1.dim; i++){
        std::cout<<v1[i]<<" ";
    }
    std::cout<<std::endl<<std::endl;

    //奇异值分解
    math::Matrix<double, 4, 5>U;
    math::Matrix<double, 5, 5> S, V;
    math::matrix_svd<double, 4, 5> (A,&U, &S, &V);
    std::cout<<"U: "<<U<<std::endl;
    std::cout<<"S: "<<S<<std::endl;
    std::cout<<"V: "<<V<<std::endl;

    return 0;
}

个人觉得代码的注释还是解释得挺清楚的，不懂的大家可以去查查资料。

Vector和Matrix的简单记录

1.矩阵Matrix和向量Vector都是由Matrix类构造的。向量是矩阵的特殊形式，只有一列（列向量）或者一行。
2.Matrix模板类有6个参数，其中前三个参数是必须的。前三个参数如下：

Matrix<typename Scalar,int RowsAtCompileTime,int ColsAtCompileTime >
//Matrix<类型名,矩阵行数, 矩阵列数>

欢迎讨论！