[论文解读] A Central Limit Theorem for Fleming-Viot Particle Systems with Hard Killing

Bernard Delyon, Frédéric Cérou+2 more|arXiv (Cornell University)|Sep 20, 2017

Stochastic processes and statistical mechanics参考文献 15被引用 6

一句话总结

本文在最小假设下建立了带硬吸收的佛朗西斯-维奥特粒子系统的中心极限定理（CLT），证明了在固定时间 T 时对生存概率和条件分布的估计量的渐近正态性。该方法依赖于鞅分解、二次变差分析以及时间一致的 L² 估计，将先前结果扩展至包含具有吸收边界条件的椭圆扩散过程。

ABSTRACT

Fleming-Viot type particle systems represent a classical way to approximate the distribution of a Markov process with killing, given that it is still alive at a final deterministic time. In this context, each particle evolves independently according to the law of the underlying Markov process until its killing, and then branches instantaneously on another randomly chosen particle. While the consistency of this algorithm in the large population limit has been recently studied in several articles, our purpose here is to prove Central Limit Theorems under very general assumptions. For this, we only suppose that the particle system does not explode in finite time, and that the jump and killing times have atomless distributions. In particular, this includes the case of elliptic diffusions with hard killing.

研究动机与目标

• 在先前研究需软吸收或有界 carré du champ 假设的限制之外，建立带硬吸收的佛朗西斯-维奥特粒子系统的中心极限定理（CLT）。
• 在一般条件下，分析对固定时间 T 的生存概率 $ p_T $ 和条件分布 $ \eta_T $ 的估计量的渐近行为。
• 在一般设定下，证明归一化粒子系统波动的收敛分布，并确定渐近方差。
• 通过假设跳跃和吸收时间分布无原子性，处理此前因缺乏有界吸收强度而被排除的、在区域边界被吸收的椭圆扩散过程的情形。
• 通过在底层过程的最小正则性条件下证明 CLT，为罕见事件估计中的序贯蒙特卡洛方法提供严格的理论基础。

提出的方法

• 通过粒子系统经验测度和生存概率估计量的鞅分解，将信号与噪声分离。
• 对粒子过程 $ p_t^N $ 使用时间一致的 $ \mathbb{L}^2 $-估计，确保在整个时间区间 $[0,T]$ 内的稳定性。
• 应用二次变差估计和涉及 $ p_t^N = (1 - 1/N)^{N\mathcal{N}_t} $ 的分部积分公式，其中 $ \mathcal{N}_t $ 表示分支次数。
• 通过二次变差过程的收敛性推导渐近方差表达式，利用跳跃和吸收时间分布无原子性的假设。
• 使用停止时间论证和独立性结构，证明跳跃事件在几乎必然意义下非同时发生，从而确保粒子系统的良好定义性。
• 在满足收敛到高斯极限的条件下，对归一化波动过程应用鞅中心极限定理。

实验结果

研究问题

• RQ1 在何种最小假设下，具有硬吸收的佛朗西斯-维奥特粒子系统满足中心极限定理？
• RQ2 当 $ N \to \infty $ 时，估计量 $ \hat{p}_T^N $ 和 $ \hat{\eta}_T^N $ 的渐近分布行为如何？其极限方差的形式是什么？
• RQ3 对于在击中边界时被吸收的椭圆扩散过程（吸收强度无界，即硬吸收），尽管缺乏有界 carré du champ 条件，是否仍可建立 CLT？
• RQ4 跳跃和吸收时间的无原子分布在此类分支事件非同时发生及 CLT 有效性的保障中起何作用？
• RQ5 如何用底层马氏过程和时间改变过程 $ Q^{T-t}(\varphi) $ 表征粒子系统的渐近方差？

主要发现

• 本文在粒子系统在有限时间内不爆炸且吸收/跳跃时间分布无原子性的最小假设下，证明了带硬吸收的佛朗西斯-维奥特粒子系统的中心极限定理。
• 归一化波动 $ p_T^N $ 的极限分布为高斯分布，其渐近方差有明确定义，且依赖于表示测试函数条件期望的时间改变过程 $ Q^{T-t}(\varphi) $。
• 渐近方差被表征为粒子系统经验测度分解中鞅分量的二次变差过程的极限。
• 该 CLT 成立于具有吸收边界条件的椭圆扩散过程（即硬吸收），此情形在早期工作中因需软吸收或有界 carré du champ 条件而被排除。
• 证明中建立了 $ p_t^N $ 的时间一致 $ \mathbb{L}^2 $-估计，这对控制整个时间区间 $[0,T]$ 内的波动至关重要。
• 附录中的分部积分公式 (4.3)–(4.5) 允许对涉及 $ p_t^N $ 的随机积分进行精确逼近，误差项有界于 $ O(1/N) $，这对渐近分析至关重要。