[논문 리뷰] Chebyshev collocation for linear, periodic ordinary and delay differential equations: a posteriori estimates

Ed Bueler|arXiv (Cornell University)|2004. 09. 23.

Numerical methods for differential equationsMathematics참고 문헌 28인용 수 18

한 줄 요약

이 논문은 정수 지연을 가진 선형 주기적 상미분방정식과 지연미분방정식(DDE)을 해결하기 위한 체비셰프 콜로케이션 방법을 제시하며, 엄밀한 후행 오차 추정을 제공한다. 힐버트 공간 위의 연산자에 대한 새로운 후행 고유값 편차 정리(증명된 고유값 편차 정리)를 도입함으로써, 단일 주기 연산자의 스펙트럼 반경에 대한 확인 가능한 경계를 도출하고, 정수 지연을 가진 DDE에 대해 수학적으로 엄밀한 안정성 인증을 가능하게 한다.

ABSTRACT

We present a Chebyshev collocation method for linear ODE and DDE problems. We first give a posteriori estimates for the accuracy of the approximate solution of a scalar ODE initial value problem. Examples of the success of the estimate are given. For linear, periodic DDEs with integer delays we define and discuss the monodromy operator U as our main goal is reliable estimation of the stability of such DDEs. We prove a theorem which gives a posteriori estimates for eigenvalues of U, our main result. This result is based on a generalization to operators on Hilbert spaces of the Bauer-Fike theorem for (matrix) eigenvalue perturbation problems. We generalize these results to systems of DDEs. A delayed, damped Mathieu equation example is given. The computation of good bounds on ODE fundamental solutions is an important technical issue; an a posteriori method for such bounds is given. Certain technical issues are also addressed, namely the evaluation of polynomials and the estimation of L^\infty norms of analytic functions. Generalization to the non-integer delays case is also considered.

연구 동기 및 목표

• 정수 지연을 가진 선형 주기적 DDE의 해와 안정성 특성에 대한 신뢰할 수 있고 확인 가능한 수치적 방법을 개발하기 위해.
• 스펙트럼 콜로케이션을 사용하여 DDE의 초기값 문제 해에 대한 후행 오차 추정을 제공하기 위해.
• 단일 주기 연산자의 고유값과 그 행렬 근사치의 고유값 사이의 거리에 대한 후행 경계를 유도하기 위해.
• 해법을 DDE 시스템으로 확장하고, 지연된 감쇠된 매쓰ieu 방정식에 적용하여 그 유용성을 입증하기 위해.
• 스펙트럼 방법에서 다항식 평가, $L^\infty$ 노름 추정, 조건수 제어와 같은 기술적 과제를 다루기 위해.

제안 방법

• Gauss-Lobatto 점에서 체비셰프 콜로케이션을 사용하여 DDE를 이산화하고 고차 정확도로 기본 해를 근사하기 위해.
• 초기값 문제에서 잔차의 $L^\infty$-노름을 경계하기 위해 후행 오차 추정 기법을 적용하기 위해.
• 서브구간에서 스펙트럼 콜로케이션을 사용하여 순차적인 ODE를 풀고, 콜로케이션 점에서 해를 보간함으로써 단일 주기 행렬 근사를 구성하기 위해.
• 힐버트 공간 위의 연산자에 대해 Bauer-Fike 정리를 일반화한 바(정리 11)를 적용하여 단일 주기 연산자의 고유값 편차를 경계하기 위해.
• 고유벡터 행렬의 조건수에 대한 후행 경계를 사용하여 이산 근사에서 고유값 정확도를 제어하기 위해.
• 지연된 매쓰ieu 방정식에 대해 방법을 검증하여, $N$ 증가와 함께 오차 경계가 지수적으로 감소하는 것을 보여주기 위해.

실험 결과

연구 질문

• RQ1 정수 지연을 가진 선형 DDE의 체비셰프 콜로케이션 해에 대해 엄밀한 후행 오차 추정을 어떻게 유도할 수 있는가?
• RQ2 단일 주기 연산자의 고유값과 그 행렬 근사치의 고유값 사이의 관계는 무엇이며, 이를 어떻게 후행적으로 경계할 수 있는가?
• RQ3 스펙트럼 콜로케이션을 사용하여 단일 주기 연산자의 스펙트럼 성질을 수학적으로 엄밀한 오차 경계로 확인할 수 있는가?
• RQ4 단일 주기 행렬 근사의 정확도와 조건수는 콜로케이션 점 수 $N$에 따라 어떻게 변화하는가?
• RQ5 이 방법은 비정수 지연, 중립 DDE 또는 함수미분방정식으로 얼마나 일반화될 수 있는가?

주요 결과

• 수치 실험을 통해 DDE 초기값 문제에 대한 후행 오차 경계가 콜로케이션 점 수 $N$에 따라 지수적으로 감소하는 것으로 확인되었다.
• 매개변수 집합 $(a,b)=(-1.1,1)$에 대해 단일 주기 행렬의 최대 고유값은 $0.9369$이며, 진정한 고유값과의 거리에 대한 후행 경계는 $r=0.0434$이다. 이는 $|\mu_1| + r < 1$이므로 안정성을 증명한다.
• 계산된 고유값 주위의 오차 원판 반경은 $N$에 따라 약 지수적으로 감소하며, $N \gtrapprox 220$일 때 약 $10^{-5}$ 수준의 정밀도에 도달한다.
• 해당 방법을 사용하여 방정식 (1)의 안정성도를 픽셀 단위로 단일 주기 행렬 고유값과 그 확인 가능한 오차 경계를 계산함으로써 성공적으로 제작하였다.
• 고유값 편차 정리(정리 14)는 단일 주기 연산자의 진정한 고유값과 그 행렬 근사치의 가장 가까운 고유값 사이의 거리에 대해 엄밀하고 계산 가능한 경계를 제공한다.
• 지연된 감쇠된 매쓰ieu 방정식에서 유사한 수렴성과 검증 성질을 보여줌으로써, 이 방법은 DDE 시스템으로 일반화 가능하다.