6.打印3*3矩阵乘法结果

编程求3*3矩阵的乘法：要求随机生成两个3*3矩阵（成员数值范围在-20到100）并打印，并输出其结果。

#include<iostream>
#include<cstdlib>
#include<ctime>
#include<cmath>
using namespace std;

//生成两个随机数组 
void grmat(int a[][3], int b[][3], int mat_size)
{
	int i,j;
	srand(time(NULL));//初始化随机种子 
	for(i=0;i<mat_size;i++)//产生a[3][3]随机矩阵 
	{
		for(j=0;j<mat_size;j++)
		{
			a[i][j]=rand()%121-20;//生成随机数 ，数值范围为 -20~100 ,总的数量为121个，所以rand()%121的结果限制在0~120之间，再减去20，限制范围就是 -20~100. 
		}
	}
	for(i=0;i<mat_size;i++)//生成b[3][3]随机矩阵 
	{
		for(j=0;j<mat_size;j++)
		{
			b[i][j]=rand()%121-20;//产生随机数  ，数值范围为 -20~100 
		}
	}
	
	cout<<"a[3][3]矩阵如下："<<endl; 
	for(i=0;i<mat_size;i++)//先打印a[3][3]矩阵 
	{
		for(j=0;j<mat_size;j++)
		{
			cout<<a[i][j]<<"\t";
		}
	cout<<endl;
	}
	
	cout<<"b[3][3]矩阵如下："<<endl; 
	for(i=0;i<mat_size;i++)//打印b[3][3]矩阵 
	{
		for(j=0;j<mat_size;j++)
		{
			cout<<b[i][j]<<"\t";
		}
	cout<<endl;
	}
}

//求矩阵乘法 
void mat_cheng(int a[][3], int b[][3], int c[][3], int mat_size)
{
	int i,j;
	for(i=0;i<mat_size;i++)//循环求矩阵中每一个元素的值 
	{
		for(j=0;j<mat_size;j++)
		{
			c[i][j]=0;//初始化 
		for(int k=0;k<mat_size;k++)
		{
			c[i][j]+=a[i][k]*b[k][j];//根据矩阵相乘的规则，左边矩阵决定行数，右边矩阵决定列数，而且左边矩阵列数等于右边矩阵行数  
		}
		}
	}
	
	cout<<"a[3][3]与b[3][3]相乘后得到矩阵c[3][3]如下："<<endl; 
	for(i=0;i<mat_size;i++)
	{
		for(j=0;j<mat_size;j++)
		{
			cout<<c[i][j]<<"\t";
		}
		cout<<endl;
	}
}

int  main()
{
	int a[3][3], b[3][3], c[3][3];//定义三个二维数组 ，即可看成矩阵 
	int mat_size = 3;
   	grmat(a,b,mat_size);
   	mat_cheng(a, b, c, mat_size);
    return 0;
}