IMEX–Crank–Nicolson方法求解Merton Jump-Diffusion PIDE（偏积分微分方程）：稳定性、收敛性与快速跳跃评估

《Axioms》：IMEX–Crank–Nicolson Methods for the Merton Jump-Diffusion PIDE: Stability, Convergence, and Fast Jump Evaluation

　　

研究人员提出了一种高效的数值框架，用于求解欧式期权定价中出现的Merton jump–diffusion（跳跃扩散）偏积分微分方程（Partial Integro-Differential Equation，PIDE）。为解决资产价格跳跃产生的非局部积分项，研究人员采用IMEX Crank–Nicolson时间步进格式，该格式保留了所得线性系统的三对角（tri-diagonal）结构。研究中引入了非均匀空间网格以及结合样条插值（spline interpolation）的快速高斯求积（Gaussian quadrature），以提高精度和计算效率。研究人员通过详细的理论分析确立了IMEX–Crank–Nicolson格式的无条件稳定性和收敛性。数值实验证实了理论结果，并说明了所提方法在代表性跳跃扩散参数下的有效性。

研究背景方面，自Merton引入泊松驱动的跳跃分量以捕捉资产价格不连续性以来，跳跃扩散动力学下的期权数值估值一直是活跃的研究课题。由此产生的偏积分微分方程（PIDE）因其非局部跳跃项带来了显著的数值挑战。早期的级数展开方法在短期到期或高跳跃强度下收敛缓慢，促使了有限差分法（FDM）等直接数值求解器的发展。隐式离散化导致稠密矩阵，计算昂贵，因此许多学者采用隐式显式（Implicit-Explicit，IMEX）格式，将扩散算子隐式处理，跳跃积分显式处理。尽管IMEX格式流行，但严格的稳定性和收敛性分析仍然有限，大多数研究强调数值性能而非理论保证。此外，均匀网格在捕捉行权价附近的非光滑性及跳跃影响区域时并非最优，非均匀网格与IMEX–CN格式的集成仍待发展。尽管跳跃积分的求积方法有所演进，但在IMEX–CN格式中集成非均匀网格及快速求积法的严格理论分析仍存在空白。为此，研究人员在《Axioms》发表该研究，旨在为Merton跳跃扩散PIDE的IMEX–CN方法提供全面的理论与数值分析，填补文献中关于严格稳定性收敛性分析、非均匀网格集成及快速求积法使用的空白。

关键技术方法方面，研究人员采用IMEX–Crank–Nicolson时间离散格式，对扩散项隐式处理而对跳跃积分项显式处理以保留三对角矩阵结构；空间离散采用非均匀网格并在均匀网格下使用二阶中心差分类比推导非均匀系数，结合快速高斯求积与三次样条插值评估对数正态跳跃的期望积分；理论分析通过傅里叶模式分析放大因子及引理证明插值算子的有界性以确立无条件稳定性；数值实验基于欧洲看涨期权标准参数集，在截断空间域上进行网格逐步加密的收敛性研究并与全隐式及经典Crank–Nicolson格式对比。

研究结果部分，各小节结论如下。

建模与PDE推导部分，研究人员通过在风险中性测度下建立资产价格的连续扩散与复合泊松跳跃过程，利用伊藤引理构建投资组合并对冲扩散风险，假设跳跃风险可分散后推导出包含非局部积分项的Merton PIDE，指出其相比Black–Scholes PDE数值求解难度更大。

离散格式与稳定性研究部分，研究人员定义了时空网格，采用后向欧拉或Crank–Nicolson进行时间离散及二阶中心差分进行空间离散，跳跃积分通过插值基求积显式处理。全隐式与Crank–Nicolson格式均保留三对角结构并通过Thomas算法求解。在稳定性分析中，研究人员证明扩散漂移算子在截断域上为退化椭圆且离散矩阵为M矩阵特征值负，Crank–Nicolson传播子无条件稳定；跳跃期望算子在最大范数下为压缩映射且有界；结合引理证明线性或三次样条插值算子有界且不破坏扩散分量无条件稳定性，从而严谨证明IMEM–CN格式对所有时间步长Δt无条件稳定，并通过傅里叶模式给出放大因子表达式。

比较结果与评估部分，研究人员分析全局收敛阶因非局部跳跃积分及显式处理限制通常在1至2阶之间约1.3至1.6；非均匀网格通过平滑映射聚类行权价附近点并提高有限差分适应性；快速高斯求积将跳跃积分转为期望并用少量Hermite节点与三次样条插值评估，降低50–70%计算成本；对比全隐式（FI）一阶光滑振荡少但引入数值扩散、经典Crank–Nicolson（CN）纯扩散稳定含跳跃可能振荡，IMEX–CN兼顾无条件稳定与二阶扩散精度且每时间步复杂度为O(N_x+N_q)。误差分析定理指出在支付Lipschitz连续、二阶空间差分、p阶插值及至少p阶求积下全局误差满足‖v?v^Δ‖≤C(Δt²+Δx^p)。

数值实验与分析部分，研究人员采用标准欧式看涨期权参数S₀=100、K=100、T=1、r=0.05、σ=0.2、λ=1、μ_J=0、σ_J=0.2及高精度参考解V_ref=15.0731。收敛研究表明立方样条插值(p=4)全局收敛率接近二阶，线性插值(p=1)全局率在1至1.6间与理论一致。可视化分析中期权价格随网格细化逼近参考曲线，空间误差分布在行权价K附近最大因支付非光滑及跳跃算子不规则，对数收敛图斜率约1.3–1.6对应线性插值结果。

讨论与未来方向部分，研究人员总结IMEX–CN结合隐式扩散与显式跳跃保留稀疏性，非均匀网格与快速求积提升分辨率，收敛率1.3–1.6符合理论，高跳跃强度与短期限增大期权价因大幅运动概率上升。未来方向包括扩展至美式期权用罚函数或前固定法、推广至Variance Gamma或CGMY等Lévy模型、引入自适应网格加密(AMR)动态优化分辨率、开发结合有限差分与FFT的混合求解器处理高维多资产问题。

结论部分翻译为，本研究开发并分析了Merton跳跃扩散PIDE的有限差分格式，尤其关注IMEX–Crank–Nicolson(IMEX–CN)方法。该框架将扩散算子隐式处理与跳跃积分显式评估相结合，保留了线性系统的三对角结构。在非均匀空间网格及结合样条插值的快速高斯求积集成下，提升了IMEX–CN框架内的数值分辨率。详细的理论分析确立了IMEX–CN格式的无条件稳定性和收敛性。数值实验支持了这些发现，并证明了不同跳跃强度和市况下的一致行为。观察到的介于1阶与2阶之间的收敛率与理论预测相符。金融解释方面，数值结果阐明了Merton模型内跳跃分量对欧式期权价值的影响。特别是较高的跳跃强度和较短的期限往往因资产大幅变动概率升高而增加期权价格。IMEX–CN格式在所提计算框架内捕捉了这些效应。

需要我帮你把这篇论文解读里的关键公式和定理（如Merton PIDE推导、稳定性证明引理）整理成结构化的笔记吗？