矩阵运算——C语言实现

1、#include "stdafx.h"#include <stdlib.h>/结构体typedef structdouble * mat;int row;int col;Matrix;void InitialMatrix(Matrix *T, int row,int col); /只分配空间不初始化；void InitialMatrixZero(Matrix *T,int row, int col); /初始化为0void InitialMatrixRand(Matrix *T,int row, int col); /初始化为50以内随机正整数void Inpu

2、tMatrix(Matrix *T);/键盘输入矩阵void DestroyMatrix(Matrix *T);/ 释放矩阵空间void PrintfMatrix(Matrix *T);/矩阵输出int AddMatrix(Matrix *A,Matrix *B,Matrix *C); / 矩阵加int MinusMatrix(Matrix *A,Matrix *B,Matrix *C); / 矩阵减int MultiMatrix(Matrix *A,Matrix *B,Matrix *C); /矩阵乘法double MeanMatrix(Matrix *T); /矩阵元素均值int SubM

3、atrix(Matrix *T1,Matrix *T2,int BeginRow,int BeginCol,int EndRow,int EndCol); /求T1的子矩阵T2;测试程序参考界面如下所示： 矩阵函数测试，请选择功能，输入对应的数字： * 1：输入一个矩阵，求矩阵均值； 2：产生一个随机数矩阵，求矩阵均值； 3：输入两个个矩阵，求矩阵和； 4：输入两个个矩阵，求矩阵差； 5：输入两个个矩阵，求矩阵积； 6：产生两个个随机数矩阵，求矩阵和； 7：产生两个个随机数矩阵，求矩阵差； 8：产生两个个随机数矩阵，求矩阵积； 9：求矩阵的子阵，如矩阵的2-4行，1-3列的子阵； 0：结束!#

4、include <stdio.h>#include <stdlib.h>/结构体typedef struct double *mat;int row;int col; Matrix;/函数声明void InitialMatrix(Matrix *T, int row, int col);void InitialMatrixZero(Matrix *T, int row, int col);void InitialMatrixRand(Matrix *T, int row, int col);void InputMatrix(Matrix *T);void Destroy

