复习数学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