C语言 求两个矩阵相乘

  矩阵相乘肯定都不陌生，那么如何用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;
}