수치해석 2010 년도 봄학기 Homework

1. 수치해석2010년도 봄학기Homework 환경공학과 20080029 김민진

2. I.1 Taylor 급수를 유도하고 절단 오차를 설명하라.

6. 절단오차(truncation error) - 수학적으로 엄밀하게 주어지는 함수f 의 값을, 유한의 사칙 연산의 반복 계산 식fa 로 근사하는 경우, (fa-f)의 오차가 생김. - 예를 들면, 테일러급수를 이용하여 삼각함수를 계산할 때, 무한급수의 계산을 유한 항까지의 계산으로 중단하기 위해서, 절단오차가 발생.

7. I.2 다음의 역행렬을 구하는 문제에 대하여1) 알고리즘을 설명하고 프로그램을 작성 실행하라.2) 책의 프로그램을 Visual Fortran, Visual C 등의 Compiler를 이용하여 실행하고 알고리즘 및 계산 결과를 설명하라. 3) 본 프로그램을 Visual Basic 프로그램으로 변환하라.

8. (1)알고리즘 • 알고리즘:알고리즘: 어떠한 문제를 해결하기 위한 여러 동작들의 유한한 모임이다. • 입력 : 외부에서 제공되는 자료가 0개 이상 존재한다. • 출력 : 적어도 1개 이상의 결과를 내어야 한다. • 명확성 : 각 명령어들은 명확하고 모호하지 않아야 한다. • 유한성 : 알고리듬의 명령어들은 유한번의 수행후에 종료되어야 한다. 이것 은 수행 시간의 현실적인 유한성을 의미한다. • 효과성 : 모든 명령어들은 원칙적으로 종이와 연필만으로 수행될 수 있는 기본적인 것이어야 한다.

9. 프래그램을 작성 실행하라.

11. Visual Fortran

15. Visual C

19. III.1 3원 연립방정식에 대하여 행렬식(Determinant)을 이용하여 역행렬을 구하는 방법과 연립방정식을 푸는 방법을 설명하라.(수학책 참조).

20. 행렬식 • 정방 행렬A의 LU분해가 A=LU인 경우, 그 행렬식|A| 는 |A|=|L| |U| • 로 표현할 수 있다. Doolittle법에 따라LU 분해에서는, L의 대각 요소는 모두 1이고, 삼각 행렬의 행렬식은 그 대각 요소의 적(積)이기 때문에 • 가 된다. 따라서, |A| 은 U의 모든 대각 요소의 적으로서 구해진다. 단,LU 분해에 있어서 Pivot의 선택이 있는 경우는, Pivot을 교환할 때마다 행렬식의 부호를 바꾸어주기 위해, 교환 회수가 m인 경우 가 된다.

21. 따라서, |A| 은 U의 모든 대각 요소의 적으로서 구해진다. 단,LU분해에 있어서 Pivot의 선택이 있는 경우는, Pivot을 교환할 때마다 행렬식의 부호를 바꾸어주기 위해, 교환 회수가 m인 경우 가 된다. • Doolittle법에 따른LU분해의 결과, 구해진U와 Gauss의 소거법에 있어서 전진 소거가 끝난 단계에서의 식 (4.14)의 계수 행렬이 같은 것은 앞서 서술하였다. 이것으로부터, Gauss 소거법에 있어서도 모든 Pivot의 적

22. 에 따라 행렬식을 구하는 것이 가능하다. 행렬식만 구하는 경우는, Gauss의 단순 소거법에 있어서 우변의 계산과 후진 대입의 계산을 생략하고, 전진 소거를 하여, 모든 Pivot의 적을 구하면 된다. 따라서, 행렬식의 계산 순서는 다음과 같다.

24. 역행렬 • 행렬 A의 역행렬 이란, 다음의 관계를 만족시키는 행렬이다. • 여기에서,I는 단위 행렬이다. 역행렬의 계산은, 식(4.4)로부터 구하는 것보다, 소거법 등을 이용하는 편이 계산 횟수가 적게 든다. 여기에서는,LU분해에 따른 계산 방법을 서술한다.

25. 은 미지이므로 X라고 치환하면, 식 (5.5)로부터 AX=I • (5.6)으로 나타내진다. Doolittle법의 LU 분해에 따라 A=LU를 구하여 LY=I의 관계로부터 행렬 를 계산하고, 계속해서 UX=Y 의 관계로부터 행렬X를 계산할 수 있다. 여기에서, 는 LY=I 의 관계로부터 대각 요소가 모두 1인 하삼각행렬이 되고, 의 순서로, 그리고 j를 1부터 n까지라고 하면

26. 으로 구해진다. 또, 는 의 순서로, 그리고 j를 n부터 1까지라고 하면 으로 구해진다. 또한, Y는 하삼각행렬이기 때문에 L의 배열 에 중복되며, X도 U의 영역에 중복되어, 이에 따라 배열 영역을 절약할 수 있다. 나아가, LU분해의 결과를 A에 중복시키면, X의 계산 결과도 A에 중복시키는 것이 가능하다. 단, LU분해에 있어서 Pivot을 교환하는 경우는, 그것에 대응하는 X의 열교환이 필요하게 된다. Pivot 교환이 있는 경우의 계산 순서를 서술하였다.

28. III.2 3원 연립방정식을 이용하여 Gauss 소거법의 알고리즘을 유도하라.

30. III.3 다음의 3원 연립방정식의 해를 Gauss의 단순 소거법 및 Pivot 선택법을 이용하여 구하라.

32. III.4 위의 3원 연립방정식의 해를 교과서에서 제시된 C 및 Fortran 프로그램을 운영하여 구하라.

36. III.5 위의 3원 연립방정식의 해를 Gauss-Jordan 법에 의하여 구하라.

37. (해)행렬식으로 표현한다면 • 제 1단계는, 계수 행렬의 제1열의 대각 요소만 1로 하고 이외는 0으로 한다. 이 때문에, 제1행을 1/2배로 한다.

38. 다음에, 제2행-제1행, 제3행-4×제1행을 계산하여 • 제 2 단계는, 계수 행렬의 제2열의 대각 요소만 1로 하고 이외는 0으로 한다. 이 때문에, 제2행을 -1/3배 한다.

39. 다음에, 제1행-2×제2행, 제3행-(-7)×제2행을 계산하여 • 제 3 단계는, 계수 행렬의 제3열의 대각 요소만 1로 하고 이외는 0으로 한다. 이 때문에, 제3행을 -3/56배 한다.

40. 다음에, 제1행-(13/3)×제3행, • 제2행-(-2/3)×제3행을 계산하여 • 이것으로부터, 의 해를 얻는다.

41. III.6 4원 연립방정식을 이용하여 LU 분해법의 알고리즘을 Doolittle 방법 및 Crout 방법에 대하여 유도하라.

47. III.7 위의 3원 연립방정식의 해를 LU 분해법으로 계산기와 Excel을 이용하여 구하라.

49. III.8 다음의 연립 1차 방정식을 초기값 로서 Jacobi법, Gauss Seidel 법, SOR 법으로 계산기, Excel, 프로그램을 이용하여 구하라.