5、Matrix(Matrix *T);/ 释放矩阵空间void PrintfMatrix(Matrix *T);/矩阵输出int AddMatrix(Matrix *A, Matrix *B, Matrix *C); / 矩阵加int MinusMatrix(Matrix *A, Matrix *B, Matrix *C); / 矩阵减int MultiMatrix(Matrix *A, Matrix *B, Matrix *C); /矩阵乘法double MeanMatrix(Matrix *T); /矩阵元素均值/int SubMatrix(Matrix *T1, Matrix *T2, i

6、nt BeginRow, int BeginCol, int EndRow, int EndCol); /求T1的子矩阵T2;void NMatrix(void);/求逆矩阵 /主函数int main()int inputevent;/int i, j;int row, col;Matrix T;Matrix A;Matrix B;Matrix C;printf(" 矩阵函数测试，请选择功能，输入对应的数字：n");printf(" *nn");printf("1：输入一个矩阵，求矩阵均值；n");printf("2：产生一

7、个随机数矩阵，求矩阵均值；n");printf("3：输入两个个矩阵，求矩阵和；n");printf("4：输入两个个矩阵，求矩阵差；n");printf("5：输入两个矩阵，求矩阵积;");printf("n6：产生两个随机数矩阵，求矩阵和；n");printf("7：产生两个随机数矩阵，求矩阵差；n");printf("8：产生两个随机数矩阵，求矩阵积;n");printf("9:求矩阵的子阵,如矩阵的2-4行D,1-3列的子阵;n");pr

8、intf("10：输入一个方阵，求其逆矩阵n");printf("0:结束!n");printf("nn选择：");scanf("%d", &inputevent);while (inputevent != 0)if (inputevent = 1) printf("您要输入的矩阵的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&T, row, col);InitialMatr

9、ixZero(&T, T.row, T.col);InputMatrix(&T);printf("n矩阵的平均值为：n%lfnn", MeanMatrix(&T);DestroyMatrix(&T);if (inputevent = 2) printf("输入您要产生的矩阵的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&T, row, col);InitialMatrixRand(&T, row, co

10、l);MeanMatrix(&T);PrintfMatrix(&T);printf("n矩阵的平均值为：n%lfnn", MeanMatrix(&T);DestroyMatrix(&T);if (inputevent = 3) printf("您要输入的矩阵A的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixZero(&A, A.row, A.col);

11、InputMatrix(&A);printf("您要输入的矩阵B的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixZero(&B, B.row, B.col);InputMatrix(&B);InitialMatrix(&C, B.row, B.col);InitialMatrixZero(&C, C.row, C.col);if (AddMatrix(&A, &a

12、mp;B, &C) = 1) printf("nn矩阵的和为：A*Bnn");PrintfMatrix(&C);else printf("nn无法计算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);if (inputevent = 4) printf("您要输入的矩阵A的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&am

13、p;A, row, col);InitialMatrixZero(&A, A.row, A.col);InputMatrix(&A);printf("您要输入的矩阵B的行数和列数e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixZero(&B, B.row, B.col);InputMatrix(&B);InitialMatrix(&C, B.row, B.col);InitialMa

14、trixZero(&C, C.row, C.col);if (MinusMatrix(&A, &B, &C) = 1) printf("n矩阵的差为：A-B=nn");PrintfMatrix(&C);else printf("nn无法计算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);if (inputevent = 5) int able;printf("您要输入的矩阵A的行数和列数e.g: 5,6: &q

15、uot;);scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixZero(&A, A.row, A.col);InputMatrix(&A);printf("您要输入的矩阵B的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixZero(&B, B.row, B.c

16、ol);InputMatrix(&B);InitialMatrix(&C, A.row, B.col);InitialMatrixZero(&C, C.row, C.col);able = MultiMatrix(&A, &B, &C);if (able = 1) printf("n积为：A*Bnn");PrintfMatrix(&C);else printf("nn无法计算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix

17、(&C);if (inputevent = 6) printf("您要产生的矩阵A的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(&A, row, col);printf("A为：nn");PrintfMatrix(&A);printf("您要产生的矩阵B的行数和列数e.g:5,6: ");scanf("%d,%d"

18、, &row, &col);InitialMatrix(&B, row, col);InitialMatrixRand(&B, row, col);printf("B为：nn");PrintfMatrix(&B);InitialMatrix(&C, A.row, A.col);InitialMatrixZero(&C, C.row, C.col);if (AddMatrix(&A, &B, &C) = 1) printf("n矩阵的和为：A+B=nn");PrintfMatr

19、ix(&C);else printf("nn无法计算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);if (inputevent = 7) printf("您要产生的矩阵A的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(&A, row, col);printf(&quo

20、t;您要产生的矩阵B的行数和列数e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixRand(&B, row, col);InitialMatrix(&C, A.row, A.col);InitialMatrixZero(&C, C.row, C.col);if (MinusMatrix(&A, &B, &C) = 1) printf("n矩阵的差为：A-B=nn")

21、;PrintfMatrix(&C);else printf("nn无法计算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);if (inputevent = 8) printf("您要产生的矩阵A的行数和列数e.g:5,6: n");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(&A, row, col)

22、;printf("A为：nn");PrintfMatrix(&A);printf("您要产生的矩阵B的行数和列数e.g:5,6: n");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixRand(&B, row, col);printf("B为：nn");PrintfMatrix(&B);InitialMatrix(&C, A.row, A.col);InitialMatr

23、ixZero(&C, C.row, C.col);if (MultiMatrix(&A, &B, &C) = 1) printf("n积为：A*B=nn");PrintfMatrix(&C);else printf("nn无法计算；nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);if (inputevent = 9) printf("对不起，该函数尚在完善中nn");if (inputevent = 1

24、0) NMatrix();if (inputevent = 0)break;printf(" 矩阵函数测试，请选择功能，输入对应的数字：n");printf(" *nn");printf("1：输入一个矩阵，求矩阵均值；n");printf("2：产生一个随机数矩阵，求矩阵均值；n");printf("3：输入两个个矩阵，求矩阵和；n");printf("4：输入两个个矩阵，求矩阵差；n");printf("5：输入两个矩阵，求矩阵积;");printf(

25、"n6：产生两个随机数矩阵，求矩阵和；n");printf("7：产生两个随机数矩阵，求矩阵差；n");printf("8：产生两个随机数矩阵，求矩阵积;n");printf("9:求矩阵的子阵,如矩阵的2-4行D,1-3列的子阵;n");printf("10：输入一个方阵，求其逆矩阵n");printf("0:结束!n");printf("nn选择：");scanf("%d", &inputevent);return 0;/其

26、他函数void InitialMatrix(Matrix *T, int row, int col)/printf("分配内存中.n");int i;int succ = 1;/T=(Matrix *)malloc(sizeof(Matrix);T->row = row;T->col = col;T->mat = (double *)malloc(T->row * sizeof(double *);if (T->mat = NULL) succ = 0;else for (i = 0; i < T->row; i+) T->m

27、ati = (double *)malloc(T->col * sizeof(double);if (T->mati = NULL) succ = 0;break;/if(succ=1)/ printf("内存分配成功|；?n");/else printf("内存分配失败；n");void InitialMatrixZero(Matrix *T, int row, int col) /printf("矩阵初始化为零中.n");int i, j;for (i = 0; i < row; i+)for (j = 0; j

28、 < col; j+)T->matij = 0;/printf("矩阵初始化为零矩阵成功；n");void InitialMatrixRand(Matrix *T, int row, int col) int i, j;for (i = 0; i < row; i+)for (j = 0; j < col; j+)(*T).matij = rand() % 50;void InputMatrix(Matrix *T) printf("输入矩阵：n");int i, j;for (i = 0; i < (*T).row; i+

29、)for (j = 0; j < (*T).col; j+)scanf("%lf", &(*T).matij);void DestroyMatrix(Matrix *T)int i;for (i = 0; i < (*T).row; i+)free(*T).mati);void PrintfMatrix(Matrix *T)int i, j;for (i = 0; i < (*T).row; i+)for (j = 0; j < (*T).col; j+)printf("%lf ", (*T).matij);printf(

30、"n");int AddMatrix(Matrix *A, Matrix *B, Matrix *C)int i, j;if (*A).row = (*B).row && (*A).col = (*B).col)for (i = 0; i < (*A).row; i+)for (j = 0; j < (*A).col; j+)(*C).matij = (*A).matij + (*B).matij;for (i = 0; i < (*A).row; i+)for (j = 0; j < (*A).col; j+)return 1;el

31、se printf("这两个矩阵不能相加！n");return 0;int MinusMatrix(Matrix *A, Matrix *B, Matrix *C)int i, j;if (*A).row = (*B).row && (*A).col = (*B).col)for (i = 0; i < (*A).row; i+)for (j = 0; j < (*A).col; j+)(*C).matij = (*A).matij - (*B).matij;return 1;elseprintf("这两个矩阵不能相减！n")

32、;return 0;int MultiMatrix(Matrix *A, Matrix *B, Matrix *C)int i=0, j=0;int k = 0;if (*A).col = (*B).row)for (i = 0; i < (*A).row; i+) for (j = 0; j < (*B).col; j+)for(k=0;k <(A->col);k+)(*C).matij += (*A).matik * (*B).matkj;return 1; elseprintf("这两个矩阵不能相乘！n");return 0;double MeanMatrix(Matrix *T) /矩阵元素均值double mean;double sum = 0;in