[论文解读] Error estimates on ergodic properties of discretized Feynman-Kac semigroups

Grégoire Ferré, Gabriel Stoltz|arXiv (Cornell University)|Dec 11, 2017

Stochastic processes and financial applications参考文献 71被引用 13

一句话总结

本文为带有加性噪声的扩散过程的费曼-卡斯半群的时间离散化提供了严格的误差估计，分析了生成元的不变测度和主导特征值中的偏差。通过塔拉伊-图博罗型分析，建立了在紧流形上离散化动力学的收敛速率，为扩散蒙特卡洛及相关领域中的数值方案提供了理论依据。

ABSTRACT

We consider the numerical analysis of the time discretization of Feynman-Kac semigroups associated with diffusion processes. These semigroups naturally appear in several fields, such as large deviation theory, Diffusion Monte Carlo or non-linear filtering. We present errors estimates a la Talay-Tubaro on their invariant measures when the underlying continuous stochastic differential equation is discretized; as well as on the leading eigenvalue of the generator of the dynamics, which corresponds to the rate of creation of probability. This provides criteria to construct efficient integration schemes of Feynman-Kac dynamics, as well as a mathematical justification of numerical results already observed in the Diffusion Monte Carlo community. Our analysis is illustrated by numerical simulations.

研究动机与目标

• 分析与扩散过程相关的费曼-卡斯半群在时间离散化下引入的偏差。
• 量化时间离散化下生成元的不变测度和主导特征值的误差。
• 为扩散蒙特卡洛和非线性滤波中使用的数值方案提供理论依据。
• 将塔拉伊-图博罗误差估计扩展至具有重加权的非线性、测度值费曼-卡斯动力学。
• 建立在紧致定义域上对费曼-卡斯动力学进行高效数值积分的收敛准则。

提出的方法

• 分析通过带加性噪声的随机微分方程和权重函数 $W$ 定义的连续时间费曼-卡斯半群。
• 使用有限时间步长 $\Delta t$ 离散化SDE，通过求积规则近似时间积分。
• 研究得到的离散时间马尔可夫链 $q^n$ 及其相关半群 $\Phi_{\Delta t,n}^W$。
• 应用塔拉伊-图博罗型误差分析，推导离散不变测度 $\nu_{W,\Delta t}$ 向连续极限 $\nu_W$ 的收敛速率。
• 利用多布拉辛的混合系数和压缩范数，控制离散半群的混合性质和收敛速度。
• 通过期望比值的统一有界性，证明主导特征值 $\lambda_{\Delta t} = \tilde{\lambda}_{\Delta t}$，从而证明特征值近似的一致性。

实验结果

研究问题

• RQ1 时间离散化如何影响费曼-卡斯半群的不变测度？
• RQ2 离散化不变测度向真实不变测度的收敛速率是多少？
• RQ3 离散化误差如何影响生成元的主导特征值，该特征值控制概率生成速率？
• RQ4 塔拉伊-图博罗框架能否被适配至非线性、测度值半群（如费曼-卡斯过程）？
• RQ5 在何种条件下，离散化半群能保持连续动力学的遍历性质？

主要发现

• 本文证明了离散化费曼-卡斯半群的不变测度 $\nu_{W,\Delta t}$ 在 $\Delta t \to 0$ 时收敛于真实不变测度 $\nu_W$，误差界通过塔拉伊-图博罗分析得出。
• 离散化生成元的主导特征值 $\lambda_{\Delta t}$ 收敛于真实特征值 $\lambda$，且通过期望比值的统一有界性，证明了 $\lambda_{\Delta t} = \tilde{\lambda}_{\Delta t}$，从而确保特征值近似的准确性。
• 离散半群的压缩范数被有界于 $1 - \alpha(Q)$，其中 $\alpha(Q)$ 为多布拉辛混合系数，确保了向不变测度的指数收敛。
• 不变测度的误差受向量场 $b$、权重函数 $W$ 以及离散化方案的弱阶的正则性控制。
• 数值模拟验证了理论误差估计，并在实际场景中确认了收敛速率。
• 该分析为在扩散蒙特卡洛中使用高阶离散化方案提供了依据，可有效降低基态能量估计中的偏差。