[논문 리뷰] Error estimates of the backward Euler-Maruyama method for multi-valued stochastic differential equations

Monika Eisenmann, Mihály Kovács외 2명|arXiv (Cornell University)|2019. 06. 27.

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

한 줄 요약

이 논문은 볼록이고 미분 가능하지 않은 포텐셜을 가진 다가우수한 확률미분방정식(MSDE)에 대해 백워드 오일러-마루야마 방법의 수렴성을 확립한다. 루트-평균-제곱 노름에서 최소 1/4의 강한 수렴성을 증명하며, 결정론적 오차 분석 기법과 다가우수한 비선형성 조건을 활용한다. 이는 확률적 경사 하강법과 반산분기적 확률적 p-라플라스 방정식에 적용 가능하다.

ABSTRACT

In this paper, we derive error estimates of the backward Euler-Maruyama method applied to multi-valued stochastic differential equations. An important example of such an equation is a stochastic gradient flow whose associated potential is not continuously differentiable, but assumed to be convex. We show that the backward Euler-Maruyama method is well-defined and convergent of order at least $1/4$ with respect to the root-mean-square norm. Our error analysis relies on techniques for deterministic problems developed in [Nochetto, Savar\'e, and Verdi, Comm.\ Pure Appl.\ Math., 2000]. We verify that our setting applies to an overdamped Langevin equation with a discontinuous gradient and to a spatially semi-discrete approximation of the stochastic $p$-Laplace equation.

연구 동기 및 목표

• 비미분 가능하고 볼록한 포텐셜을 가진 다가우수한 확률미분방정식(MSDE)에 백워드 오일러-마루야마 방법을 적용한 강한 수렴 오차 추정을 유도하는 것.
• 전진 오일러-마루야마와 같은 명시적 방법은 초선형 성장 또는 비연속성 드리프트에서 발산하는 한계를 해결하는 것.
• 드리프트가 볼록 포텐셜의 준미분인 다가우수한 MSDE로 수치 해석을 확장하는 것, 예를 들어 비미분 가능 포텐셜을 가진 확률적 경사 하강법에서의 적용.
• 과다진동 랑주르 방정식의 비연속 기울기나 공간적으로 반산분기적 확률적 p-라플라스 방정식과 같은 구체적 문제에 적용 가능성을 검증하는 것.
• 정확한 해의 시간 정규성 조건을 요구하지 않고 수렴성을 확립하는 것 — 이는 준선형 SPDE에서 흔한 장애 요소이다.

제안 방법

• MSDE를 dX(t) + f(X(t)) dt ∋ b(X(t)) dt + g(X(t)) dW(t) 형태로 설정하며, f는 최대 단조성 연산자(예: 볼록 포텐셜의 준미분)이다.
• 백워드 오일러-마루야마 스킴을 적용: Xn ∈ Xn−1 − k f(Xn) + k b(Xn) + g(Xn−1) ΔWn, 등간격 시간 스텝을 사용한다.
• 일반화된 단조성 조건을 활용: 모든 v, w, z ∈ D(f)에 대해 ⟨fv − fz, z − w⟩ ≤ γ ⟨fv − fw, v − w⟩, 여기서 fv ∈ f(v) 등. 이 조건은 볼록 함수의 준미분에 대해 성립한다.
• 결정론적 오차 분석 기법(Nochetto 등, 2000)을 사용하여 이산 오차를 유계로 제한하며, 그론발라 타입의 추론을 피한다.
• b와 g가 전역 리프시츠 조건을 만족하고 f가 단조성 및 강력성 조건을 만족할 경우, 연속적 MSDE와 이산 스킴의 잘 정의됨을 증명한다.
• 반산분기 유한요소 근사에 의한 확률적 p-라플라스 방정식에 이 방법을 적용하여 문제를 유한차원 SDE로 매핑한다.

실험 결과

연구 질문

• RQ1 클래식한 미분 가능성 가정이 실패하는 비미분 가능하고 볼록한 드리프트를 가진 MSDE에 대해 백워드 오일러-마루야마 방법을 엄밀히 분석할 수 있는가?
• RQ2 드리프트가 비연속적이거나 비리프시츠일 경우, 이러한 MSDE에 대해 백워드 오일러-마루야마 방법의 수렴 속도는 어떠한가?
• RQ3 비미분 가능 포텐셜을 가진 확률적 경사 하강법에 적용했을 때, 이 방법은 여전히 잘 정의되고 수렴하는가?
• RQ4 시간 정규성이 해에 요구되지 않는 무한차원 문제, 예를 들어 반산분기 SPDE에 대해 오차 분석을 확장할 수 있는가?
• RQ5 f에 대한 일반화된 단조성 조건은 클래식한 리프시츠 연속성 조건 없이 어떻게 수렴 분석을 가능하게 하는가?

주요 결과

• 백워드 오일러-마루야마 방법은 볼록하고 비미분 가능한 드리프트를 가진 MSDE에 대해 L²(Ω; Rd) 노름에서 최소 1/4의 강한 수렴 차수를 보이며 잘 정의되고 수렴한다.
• 수렴 결과는 최소한의 가정 조건에서 성립한다: b와 g의 전역 리프시츠 연속성, f가 최대 단조성 연산자(예: 볼록 포텐셜의 준미분)임.
• 이 방법은 Φ(x) = |x|와 같이 0에서 다가우수한 기울기를 가진 과다진동 랑주르 방정식에 적용 가능하다.
• 이 방법은 p ∈ [1, 2)에 대해 공간적으로 반산분기된 확률적 p-라플라스 방정식에 대해서도 적용 가능하며, 이 경우 드리프트가 비연속적이다.
• 오차 유계는 해의 시간 정규성에 영향을 받지 않으며, 이는 이전 방법이 이러한 가정을 요구하는 것과 비교해 큰 이점이다.
• 오차 유계 Ck^{1/4}의 상수는 시간 스텝 k에 의존하지 않지만, 유한요소 차원 d에 따라 달라질 수 있어, 전체 시공간 이산화에서의 추가 분석이 필요하다.