[논문 리뷰] Linear Multistep Numerical Methods for Ordinary Differential Equations

Nikesh S. Dattani|ArXiv.org|2008. 10. 28.

Numerical methods for differential equations인용 수 4

한 줄 요약

이 논문은 비선형 미분방정식을 해결하기 위한 선형 다단계 방법에 대한 종합적인 리뷰를 제공하며, 아담스-바슈포스팅(해당 방법은 명시적), 아담스-مول턴(해당 방법은 암시적), 후진 미분 공식(BDF) 방법에 초점을 맞춘다. 이는 우도-바이어슈트라스 정리로부터 우측항을 다항식 근사화함으로써 이러한 방법들을 유도하며, 그들의 순서, 일致성 및 안정성(0-안정성, A-안정성, L-안정성 포함)을 분석한다. 또한 명시적 방법은 A-안정성이 불가능하며, 2보다 큰 순서를 가진 A-안정성 방법은 존재하지 않는다는 것을 입증한다.

ABSTRACT

A review of the most popular Linear Multistep (LM) Methods for solving Ordinary Differential Equations numerically is presented. These methods are first derived from first principles, and are discussed in terms of their order, consistency, and various types of stability. Particular varieties of stability that may not be familiar, are briefly defined first. The methods that are included are the Adams-Bashforth Methods, Adams-Moulton Methods, and Backwards Differentiation Formulas. Advantages and disadvantages of these methods are also described. Not much prior knowledge of numerical methods or ordinary differential equations is required, although knowledge of basic topics from calculus is assumed.

연구 동기 및 목표

• 비스터프 및 스터프 ODE를 해결하기 위해 가장 널리 사용되는 선형 다단계 방법을 체계적으로 유도하고 분석하는 것.
• 선형 다단계 스킴의 맥락에서 순서, 일치성 및 다양한 안정성 개념(0-안정성, A-안정성, L-안정성)의 이론적 기초를 확립하는 것.
• 정확도 및 안정성 특성 측면에서 명시적(아담스-바슈포스팅), 암시적(아담스-몰턴), BDF 방법의 장단점을 비교하는 것.
• 선형 다단계 방법의 한계를 명확히 하는 것, 예를 들어 명시적 방법에서는 A-안정성이 불가능하고, A-안정성 방법의 순서 장벽이 존재한다는 점

제안 방법

• k개 이전 점에서 f(t,y)를 다항식 차수 ≤k−1로 다항식 보간하여 [t_{n-1}, t_n] 간격에서 다항식 근사치를 통합함으로써 아담스-바슈포스팅 방법을 유도한다.
• f_n를 포함하는 암시적 다항식 보간을 사용하여 아담스-몰턴 방법을 유도하며, 이는 k-단계 방법에서 순서 k+1의 더 높은 정확도를 제공한다.
• 절대 안정성 분석을 위해 시험 방정식 y' = λy를 적용하고 안정성 다항식 ζ(ξ) = ρ(ξ) − hλσ(ξ)를 정의하며, 안정성 영역을 |ξ| ≤ 1을 만족하는 z = λh의 집합으로 정의한다.
• 특성 다항식 ρ(ξ)의 근 조건을 통해 0-안정성을 정의하며, |ξ_i| ≤ 1 이며 단순 근이 단위 원 위에 있는 조건을 요구한다.
• A-안정성을 안정성 영역가 전체 왼쪽 반평면(Re(λh) < 0)를 포함하는 성질로 정의하고, L-안정성을 Re(z) → −∞ 일 때 |R(z)| → 0 가 되는 조건으로 정의한다.
• 국소 절단 오차를 정량화하기 위해 오차 상수 C_{p+1}를 사용하며, 순서 p의 방법에 대해 오차가 O(h^{p+1})임을 입증한다.

실험 결과

연구 질문

• RQ1 아담스-바슈포스팅 및 아담스-몰턴 방법은 f(t,y)의 이전 단계에서의 다항식 보간으로 어떻게 유도되는가?
• RQ2 선형 다단계 방법에서 계수 α_j 및 β_j에 어떤 조건이 성립하면 일치성과 순서 p를 확보할 수 있는가?
• RQ3 왜 명시적 선형 다단계 방법은 A-안정성을 달성할 수 없는가?
• RQ4 L-안정성과 스터프 감쇠 사이의 관계는 무엇이며, 왜 낮은 순서임에도 불구하고 뒷단 미분법(BDF)이 트라페지oidal 법칙보다 더 우수한 성능을 보이는가?
• RQ5 A-안정성 선형 다단계 방법의 최대 순서는 얼마이며, 이 한계가 존재하는 이유는 무엇인가?

주요 결과

• 아담스-바슈포스팅 방법은 명시적이고 k-단계 스킴에서 순서 k를 가지며, k개 이전 점에서 f를 다항식으로 보간하고 그 다항식 근사치를 통합함으로써 유도된다.
• 아담스-몰턴 방법은 암시적이고 k-단계 스킴에서 순서 k+1를 가지며, 동일한 수의 단계를 사용할 때 아담스-바슈포스팅보다 더 높은 정확도를 제공한다.
• BDF 방법은 k-단계 스킴에서 순서 k이며, k ≤ 2일 경우 A-안정성이지만 k > 2일 경우 그렇지 않다.
• A-안정성 선형 다단계 방법은 순서가 2를 초과할 수 없으며, 순서 2인 트라페지oidal 법칙(오차 상수 C_3 = 1/12)은 가장 정확한 A-안정성 방법이다.
• L-안정성은 Re(z) → −∞ 일 때 |R(z)| → 0 가 되는 조건이며, 이는 뒷단 미분법(L-안정성)이 트라페지oidal 법칙(A-안정성이지만 L-안정성 아님)보다 스터프 문제에서 더 우수한 성능을 보이는 이유를 설명한다.
• 순서 p의 방법에 대해 국소 절단 오차는 O(h^{p+1})이며, 오차 상수 C_{p+1}가 주로 계수를 결정한다.