[논문 리뷰] Higher Order Methods for Differential Inclusions

Sanja Gonzalez Zivanovic, Pieter Collins|arXiv (Cornell University)|2012. 06. 28.

Numerical methods for differential equations인용 수 5

한 줄 요약

이 논문은 Fliess 유사 전개와 로그 노름 기반 오차 추정을 사용하여 미분 포함의 달성 가능한 집합에 대한 엄밀한 초과 근사값을 계산하기 위한 고차수 수치 방법을 제안한다. 이 방법은 1차 수준의 방법보다 향상된 국소 오차 한계를 달성하여 정확한 장시간 통합을 가능하게 하며, 불확실하거나 노이즈가 있는 입력을 가진 시스템에 대해 구성요소별 및 균일한 오차 제어를 유지한다.

ABSTRACT

We present a numerical method for rigorous over-approximation of a reachable set of differential inclusions. The method gives high-order error bounds for single step approximations and a uniform bound on the error over the finite time interval. We provide formulas for the local error based on Lipschitz constants and bounds on higher-order derivatives. The method is based on a Fliess-like expansion, and extends previous results by providing error estimates which are valid for all possible inputs.

연구 동기 및 목표

• 미분 포함에서 달성 가능한 집합에 대한 엄밀한 초과 근사값을 계산하기 위한 고차수 수치 방법을 개발한다.
• 제어 및 노이즈가 있는 시스템에서 가능한 모든 입력에 대해 유효한 균일한 고차수 오차 한계를 제공함으로써 1차 수준 방법을 확장한다.
• 우수한 볼륨과 볼륨이 유한한 우변을 가진 미분 포함의 장시간 통합에서 정확성과 효율성을 향상시킨다.
• 불확실성, 비연속적인 역학 또는 모델 단순화 오차가 있는 시스템의 신뢰성 있는 검증을 가능하게 한다.
• 달성 가능한 집합 계산에서 국소 오차 정확도와 계산 효율성 사이의 상호 교환 관계를 조사한다.

제안 방법

• 시간 단계 동안 해의 진화를 모델링하기 위해 Fliess 유사 거듭제곱 급수 전개를 사용한다.
• 리프시츠 상수보다 더 날카운 오차 한계를 도출하기 위해 로그 노름을 적용하여 안정적이고 불안정한 역학을 구분할 수 있다.
• 고차수 도함수를 제한하고 균일한 입력 제한을 사용하여 국소 오차 추정치를 2차 및 3차로 계산한다.
• 각 시간 간격 동안 근사 해와 균일한 오차 한계를 조합하여 달성 가능한 집합의 초과 근사값을 구성한다.
• 그리드 기반 방법을 피하기 위해 적응형 시간 간격을 사용하는 단일 단계 통합 기법을 채택하여 효율성을 향상시킨다.
• 함수 해석학과 자동 오차 전파를 사용하여 상태 공간 내 달성 가능한 집합의 봉투를 계산한다.

실험 결과

연구 질문

• RQ1 다중값 역학을 가진 미분 포함에 대해 고차수 오차 추정이 효과적으로 적용될 수 있는가?
• RQ2 오차 한계의 날카움 정도와 적용 가능성 측면에서 로그 노름 기반 오차 한계는 리프시츠 기반 오차 한계보다 어떻게 비교되는가?
• RQ3 고차수 달성 가능한 집합 계산에서 국소 오차 정확도와 계산 비용 사이의 상호 교환 관계는 무엇인가?
• RQ4 이 방법은 노이즈가 있는 입력이나 모델 단순화 오차가 있는 시스템에 대해 신뢰할 수 있는 초과 근사값을 제공할 수 있는가?
• RQ5 불확실한 매개변수를 가진 로슬러 방정식과 같은 복잡한 시스템에서 이 방법의 성능은 어떠한가?

주요 결과

• 편향된 조화 진동자에 대해 3차 오차 추정치는 h=0.01일 때 달성 가능한 집합 반경 0.00100759를 도출하였으며, 이는 [17]에서 2-노름 로그 노름을 사용하여 도출한 0.00100502와 유사하다.
• h=0.25일 때 이 방법의 2차 추정치는 반경 0.0420586을 제공하였고, [17]의 0.0284025보다 약간 열등하지만, h가 더 작아질수록 결과가 크게 향상된다.
• a=5.7이고 노이즈 v_i ∈ [-10⁻⁴, 10⁻⁴]인 로슬러 방정식에 대해, 이 방법은 h=0.005, T=11.1일 때 R(X₀,T) ≈ {0}×[-0.15572, 0.15391]×[-3.75926, -3.41772]×[0.03139, 0.03398]를 계산한다.
• 해석적 오차 한계는 e = 8.586×10⁻⁸이며, 이는 봉투의 높은 정밀도를 보여준다.
• 수동 계산된 한계가 아닌 자동 오차 추정을 사용했음에도 불구하고, 이 방법의 결과는 [17]의 결과와 품질 면에서 유사하다.
• 이 방법은 불확실한 입력을 가진 복잡하고 혼돈적인 시스템에 대해 고차수 오차 추정이 실현 가능하고 효과적임을 보여주며, 엄밀한 초과 근사값을 유지한다.