[论文解读] First order strong approximations of scalar SDEs with values in a domain

Andreas Neuenkirch, Łukasz Szpruch|arXiv (Cornell University)|Sep 3, 2012

Stochastic processes and financial applications参考文献 18被引用 3

一句话总结

本文提出了一种基于Lamperti变换的后向欧拉-马里亚米方法，用于强逼近具有非利普希茨系数且保持在某一区域内的单变量SDE。通过将SDE变换为具有恒定扩散系数的形式，应用无漂移隐式欧拉格式，并反向变换，该方法保持了区域边界并实现了强收敛阶为一的逼近，适用于金融与生物学中的关键模型，如CIR、CEV和Ait-Sahalia模型。

ABSTRACT

We are interested in strong approximations of one-dimensional SDEs which have non-Lipschitz coefficients and which take values in a domain. Under a set of general assumptions we derive an implicit scheme that preserves the domain of the SDEs and is strongly convergent with rate one. Moreover, we show that this general result can be applied to many SDEs we encounter in mathematical finance and bio-mathematics. We will demonstrate flexibility of our approach by analysing classical examples of SDEs with sublinear coefficients (CIR, CEV models and Wright-Fisher diffusion) and also with superlinear coefficients (3/2-volatility, Ait-Sahalia model). Our goal is to justify an efficient Multi-Level Monte Carlo (MLMC) method for a rich family of SDEs, which relies on good strong convergence properties.

研究动机与目标

• 开发一种能够保持标量SDE在非利普希茨漂移或扩散系数下区域边界的数值格式。
• 实现此类SDE的强收敛阶一，克服经典显式格式的局限性。
• 提供一个适用于数学金融与生物数学中具有次线性和超线性系数的SDE的通用框架。
• 通过确保高阶强收敛，为多级蒙特卡洛方法的效率提供理论依据。
• 建立变换后后向欧拉格式与原始SDE上无漂移隐式米尔斯坦格式之间的联系。

提出的方法

• 应用Lamperti变换，将原始具有状态依赖扩散系数的SDE转化为具有恒定扩散系数的SDE。
• 使用后向欧拉-马里亚米（BEM）格式近似变换后的SDE，以确保稳定性与区域保持。
• 通过逆Lamperti变换将BEM格式的解反变换，得到原始SDE的数值格式。
• 采用无漂移隐式格式以处理非利普希茨系数，保持正性与单调性。
• 通过先验估计与矩界建立收敛性，利用变换后方程中漂移项的单边利普希茨性质。
• 通过逆变换的渐近展开，将变换后SDE上的BEM格式与原始SDE上的无漂移隐式米尔斯坦格式联系起来。

实验结果

研究问题

• RQ1 能否构造一种数值格式，使其在保持标量SDE非利普希茨系数下区域边界的前提下，实现强收敛阶一？
• RQ2 Lamperti变换如何与隐式格式有效结合，以处理次线性和超线性扩散系数？
• RQ3 原始SDE系数需满足何种条件，才能保证所提格式的强收敛性？
• RQ4 变换后SDE上的后向欧拉-马里亚米格式是否等价于原始SDE上的已知无漂移隐式米尔斯坦格式？
• RQ5 该框架能否推广以支持金融与生物学中SDE的高效多级蒙特卡洛方法？

主要发现

• 在温和假设下，所提出的Lamperti-Euler格式在离散点上的最大$L^p$-距离下实现了强收敛阶一。
• 对于CIR过程，当满足$\frac{\kappa\theta}{\sigma^2} > \frac{3}{2}$时，该格式实现了强收敛阶一，优于先前结果中$\frac{1}{2}$阶收敛的结论。
• 该方法保持了对初值的正性与单调性，而标准显式格式则不具备此性质。
• 变换后SDE上的后向欧拉格式等价于原始SDE上的无漂移隐式米尔斯坦格式，提供了新的解释与验证。
• 该框架适用于广泛的SDE类，包括CEV模型、Wright-Fisher扩散、3/2波动率模型以及Ait-Sahalia模型，均实现了强收敛阶一。
• 针对无漂移隐式米尔斯坦格式的分段线性插值，推导出一个精确的误差估计，表明在相同条件下有$\mathbb{E}\left[\max_{t\in[0,T]}|y(t)-\overline{Z}_t|^2\right]^{1/2} \leq C \sqrt{\log(\Delta t)} \sqrt{\Delta t}$。