수다닷컴

질문 제목 : 가우스 조던 역행렬 구하기. 마지막 a*a-1 어떻게 하죠~?일단 코딩을 어떻게든 다 짜놨습니다. 근데.. 마지막 a*a-1 을 하려니.. 뭔가 처음부터 다시 해야하나 싶기도하고..질문 내용 : #include stdio.h
#include math.h

void main()
{
float a[11][21]={0,};/////////////범용으로 사용하기위해 충분히 큰 배열 배정
/*float templine[102], tempmax;/////partial pivoting
int pivotmaxline[101];////////각 대각성분별선정이벤트별로 최대성분이 들어간 라인 추출
*/

float aik;////////gauss중 upper trigonal matrix의 0되기전 a[i][k]값
int i,j,n,k/*,l*/;/////////////////////좌표사용을 위한 계수들

printf([n*n]행렬 크기(n)를 입력하세요. \n);
scanf(%d,&n);
printf(\n);

printf(각 행렬 성분을 입력하시오.\n);
for(i=1;in+1;i++)
{
for(j=1;jn+1;j++)
{
printf(a[%d][%d] = , i, j);
scanf(%f,&a[i][j]);
}
}////////////////행렬 성분들 입력
printf(\n);for(i=1;in+1;i++)
{
a[i][i+n]=1;
}///////////////역행렬 자리의 단위행렬화
/*
for(k=1;kn;k++)
{
tempmax=0;///////////tempmax는 각 k단계 주대각선상에 놓일 최대값을 찾는다.
for(i=k;in+1;i++)
{
if(fabs(a[i][k])fabs(tempmax))
{
pivotmaxline[k]=i;
tempmax=a[i][k];
}
}
l=pivotmaxline[k];
for(j=1;j2*n+1;j++)
{

templine[j]=a[l][j];
a[l][j]=a[k][j];
a[k][j]=templine[j];
}
}/////////////////partial pivoting
*/

printf([a | i] 행렬.\n);
for(i=1;in+1;i++)
{
for(j=1;j2*n+1;j++)
{
printf(%f ,a[i][j]);
}
printf(\n);
}//////////////입력 출력
printf(\n);
for(k=1;kn;k++)
{
for(i=k+1;in+1;i++)
{
aik=a[i][k];
for(j=1;j2*n+1;j++)
{
a[i][j]=a[i][j]*a[k][k]-a[k][j]*aik;
}

for(j=k;j2*n+1;j++)
{
a[i][j]=a[i][j]/a[k][k];
}////////가우스 소거중에 계수가 매우 커지는 걸 막는다.
}
}/////////////가우스 소거법upper trigonal matrix
printf(upper triangle 행렬. \n);
for(i=1;in+1;i++)
{
for(j=1;j2*n+1;j++)
{
printf(%f ,a[i][j]);
}
printf(\n);
}//////////////결과 출력
printf(\n);for(k=1;kn+1;k++)
{
a[n][n+k]=a[n][n+k]/a[n][n];
}
a[n][n]=1;//////////////마지막행은 알고리즘상 미리 계산한다.
for(k=n;k1;k--)
{
for(i=k-1;i0;i--)
{
aik=a[i][k];
for(j=2*n+1;j0;j--)
{
a[i][j]=a[i][j]*a[k][k]-a[k][j]*aik;
}

for(j=2*n+1;j0;j--)
{
a[i][j]=a[i][j]/a[i][i];
}////////가우스 소거중에 대각성분값을 1로
}
}/////////////가우스 소거법lower trigonal matrixprintf(최종 [n*2n] 역행렬. \n);
for(i=1;in+1;i++)
{
for(j=1;j2*n+1;j++)
{
printf(%f ,a[i][j]);
}
printf(\n);
}//////////////결과 출력
printf(\n);}
--------------------------------------------------------------------------------------------

이렇게 해서. 가우스 조던 역행렬 구하는 것의求?것을 했습니다.

근데 마지막에 a*a-1 = i 임을 출력을 해야하는데...

처음 a 행렬의 값을 다시 되돌리는 것과.. a-1을 다시 곱해야하는 것..을 어떻게 해야하는지 막막하네요..

도와주세요.. ㅜㅜ

신청하기