矩陣運(yùn)算——C語言實現(xiàn)

1、#include "stdafx.h"#include <stdlib.h>/ 結(jié)構(gòu)體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 以內(nèi)隨機(jī)正整數(shù)void

2、InputMatrix(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 Sub

3、Matrix(Matrix *T1,Matrix *T2,int BeginRow,int BeginCol,int EndRow,int EndCol); / 求T1 的子矩陣T2;測試程序參考界面如下所示：矩陣函數(shù)測試，請選擇功能，輸入對應(yīng)的數(shù)字：*1：輸入一個矩陣，求矩陣均值；2：產(chǎn)生一個隨機(jī)數(shù)矩陣，求矩陣均值；3：輸入兩個個矩陣，求矩陣和；4：輸入兩個個矩陣，求矩陣差；5：輸入兩個個矩陣，求矩陣積；6：產(chǎn)生兩個個隨機(jī)數(shù)矩陣，求矩陣和；7：產(chǎn)生兩個個隨機(jī)數(shù)矩陣，求矩陣差；8：產(chǎn)生兩個個隨機(jī)數(shù)矩陣，求矩陣積；9：求矩陣的子陣，如矩陣的2-4 行， 1-3 列的子陣；0：結(jié)束!#includ

4、e <stdio.h>#include <stdlib.h>/ 結(jié)構(gòu)體typedef struct double *mat;int row;int col; Matrix;/ 函數(shù)聲明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 DestroyMatr

5、ix(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, int B

6、eginRow, int BeginCol, int EndRow, int EndCol); / 求T1 的子矩陣T2;void NMatrix(void);/ 求逆矩陣/ 主函數(shù)int main()int inputevent;/int i, j;int row, col;Matrix T;Matrix A;Matrix B;Matrix C;printf（" 矩陣函數(shù)測試，請選擇功能，輸入對應(yīng)的數(shù)字：n"）;printf(" *printf（"1 ：輸入一個矩陣，求矩陣均值；n"）;printf（"2 ：產(chǎn)生一個隨機(jī)數(shù)矩陣，求矩

7、陣均值；n"）;printf（"3 ：輸入兩個個矩陣，求矩陣和；n"）;printf（"4 ：輸入兩個個矩陣，求矩陣差；n"）;printf（"5 ：輸入兩個矩陣，求矩陣積;"）;printf（"n6 ：產(chǎn)生兩個隨機(jī)數(shù)矩陣，求矩陣和；n"）;printf（"7 ：產(chǎn)生兩個隨機(jī)數(shù)矩陣，求矩陣差；n"）;printf（"8 ：產(chǎn)生兩個隨機(jī)數(shù)矩陣，求矩陣積;n"）;printf（"9: 求矩陣的子陣,如矩陣的2-4 行 D,1-3 列的子陣;n"）;p

8、rintf（"10 ：輸入一個方陣，求其逆矩陣n"）;printf（"0: 結(jié)束 !n"）;printf（"nn 選擇： "）;scanf("%d", &inputevent);while (inputevent != 0)if (inputevent = 1) printf(" 您要輸入的矩陣的行數(shù)和列數(shù)e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&T, row, col);Init

9、ialMatrixZero(&T, T.row, T.col);InputMatrix(&T);printf("n 矩陣的平均值為：n%lfnn", MeanMatrix(&T);DestroyMatrix(&T);if (inputevent = 2) printf(" 輸入您要產(chǎn)生的矩陣的行數(shù)和列數(shù)e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&T, row, col);InitialMatrixRand(&T

10、, row, col);MeanMatrix(&T);PrintfMatrix(&T);printf("n 矩陣的平均值為：n%lfnn", MeanMatrix(&T);DestroyMatrix(&T);if (inputevent = 3) printf("您要輸入的矩陣A的行數(shù)和列數(shù)e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixZero(&A, A.ro

11、w, A.col);InputMatrix(&A);printf(" 您要輸入的矩陣B 的行數(shù)和列數(shù)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 (AddMatri

12、x(&A, &B, &C) = 1) printf("nn 矩陣的和為：A*Bnn");PrintfMatrix(&C);else printf("nn 無法計算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);if (inputevent = 4) printf("您要輸入的矩陣A的行數(shù)和列數(shù)e.g:5,6:");scanf("%d,%d", &row, &col);Ini

13、tialMatrix(&A, row, col);InitialMatrixZero(&A, A.row, A.col);InputMatrix(&A);printf(" 您要輸入的矩陣B 的行數(shù)和列數(shù)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,

14、B.col);InitialMatrixZero(&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("您要輸入的矩陣

15、A的行數(shù)和列數(shù)e.g: 5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixZero(&A, A.row, A.col);InputMatrix(&A);printf(" 您要輸入的矩陣B 的行數(shù)和列數(shù)e.g:5,6: ");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixZer

16、o(&B, B.row, B.col);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);DestroyMatri

17、x(&B);DestroyMatrix(&C);if (inputevent = 6) printf("您要產(chǎn)生的矩陣A的行數(shù)和列數(shù)e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(&A, row, col);printf("A 為： nn");PrintfMatrix(&A);printf(" 您要產(chǎn)生的矩陣B 的行數(shù)和列數(shù)e.g:5,6: "

18、;);scanf("%d,%d", &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 矩

19、陣的和為：A+B=nn");PrintfMatrix(&C);else printf("nn 無法計算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C);if (inputevent = 7) printf("您要產(chǎn)生的矩陣A的行數(shù)和列數(shù)e.g:5,6:");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);InitialMatrixRand(

20、&A, row, col);printf(" 您要產(chǎn)生的矩陣B 的行數(shù)和列數(shù)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) pr

21、intf("n 矩陣的差為：A-B=nn");PrintfMatrix(&C);else printf("nn 無法計算！nn");DestroyMatrix(&A);DestroyMatrix(&B);DestroyMatrix(&C); if (inputevent = 8) printf("您要產(chǎn)生的矩陣A的行數(shù)和列數(shù)e.g:5,6: n");scanf("%d,%d", &row, &col);InitialMatrix(&A, row, col);I

22、nitialMatrixRand(&A, row, col);printf("A 為： nn");PrintfMatrix(&A);printf(" 您要產(chǎn)生的矩陣B 的行數(shù)和列數(shù)e.g:5,6: n");scanf("%d,%d", &row, &col);InitialMatrix(&B, row, col);InitialMatrixRand(&B, row, col);printf("B 為： nn");PrintfMatrix(&B);Initial

23、Matrix(&C, A.row, A.col);InitialMatrixZero(&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(&qu

24、ot; 對不起，該函數(shù)尚在完善中nn");if (inputevent = 10) NMatrix();if (inputevent = 0) break;n");printf（" 矩陣函數(shù)測試，請選擇功能，輸入對應(yīng)的數(shù)字：printf(n *nn");printf（"1 ：輸入一個矩陣，求矩陣均值；n"）;printf("2 ：產(chǎn)生一個隨機(jī)數(shù)矩陣，求矩陣均值；產(chǎn)一個隨機(jī)數(shù)矩陣，求矩陣均值；n");printf("3 ：輸入兩個個矩陣，求矩陣和；n");printf("4 ：輸入兩個個

25、矩陣，求矩陣差；n");printf("5 ：輸入兩個矩陣，求矩陣積;");printf("n6：產(chǎn)生兩個隨機(jī)數(shù)矩陣，求矩陣和；n");printf("7 ：產(chǎn)生兩個隨機(jī)數(shù)矩陣，求矩陣差；產(chǎn)生兩個隨機(jī)數(shù)矩陣，求矩陣差；n");printf("8 ：產(chǎn)生兩個隨機(jī)數(shù)矩陣，求矩陣積;n");printf("9: 求矩陣的子陣,如矩陣的2-4 行 D,1-3 列的子陣;n");printf("10 ：輸入一個方陣，求其逆矩陣n");printf("0: 結(jié)束 !n

26、");printf("nn 選擇： ");scanf("%d", &inputevent);return 0;/ 其他函數(shù)void InitialMatrix(Matrix *T, int row, int col) /printf(" 分配內(nèi)存中n");int i;int succ = 1;/T=(Matrix *)malloc(sizeof(Matrix);T->row = row; T->col = col;T->mat = (double *)malloc(T->row * sizeo

27、f(double *);if (T->mat = NULL) succ = 0; else for (i = 0; i < T->row; i+) T->mati = (double *)malloc(T->col * sizeof(double); if (T->mati = NULL) succ = 0; break;/if(succ=1)/ printf(" 內(nèi)存分配成功| ； ?n");/else printf(" 內(nèi)存分配失敗；n");void InitialMatrixZero(Matrix *T, int

28、 row, int col) /printf(" 矩陣初始化為零中n");int i, j;for (i = 0; i < row; i+)for (j = 0; j < 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() %

29、50; void InputMatrix(Matrix *T) printf(" 輸入矩陣：n");int i, j;for (i = 0; i < (*T).row; i+)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

30、 = 0; i < (*T).row; i+) for (j = 0; j < (*T).col; j+) printf("%lf ", (*T).matij);printf("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).ma

31、tij = (*A).matij + (*B).matij;for (i = 0; i < (*A).row; i+)for (j = 0; j < (*A).col; j+)return 1;else 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 (

32、j = 0; j < (*A).col; j+)(*C).matij = (*A).matij - (*B).matij;return 1;elseprintf(" 這兩個矩陣不能相減！n");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)/ 矩陣