[논문 리뷰] Implementation of Strong Numerical Methods of Orders 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 for Ito SDEs with Non-Commutative Noise Based on the Unified Taylor-Ito and Taylor-Stratonovich Expansions and Multiple Fourier-Legendre Series

Mikhail Kuznetsov, Dmitriy F. Kuznetsov|arXiv (Cornell University)|2020. 01. 01.

Stochastic processes and financial applications참고 문헌 70인용 수 3

한 줄 요약

이 논문은 다차원 비가환성 노이즈를 가진 이토 스토크래틱 미분 방정식(SDEs)을 위한 고차수 강한 수치적 방법을 제시한다. 이는 통합 타일러-이토 및 타일러-스트라토니비치 전개를 사용하고, 반복 스토크래틱 적분의 정확한 근사화를 위해 다중 푸리에-레지온드르 급수를 적용한다. 주요 기여는 평균 제곱수렴을 보장하는 0.5에서 3.0까지의 차수를 갖는 수치적 스킴을 구현한 파이썬 기반의 SDE-MATH 소프트웨어 패키지이다. 이는 비선형 및 선형 SDE 시스템, 특히 태양 활동 모델을 포함한 검증을 통해 이루어졌다.

ABSTRACT

The article is devoted to the implementation of strong numerical methods with convergence orders $0.5,$ $1.0,$ $1.5,$ $2.0,$ $2.5,$ and $3.0$ for Ito stochastic differential equations with multidimensional non-commutative noise based on the unified Taylor--Ito and Taylor-Stratonovich expansions and multiple Fourier-Legendre series. Algorithms for the implementation of these methods are constructed and a package of programs in the Python programming language is presented. An important part of this software package, concerning the mean-square approximation of iterated Ito and Stratonovich stochastic integrals of multiplicities 1 to 6 with respect to components of the multidimensional Wiener process is based on the method of generalized multiple Fourier series. More precisely, we used the multiple Fourier-Legendre series converging in the sense of norm in Hilbert space $L_2([t, T]^k)$ $(k=1,\ldots,6)$ for the mean-square approximation of iterated Ito and Stratonovich stochastic integrals.

연구 동기 및 목표

• 비가환성 노이즈를 가진 다차원 이토 SDEs를 위한 고차수 강한 수치적 스킴(0.5에서 3.0까지의 차수)을 개발한다.
• 다중성 1~6인 반복 이토 및 스트라토니비치 스토크래틱 적분의 근사화에 대한 계산적 과제를 해결한다.
• 힐베르트 공간에서 수렴하는 일반화된 다중 푸리에-레지온드르 급수를 사용하여 효율적인 평균 제곱수 근사화를 구현한다.
• 복잡한 SDEs의 정확하고 효율적인 시뮬레이션을 가능하게 하는 소프트웨어 패키지를 구축한다.
• 비선형 및 선형 SDE 시스템, 특히 태양 활동 모델을 포함한 검증을 통해 수렴성과 정확성을 입증한다.

제안 방법

• 이토 SDEs를 위한 고차수 강한 수치적 스킴를 유도하기 위해 통합 타일러-이토 및 타일러-스트라토니비치 전개를 사용한다.
• 다중 푸리에-레지온드르 급수를 사용하여 다중성 6까지의 반복 스토크래틱 적분의 평균 제곱수 근사화를 수행한다.
• 힐베르트 공간 노름에서 수렴하는 일반화된 다중 푸리에 급수를 적용하여 근사화의 안정성과 정확성을 보장한다.
• 기호 계산을 위해 SymPy를 사용하여 미분 연산자 및 계수의 해석적 유도를 수행한다.
• 고성능 수치 계산을 위해 NumPy를 사용하고, 시뮬레이션 결과의 시각화를 위해 Matplotlib를 사용한다.
• 푸리에-레지온드르 계수와 시뮬레이션 메타데이터를 저장 및 관리하기 위해 SQLite 데이터베이스를 통합한다.

실험 결과

연구 질문

• RQ1 비가환성 노이즈를 가진 이토 SDEs를 위한 고차수 강한 수치적 스킴(최대 차수 3.0까지)을 체계적으로 유도할 수 있는 방법은 무엇인가?
• RQ2 다중성 1에서 6까지의 반복 이토 및 스트라토니비치 스토크래틱 적분을 근사하는 데 가장 정확하고 효율적인 방법은 무엇인가?
• RQ3 다중 푸리에-레지온드르 급수는 어떻게 적용되어 스토크래틱 적분 근사에서 평균 제곱수 수렴을 보장할 수 있는가?
• RQ4 SDE-MATH 소프트웨어 패키지는 비선형 및 선형 SDE 시스템의 시뮬레이션에서 높은 정확도와 성능를 얼마나 달성할 수 있는가?
• RQ5 제안된 방법은 태양 활동 시스템과 같은 실제 모델에 효과적으로 적용될 수 있는가?

주요 결과

• SDE-MATH 소프트웨어 패키지는 비가환성 노이즈를 가진 이토 SDEs를 위한 강한 수치적 스킴(0.5, 1.0, 1.5, 2.0, 2.5, 3.0 차수)을 성공적으로 구현하였다.
• 다중 푸리에-레지온드르 급수를 사용한 반복 스토크래틱 적분의 평균 제곱수 근사화는 수렴 분석을 통해 높은 정확도를 확보하였다.
• 패키지는 비선형 이토 SDE 시스템의 정확한 시뮬레이션을 수행하였으며, 결과는 시각화 및 수치적으로 검증되었다.
• 선형 태양 활동 모델을 위한 시스템에 대해 방법이 성공적으로 적용되어 일관된 수렴성과 신뢰성을 보였다.
• 기호 계산(SymPy), 수치 라이브러리(NumPy), 그리고 데이터 관리용 SQLite의 통합은 효율적이고 확장 가능한 SDE 시뮬레이션을 가능하게 하였다.
• 소프트웨어 패키지는 완전한 그래픽 사용자 인터페이스를 포함하고 있으며, 비선형 및 선형 SDE 시스템을 모두 지원하며, 소스 코드는 공개되어 있다.