런게-쿠타 법에 대하여
이 도구는 수치 런게-쿠타 법을 탐색하기 위한 대화형 도구입니다. 이 방법은 미분 방정식의 해를 근사하는 데 사용되며 과학 및 공학의 광범위한 문제를 해결하는 데 매우 강력합니다.
다음 미분 방정식을 고려해 보세요:
\( y' = f(x, y) \) 이고 초기값 \( y(x_0) = K \)
구간 \( [a, b] \)에서 해를 근사하고자 합니다. 이 구간을 크기가 \( h \)인 \( n \)개의 더 작은 구간으로 나눕니다. 4차 런게-쿠타 법은 다음과 같이 근사값을 제공합니다:
\( y_0 = K \) (시작 값)이라 할 때
\( i = 0, 1, \dots, n-1 \)에 대하여 \( y_{i+1} = y_i + \frac{1}{6}[k_1 + 2k_2 + 2k_3 + k_4] \)
여기서:
\( k_1 = h \cdot f(x_i, y_i) \)
\( k_2 = h \cdot f(x_i + \frac{h}{2}, y_i + \frac{k_1}{2}) \)
\( k_3 = h \cdot f(x_i + \frac{h}{2}, y_i + \frac{k_2}{2}) \)
\( k_4 = h \cdot f(x_i + h, y_i + k_3) \)
\( i = 0, 1, \dots, n-1 \)에 대하여 \( y_{i+1} = y_i + \frac{1}{6}[k_1 + 2k_2 + 2k_3 + k_4] \)
여기서:
\( k_1 = h \cdot f(x_i, y_i) \)
\( k_2 = h \cdot f(x_i + \frac{h}{2}, y_i + \frac{k_1}{2}) \)
\( k_3 = h \cdot f(x_i + \frac{h}{2}, y_i + \frac{k_2}{2}) \)
\( k_4 = h \cdot f(x_i + h, y_i + k_3) \)
국소 절단 오차는 \( O(h^5) \) 차수이며 원칙적으로 \( h \)가 감소함에 따라 감소합니다.
이 데모에서 사용된 방정식
- \( y' = x^2 \), \( y(0) = 1 \), 정확한 해: \( y = \frac{x^3}{3} + 1 \)
- \( y' = x^4 \), \( y(0) = 1 \), 정확한 해: \( y = \frac{x^5}{5} + 1 \)
- \( y' = e^x \), \( y(0) = 1 \), 정확한 해: \( y = e^x \)
- \( y' = \frac{1}{1+x} \), \( y(0) = 1 \), 정확한 해: \( y = \ln(1+x) + 1 \)
튜토리얼
이 데모의 모든 미분 방정식은 동일한 초기값 \( y(0) = 1 \)과 비교를 위한 정확한 해를 가집니다.
- 첫 번째 미분 방정식 \( y' = x^2 \)를 선택합니다. 시작할 때 \( h = 1.25 \)이고 \( n = 8 \)입니다.
- 그래프에서 정확한 해(파란색 곡선)와 근사치(빨간색 점)를 검토합니다.
- \( n \)을 늘리거나(또는 \( h \)를 직접 줄여서) \( h \)를 줄이고 근사치가 어떻게 향상되는지 관찰합니다.
- 다른 미분 방정식을 선택하고 결과를 분석합니다.
- 그래프 아래의 표에서 정확한 값과 근사값을 비교합니다.
런게-쿠타 계산기
결과
정확한 해와 런게-쿠타 근사치 사이의 비교:
| 단계 (i) | xᵢ | 근사치 yᵢ | 정확한 값 y(xᵢ) | 오차 | k₁ | k₂ | k₃ | k₄ |
|---|
추가 자료
미분 방정식에 대한 추가 참고 자료: