C语言综合实验矩阵运算

/*
  综合实验一	
    简易矩阵计算程序
	*/

#include<stdio.h>
#include<malloc.h>
#include<stdlib.h>
#include<time.h>
#include<windows.h>
#include<graphics.h>

typedef struct
{
	double** mat;
	int row,col;
}Matrix;

void wait(int ms)            //延时函数，延长输入时长的时间
{
	clock_t endwait;
	endwait = clock() + ms * CLOCKS_PER_SEC / 1000;
	while (clock() < endwait) {}
}

void InitialMatrix(Matrix *T,int row, int col)
{
	int i = 0;
	T->row = row;
	T->col = col;
	T->mat= (double**)malloc(row * sizeof(double *));
	for (i = 0; i < row; i++)
	{
		T->mat[i] = (double*)malloc(col * sizeof(double));		
	}
}

void InitialMatrixZero(Matrix* T, int row, int col)
{
	int i, j;
	InitialMatrix(T, row, col);
	for (i = 0; i < row; i++)
	{
		for (j = 0; j < col; j++)
		{
			T->mat[i][j] = 0;
		}
	}
}

void InitialMatrixRand(Matrix* T, int row, int col)
{
	int i, j;
	InitialMatrix(T, row, col);
	srand((unsigned)time(NULL));
	for (i = 0; i < row; i++)
	{
		for (j = 0; j < col; j++)
		{
			T->mat[i][j] = rand()%50+1;
		}
	}
}

void InputMatrix(Matrix* T)
{
	int i, j,row,col;
	printf("输入该矩阵的行数与列数：");
	scanf_s("%d %d", &row, &col);
	InitialMatrix(T, row, col);		//输入一个矩阵（创建一个新的矩阵之前），要先将其初始化为其分配内存空间。
	printf("请输入一个矩阵（两个元素间键入空格，换行时键入回车）：\n");
	for (i = 0; i < T->row; i++)
	{
		for (j = 0; j < T->col; j++)
		{
			scanf_s("%lf", &T->mat[i][j]);
		}
	}
}

void DestroyMatrix(Matrix* T)
{
	int i;
	for (i = 0; i < T->row; i++)
	{
		free(T->mat[i]);
	}
	free(T->mat);
}

void PrintfMatrix(Matrix* T)
{
	int i, j;
	for (i = 0; i < T->row; i++)
	{
		for (j = 0; j < T->col; j++)
		{
			printf("%d   ", (int)T->mat[i][j]);
			if (j == T->col-1)
				printf("\n");
		}
	}
}

void AddMatrix(Matrix* A, Matrix* B, Matrix* C)
{
	int i, j;
	for (i = 0; i < A->row; i++)
	{
		for (j = 0; j < B->col; j++)
		{
			C->mat[i][j] = A->mat[i][j] + B->mat[i][j];
		}
	}
}

void MinusMatrix(Matrix* A, Matrix* B, Matrix* C)
{
	int i, j;
	for (i = 0; i < A->row; i++)
	{
		for (j = 0; j < B->col; j++)
		{
			C->mat[i][j] = A->mat[i][j] - B->mat[i][j];
		}
	}
}

void MultiMatrix(Matrix* A, Matrix* B, Matrix* C)
{
	int i, j, z;
	double sum = 0;
	for (i = 0; i < A->row; i++)
	{
		for (z = 0; z < B->col; z++)
		{
			for (j = 0; j < A->col; j++)
			{
				sum = (A->mat[i][j] * B->mat[j][z]);
				C->mat[i][z] += sum;
			}
		}
	}
	C->row = A->row;
	C->col = B->col;
}

double MeanMatrix(Matrix* T)
{
	int i, j;
	double sum = 0;
	double avr = 0;
	for (i = 0; i < T->row; i++)
	{
		for (j = 0; j < T->col; j++)
		{
			sum += T->mat[i][j];
		}
	}
	avr = sum / ((double)T->row * (double)T->col);
	return avr;
}

void SubMatrix(Matrix* T1, Matrix* T2, int BeginRow, int BeginCol, int EndRow, int EndCol)
{
	int i, j,m,n;
	for (i = BeginRow - 1,m=0; i < EndRow; i++, m++)
	{
		for (j = BeginCol - 1, n = 0; j < EndCol; j++, n++)
		{
			T2->mat[m][n] = T1->mat[i][j];
		}
	}
}

