특이 핵을 가진 확률적 Volterra 방정식에 대한 복합 Poisson 근사

기본 정보

• 논문 ID: 2510.27323
• 제목: Compound Poisson Approximation for Stochastic Volterra Equations with Singular Kernels
• 저자: Xicheng Zhang (베이징 이공대학교), Yuanlong Zhao (우한대학교)
• 분류: math.PR (확률론)
• 제출 시간: 2025년 10월 31일
• 논문 링크: https://arxiv.org/abs/2510.27323

초록

본 논문은 확률미분방정식(SDE)과 Volterra형 SDE에 대한 복합 Poisson 과정 근사 하에서의 강 수렴성을 확립하고 명시적 수렴 속도를 도출한다. 고전적인 Euler-Maruyama 방법과 비교하여, 복합 Poisson 근사의 핵심 장점은 시간 변수에 대한 표류 계수의 연속성을 요구하지 않으며, 심지어 특이성 문제를 처리할 수 있다는 것이다. 수치 실험은 본 방법의 안정성을 검증한다.

연구 배경 및 동기

1. 핵심 문제

본 논문은 확률미분방정식과 확률적 Volterra 방정식의 수치 근사 문제를 연구하며, 특히 계수가 시간 특이성을 갖는 경우에 초점을 맞춘다.

2. 문제의 중요성

• 이론적 의의: 확률적 Volterra 방정식은 확률 분석에서 중요한 위치를 차지하며, 특히 기억 효과를 갖는 확률 시스템을 기술할 때 그러하다
• 응용 가치: 분수 Brown 운동이 구동하는 확률 방정식, 확률적 Navier-Stokes 방정식 등 중요한 모델들이 이러한 유형의 문제로 귀결된다
• 수치적 도전: 핵 함수가 특이성을 가질 때, 전통적인 수치 방법은 수렴성과 안정성 문제에 직면한다

3. 기존 방법의 한계

• Euler-Maruyama 방법: 시간 변수에 대한 Hölder 연속성을 요구하며, 불연속 또는 특이한 경우를 처리할 수 없다
• 개선된 EM 방법: 일부 조건(예: one-sided Lipschitz 조건)을 완화하지만, 여전히 특정 시간 정규성이 필요하다
• 특이 핵에 대한 방법: 기존 문헌은 주로 비특이 핵에 초점을 맞추며, 특이 핵의 처리 이론이 불완전하다

4. 연구 동기

• 시간 특이성을 처리할 수 있는 수치 방법 개발
• 명시적 수렴 속도 추정 확립
• 분수 Brown 운동 등 중요한 확률 과정에 대한 신뢰할 수 있는 수치 도구 제공

핵심 기여

1. 복합 Poisson 근사의 강 수렴 이론 확립: SDE와 확률적 Volterra 방정식에 대한 복합 Poisson 근사 하에서의 강 수렴성을 처음으로 체계적으로 증명
2. 명시적 수렴 속도 도출:
   • 표준 SDE의 경우: 수렴 속도는 $\varepsilon^{\gamma \wedge \frac{\beta}{2}}$
   • 확률적 Volterra 방정식의 경우: 수렴 속도는 $\varepsilon^{\gamma/(2(2+\gamma))}$
3. 특이 계수의 획기적 처리: 본 방법은 표류 계수의 시간 연속성을 요구하지 않으며 시간 특이성을 처리할 수 있다
4. 완전한 이론 체계 제공: 모멘트 추정, Gronwall 부등식 적용, 기술적 보조정리 등 완전한 증명 체계 포함
5. 구체적 응용 사례 제공: 분수 Brown 운동(Hurst 매개변수 $H \in (0,1/2) \cup (1/2,1)$) 경우에 대한 가정 조건의 상세한 검증
6. 수치 검증: 두 개의 수치 실험을 통해 방법의 유효성과 안정성 검증

방법 상세 설명

작업 정의

문제 1: 표준 SDE$\mathbb{R}^d$ 위의 확률미분방정식을 고려: $X_t = X_0 + \int_0^t \sigma(s, X_s) dW_s + \int_0^t b(s, X_s) ds$

여기서 $\sigma: \mathbb{R}_+ \times \mathbb{R}^d \to \mathbb{R}^d \otimes \mathbb{R}^m$, $b: \mathbb{R}_+ \times \mathbb{R}^d \to \mathbb{R}^d$, $W_t$는 $m$차원 표준 Brown 운동이다.

문제 2: 확률적 Volterra 방정식$Y_t = Y_0 + \int_0^t \sigma(t, s, Y_s) dW_s + \int_0^t b(t, s, Y_s) ds$

핵심 특징은 계수가 현재 시간 $t$와 적분 변수 $s$에 의존하며, 기억 효과를 나타낸다.

모델 구조

1. 복합 Poisson 과정 구성

기본 Poisson 과정: $(T_k)_{k \in \mathbb{N}}$을 매개변수 1인 독립 동일 분포 지수 확률변수라 하자: $P(T_k \geq t) = e^{-t}, \quad t \geq 0$

