[논문 리뷰] $V$-Integrability, Asymptotic Stability And Comparison Theorem of Explicit Numerical Schemes for SDEs

Łukasz Szpruch, Xīlíng Zhāng|arXiv (Cornell University)|2013. 10. 02.

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

한 줄 요약

이 논문은 비글로벌 리프시츠 계수를 가진 확률미분방정식(SDEs)에 대한 수정된 명시적 수치적 스킴의 가족을 제안하며, 케이스미스키의 리아푸노프 함수 기법을 적용하여 이산 시간에서 $V$-적분 가능성, 점근적 안정성, 비교 정리 등을 확립한다. 주요 기여는 표준 오일러 스킴이 실패할 수 있는 상황에서도 연속 SDE가 지닌 정성적 성질(예: 기대값에서의 유계성, 균형점으로의 거의 확실한 수렴)을 스킴이 유지할 수 있도록 하는 것이다.

ABSTRACT

Khasminski's \cite{chas1980stochastic} showed that many of the asymptotic stability and the integrability properties of the solutions to the Stochastic Differential Equations (SDEs) can be obtained using Lyapunov functions techniques. These properties are rarely inherited by standard numerical integrators. In this article we introduce a family of explicit numerical approximations for the SDEs and derive conditions that allow to use Khasminski's techniques in the context of numerical approximations, particularly in the case where SDEs have non globally Lipschitz coefficients. Consequently, we show that it is possible to construct a numerical scheme, that is bounded in expectation with respect to a Lyapunov function, and/or inherit the asymptotic stability property from the SDEs. Finally we show that using suitable schemes it is possible to recover comparison theorem for scalar SDEs.

연구 동기 및 목표

• 기존에 연속 시간 SDEs에만 적용 가능했던 케이스미스키의 리아푸노프 함수 기법을 이산 시간 수치 스킴으로 확장한다.
• 계수들이 비글로벌 리프시츠 조건을 만족할 경우 표준 명시적 오일러 스킴이 적분 가능성과 안정성을 유지하지 못하는 문제를 다룬다.
• 기본 SDE에서 유도된 $V$-적분 가능성과 점근적 안정성을 유지하는 명시적 수치 스킴을 구축한다. 이는 약한 모멘트 조건 하에서도 가능하다.
• 적절한 스킴 수정을 통해 스칼라 SDE의 비교 정리를 수치적 설정에서 복원한다.
• 일반 리아푸노프 함수 $V(x) = \sum_{i=1}^d c_i x_i^{p_i}$를 사용하여 $L^p$-노름을 초월한 수치 해의 정성적 성질 분석 프레임워크를 제공한다.

제안 방법

• 성장 제어를 위해 억제된 드리프트와 확산 계수를 가진 수정된 명시적 스킴 $\bar{X}_{k+1} = \bar{X}_k + b^h(t_k, \bar{X}_k)h + \sigma^h(t_k, \bar{X}_k)\Delta W_{k+1}$를 도입한다.
• 증분 항에 대해 투영 또는 잘라내기 메커니즘을 적용하여, $p_i$가 서로 다를 수 있는 경우에도 리아푸노프 함수 $V(x) = \sum_{i=1}^d c_i x_i^{p_i}$ 하에서 기대값에서의 유계성을 보장한다.
• 연속 시간 리아푸노프 분석을 모방하기 위해 이산 시간 버전의 무한소 생성자 $\mathcal{L}V$의 유사체를 유도한다. 이를 통해 안정성 및 적분 가능성 증명이 가능해진다.
• 모멘트 추정과 코시-슈바르츠 및 체비셰프 부등식을 사용하여, 특히 $L^{2+\epsilon}$ 모멘트 조건 하에서 스킴의 약한 및 강한 오차를 유계화한다.
• 드라이프 조건에 따라 지수 감쇠 또는 유계성을 보장하기 위해 $\mathds{E}V(\bar{X}_k)$의 기대값에 그론월라의 유형의 추론을 적용한다.
• 오차 추정을 통해 수렴성을 검증하며, $\mathds{E}|\bar{X}_{t,x} - \tilde{X}_{t,x}|^2 \leq K(1 + |x|^{2+2/r})h^2$ 및 $\mathds{E}|\bar{X}_{t,x} - \mathds{E}\tilde{X}_{t,x}| \leq K(1 + |x|^{1+3/(2r)})h^{3/2}$ 를 도출한다.

실험 결과

연구 질문

• RQ1 비글로벌 리프시츠 계수를 가진 SDEs에 대한 명시적 수치 스킴에서 리아푸노프 함수 기법을 어떻게 조정하여 $V$-적분 가능성을 확보할 수 있는가?
• RQ2 표준 오일러 스킴이 실패할 경우, 어떤 조건에서 명시적 스킴이 기저 SDE로부터 점근적 안정성을 상속할 수 있는가?
• RQ3 스칼라 SDE의 비교 정리를 수정된 스킴을 통해 수치 근사 설정에서 복원할 수 있는가?
• RQ4 일반 리아푸노프 함수 $V(x) = \sum c_i x_i^{p_i}$는 $L^p$-노름을 초월한 수치 스킴의 정성적 성질 분석에 어떻게 활용될 수 있는가?
• RQ5 $b^h$와 $\sigma^h$에 대한 어떤 모멘트 및 성장 조건이 $V$-리아푸노프 함수 하에서 강한 수렴성과 기대값에서의 유계성을 보장하는가?

주요 결과

• 제안된 억제된 명시적 스킴은 조건 $\mathcal{L}V(x) \leq \rho V(x)$ 하에서 $\mathds{E}V(\bar{X}_k) \leq e^{\rho T} \mathds{E}V(\bar{X}_0)$ 를 만족하며, 이는 연속 시간의 적분 가능성 결과를 모방한다.
• 점근적 안정성의 경우, $\mathcal{L}V(x) \leq -z(x)$ 이고 $z(x) \geq \rho V(x)$ 라면 $\mathds{E}V(\bar{X}_k) \leq e^{-\rho k h} \mathds{E}V(\bar{X}_0)$ 를 도출하여 기대값에서의 지수 감쇠를 증명한다.
• 스킴은 $L^{4+4/r}$ 모멘트 조건 하에서 $L^2$-수렴성을 확보하며 오차 유계 $\mathds{E}|\bar{X}_{t,x} - \tilde{X}_{t,x}|^2 \leq K(1 + |x|^{2+2/r})h^2$ 를 만족한다.
• 이 방법은 스킴이 해의 단조성과 순서 성질을 유지하도록 함으로써 스칼라 SDE의 비교 정리를 복원할 수 있다.
• 스킴은 수렴을 위해 $4\kappa$ 모멘트만 필요하며, 균형 스킴과 유사하지만, 투영 메커니즘 덕분에 실질적으로 더 적은 모멘트 요구 조건을 가진다.
• 기대값에서 $X_t^-$에 그론월라 추론를 적용하면 $\mathds{E}X_t^- = 0$ 를 도출할 수 있으며, 이는 적절한 드라이프 구조 하에서 $X_t \geq 0$ 거의 확실하게 성립함을 의미한다.