矩阵运算——C语言实现

} } //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 < col; j++) T->mat[i][j] = 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).mat[i][j] = rand() % 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(\}

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

void PrintfMatrix(Matrix *T) {

int i, j; for (i = 0; i < (*T).row; i++) { for (j = 0; j < (*T).col; j++) printf(\ \ printf(\ } }

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).mat[i][j] = (*A).mat[i][j] + (*B).mat[i][j]; 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 (j = 0; j < (*A).col; j++)

(*C).mat[i][j] = (*A).mat[i][j] - (*B).mat[i][j]; return 1; } else

printf(\这两个矩阵不能相减！\\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).mat[i][j] += (*A).mat[i][k] * (*B).mat[k][j]; } return 1; } else printf(\这两个矩阵不能相乘！\\n\ return 0; }

double MeanMatrix(Matrix *T) //矩阵元素均值 { double mean; double sum = 0; int i, j; for (i = 0; i < (*T).row; i++) for (j = 0; j < (*T).col; j++) sum += (*T).mat[i][j]; mean = sum / ((*T).row*(*T).col); return mean; }

void NMatrix(void) {

#define M 20 #define N 20

int i,j,m,n; float y=1.0; float a[M][2 * M]; float b[N][2 * N]; float t, x; int k, T;

printf(\输入方阵的维数：\\n\请输入方阵，即行和列相等的矩阵。 scanf(\请输入矩阵:\\n\ for(i=0;i

for(i=0;i

}

for(j=0;j<(2*T);j++) a[i][j]=a[i][j]/t; }

for(i=0;i

printf(\逆矩阵为：\\n\ for (i=0;i