[논문 리뷰] Pathwise approximation of SDEs by coupling piecewise abelian rough paths

Guy Flint, Terry Lyons|arXiv (Cornell University)|2015. 05. 06.

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

한 줄 요약

이 논문은 다차원 브라운 운동에 의해 구동되는 확률미분방정식(SDEs)에 대한 새로운 경로 기반 근사 방법을 제안한다. 이 방법은 조각별 아벨리안 거친 경로와 가우시안 근사법을 사용하여 레비 영역 증분을 근사한다. 이 방법은 γ-립시츠 벡터장(γ > 2)에 대해 워샤어스타인 수렴 속도 $ O(h^{1-2/\gamma - \varepsilon}) $ 를 달성하며, 표준 오일러 방법의 $ O(\sqrt{h}) $ 속도를 뛰어넘는다. 레비 영역의 명시적 시뮬레이션을 요구하지 않아도 된다.

ABSTRACT

We present a new pathwise approximation scheme for stochastic differential equations driven by multidimensional Brownian motion which does not require the simulation of Lévy area and has a Wasserstein convergence rate better than the Euler scheme's strong error rate of $O(\sqrt{h})$, where $h$ is the step-size. By using rough path theory we avoid imposing any non-degenerate Hörmander or ellipticity assumptions on the vector fields of the SDE, in contrast to the similar papers of Alfonsi, Davie, and Malliavin et al. The scheme is based on the log-ODE method with the Lévy area increments replaced by Gaussian approximations with the same covariance structure. The Wasserstein coupling is achieved by making small changes to the argument of Davie, the latter being an extension of the Komlós-Major-Tusnády Theorem. We prove that the convergence of the scheme in the Wasserstein metric is of the order $O(h^{1-2/γ-\varepsilon})$ when the vector fields are $γ$-Lipschitz in the sense of Stein.

연구 동기 및 목표

• d ≥ 2 인 경우 계산적으로 비현실적인 레비 영역 증분의 시뮬레이션을 피하는 경로 기반 SDE 근사 방법을 개발하는 것.
• 특히 고차원 또는 비퇴화된 설정에서 표준 오일러 방법의 $ O(\sqrt{h}) $ 수렴 속도를 뛰어넘는 워샤어스타인 수렴 속도를 달성하는 것.
• 벡터장에 대한 비퇴화 허르만더 또는 타원성 조건이 필요 없도록 제거하여 더 넓은 범위의 SDE에 적용 가능하게 하는 것.
• 거친 경로 이론과 로그-ODE 방법을 활용하여 진짜 SDE 해와 근사값 사이의 쌍용을 구성하고 워샤어스타인 오차를 통제하는 것.
• Komlós-Major-Tusnády 정리의 프레임워크를 SDE 맥락에서 다차원 브라운 운동에 의해 구동되는 확률과정으로 확장하여 워샤어스타인 쌍용을 활용하는 것.

제안 방법

• 레비 영역 증분을 근사하기 위해 독립적인 가우시안 난수 $ W^{(j)} \sim N(0, hI_d) $, $ z^{(j)} \sim N(0, 12^{-1}hI_d) $, $ \lambda^{(j)} \sim N(0, 12^{-1}h^2I_{d(d-1)/2}) $ 를 사용한다.
• 근사값 $ x_j^h $ 는 지수 사상에 의해 반복적으로 갱신된다: $ x_{j+1}^h = \exp\left(hV_0 + \sum_k W_k^{(j)} V_k + \sum_{k<l} B_{kl}^{(j)} [V_k, V_l] \right)(x_j^h) $, 여기서 $ B_{kl}^{(j)} $ 는 반복 적분의 가우시안 근사이다.
• 이 방법은 로그-ODE 접근법을 활용하여 비예측적이고 경로 기반의 근사를 다룰 수 있도록 SDE 해를 거친 경로 프레임워크에 통합한다.
• 드라이빙 브라운 운동과 근사 노이즈를 동일한 확률공간 위에 정렬함으로써 워샤어스타인 쌍용을 구성하여 경로 기반 오차를 공동으로 통제한다.
• 수렴 분석은 Komlós-Major-Tusnády 정리와 그 거친 경로 설정으로의 확장에 기반하며, Baker–Campbell–Hausdorff 공식의 리 대수 전개를 사용한다.
• 핵심 기술적 혁신은 리 대수 수준의 변화에 의해 BCH 전개의 차이를 유한하게 제어하는 것으로, Lemma 11.3을 통해 전역 군 수준 오차를 통제한다.

실험 결과

연구 질문

• RQ1 레비 영역의 시뮬레이션 없이도 $ O(\sqrt{h}) $ 수렴 속도를 뛰어넘는 경로 기반 SDE 근사 방법을 구축할 수 있는가?
• RQ2 비퇴화되거나 타원성이 없는 벡터장을 갖는 SDE에 대해, 워샤어스타인 거리에서 오일러 방법의 $ O(\sqrt{h}) $ 수렴 속도를 초월할 수 있는가?
• RQ3 단지 가우시안 난수를 사용하고 레비 영역의 명시적 계산 없이도 워샤어스타인 거리에서 고차수 수렴을 달성할 수 있는가?
• RQ4 Komlós-Major-Tusnády 정리는 SDE 맥락에서 다차원 브라운 운동에 의해 구동되는 확률과정으로 어떻게 일반화될 수 있는가?
• RQ5 거친 경로 기반 근사 방법을 사용할 때, γ-립시츠 벡터장을 갖는 SDE에 대해 워샤어스타인 거리에서 최적의 수렴 속도는 무엇인가?

주요 결과

• 제안된 방법은 $ \gamma > 2 $ 인 벡터장의 허더-립시츠 정규성에 대해 임의의 $ \varepsilon > 0 $ 에 대해 워샤어스타인 수렴 속도 $ O(h^{1-2/\gamma - \varepsilon}) $ 를 달성한다.
• 레비 영역의 시뮬레이션을 완전히 피하고, 브라운 운동의 반복 적분과 일치하는 공분산 구조를 갖는 가우시안 근사에 의존한다.
• 표준 오일러 방법의 $ O(\sqrt{h}) $ 오차 속도를 뛰어넘으며, 특히 $ \gamma > 4 $ 인 경우 수렴 속도가 $ O(h^{1-\varepsilon}) $ 로 향상된다.
• 최소한의 가정 하에 유효하다: $ V_0 \in \text{Lip}^1(\mathbb{R}^q) $, $ V \in \text{Lip}^\gamma(\mathbb{R}^q) $ 이며 $ \gamma > 2 $ 이고, 비퇴화성 또는 허르만더 유형 조건이 필요하지 않다.
• 워샤어스타인 쌍용은 동일한 확률공간 위에 명시적으로 구성되어 있어 진짜 해와 근사값 간의 공동 분포 제어를 보장한다.
• 이론적 기초는 리 대수에서의 Baker–Campbell–Hausdorff 공식의 정밀한 분석에 기반하며, 리 대수 수준의 변화에 의해 군 원소 간의 차이를 제어하는 새로운 경계(논문 내 Lemma 11.3)를 제공한다.