|
0, 0; 0, 1, 0; 0, 0, 1]; %역행렬구하기
D=L^(-1)*E;
fprintf('\n [d]행렬 \n');
disp(D);
fprintf('\n A함수의 역행렬 \n');
Ainv=U^(-1)*D;
disp(Ainv); 1.크래머규칙
2.가우스소거법
3.LU분해법
4.LU분해법으로 역행렬구하기
5.소스코드
|
|
|
LU분해를 이용하자.
역행렬을 LU분해를 이용해서 구하는 방법은 다음과 같다.
각 단계를 계산해보자.
따라서 이다.
물론 자명하게 이다.
11.12 a) system은 다음과 같다.
가우스 소거법을 이용해서 역행렬을 구하면 다음과 같다.
해를 구하면 다음과
|
|
|
소거법
<n의 값으로 4를 입력하고 원소들을 입력함>
Ⅰ-3. 결 과 - 가우스 조르단 법
- n값은 4로 하고 행렬의 원소들은 위의 가우스 소거법에서 사용한 원소들을 그대로 데이터파일에 저장함.
<데이터파일>
<실행결과>
- 역행렬을 구
|
|
|
법에 비해 그 과정이 오래 걸리므로 이것은 알고자 하는 값이 적을 경우만 사용하여야 할 것이다. LU분해법은 때에 따라 전혀 다른 값을 보여주고 있다. 이는 가우스 소거법이나 가우스 조던법에 비하여 좀더 많은 마무리 오차가 발생하는 것
|
|
|
법에 비해 그 과정이 오래 걸리므로 이것은 알고자 하는 값이 적을 경우만 사용하여야 할 것이다. LU분해법은 때에 따라 전혀 다른 값을 보여주고 있다. 이는 가우스 소거법이나 가우스 조던법에 비하여 좀더 많은 마무리 오차가 발생하는 것
|
|
|
3X3 행렬에서 1. 역행렬 구하기
2. 가우스 소거법으로 해 계산하기 3X3 행렬에서 1. 역행렬 구하기
2. 가우스 소거법으로 해 계산하기
|
|
|
Gauss Jordan 실행 후 행렬 출력하기
printf("<-------- 가우스 조단 소거법 -------->");
printf("\n\n");
for(i = 0; i < row_num ; i++)
{
for(j = 0 ; j < col_num; j++)
{
printf("%.2lf",matrix[i][j]);
printf("\t\t");
}
printf("\n");
}
printf("\n\n");
//역행렬 출력하기
printf("<------
|
|
|
or( row = 0 ; row < matrix.matrix_r_size ; row++) // 주대각 값이 0인 경우 역행렬이 없다고 에러 메세지 검출
{
for( col = 0 ; col < matrix.matrix_c_size ; col++)
{
if(matrix.element[row][row]==0) //주대각 값이 0인 경우 강제 종료
{
printf("\nNot Inverse Matrix!!!\n");
exit(0);
}
}
|
|
|
0.000095646
0.809308
-2.96388
-0.3196
0.000034588
0.809279
-2.96394
-0.3196
0.000037162
0.809267
-2.96393
-0.31958
0.000013272
0.809279
-2.96391
-0.31958
0.000013995
0.809283
-2.96391
-0.31959
|
|
|
가우스소거법
1) 선형연립방정식
2) 가우스소거법
①한 행에 k배하여 다른 행에 더함
②두 행을 서로 바꿈
③한 행에 0이 아닌 상수를 곱한다
①+②+③을 통해 A를 ‘삼각형태’ 행렬(에쉴론 행렬)로 변환
Ⅴ. 수학 방정식의 수치 해법
1. 방정식
|
|
메뉴펼치기
