수다닷컴

질문 제목 : 역행렬 구하는 프로그램 도와주세요file(input.txt)를 읽어들여서 역행렬을 구하는 프로그램을 만드려고 하는데 안되네요질문 내용 :
input.txt라는 파일을 가지고 행렬의 크기와 성분을 입력해주려구요

input.txt는 아래같은 입력방식을 취해줍니다
2
1 0
0 1아래와 같이 하면 되긴하는데 일일히 행렬의 성분을 입력하는게 뭔가 아니다 싶어서요
#include stdio.h
#include stdlib.h
#include math.h
typedef 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\n);
printf(행렬의 크기 입력 : );
scanf(%d, &n);
initmatrix(&matrix, n);
initmatrix(&inverse, n);
inputmatrix(&matrix);
det = determinant(matrix);
inverse = inversematrix(matrix);
if(det == 0)
{
printf(\n행렬식의 값이 0이므로 역행렬이 존재하지 않습니다.\n);
}
else
{
printf(\n\n 입력한 행렬의 행렬식 값\n);
printf( det(a) = %.3lf\n\n, det);
printf(역행렬 = \n);
printmatrix(inverse);
printf(\n);
}
deletematrix(&matrix);
deletematrix(&inverse);
return 0;
}
void initmatrix(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 * %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 i, j;
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;
}
matrix 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 =factor = cofactormatrix(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] = matrix.m_data[0][0];
temp.m_data[0][1] = -(matrix.m_data[1][0]);
temp.m_data[1][0] = -(matrix.m_data[0][1]);
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;
}

신청하기