[论文解读] Taming singular stochastic differential equations: A numerical method

Khoa Lê, Chengcheng Ling|arXiv (Cornell University)|Oct 4, 2021

Stochastic processes and financial applications参考文献 36被引用 15

一句话总结

本文提出了一种针对具有可积漂移项且满足Ladyzhenskaya–Prodi–Serrin条件的多维时变非齐次SDEs的驯服伊藤-丸柳马 scheme。该方法在漂移项适当逼近的条件下，建立了具有对数因子的 $1/2$ 阶强收敛率，从而实现了对奇异SDEs和随机输运方程的数值求解。

ABSTRACT

We consider a generic and explicit tamed Euler--Maruyama scheme for multidimensional time-inhomogeneous stochastic differential equations with multiplicative Brownian noise. The diffusive coefficient is uniformly elliptic, Hölder continuous and weakly differentiable in the spatial variables while the drift satisfies the strict Ladyzhenskaya--Prodi--Serrin condition, as considered by Krylov and Röckner (2005). In the discrete scheme, the drift is tamed by replacing it by an approximation. A strong rate of convergence of the scheme is provided in terms of the approximation error of the drift in a suitable and possibly very weak topology. A few examples of approximating drifts are discussed in detail. The parameters of the approximating drifts can vary and -- under suitable conditions -- be fine-tuned to achieve a strong convergence rate which is arbitrarily close to the benchmark $0.5$ rate. The result is then applied to provide numerical solutions for stochastic transport equations with singular vector fields satisfying the aforementioned condition.

研究动机与目标

• 开发一种针对奇异（可积）漂移项且不适用于标准伊藤-丸柳马 scheme 的SDEs 的数值稳定且收敛的格式。
• 解决在一般一致椭圆扩散系数下，SDEs 在 Ladyzhenskaya–Prodi–Serrin 条件下的强收敛率这一开放问题。
• 为由奇异向量场驱动的随机输运方程的数值逼近提供一个框架。
• 在弱拓扑意义下，以漂移项逼近误差为基准建立收敛速率，从而支持灵活且精细调节的数值格式。

提出的方法

• 提出一种驯服伊藤-丸柳马 scheme，其中漂移项被近似 $b^n$ 替代，以避免模拟过程中的奇异性。
• 采用时间离散格式，使用分段常数系数：$dX_t^n = b^n(t, X_{k_n(t)}^n)dt + \sigma(t, X_{k_n(t)}^n)dB_t$，确保稳定性。
• 采用逼近技术，如截断 $b^n_r(x) = b_r(x)\mathbf{1}_{(|b_r(x)| \leq Cn^\chi \|b_r\|_{L_p})}$ 和与高斯核的卷积 $p_{1/n^\chi} * b_r(x)$。
• 通过噪声正则化和Zvonkin变换控制分析中奇异漂移的影响。
• 利用随机缝合法与定量Khasminskii引理推导路径样本估计，并控制真实解与近似解之间的误差。
• 通过在弱 $L^p$-型拓扑中以漂移项逼近误差为基准，建立 $L^2$-范数意义下的强收敛性。

实验结果

研究问题

• RQ1 能否为仅可积而非局部有界的漂移项SDEs 构造出稳定且收敛的数值格式？
• RQ2 在 Ladyzhenskaya–Prodi–Serrin 条件下，驯服伊藤-丸柳马 scheme 的最优强收敛率是多少？
• RQ3 漂移项逼近方式（截断、卷积）的选择如何影响格式的收敛速率与稳定性？
• RQ4 该格式能否推广用于求解具有奇异向量场的随机输运方程？
• RQ5 漂移项逼近的弱拓扑在决定数值格式收敛速率方面起到何种作用？

主要发现

• 在适当调节逼近参数的条件下，驯服伊藤-丸柳马 scheme 在具有对数因子的条件下实现了 $1/2$ 阶强收敛率。
• 收敛速率以漂移项在弱 $L^p$-型拓扑中的逼近误差为基准推导，从而支持多种逼近方案的应用。
• 当参数精细调节时，截断与正则化漂移逼近（如与高斯核卷积）均能实现相同的收敛速率。
• 该方法适用于满足 Ladyzhenskaya–Prodi–Serrin 条件的奇异向量场驱动的随机输运方程。
• 分析依赖于高级工具，如随机缝合法、Zvonkin变换与定量Khasminskii引理，以控制奇异漂移的影响。
• 该结果解决了可积漂移项SDE数值分析中的一个开放问题，为多水平蒙特卡洛方法提供了计算上可行的方案。