[论文解读] Expected Number and Height Distribution of Critical Points of Smooth Isotropic Gaussian Random Fields

Dan Cheng, Armin Schwartzman|arXiv (Cornell University)|Nov 21, 2015

Geometry and complex manifolds参考文献 36被引用 4

一句话总结

本文推导了欧氏空间和任意维度球面上平滑各向同性高斯随机场中临界点的期望数量与高度分布的显式公式，消除了对协方差函数的先前限制。通过使用高斯正交不变（GOI）矩阵——其为高斯正交系综（GOE）的推广——来表示海森矩阵，作者获得了仅依赖于一个或两个参数的可计算表达式，适用于所有各向同性场，包括参数空间边界处的随机拉普拉斯特征函数。

ABSTRACT

We obtain formulae for the expected number and height distribution of critical points of smooth isotropic Gaussian random fields parameterized on Euclidean space or spheres of arbitrary dimension. The results hold in general in the sense that there are no restrictions on the covariance function of the field except for smoothness and isotropy. The results are based on a characterization of the distribution of the Hessian of the Gaussian field by means of the family of Gaussian orthogonally invariant (GOI) matrices, of which the Gaussian orthogonal ensemble (GOE) is a special case. The obtained formulae depend on the covariance function only through a single parameter (Euclidean space) or two parameters (sphere), and include the special boundary case of random Laplacian eigenfunctions.

研究动机与目标

• 推导平滑各向同性高斯随机场在 ℝᴺ 和 𝕊ᴺ 上的临界点期望数量，且不施加对协方差函数的限制。
• 将费约多罗夫对临界点期望的显式计算结果推广至协方差函数并非对所有维度 N 都有效的场类之外。
• 通过高斯正交不变（GOI）矩阵表征场的海森矩阵，这是高斯正交系综（GOE）的推广。
• 推导仅依赖于协方差函数在零点的一阶与二阶导数的临界点期望数量与高度分布的可计算公式。
• 将随机拉普拉斯特征函数作为参数空间边界处的特殊情况处理，其满足亥姆霍兹方程，且其海森矩阵与场之间呈现退化相关性。

提出的方法

• 以 Kac-Rice 公式作为计算临界点期望数量的基础工具。
• 将高斯场的海森矩阵表示为 GOI 矩阵，其通过协方差参数 c ≥ –1/N 允许对角元素之间存在依赖，从而推广了 GOE。
• 利用 GOI 矩阵分布推导海森矩阵特征值的有序联合密度，从而实现对临界点指标的积分。
• 通过 GOI 特征值密度上的积分表达临界点指数 i 且高度超过水平 u 的期望数量，得到隐式但可计算的公式。
• 通过变量替换与辅助高斯积分，建立 GOI 矩阵与大小为 N+1 的 GOE 之间的联系，尤其当 c > 0 时。
• 将推导出的公式应用于 ℝ² 和 𝕊² 上的显式情形，验证了该方法的普适性，并在费约多罗夫的限制条件下与已有结果保持一致。

实验结果

研究问题

• RQ1 在平滑各向同性高斯场中，是否可以不施加协方差函数在所有维度 N 上有效的限制，显式计算临界点的期望数量？
• RQ2 如何利用随机矩阵理论，超越 GOE，表征一般各向同性高斯场的海森矩阵？
• RQ3 在 GOI 矩阵表示中，协方差参数 c 的作用是什么？它如何影响临界点的分布？
• RQ4 结果如何在先前限制下恢复或推广费约多罗夫的临界点期望公式？
• RQ5 在参数空间边界处，临界点分布的行为如何，特别是对于随机拉普拉斯特征函数？

主要发现

• 指数 i 且高度超过水平 u 的临界点的期望数量由仅依赖于协方差函数在零点的一阶与二阶导数的可计算公式给出，适用于所有平滑各向同性场。
• 场的海森矩阵以 GOI 矩阵表示，其推广了 GOE，使得所有各向同性场（包括费约多罗夫限制所排除的场）均可被处理。
• 协方差结构的参数空间被扩展至超越费约多罗夫条件的范围：对于 ℝᴺ，其范围为 0 < κ² ≤ (N+2)/N；对于 𝕊ᴺ，其范围为 κ² – η² ≤ (N+2)/N。
• 随机拉普拉斯特征函数的特殊情况对应于参数空间的边界（在 ℝᴺ 上为 κ² = (N+2)/N，在 𝕊ᴺ 上为 κ² – η² = (N+2)/N），此时场与其海森矩阵之间呈现退化相关性。
• 当 c > 0 时，临界点期望数量的 Kac-Rice 积分可改写为对 GOE(N+1) 的期望形式，从而将费约多罗夫的结果作为特例恢复。
• 推导并计算了 ℝ² 和 𝕊² 上各向同性高斯场的显式公式，验证了该方法在低维情形下的适用性与正确性。