【线性代数｜C++】同型矩阵及矩阵加法

1 矩阵的定义

1.1 定义

• 由$m\times n$个数$a_{ij}(i=1,2,\cdots ,m;j=1,2,\cdots ,n)$排成的$m$行$n$列的数表$$\begin{array}{cccc}{a}_{11}& a_{12}& \cdots & a_{1n}\\ a_{21}& a_{22}& \cdots & a_{2n}\\ ⋮& ⋮& & ⋮\\ a_{m1}& a_{m2}& \cdots & a_{mn}\end{array}$$称为$m$行$n$列矩阵，简称$m\times n$矩阵,记作$$A=\left[\begin{array}{cccc}{a}_{11}& a_{12}& \cdots & a_{1n}\\ a_{21}& a_{22}& \cdots & a_{2n}\\ ⋮& ⋮& & ⋮\\ a_{m1}& a_{m2}& \cdots & a_{mn}\end{array}\right].$$

1.2 C++代码

//新建CMatrix.h文件，在其中定义矩阵类CMatrix，用于管理矩阵算法
#ifndef CMATRIX_H_
#define CMATRIX_H_

#include <vector>
#include <cmath>

using namespace std;

class CMatrix
{
public:
  CMatrix();  //构造函数
  ~CMatrix(); //析构函数   
};
#endif

2 同型矩阵

2.1 定义

• 两个矩阵的行数相等、列数也相等时，称其为同型矩阵.

2.2 C++代码

//==========CMatrix.h==========
//成员函数（静态）名：IsHomomorphicMatrix
//功能：判断两矩阵是否为同型矩阵
//形参vvMatrixInput1：被判断矩阵之一
//形参vvMatrixInput2：被判断矩阵之二
static bool IsHomomorphicMatrix
(
  const vector<vector<int>> &vvMatrixInput1,
  const vector<vector<int>> &vvMatrixInput2
);


//==========CMatrix.cpp==========
//IsHomomorphicMatrix函数定义
bool CMatrix::IsHomomorphicMatrix
(
  const vector<vector<int>> &vvMatrixInput1,
  const vector<vector<int>> &vvMatrixInput2
)
{
  //判断两矩阵（vector变量）的行数是否相等
  if(vvMatrixInput1.size() != vvMatrixInput2.size())
  {
    return false;
  }

  //判断两vector变量所有行中列数是否相等，也实现了对两矩阵列数是否相等的判断
  for(unsigned int i = 0; i < vvMatrixInput1.size(); i++)
  {
    if(vvMatrixInput1[0].size() != vvMatrixInput1[i].size())
    {
      return false;
    }

    if(vvMatrixInput1[0].size() != vvMatrixInput2[i].size())
    {
      return false;
    }
  }

  return true;
}

3 矩阵加法

3.1 规则

• 设有两个$m\times n$矩阵$A=(a_{ij})$和$B=(b_{ij})$，那么矩阵$A$和$B$的和记作$A+B$，规定为$$A+B=\left[\begin{array}{cccc}{a}_{11}+b_{11}& a_{12}+b_{12}& \cdots & a_{1n}+b_{1n}\\ a_{21}+b_{21}& a_{22}+b_{22}& \cdots & a_{2n}+b_{2n}\\ ⋮& ⋮& & ⋮\\ a_{m1}+b_{m1}& a_{m2}+b_{m2}& \cdots & a_{mn}+b_{mn}\end{array}\right].$$

3.2 C++代码

//==========CMatrix.h==========
//成员函数（静态）名：MatrixAddition
//功能：将两同型矩阵相加
//形参vvMatrixInput1：加数（矩阵）1
//形参vvMatrixInput2：加数（矩阵）2
//形参vvMatrixOutput：和（矩阵）
static bool MatrixAddition
(
  const vector<vector<int>> &vvMatrixInput1,
  const vector<vector<int>> &vvMatrixInput2,
  vector<vector<int>> &vvMatrixOutput
);

//==========CMatrix.cpp==========
//MatrixAddition函数定义
bool CMatrix::MatrixAddition
(
  const vector<vector<int>> &vvMatrixInput1,
  const vector<vector<int>> &vvMatrixInput2,
  vector<vector<int>> &vvMatrixOutput
)
{
  //首先判断输入矩阵是否为同型矩阵
  if (!IsHomomorphicMatrix(vvMatrixInput1, vvMatrixInput2))
  {
    return false;
  }
  
  vvMatrixOutput.clear();
  vector<int> vRow;
  for (unsigned int i = 0; i < vvMatrixInput1.size(); i++)
  {
    vRow.clear();
    for (unsigned int j = 0; j < vvMatrixInput1[i].size(); j++)
    {
      vRow.push_back(vvMatrixInput1[i][j] + vvMatrixInput2[i][j]);
    }
    vvMatrixOutput.push_back(vRow);
  }    

  return true;
}

//==========test.cpp（测试）==========
#include <iostream>
#include <iomanip>
#include <vector>

#include "CMatrix.h"

using namespace std;

int main()
{
  vector<vector<int>> vvMatrix1{{3,4,5},{6,7,9}};
  vector<vector<int>> vvMatrix2{{2,5,6},{7,9,4}};
  vector<vector<int>> vvRet;

  if(!CMatrix::MatrixAddition(vvMatrix1, vvMatrix2, vvRet))
  {
    cout << "两矩阵非同型矩阵，不能相加" << endl;
  }
  else
  {
    for (unsigned int i = 0; i < vvRet.size(); i++)
    {
      for (unsigned int j = 0; j < vvRet[i].size(); j++)
      {
        cout << setw(3) << vvRet[i][j] << "  ";
      }
      cout << endl;
    }    
  }

  return 0;
}

---

1. 引用文献：《工程数学 线性代数（第五版）》同济大学数学系编，高等教育出版社
2. 以上内容均为个人学习、练习的记录，如有错误，欢迎指出。