|
Gauss 소거법을 이용한 선형방정식의 풀이~!! >> \n\n");
printf("\n본래 행렬 값 \n");//본래 행렬 값을 표현한다.
for(j=0;j<4;j++)
{
for(k=0;k<4;k++)
{
printf("%lf ",A[j][k]);
}
printf("\n");
}
printf("\n");
for(i=1;i<4;i++) // 가우스 소거법을 실행한다
{
m[0]=-A[i][0]
|
|
|
가우스 소거법을 수행하는 함수
int BS (E_TYPE *mat, int phase, int size); //후진대입법 적용
E_TYPE *factor(int size) //계수행렬을 만드는 함수
{
int i,j;
E_TYPE *mat;
mat=(E_TYPE*)malloc(sizeof(E_TYPE)*(size)*(size));
printf("
|
|
|
gauss
Ax=B의 형태의 방정식을 가우스소거법으로 풀기
A행렬을 입력하시오: A=[70 1 0;60 -1 1;40 0 -1]
B행렬을 입력하시오: B=[636;518;307]
연산을 시작합니다.
연산을 종료합니다.
ans =
8.59411764705882
34.41176470588233
36.76470588235292
전체 pivoting을 한 횟수
==>
|
|
|
법은 다음과 같다.
를 만족하므로 가 되고
일 때
가 되어야 하므로 이 식에서 각 를 구한 후
각 에 대응하는 Gauss소거법 등을 이용하여 를 구한다.
여기서 가 가 특성다항식 이다.
접근 2. 특성다항식을 이용하여 나의 방정식의 고유치, 고유벡
|
|
|
[j][i])
{
pivot = j;
con = MatrixA[j][i];
}
}
for(j=0;j<column;j++)
{
Sort[j] = MatrixA[pivot][j];
MatrixA[pivot][j] = MatrixA[i][j];
MatrixA[i][j]=Sort[j];
}
con = 0.0;
}
/////////////////// Gauss Elimination - forward elimination /////////////////////////////////
int row_in = row-1; //연산 상에
|
|
메뉴펼치기
