复习数学1的线性代数，矩阵相乘这块有点晕，想编个C++对象化的矩阵相乘小程序。
相乘部分
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
|
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维）的两向量相乘公式而已【结果为数，是新矩阵的一个元素】
可运行代码：
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
#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矩阵的实际元素值）
|
1
2
|
void sum(juzhen a, juzhen b, juzhen &c)
//矩阵相乘公式所在。。。【要改实参值的要用&引用】
|
以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持服务器之家。
原文链接：https://blog.csdn.net/sinat_27382047/article/details/76973054