亲宝软件园·资讯

展开

C++矩阵相乘小程序

超自然祈祷 人气:0

复习数学1的线性代数,矩阵相乘这块有点晕,想编个C++对象化的矩阵相乘小程序。

相乘部分

void sum(juzhen a, juzhen b, juzhen &c)
{
 int s=0;
 for (int i = 1; i <= a.m1(); i++)//A矩阵的M
  for (int j = 1; j <= b.n1(); j++)//B矩阵的S
  {
   for (k0 = 1; k0 <= a.n1(); k0++)//a.n1也就是b.m1(a的n,b的n)【行向量*列向量】
   {
    s += a.read(i,k0)*b.read(k0,j);
   }
    c.write(i, j, s);
    s = 0;
  }
}

公式:

代码中的头两个for循环就是i,j的。公式中的k从1到p求和就是里面的k0的for循环。

容易出现误解的就是公式中只是表示第“[i][j]”元素,而不是整个矩阵,整个矩阵的结果需要外面的两个for循环。

本质:这就是个p维向量(高中就记2维)的两向量相乘公式而已【结果为数,是新矩阵的一个元素】

可运行代码:

#include<iostream>
#include <string>
using namespace std;
 
class juzhen
{
private:
 int m,n;//长宽
 int num[10][10] = {0};
 string name;
 
public:
  void size(int a,int b)
  { m = a;
  n = b; }
  void set()
  {
   cout << "此矩阵规模:" << this->m <<","<< this->n << endl;//=====?
   for (int i = 1; i <= this->m; i++)
    for (int j = 1; j <= this->n; j++)
    {
     cin >> this->num[i][j];
    }
   cout << "输入完成"<< endl;
  }
  void display()
  {
   for (int i = 1; i <= this->m; i++)//===i为行号(第几行),j为列号
    for (int j = 1; j <= this->n; j++)
    {
     cout << this->num[i][j] << " ";
    if (j == this->n) cout << endl;//先输出再换行
    }
  }
  int read(int a, int b) { return num[a][b]; }//调用此函数,得[m][n]元素的值
  void write(int a, int b,int count) {  num[a][b]=count; }//第三个参数的值,传递给[a][b]元素
  int m1() { return m; }//调用得到矩阵的m
  int n1() { return n; }//调用得到矩阵的n
};
 
int m0, n0, s0, k0;//矩阵规模(容易搞混的东西)
 
void sum(juzhen a, juzhen b, juzhen &c)//矩阵相乘公式所在。。。【要改实参值的要用&引用】
{
 for (int i = 1; i <= a.m1(); i++)//A矩阵的M
  for (int j = 1; j <= b.n1(); j++)//B矩阵的S
  {
 
   int s = 0;
   for (k0 = 1; k0 <= a.n1(); k0++)//a.n1也就是b.m1(a的n,b的n)【行向量*列向量】
   {
    s += a.read(i,k0)*b.read(k0,j);//因为用了C++,所以没那么直观,就是a[i][k]*b[k][j],套个for循环求累加和(就是高中时向量的点乘)
   }
    c.write(i, j, s);
  }
}
int main()
{
 juzhen A,B,C;
 cout << "设定m,s,n。A的m*s,B的s*n(横条数*纵条数)"<<endl;
 cin >> m0>> s0>> n0;
 A.size(m0, s0);
 B.size(s0, n0);
 
 C.size(m0, n0);
 
 A.set();
 B.set();
 
 sum(A, B, C);//C=A*B
 cout << "结果C的m*n:" << endl;
 C.display();
 
 return 0;
}

样例输入及输出:

1 0          1 0           1 0
0 1      *   0 1      =   0 1

如图所示两个矩阵

懒得打了……就如图所示两个矩阵

PS:函数中形参引用真好用,过去一直不知道,省得用指针了。(不然改不了C矩阵的实际元素值)

void sum(juzhen a, juzhen b, juzhen &c) 
//矩阵相乘公式所在。。。【要改实参值的要用&引用】

加载全部内容

相关教程
猜你喜欢
用户评论