[논문 리뷰] Modify the Improved Euler scheme to integrate stochastic differential equations

A. J. Roberts|arXiv (Cornell University)|2012. 10. 02.

Stochastic processes and financial applications참고 문헌 9인용 수 26

한 줄 요약

이 논문은 확률미분방정식(SDE)을 통합하기 위한 수정된 룬게쿠타 방법을 제안한다. 이 방법은 결정론적 극한에서 잘 알려진 개선된 오일러(헤운) 방법으로 수렴한다. 랜덤 부호 $ S_k = \pm 1 $을 소음 증분에 포함시켜 이토 및 스트라토노비치 SDE 모두에 대해 일阶 강한 수렴을 달성하며, 이론적 증명과 수치 예제를 통해 전역 오차가 $ \mathcal{O}(h) $임을 확인하였다.

ABSTRACT

A practical and new Runge--Kutta numerical scheme for stochastic differential equations is explored. Numerical examples demonstrate the strong convergence of the method. The first order strong convergence is then proved using Ito integrals for both Ito and Stratonovich interpretations. As a straightforward modification of the deterministic Improved Euler/Heun method, the method is a good entry level scheme for stochastic differential equations, especially in conjunction with Higham's introduction [SIAM Review, 43:525--546, 2001].

연구 동기 및 목표

• 확률미분방정식(SDE)을 통합하기 위한 실용적이고 직관적인 수치적 방법을 개발하여 기존의 결정론적 방법과의 명확한 연결을 유지하고자 한다.
• 이 방법이 이토 및 스트라토노비치 SDE의 해석 모두에 대해 일阶 강한 수렴을 달성하도록 보장하고자 한다.
• 스토케스틱 역학의 교육용 입문 수준 방법을 제공하고자 하며, 특히 하이엄의 SDE 시뮬레이션에 대한 소개적 리뷰와 함께 사용하기에 적합하다.
• 다중 소음원을 다룰 수 있도록 방법을 일반화하고자 하지만, 이는 아직 해결되지 않은 과제이다.

제안 방법

• 이 방법은 결정론적 개선된 오일러 방법을 수정하여, 소음 항에 확률적 부호 $ S_k = \pm 1 $을 도입함으로써 확률성을 반영한다. 각각의 확률은 $ 1/2 $이다.
• 각 시간 단계 $ h $ 에서 두 개의 중간 단계를 계산한다: 현재 상태를 사용한 $ \vec{K}_1 $ 과 예측된 다음 상태를 사용한 $ \vec{K}_2 $ 로, 둘 다 $ \Delta W_k = \sqrt{h} Z_k $ 를 통해 소음을 포함한다.
• 최종 갱신은 $ \vec{X}_{k+1} = \vec{X}_k + \frac{1}{2}(\vec{K}_1 + \vec{K}_2) $ 로 주어지며, 결정론적 헤운 방법의 구조를 유지한다.
• 스트라토노비치 SDE의 경우, $ S_k = 0 $ 으로 설정하면 일阶 강한 수렴이 달성되며, 이는 방법이 스트라토노비치 해석과 일치함을 의미한다.
• 이 방법은 이토 미적분과 테일러 전개를 사용하여 局부 및 전역 오차를 분석하며, 부드러운 조건 하에 일반적인 SDE에 대해 $ \mathcal{O}(h) $ 수렴이 엄밀히 증명된다.
• 이 방법은 해석적 해가 알려진 테스트 SDE에 대해 수치적으로 검증되었으며, 기대한 $ \mathcal{O}(h) $ 경로 오차가 관찰되었다.

실험 결과

연구 질문

• RQ1 소음이 없는 경우 정확히 결정론적 개선된 오일러 방법으로 수렴하는, SDE에 대한 룬게쿠타 유사 방법을 구성할 수 있는가?
• RQ2 이 수정된 방법이 이토 SDE에 대해 강한 수렴의 차수는 무엇이며, 엄밀한 증명이 가능한가?
• RQ3 이 방법은 스트라토노비치 SDE에 대해 어떻게 작동하는가? 예를 들어 $ S_k = 0 $ 으로 설정하는 수정이 일관된 수렴을 가능하게 하는가?
• RQ4 어떤 조건에서 이 방법이 일阶 초과 수렴, 예를 들어 $ \mathcal{O}(h^2) $ 를 달성하는가?
• RQ5 이 방법을 다중 독립 브라운 운동으로 자연스럽게 일반화할 수 있는가?

주요 결과

• 제안된 방법은 이토 SDE에 대해 일阶 강한 수렴을 달성하며, 전역 오차가 $ \mathcal{O}(h) $ 임을 이토 적분과 테일러 전개 기법을 통해 증명하였다.
• 스트라토노비치 SDE의 경우, $ S_k = 0 $ 으로 설정하면 동일한 $ \mathcal{O}(h) $ 전역 오차를 얻으며, 이는 스트라토노비치 해석과의 일致성을 확인한다.
• 수치 예제를 통해 $ \mathcal{O}(h) $ 오차 행동이 확인되었으며, 한 예에서는 유리한 드리프트 및 변동성 구조로 인해 $ \mathcal{O}(h^2) $ 수렴이 관찰되었다.
• 소음이 없는 경우(즉, $ \vec{b} = \vec{0} $)에 표준 개선된 오일러/헤운 방법으로 수렴하므로 직관적 일致성이 확보된다.
• 선형 드리프트와 변동성을 갖는 경우, $ \dot{\vec{b}} = \beta \vec{b} $ 를 만족할 때 국소 오차는 $ \mathcal{O}(h^3, h^5) $ 로 향상되며, 이로 인해 전역 오차는 $ \mathcal{O}(h^2, h^4) $ 가 된다. 이는 수치적으로 확인되었다.
• 이 방법은 교육용 및 다중 스케일 시뮬레이션에 적합하며, 특히 소음 강도가 시간 단계에 비례할 경우 유용하지만, 다중 소음원으로의 확장은 아직 해결되지 않은 문제이다.