[논문 리뷰] AN OVERVIEW OF NUMERICAL AND ANALYTICAL METHODS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS

Byakatonda Denis|arXiv (Cornell University)|2020. 01. 01.

Numerical methods for differential equations참고 문헌 18인용 수 3

한 줄 요약

이 논문은 상미분방정식(OEDs)을 해결하기 위한 분석적 및 수치적 방법에 대한 종합적인 개요를 제공하며, MATLAB에서 오일러 방법과 4차 룬게-쿠타 방법을 비교한다. 이는 룬게-쿠타 방법이 오일러 방법보다 훨씬 높은 정확도를 보임을 보여주며, 특히 작은 스텝 크기에서 두드러진다. 또한 향후 ODE와 강성 시스템에 대한 연구를 위해 고급 수치 해법 및 오픈소스 도구를 권장한다.

ABSTRACT

Differential Equations are among the most important Mathematical tools used in creating models in the science, engineering, economics, mathematics, physics, aeronautics, astronomy, dynamics, biology, chemistry, medicine, environmental sciences, social sciences, banking and many other areas [7]. A differential equation that has only one independent variable is called an Ordinary Differential Equation (ODE), and all derivatives in it are taken with respect to that variable. Most often, the variable is time, t; although, I will use x in this paper as the independent variable. The differential equation where the unknown function depends on two or more variables is referred to as Partial Differential Equations (PDE). Ordinary differential equations can be solved by a variety of methods, analytical and numerical. Although there are many analytic methods for finding the solution of differential equations, there exist quite a number of differential equations that cannot be solved analytically [8]. This means that the solution cannot be expressed as the sum of a finite number of elementary functions (polynomials, exponentials, trigonometric, and hyperbolic functions). For simple differential equations, it is possible to find closed form solutions [9]. But many differential equations arising in applications are so complicated that it is sometimes impractical to have solution formulas; or at least if a solution formula is available, it may involve integrals that can be calculated only by using a numerical quadrature formula. In either case, numerical methods provide a powerful alternative tool for solving the differential equations under the prescribed initial condition or conditions [9]. In this paper, I present the basic and commonly used numerical and analytical methods of solving ordinary differential equations.

연구 동기 및 목표

• 상미분방정식을 해결하기 위한 분석적 및 수치적 방법을 조사하고 비교하기.
• 오일러 방법과 4차 룬게-쿠타 방법의 정확도 및 수렴성을 평가하기.
• 스텝 크기가 수치적 해의 정확도에 미치는 영향을 평가하기.
• 특히 강성 방정식과 시스템에 대해 ODE를 해결하기 위한 계산 도구 및 향후 연구 방향을 제안하기.

제안 방법

• 연구는 상수 계수와 변수 계수를 가진 일계 및 이계 선형 ODE에 대해 분석 기법을 적용한다.
• 수치적 해는 다양한 스텝 크기(h = 0.1, 0.05, 0.01)로 오일러 방법과 4차 룬게-쿠타 방법을 사용하여 계산한다.
• MATLAB을 사용하여 두 수치 알고리즘을 구현하고 비교 분석을 수행한다.
• 정확한 분석적 해와의 비교를 통해 오차 분석을 수행한다.
• 다양한 스텝 크기에서의 수렴성과 정확도를 시각화하기 위해 그래픽적 비교를 생성한다.
• 이산 점에서 절대 오차 지표를 사용하여 수치 방법의 성능을 평가한다.

실험 결과

연구 질문

• RQ1 같은 ODE를 동일한 스텝 크기로 풀었을 때 오일러 방법과 4차 룬게-쿠타 방법은 정확도에서 어떻게 비교되는가?
• RQ2 스텝 크기를 줄임으로써 ODE의 수치적 해의 정확도는 어느 정도 향상되는가?
• RQ3 스텝 크기가 0에 가까워질수록 수치적 해는 정확한 분석적 해로 수렴하는가?
• RQ4 초기값 문제에 대해 오일러 또는 룬게-쿠타 방법 중 어느 것이 더 빠른 수렴성과 낮은 오차를 제공하는가?
• RQ5 분석적 해가 존재하지 않을 경우 ODE를 해결하기 위해 가장 효과적인 계산 도구와 전략은 무엇인가?

주요 결과

• 같은 스텝 크기에서 4차 룬게-쿠타 방법은 오일러 방법보다 훨씬 작은 오차를 생성하며, h = 0.01일 때 오차가 최대 1.69 × 10⁻⁸까지 낮아진다.
• h = 0.01일 때, 룬게-쿠타 방법의 최대 오차는 1.69 × 10⁻⁸이었고, 오일러 방법은 x = 0.1에서 최대 오차 0.010196522를 기록했다.
• 더 작은 스텝 크기일수록 수치 오차가 지속적으로 감소하였으며, 두 방법 모두 h → 0에 가까워질수록 정확한 해로 수렴하는 경향을 보였다.
• h < 0.05일 경우, 룬게-쿠타 방법의 해 곡선은 정확한 해와 거의 구분되지 않아, 뛰어난 안정성과 정확도를 보였다.
• 오일러 방법은 작은 스텝 크기일지라도 큰 오차를 보이며, 실용적 응용에서의 한계를 드러냈다.
• 이 연구는 룬게-쿠타 방법이 ODE의 초기값 문제를 해결하는 데 오일러 방법보다 더 효과적이고 효율적임을 확인한다.