int main()
{
	Matrix T, T1, T2, A, B, C;
	int n, row, col, BeginRow, EndRow, BeginCol, EndCol,check=1;
	double mean;
	system("title 简易矩阵计算程序");        //设置控制台窗口的标题
	system("color 8F");                    //设置控制台窗口的颜色
	printf("    [矩阵函数测试，请选择功能，输入对应的数字]:\n");
	printf("\n");
	printf("**************************************************\n");
	printf("\n");
	printf("	1：输入一个矩阵，求矩阵均值；\n");
	printf("	2：产生一个随机数矩阵，求矩阵均值；\n");
	printf("	3：输入两个矩阵，求矩阵和；\n");
	printf("	4：输入两个矩阵，求矩阵差；\n");
	printf("	5：输入两个矩阵，求矩阵积；\n");
	printf("	6：产生两个随机数矩阵，求矩阵和；\n");
	printf("	7：产生两个随机数矩阵，求矩阵差；\n");
	printf("	8：产生两个随机数矩阵，求矩阵积；\n");
	printf("	9：求矩阵的子阵，如矩阵的2～4行，1～3列的子阵；\n");
	printf("	0：结束程序\n");
	do 
	{
		printf("\n");
		printf("请输入对应的数字:");
		scanf_s("%d", &n);
		switch (n)
		{
		case 1:
			InputMatrix(&T);
			printf("矩阵元素平均值为\n");
			mean = MeanMatrix(&T);
			printf("%lf\n", mean);
			break;
		case 2:
			printf("下面进行矩阵的均值运算：\n");
			printf("请输入矩阵的行和列：");
			scanf_s("%d %d", &row, &col);
			InitialMatrixRand(&T, row, col);
			printf("您获得了一个随机数矩阵：\n");
			PrintfMatrix(&T);
			printf("矩阵的均值为：");
			mean = MeanMatrix(&T);
			printf("%lf\n", mean);
			break;
		case 3:
			printf("下面进行矩阵的加法运算：\n");
			printf("请输入第一个矩阵\n");
			InputMatrix(&A);
			printf("请输入第二个矩阵\n");
			InputMatrix(&B);
			InitialMatrixZero(&C, A.row, A.col);
			AddMatrix(&A, &B, &C);
			printf("相加后的矩阵如下：\n");
			PrintfMatrix(&C);
			break;
		case 4:
			printf("下面进行矩阵的减法运算：\n");
			printf("请输入第一个矩阵\n");
			InputMatrix(&A);
			printf("请输入第二个矩阵\n");
			InputMatrix(&B);
			InitialMatrixZero(&C, A.row, A.col);
			MinusMatrix(&A,&B,&C);
			printf("相减后的矩阵如下：\n");
			PrintfMatrix(&C);
			break;
		case 5:
			printf("下面进行矩阵的乘法运算：\n");
			printf("请输入第一个矩阵\n");
			InputMatrix(&A);
			printf("请输入第二个矩阵\n");
			InputMatrix(&B);
			InitialMatrixZero(&C, A.row, A.col);
			MultiMatrix(&A, &B, &C);
			printf("相乘后的矩阵如下：\n");
			PrintfMatrix(&C);
			break;
		case 6:
			printf("下面进行随机数矩阵的加法运算：\n");
			printf("请输入随机数矩阵的行和列：");
			scanf_s("%d %d", &row, &col);
			printf("产生的第一个随机数矩阵如下：\n");
			InitialMatrixRand(&T1, row, col);
			PrintfMatrix(&T1);
			wait(800);
			printf("产生的第二个随机数矩阵如下：\n");
			InitialMatrixRand(&T2, row, col);
			PrintfMatrix(&T2);
			printf("相加后的随机数矩阵如下：\n");
			InitialMatrix(&T, row, col);
			AddMatrix(&T1, &T2, &T);
			PrintfMatrix(&T);
			break;
		case 7:
			printf("下面进行随机数矩阵的减法运算：\n");
			printf("请输入随机数矩阵的行和列：");
			scanf_s("%d %d", &row, &col);
			printf("产生的第一个随机数矩阵如下：\n");
			InitialMatrixRand(&T1, row, col);
			PrintfMatrix(&T1);
			wait(800);            //产生时间间隔，使随机数种子改变，达到真正随机的效果
			printf("产生的第二个随机数矩阵如下：\n");
			InitialMatrixRand(&T2, row, col);
			PrintfMatrix(&T2);
			printf("相减后的随机数矩阵如下：\n");
			InitialMatrix(&T, row, col);
			MinusMatrix(&T1, &T2, &T);
			PrintfMatrix(&T);
			break;
		case 8:
			printf("下面进行随机数矩阵的乘法运算：\n");
			printf("请输入随机数矩阵的行和列：");
			scanf_s("%d %d", &row, &col);
			printf("产生的第一个随机数矩阵如下：\n");
			InitialMatrixRand(&T1, row, col);
			PrintfMatrix(&T1);
			wait(800);
			printf("产生的第二个随机数矩阵如下：\n");
			InitialMatrixRand(&T2, row, col);
			PrintfMatrix(&T2);
			InitialMatrixZero(&T, row, col);
			MultiMatrix(&T1, &T2, &T);
			printf("相乘后的随机数矩阵如下：\n");
			PrintfMatrix(&T);
			break;
		case 9:
			printf("下面进行随机数矩阵的子矩阵运算：\n");
			printf("请输入矩阵的行和列：");
			scanf_s("%d %d", &row, &col);
			InitialMatrixRand(&T1,row,col);
			printf("您获得了一个如下的随机数矩阵：\n");
			PrintfMatrix(&T1);
			printf("请输入子矩阵的起始行，结束行：");
			scanf_s("%d %d", &BeginRow, &EndRow);
			printf("请输入子矩阵的起始列，结束列：");
			scanf_s("%d %d", &BeginCol, &EndCol);
			InitialMatrix(&T2, EndRow - BeginRow + 1, EndCol - BeginCol + 1);
			SubMatrix(&T1, &T2, BeginRow, BeginCol, EndRow, EndCol);
			printf("按要求所求得的子矩阵如下：\n");
			PrintfMatrix(&T2);
			break;
		case 0:
			printf("\n");
			printf("	[!*****本程序结束*****!]");
			printf("\n");
			check = 0;
		}
	} while (check);
	return 0;
}