점프 시간 정의: $S_0 = 0$, $S_n = S_{n-1} + T_n$

표준 Poisson 과정: $N_t = \max\{n: S_n \leq t\}$

스케일된 Poisson 과정: 주어진 $\varepsilon > 0$에 대해 정의: $N_t^\varepsilon := \varepsilon N_{t/\varepsilon}, \quad \tilde{N}_t^\varepsilon := N_t^\varepsilon - t$

여기서 $N_t^\varepsilon$는 점프 크기가 $\varepsilon$, 강도가 $1/\varepsilon$인 Poisson 과정이다.

복합 Poisson 과정: $W_{N_t^\varepsilon} = \sum_{k=1}^{N_{t/\varepsilon}} (W_{k\varepsilon} - W_{(k-1)\varepsilon})$

이는 강도 측도가 $\nu_\varepsilon \sim N(0, \varepsilon I_d)$인 복합 Poisson 과정이다.

2. 수치 격식

SDE의 복합 Poisson 근사: $X_t^\varepsilon = X_0 + \int_0^t \sigma(s, X_{s-}^\varepsilon) dW_{N_s^\varepsilon} + \int_0^t b(s, X_{s-}^\varepsilon) dN_s^\varepsilon$

완전 이산화 형식: $X_t^\varepsilon = X_0 + \sum_{k=1}^{N_{t/\varepsilon}} [\sigma(S_k^\varepsilon, X_{S_{k-1}^\varepsilon}^\varepsilon)(W_{k\varepsilon} - W_{(k-1)\varepsilon}) + \varepsilon b(S_k^\varepsilon, X_{S_{k-1}^\varepsilon}^\varepsilon)]$

확률적 Volterra 방정식의 복합 Poisson 근사: $Y_t^\varepsilon = Y_0 + \int_0^t \sigma(t, s, Y_{s-}^\varepsilon) dW_{N_s^\varepsilon} + \int_0^t b(t, s, Y_{s-}^\varepsilon) dN_s^\varepsilon$

완전 이산화 형식: $Y_t^\varepsilon = Y_0 + \sum_{k=1}^{N_{t/\varepsilon}} [\sigma(t, S_k^\varepsilon, Y_{S_{k-1}^\varepsilon}^\varepsilon)(W_{k\varepsilon} - W_{(k-1)\varepsilon}) + \varepsilon b(t, S_k^\varepsilon, Y_{S_{k-1}^\varepsilon}^\varepsilon)]$

기술적 혁신점

1. 확률적 시간 이산화

확정적 시간 격자와 달리, Poisson 과정의 점프 시간을 이산화 점으로 사용하며, 이는 확률적 이산화 전략이고 핵심 장점은:

• 자연스럽게 확률성을 도입하며 Brown 운동의 증분과 독립적
• 계수의 시간 연속성이 필요하지 않음
• 시간 특이성을 처리할 수 있음

2. 오차 분해 전략

SDE의 경우, 오차를 네 항으로 분해: $X_t - X_t^\varepsilon = I_1(t) + I_2(t) + I_3(t) + I_4(t)$

여기서:

• $I_1(t)$: 실제 경로와 확률적 시간 점에서의 값의 차이
• $I_2(t)$: Brown 운동과 복합 Poisson 과정의 차이 (가장 핵심 항)
• $I_3(t)$: 서로 다른 경로에서 계수의 차이
• $I_4(t)$: 표류항의 이산화 오차

3. 핵심 기술 보조정리

보조정리 2.1 (Poisson 점프 시간의 모멘트 추정): $\alpha \in (0,1]$, $\beta \geq 0$에 대해: $\mathbb{E}|r^\alpha - (S_k^\varepsilon)^\alpha|^\beta \lesssim (k\varepsilon)^{\alpha\beta} k^{-\beta/2}, \quad r \in [(k-1)\varepsilon, k\varepsilon]$

이 추정은 직접적으로 $|a^\alpha - b^\alpha| \leq |a-b|^\alpha$를 사용하는 것보다 더 정밀하다.

보조정리 2.3 (Poisson 과정 편차 추정): 임의의 $p > 0$에 대해: $\mathbb{E}|N_t^\varepsilon - t|^p \leq C\varepsilon^{p/2}(t^{p/2} \vee t)$

4. $I_2(t)$ 항의 정밀한 처리

이는 가장 어려운 항이며, 다음을 활용해야 함:

• 보정된 Poisson 확률 측도 표현
• Doob 극대 부등식
• 시간 변수의 Hölder 연속성 가정 $(H_\sigma^t)$
• 정밀한 합 추정

Volterra 방정식의 경우, $I_2$의 처리는 더욱 복잡하며 중간 이산화 매개변수 $\delta$를 도입하고 네 개의 부분항으로 분해해야 한다.

실험 설정

데이터 집합 (수치 실험 시나리오)

실험 1: 특이 시간 의존 표류를 갖는 선형 SDE

$dX_t = \mu(t)X_t dt + \sigma(t)X_t dW_t$

여기서: