矩阵相乘肯定都不陌生,那么如何用C语言实现矩阵相乘呢?先看下面两个矩阵相乘:
矩阵A和B都是两行两列,矩阵相乘时,A的行元素乘以B的列元素,然后相加。即
20*2+5*1=45作为C11(矩阵C的第一行第一列)
20*1+5*4=40作为C12(矩阵C的第一行第二列)
15*2+10*1=40作为C21(矩阵C的第二行第一列)
15*1+10*4=55作为C22(矩阵C的第二行第二列)
因此我们可以用两个for循环来求解。第一个for循环指定行,第二个for循环指定列。此时还没有结束,因为是矩阵A行的各元素与矩阵B的列的各元素相乘,再相加。必须要定义一个变量k来标识A的行中个元素以及标识B的列中个元素。通过k的增加来遍历A的行中的元素和遍历B的列中的元素。
下面是代码的实现:
#include<stdio.h>
#include<iostream>
using namespace std;
#define X 3
#define Y 2
#define Z 4
int main()
{
int i, j, k;
int a[X][Y], b[Y][Z], c[X][Z];
for (i = 0; i < X; i++) //输入矩阵a,3行2列
for (j = 0; j < Y; j++)
scanf_s("%d", &a[i][j]);
for (i = 0; i < Y; i++) //输入矩阵b,2行4列
for (j = 0; j < Z; j++)
scanf_s("%d", &b[i][j]);
for (i = 0; i < X; i++) //初始化矩阵c,3行4列
for (j = 0; j < Z; j++)
c[i][j] = 0;
for (i = 0; i<X; i++) //矩阵a与矩阵b相乘
for (j = 0; j < Z; j++)
{
c[i][j] = 0;
for (k = 0; k < Y; k++)
c[i][j] += a[i][k] * b[k][j];
}
for (i = 0; i < X; i++) //输出矩阵c
for (j = 0; j < Z; j++)
printf("%3d", c[i][j]);
system("pause");
return 0;
}