[论文解读] Two-dimensional Kac-Rice formula. Application to shot noise processes excursions

Raphaël Lachièze-Rey|arXiv (Cornell University)|Jul 19, 2016

Simulation Techniques and Applications被引用 4

一句话总结

本文推导出一个二维 Kac-Rice 公式，将 $\mathbb{R}^2$ 上一个确定性 $ C^2 $ 函数 $ f: \mathbb{R}^2 \to \mathbb{R} $ 的上穿集的欧拉示性数表示为 $ \mathbb{R}^2 $ 上的积分，该积分涉及 $ f $、其梯度 $ \nabla f $ 以及其黑塞矩阵。该公式使得在不假设场或其导数具有密度的条件下，能够计算包括散粒噪声过程在内的随机场的期望欧拉示性数，并在较弱的正则性条件下保证连续性与矩的有限性。

ABSTRACT

Given a deterministic function f:R^2->R atisfying suitable assump- tions, we show that for h smooth with compact support, the integral of the Euler characteristic of the excursion set of f above some level u against a test function h corresponds to the Lebesgue integral on R^2 of a bounded quantity depending on grad(f)(x),h(f(x)),h'(f(x)) and \\partial_{ii}f(x),i = 1,2. This formula can be seen as a 2-dimensional analogue of Kac-Rice formula. It yields in particular that the left hand member is continuous in the argument f, for an appropriate norm on the space of C2 functions. If f is a random field, the expectation can be passed under integrals in this identity under minimal requirements, not involving any density assumptions on the marginals of f or his derivatives. We apply these results to give a weak expression of the mean Euler characteristic of a shot noise process, and the finiteness of its moments.

研究动机与目标

• 推导二维上穿集欧拉示性数的确定性积分公式，扩展经典 Kac-Rice 方法。
• 在 $ C^2 $ 拓扑下建立欧拉示性数泛函的连续性，从而实现对随机场扰动的稳健分析。
• 在最小概率假设下，为随机 $ C^2 $ 场（特别是散粒噪声过程）提供均值欧拉示性数的弱公式化表达。
• 避免对场或其导数依赖密度条件，从而扩大其在非高斯与非正则场中的适用范围。
• 通过所推导的公式，证明散粒噪声过程的欧拉示性数泛函的矩是有限的。

提出的方法

• 通过应用面积分与变分技术，推导出二维 Kac-Rice 公式，将上穿集的欧拉示性数表示为 $ \mathbb{R}^2 $ 上的积分。
• 引入泛函 $ \chi_f(h) = \int_\mathbb{R} h(u) \chi(\{f \geq u\}) du $，即欧拉原函数，并证明其等于涉及 $ h(f(x)) $、$ h'(f(x)) $、$ \partial_i f(x) $ 与 $ \partial_{ii}^2 f(x) $ 的勒贝格积分。
• 基于梯度分量的符号，定义平面区域 $ Q_1, Q_2 \subset \mathbb{R}^2 $，以局部化对欧拉示性数的贡献。
• 采用 $ \mathcal{C}^{1,1} $ 正则性（几乎处处二次可微）而非莫尔斯条件，从而将适用范围扩展至临界点计数之外。
• 通过利用平稳性与各向同性，将公式应用于散粒噪声过程，以计算特征函数并推导矩的有限性。
• 利用 $ \int_0^\infty \frac{\sin(ux)}{u} du = \frac{\pi}{2} \text{sign}(u) $ 对符号函数的非正常积分表示，处理涉及导数符号函数乘积的期望表达式。

实验结果

研究问题

• RQ1 是否可以不依赖临界点计数，将二维上穿集的欧拉示性数表示为定义域上的光滑积分？
• RQ2 在不假设莫尔斯条件或密度条件的情况下，欧拉示性数泛函在 $ C^2 $ 拓扑收敛下是否仍保持连续性？
• RQ3 是否可以在不假设其边际或导数密度的条件下，计算随机场的均值欧拉示性数？
• RQ4 在何种条件下可保证散粒噪声过程的欧拉示性数矩的有限性？
• RQ5 如何通过弱公式化方法，将 Kac-Rice 公式推广至非高斯、非莫尔斯场？

主要发现

• 上穿集 $ \{f \geq u\} $ 的欧拉示性数被表示为 $ \mathbb{R}^2 $ 上的勒贝格积分，涉及 $ h'(f(x)) \partial_i f(x)^2 $ 与 $ h(f(x)) \partial_{ii}^2 f(x) $，并由梯度区域 $ Q_i $ 上的指示函数加权。
• 该公式在 $ \mathcal{C}^{1,1} $ 正则性下成立，使其可应用于非莫尔斯函数，并确保欧拉原函数在 $ C^2 $ 范数下的连续性。
• 对于随机 $ C^2 $ 场，其欧拉原函数的期望可在不假设 $ f(x) $ 或其导数具有密度的条件下计算，从而实现对非高斯场的分析。
• 散粒噪声过程的均值欧拉示性数被证明为有限，且在场及其导数的温和矩条件下，其矩也被证明为有限。
• 对于各向同性、平稳的散粒噪声场，欧拉原函数的特征函数可通过 $ (f(0), \nabla f(0), \partial_{ii}^2 f(0)) $ 的特征函数的导数表示，从而在特定情形下实现显式计算。
• 该公式允许通过弱收敛与 $ o(|W_n|) $ 误差控制，推导出大区域中欧拉示性数的渐近行为。