수다닷컴

#include stdio.h
#include stdlib.h
#include cmathtypedef struct _MATRIX
{
double **m_data;
int m_size;
}MATRIX;void initMatrix(MATRIX *A, int n);//동적배열 생성함수
void deleteMatrix(MATRIX *A);//동적배열 해제함수
void inputMatrix(MATRIX *A);//행렬 값 입력함수
void printMatrix(MATRIX matrix);//행렬 값 출력함수
double determinant(MATRIX matrix);//행렬식 계산함수
MATRIX transpose(MATRIX matrix); //전치행렬 계산함수
MATRIX minorMatrix(MATRIX matrix, int col, int row);//소행렬식 계산함수
MATRIX cofactorMatrix(MATRIX matrix);//여인자행렬 계산함수
MATRIX adjoint(MATRIX matrix);//수반행렬 계산함수
MATRIX inverseMatrix(MATRIX matrix);//역행렬 계산함수int main(void)
{
MATRIX matrix;
MATRIX inverse;int n;
double det = 0;printf(******************************************************************\n);
printf(수반행렬, 여인수행렬,역행렬,행렬식의값 \n);
printf(******************************************************************\n\n);
printf( 행렬의 크기입력: );
scanf(%d, &n);initMatrix(&matrix, n);
initMatrix(&inverse, n);inputMatrix(&matrix);det = determinant(matrix);inverse = inverseMatrix(matrix);printf(\n입력한 행렬의 행렬식 값 \n);
printf( Det(A) = %.3lf\n\n, det);printf(역행렬 =\n);
printMatrix(inverse);
printf(\n);deleteMatrix(&matrix);
deleteMatrix(&inverse);return 0;
}void initiMatrix (MATRIX *A, int n)
{
int i = 0;A-m_data = (double **)malloc(sizeof(double*) * n);for( i = 0; i n; i++)
{
A-m_data[i] = (double*)malloc(sizeof(double) *n);
}A-m_size = n;
}void deleteMatrix (MATRIX *A)
{
int i = 0;for( i = 0 ; i A-m_size ; i++ )
{
free(A-m_data[i]);
}
free (A-m_data);
}void inputMatrix(MATRIX* A)
{
int i = 0, j = 0;
double input = 0;printf(\n);
printf(행렬의 값을 입력하세요.\n);for( i = 0; i A-m_size; i++)
{
for( j = 0; j A-m_size; j++ )
{
fflush(stdin);
printf(%d X %d 행렬의 값을 입력하세요: , i+1, j+1);
scanf(%lf, &input);
A-m_data[i][j] = input;
}
}
}void printMatrix(MATRIX matrix)
{
int i = 0, j = 0;for( i = 0; i matrix.m_size; i++)
{
printf(\t|\t);
for( j = 0; j matrix.m_size; j++)
{
printf(%.3lf \t, matrix.m_data[i][j]);
}
printf(|\n);
}
}double determinant(MATRIX matrix)
{
int i = 0;
double det = 0;
int sign = 1;if( matrix.m_size == 2)
{
det = matrix.m_data[0][0] * matrix.m_data[1][1] -matrix.m_data[1][0] * matrix.m_data[0][1];return det;
}
for ( i = 0; i matrix.m_size ; i++)
{
MATRIX minor;
initMatrix(&minor, matrix.m_size);minor = minorMatrix(matrix, 0, i);
det = det + sign * matrix.m_data[i][0] * determinant(minor);
sign = sign * -1;
}
return det;
}MATRIX transpose(MATRIX matrix)
{
MATRIX Result;
int i, j;initMatrix( &Result, matrix.m_size);for( i = 0; i matrix.m_size ; i++)
{
for( j = 0; j matrix.m_size ; j++)
{
Result.m_data[i][j] = matrix.m_data[j][i];
}
}
return Result;
}MATRIX minorMatrix(MATRIX matrix, int col, int row)
{
MATRIX Result;
int rowindex = 0;
int colindex = 0;initMatrix( &Result, matrix.m_size - 1);for( i = 0; i matrix.m_size ; i++)
{
for( j = 0; j matrix.m_size ; j++)
{
if( i != row && j != col)
{
Result.m_data[rowindex][colindex] = matrix.m_data[i][j];
colindex++;
}
}
if( i != row && j != col)
{
colindex = 0;
rowindex++;
}
}
return Result;
}
MATRIX cofactorMatrix(MATRIX matrix)
{
MATRIX Result;
int i, j;initMatrix(&Result, matrix.m_size);for( i = 0; i matrix.m_size ; i++)
{
for( j = 0; j matrix.m_size ; j++)
{
Result.m_data[j[i] = determinant(minorMatrix(matrix , i , j));
}
}
return Result;
}rm_MainMATRIX adjoint(MATRIX matrix)
{
MATRIX Result;
MATRIX confactor;
MATRIX transposed;
int i,j;
double ipow = 1;initMatrix(&Result , matrix.m_size);
initMatrix(&confactor, matrix.m_size);
initMatrix(&transposed, matrix.m_size);confactor = confactorMatrix(matrix);
transposed = transpose(confactor);for( i = 0 ; i matrix.m_size; i++)
{
for( j = 0; j matrix.m_size ; j++)
{
ipow = pow(-1, (i+j));
Result.m_data[i][j] = ipow * transposed.m_data[i][j];
}
}
deleteMatrix(&confactor);
deleteMatrix(&transposed);return Result;
}MATRIX inverseMatrix(MATRIX matrix)
{
MATRIX Result;
MATRIX temp;
double det;
int i, j;det = determinant(matrix);initMatrix(&temp, matrix.m_size);
initMatrix(&Result, matrix.m_size);if( matrix.m_size ==2)
{
temp.m_data[0][0] = matrix.m_data[1][1];
temp.m_data[1][1] = matirx.m_data[0][0];
temp.m_data[0][1] = -(matrix.m_data[0][1]);
temp.m_data[1][0] = -(matrix.m_data[1][0]);for( i = 0; i matrix.m_size; i++)
{
for( j = 0; j matrix.m_size; j++)
{
Result.m_data[i][j] = (double)((double)1/det)*temp.m_data[i][j];
}
}
deleteMatrix(&temp);return Result;
}
temp = adjoint(matrix);for ( i = 0; i matrix.m_size ; i++)
{
for( j = 0; j matrix.m_size ; j++)
{
Result.m_data[i][j] = (double)((double)1/det)*temp.m_data[i][j];
}
}
deleteMatrix(&temp);return Result;
}이런 소스가 있는데요.. 수반행렬을 이용해서 역행렬이 구해지는건데요;여기서도 에러가 몇개 있는데 .. .수정을잘 ;;그리고 수반행렬 여인수행렬, 행렬식 값이 중간중간에 구해지는데그걸 표현하는 방법을 모르겠어요 ;;텍스트 파일에 행렬을 입력하면 수반행렬, 여인수행렬,역행렬,행렬식의값 이렇게 네가지가 구해지는 소스로좀 고쳐주세요 ;;제 능력으로 못하겠네요 ;;도와주세요~~

신청하기