[논문 리뷰] A comparison of a few numerical schemes for the integration of stochastic differential equations in the Stratonovich interpretation

David García-Álvarez|arXiv (Cornell University)|2011. 02. 22.

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

한 줄 요약

이 논문은 스트라토노비치 해석에서 확률미분방정식을 통합하기 위해 헤운, 밀스타인, 도함수 없는 밀스타인의 세 가지 수치적 방법을 평가한다. 강한 수렴 분석과 계산 시간 측정을 통해 헤운 방법이 더 낮은 시간 간격에서 정확도와 효율성 면에서 다른 방법들을 능가함을 발견하였으며, 이는 두 번째 차수 예측-수정 구조와 낮은 계산 오버헤드 덕분이다. 비록 단계당 비용이 약간 더 높지만, 이는 헤운 방법이 더 높은 정확도를 달성함으로써 상쇄된다.

ABSTRACT

Three schemes, whose expressions are not too complex, are selected for the numerical integration of a system of stochastic differential equations in the Stratonovich interpretation: the integration methods of Heun, Milstein, and derivative-free Milstein. The strong (path-wise) convergence is studied for each method by comparing the final points after integrating with $2^n$ and $2^{n-1}$ time steps. We also compare the time that the computer takes to carry out the integration with each scheme. Putting both things together, we conclude that, at least for our system, the Heun method is by far the best performing one.

연구 동기 및 목표

• 스트라토노비치 해석 하에서 세 가지 수치적 통합 방법의 성능을 평가하고 비교하는 것.
• 다양한 시간 간격 해상도에서 결과를 비교하여 각 방법의 강한(경로 기반) 수렴성을 평가하는 것.
• 각 방법의 계산 비용을 분석하여 효율성의 상충 관계를 규명하는 것.
• 가산성 노이즈가 있는 시스템에서 실용적 구현에 가장 효과적인 방법을 특정하는 것.
• 간단한 구조임에도 불구하고 헤운 방법이 정확도 대 비용 비율에서 뛰어나다는 것을 검증하는 것.

제안 방법

• 헤운 방법을 사용하며, 이는 두 번째 차수 예측-수정 방법으로, 예측자로 오일러 단계를 사용하고 보정에는 평균을 취한다.
• 밀스타인 및 도함수 없는 밀스타인 방법을 적용하며, 이들은 확산 계수의 도함수를 기반으로 한 두 번째 차수 보정 항을 포함한다.
• 각 시간 단계에서 평균 0, 분산 h인 정규 분포를 따르는 독립적인 정규 변수 ∼N(0, h)를 사용하여 이산화된 위너 과정을 시뮬레이션한다.
• 2^n 및 2^(n-1) 시간 단계로 통합하여 수렴 차수를 평가하기 위해 수치 결과를 비교한다.
• 각 방법의 계산 시간을 측정하여 효율성의 상충 관계를 평가한다.
• 다수의 실행 동안 평균, 두 번째부터 네 번째 중심 모멘트까지의 통계적 모멘트를 분석하여 안정성과 정확도를 평가한다.

실험 결과

연구 질문

• RQ1 헤운, 밀스타인, 도함수 없는 밀스타인 방법은 스트라토노비치 SDE에서 강한 수렴 면에서 어떻게 비교될 수 있는가?
• RQ2 각 방법의 계산 비용은 얼마이며, 전체 성능에 어떻게 영향을 미치는가?
• RQ3 동일한 계산 노력에서 헤운 방법은 밀스타인 방법들보다 더 높은 정확도를 달성하는가?
• RQ4 왜 밀스타인 방법들은 더 넓은 시간 간격에서 발산 현상이 나타나며, 이는 헤운 방법과 어떻게 비교되는가?
• RQ5 가산성 노이즈가 존재할 경우, 헤운 방법은 고차수 방법들보다 더 뛰어난 복원성과 효율성을 보이는가?

주요 결과

• 헤운 방법은 모든 테스트 시간 간격에서 밀스타인 방법들보다 평균 및 중심 모멘트의 절대값이 훨씬 작게 나타나며, 최대 한 계단 차수까지 작다.
• 밀스타인 방법들은 2^10 단계에서 발산을 보이지만, 헤운 방법은 2^9 단계까지 안정성을 유지함으로써 더 뛰어난 수치적 안정성을 보여준다.
• 단계당 계산 시간이 약간 증가하지만, 헤운 방법은 시간 간격을 두 배로 늘인 밀스타인 방법들과 비교해도 정확도가 유사하게 유지된다.
• 헤운 방법의 두 번째 차수 정확도와 단순성 덕분에 전체적으로 더 효율적이며, 이는 약간의 계산 비용 증가를 상쇄할 정도로 더 높은 정확도를 제공한다.
• 가산성 노이즈의 경우, 밀스타인 방법은 오일러 방법으로 간소화되므로, 일관된 두 번째 차수 수렴성을 보이는 헤운 방법이 더 유리하다.
• 정확도, 안정성, 구현 용이성의 균형을 고려할 때 헤운 방법을 가장 우수한 성능을 보이는 방법으로 권장한다.