【线性代数|C++】矩阵乘法

1 数与矩阵相乘

1.1 定义

• 数$\lambda$与矩阵$A$的乘积记作$\lambda A$或$A\lambda$，规定$$\lambda A=A\lambda =\left[\begin{array}{cccc}\lambda a_{11}& \lambda a_{12}& \cdots & \lambda a_{1n}\\ \lambda a_{21}& \lambda a_{22}& \cdots & \lambda a_{2n}\\ ⋮& ⋮& & ⋮\\ \lambda a_{m1}& \lambda a_{m2}& \cdots & \lambda a_{mn}\end{array}\right].$$

1.2 C++代码

1.2.1 增加成员函数MatMulti()

• 在CMatrix.h中声明函数

  //成员函数（静态）名：MatMulti
  //功能：求数与矩阵相乘的乘积
  //形参fNum：数
  //形参vvMatInput：参与计算的矩阵
  //形参vvMatOutput：结果矩阵
  static bool MatMulti
  (
    double fNum,
    const vector<vector<double>> &vvMatInput,
    vector<vector<double>> &vvMatOuput 
  );

• 在CMatrix.cpp中定义函数

bool CMatrix::MatMulti
(
  double fNum,
  const vector<vector<double>> &vvMatInput,
  vector<vector<double>> &vvMatOuput 
)
{
  vvMatOuput.clear();
  vvMatOuput = vvMatInput;
  for (unsigned int i = 0; i < vvMatOuput.size(); i++)
  {
    for (unsigned int j = 0; j < vvMatOuput[i].size(); j++)
    {
      vvMatOuput[i][j] *= fNum; 
    }
  }  
  return true;
}

• 在test.cpp中测试函数

#include <iostream>
#include <iomanip>
#include <vector>

#include "CMatrix.h"

using namespace std;

int main()
{
  vector<vector<double>> vvMat1{{ 2, 4},
                                {-3,-6}};

  vector<vector<double>> vvMatRet;

  CMatrix::MatMulti(5.0, vvMat1, vvMatRet);

  for (unsigned int i = 0; i < vvMatRet.size(); i++)
  {
    for (unsigned int j = 0; j < vvMatRet[i].size(); j++)
    {
        cout << setw(5) << vvMatRet[i][j];
    }
    cout << endl;
  }

  cout << "DONE" << endl; 

  return 0;
}

2 矩阵与矩阵相乘

2.1 定义

• 设$A=(a_{ij})$是一个$m\times s$矩阵，$B=(b_{ij})$是一个$s\times n$矩阵，那么规定矩阵$A$与矩阵$B$的乘积是一个$m\times n$矩阵$C=(c_{ij})$，其中$$c_{ij}=a_{i1}{b}_{1j}+a_{i2}{b}_{2j}+\cdots +a_{is}{b}_{sj}=\sum _{k=1}^s{a}_{ik}{b}_{kj}$$$$(i=1,2,\cdots ,m;j=1,2,\cdots ,n)$$并把此乘积记作$C=AB.$
• 注意：只有当第一个矩阵（左矩阵）的列数等于第二个矩阵（右矩阵）的行数时，两个矩阵才能相乘.

2.2 C++代码

2.2.1 增加成员函数GetMatRowsCols()

矩阵与矩阵相乘前，需首先判断左矩阵的列数是否等于右矩阵的行数，为此在类CMatrix中增加一个读取矩阵行数和列数的成员函数GetMatRowsCols().

• 在CMatrix.h中声明函数

//成员函数（静态）名：GetMatRowsCols
//功能：读取矩阵的行数和列数，如果非矩阵结构（如：各行元素数量不等）则返回false
//形参vvMatInput：待读取行数和列数的矩阵
//形参iRowsNum：矩阵的行数
//形参iColsNum：矩阵的列数
static bool GetMatRowsCols
(
  const vector<vector<double>> &vvMatInput,
  int &iRowsNum,
  int &iColsNum
);

• 在CMatrix.cpp中定义函数

bool CMatrix::GetMatRowsCols
(
  const vector<vector<double>> &vvMatInput,
  int &iRowsCount,
  int &iColsCount
)
{
  iRowsCount = vvMatInput.size();
  if (0 == iRowsCount)
  {
    iColsCount = 0;
    return false;
  }    

  iColsCount = vvMatInput[0].size();
  if (0 == iColsCount)
  {
    iRowsCount = 0;
    return false;
  }    
  
  for (unsigned int i = 1; i < vvMatInput.size(); i++)
  {
    if (iColsCount != vvMatInput[i].size())
    {
      iRowsCount = 0;
      iColsCount = 0;
      return false;
    }
  } 
  
  return true;
}

2.2.2 增加成员函数MatrixMulti()

矩阵乘矩阵函数与数乘矩阵函数为重载函数，通过形参进行区分。

• 在CMatrix.h中声明函数

//成员函数（静态）名：MatMulti
//功能：求矩阵与矩阵相乘的乘积
//形参vvLMat：左矩阵
//形参vvRMat：右矩阵
//形参vvMatRet：计算结果
static bool MatrixMulti
(
  const vector<vector<double>> &vvLMat,
  const vector<vector<double>> &vvRMat,
  vector<vector<double>> &vvMatRet
); 

• 在CMatrix.cpp中定义函数

bool CMatrix::MatrixMulti
(
  const vector<vector<double>> &vvLMat,
  const vector<vector<double>> &vvRMat,
  vector<vector<double>> &vvMatRet
)
{
  //首先判断左矩阵的列数是否与右矩阵的行数相等
  int iLMatRowsCount, iLMatColsCount; 
  int iRMatRowsCount, iRMatColsCount;
  if (!GetMatRowsCols(vvLMat, iLMatRowsCount, iLMatColsCount))
    return false;
  
  if (!GetMatRowsCols(vvRMat, iRMatRowsCount, iRMatColsCount))
    return false;

  if (iLMatColsCount != iRMatRowsCount)
    return false;
  
  vvMatRet.clear();
  vector<double> v0;  //用于临时存储矩阵行元素
  double f0;//用于临时存储矩阵元素值
  for (int i = 0; i < iLMatRowsCount; i++)
  {
    v0.clear();
    for (int j = 0; j < iRMatColsCount; j++)
    {
      f0 = 0;
      for (int n = 0; n < iLMatColsCount; n++)
      {
        f0 += (vvLMat[i][n] * vvRMat[n][j]);
      }
      v0.push_back(f0);
    } 
    vvMatRet.push_back(v0) ; 
  }  
  
  return true;
}

• 在text.cpp中测试函数

#include <iostream>
#include <iomanip>
#include <vector>

#include "CMatrix.h"

using namespace std;

int main()
{
  vector<vector<double>> vvMat1{{ 2, 4},
                                {-3,-6}};

  vector<vector<double>> vvMat2{{-2, 4},
                                { 1,-2}};
  vector<vector<double>> vvMatRet;

  if(!CMatrix::MatrixMulti(vvMat1, vvMat2, vvMatRet))
  {
    cout << "输入矩阵不具备相乘的条件" << endl;
  }
  else
  {
    for (unsigned int i = 0; i < vvMatRet.size(); i++)
    {
      for (unsigned int j = 0; j < vvMatRet[i].size(); j++)
      {
         cout << setw(3) << vvMatRet[i][j];
      }
      cout << endl;
    }    
  }

  cout << "DONE" << endl; 

  return 0;
}

---

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