【C++】矩阵的乘法

先复习一下矩阵的乘法 。

已知

求AB。 

因为矩阵A是2×3矩阵，矩阵B是3×3矩阵，A的列数等于B的行数，所以矩阵A与B可以相乘，乘积AB是一个2×3矩阵。

矩阵相乘时需要注意两点，一点是矩阵1的列数要等与矩阵2的行数，一点是矩阵相乘后的矩阵

c[i][j] = a[i][k]*b[k][j] 

 由矩阵相乘的性质可以用三层循环求出结果，注意i的循环需要放在最外层，j的循环需要放在中间层，k的循环需要放在最低层。

for(int i=0;i<ROW1;i++)
            for(int j=0;j<COL2;j++)
                for(int k=0;k<COL1;k++)
                    arr3[i][j] += arr1[i][k]*arr2[k][j];

 完整代码：

#include<iostream>
using namespace std;

//第i个矩阵的行数和列数 
#define ROW1 3
#define COL1 3
#define ROW2 3
#define COL2 3

/*
矩阵乘法 
*/
void matrix_mutiplicate(double(*arr1)[COL1], double (*arr2)[COL2], double(*arr3)[COL2]) {
	if(COL1==ROW2) {
		//矩阵相乘的性质c[i][j] += a[i][k]*b[k][j] 
		for(int i=0;i<ROW1;i++)
			for(int j=0;j<COL2;j++)
				for(int k=0;k<COL1;k++)
					//	*(* (arr3 + i)+j) += *(*(arr1 + i) + k) * *(*(arr2 +k) + j);
					arr3[i][j] += arr1[i][k]*arr2[k][j];
	}
	return;
}


/*
打印矩阵 
*/
void matrix_print(double(*arr)[COL2]) {
	for(int i=0;i<ROW1;i++) {
		for(int j=0;j<COL2;j++)
			cout<<arr[i][j]<<" ";
		cout<<endl;
	}
}

int main()
{
    double arr1[ROW1][COL1] = {1,2,3,1,2,3,1,2,3};
    double arr2[ROW2][COL2] = {1,1,1,1,1,1,1,1,1};
    double arr3[ROW1][COL2] = {0};
    matrix_mutiplicate(arr1,arr2,arr3);
    matrix_print(arr3);
	return 0;   	  		   		 
}

输出结果：

手算一下AB，不难发现结果是正确的。

总结 

我的想法是把矩阵相乘封装成一个函数，但是二维数组作为参数时列数必须是常量，这就限制了函数的调用。其实c99标准支持变量，有兴趣的可以尝试。