Compound Poisson Approximation for Stochastic Volterra Equations with Singular Kernels

基本信息

• 论文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

摘要

本文建立了随机微分方程(SDEs)和Volterra型SDEs在复合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$

其中: