[논문 리뷰] Approximation of solutions of DDEs under nonstandard assumptions via Euler scheme

Natalia Czyżewska, Paweł Morkisz외 1명|arXiv (Cornell University)|2021. 06. 07.

Numerical methods for differential equations참고 문헌 22인용 수 2

한 줄 요약

이 논문은 오른쪽 변이 함수의 국소적 허더 연속성, 전역적인 한계선형 성장 및 전역적인 한쪽 립시츠 연속성 조건을 만족하는 지연 미분 방정식(DDEs)에 대한 오일러 스킴의 엄밀한 오차 분석을 제시한다. 주요 기여는 허더 지수에 명시적으로 의존하는 점별 오차에 대한 이론적 상한을 제공하는 것으로, 비선형성이 립시츠 조건을 만족하지 않는 다차원 DDE에 대한 수치 실험을 통해 검증된다.

ABSTRACT

We deal with approximation of solutions of delay differential equations (DDEs) via the classical Euler algorithm. We investigate the pointwise error of the Euler scheme under nonstandard assumptions imposed on the right-hand side function $f$. Namely, we assume that $f$ is globally of at most linear growth, satisfies globally one-side Lipschitz condition but it is only locally H\"older continuous. We provide a detailed error analysis of the Euler algorithm under such nonstandard regularity conditions. Moreover, we report results of numerical experiments.

연구 동기 및 목표

• 오른쪽 변이 함수 f가 전역 립시츠 조건을 만족하지 않지만 더 약한 정규성 조건을 만족할 때 고전적 오일러 스킴의 수렴성을 분석하는 것.
• 전역 선형 성장, 한쪽 립시츠 조건, 국소적 허더 연속성이라는 비표준 가정 하에서 오일러 방법에 대한 엄밀한 점별 오차 상한을 수립하는 것.
• 일般적으로 전역 립시츠 연속성을 가정하는 기존 오차 분석 프레임워크를 물리학적 및 공학적 모델에서 흔히 볼 수 있는 더 현실적이고 덜 정규적인 설정으로 확장하는 것.
• 비선형성이 립시츠 조건을 만족하지 않는 대표적인 DDE에 대한 수치 실험을 통해 이론적 결과를 검증하는 것.
• 유사한 약한 정규성 가정 하에서 고차수 수치 스킴로의 분석 확장을 위한 기초를 마련하는 것.

제안 방법

• 균일한 스텝 크기 h = τ/N를 사용하여 시간 간격 [jτ, (j+1)τ]에 대해 반복적으로 오일러 스킴을 적용하며, 이전 간격의 해를 다음 간격의 초기 조건으로 사용한다.
• 이론적 분석은 허더 연속성과 한쪽 립시츠 조건을 통해 비연속성의 영향을 반영한 그론발라 타입의 점진적 오차 부등식에 기반한다.
• 핵심 단계로는 동일한 약한 정규성 조건 하에서 ODE 해의 거동을 분석하는 보조 보조정리 증명이 포함되며, 특히 허더 지수 β_i에 대해 국소 오차를 h^α 및 h^β_i로 유계화하는 것에 중점을 둔다.
• 국소 절단 오차 추정과 오차의 전역적 전파를 결합하여 오차 상한을 유도하며, 한쪽 립시츠 조건을 활용해 오차 성장을 통제한다.
• 분석 과정에서 함수 f의 허더 지수 β1, ..., βp에 따른 오차 의존성은 명시적으로 추적되며, 이는 새로운 기여이다.
• 비선형성이 립시츠 조건을 만족하지 않는 세 가지 테스트 DDE에 대해 수치 실험을 수행하여 이론적 수렴률과 허더 지수의 오차에 대한 영향을 확인한다.

실험 결과

연구 질문

• RQ1 오른쪽 변이 함수 f가 전역적으로 립시츠 연속성이 아니고 국소적 허더 연속성과 한쪽 립시츠 조건만 만족할 경우, DDE에 대한 오일러 스킴의 수렴 속도는 어떻게 되는가?
• RQ2 이러한 비표준 가정 하에서 오일러 스킴의 점별 오차는 함수 f의 허더 지수에 어떻게 의존하는가?
• RQ3 전역 립시츠 연속성 없이 선형 성장과 한쪽 립시츠 조건에 의존하여 오일러 방법에 대한 엄밀한 오차 상한을 수립할 수 있는가?
• RQ4 수치 실험은 이론적 오차 추정, 특히 허더 지수에 대한 의존성에 대해 어느 정도 검증하는가?
• RQ5 제안된 분석 프레임워크는 유사한 정규성 조건을 갖는 DDE에 대한 고차수 수치 스킴로 확장 가능한가?

주요 결과

• 논문은 오른쪽 변이 함수 f의 허더 지수 β1, ..., βp에 명시적으로 의존하는 오일러 스킴의 점별 오차에 대한 이론적 상한을 수립한다.
• 오차 상한은 O(h^α + ∑h^β_i) 형태이며, α는 시간에 대한 허더 지수이고 β_i는 공간에 대한 허더 지수이다. 이는 정규성이 낮아질수록 수렴 속도가 떨어짐을 보여준다.
• 한쪽 립시츠 상수 H+가 0인 경우 오차 상한은 ∆ + C(1 + ∥ξ∥)(1 + ∆)(b−a)(h^α + ∑h^β_i) 형태이며, 유리한 조건 하에서 h에 대한 선형 수렴을 나타낸다.
• H+ > 0인 경우 오차 상한에는 지수 인자 e^{H+(b−a)}와 H+ 및 문제 매개변수에 의존하는 곱셈 상수 항이 포함되어 있어, 한쪽 립시츠 상수가 영향을 미침을 반영한다.
• 수치 실험은 이론적 수렴 속도를 확인하고, 허더 지수가 감소할수록 오차가 증가함을 보여주며, 정규성에 대한 이론적 의존성이 검증된다.
• 다양한 테스트 케이스에 걸쳐 결과가 안정적이며, 비선형성이 립시츠 조건을 만족하지 않는 다차원 DDE에 대해서도 적용 가능함을 보여주며, 재료 과학 및 상전이 모델에 대한 실제 적용 가능성을 입증한다.