[论文解读] On the size-distribution of Poisson Voronoi cells

F. Jarai-Szabo, Zoltán Néda|ArXiv.org|Jun 4, 2004

Computational Geometry and Mesh Generation被引用 8

一句话总结

本文提出了一种用于二维和三维泊松Voronoi胞的归一化大小分布的紧凑解析近似，采用具有单参数 $ a = (3d+1)/2 $ 的广义伽马形式，其中 $ d $ 为维度。公式 $ f_d(y) \propto y^{(3d-1)/2} \exp\left(-\frac{3d+1}{2}y\right) $ 提供了一个简单且可解析处理的模型，其与大规模蒙特卡洛模拟结果高度吻合，尽管其精度低于多参数拟合，但为实验工作者提供了一项实用工具。

ABSTRACT

Poisson Voronoi diagrams are useful for modeling and describing various natural patterns and for generating random lattices. Although this particular space tessellation is intensively studied by mathematicians, in two- and three dimensional spaces there is no exact result known for the size-distribution of Voronoi cells. Motivated by the simple form of the distribution function in the one-dimensional case, a simple and compact analytical formula is proposed for approximating the Voronoi cell's size distribution function in the practically important two- and three dimensional cases as well. Denoting the dimensionality of the space by d (d=1,2,3) the $f(y)=Const*y^{(3d-1)/2}exp(-(3d+1)y/2)$ compact form is suggested for the normalized cell-size distribution function. By using large-scale computer simulations the validity of the proposed distribution function is studied and critically discussed.

研究动机与目标

• 解决在二维和三维中泊松Voronoi胞大小分布函数缺乏精确或简单解析表达式的问题，尽管其在自然与社会现象建模中被广泛应用。
• 为实验科学家提供一种实用且可解析使用的近似方法，使其能以有限的统计数据表征复杂模式。
• 将一维中已知的精确结果 $ f_1(y) \propto y \exp(-2y) $ 扩展至更高维度，采用统一的函数形式。
• 通过大规模蒙特卡洛模拟验证所提近似，确保其在所有细胞尺寸范围内的准确性。

提出的方法

• 提出形式为 $ f_d(y) = C \cdot y^{(3d-1)/2} \exp\left(-\frac{3d+1}{2}y\right) $ 的归一化大小分布函数，其中 $ y = S/\langle S \rangle $ 为归一化细胞尺寸，$ d $ 为空间维度。
• 利用伽马函数推导归一化常数 $ C $，确保 $ \int_0^\infty f_d(y) dy = 1 $。
• 利用大规模计算机模拟（最多 $ 10^7 $ 个细胞）生成二维和三维的实证细胞大小分布以供验证。
• 将所提公式与三参数伽马函数拟合及模拟数据进行比较，评估其在尾部和峰值区域的精度。
• 解析计算矩（均值、方差、偏度），并与模拟结果和理论值比较，以验证模型。
• 通过尾部分布的双对数图和对数正态图评估该近似在极端尺寸区域的表现。

实验结果

研究问题

• RQ1 是否能够使用一种简单紧凑的解析形式，准确近似二维和三维泊松Voronoi胞的大小分布，尽管目前尚无精确解？
• RQ2 与更复杂的三参数拟合相比，所提出的单参数形式在精度和实际可用性方面表现如何？
• RQ3 函数形式 $ f_d(y) \propto y^{(3d-1)/2} \exp\left(-\frac{3d+1}{2}y\right) $，其中 $ a = (3d+1)/2 $，是否能在所有细胞尺寸范围内良好拟合大规模模拟数据？
• RQ4 所提分布的关键统计矩（均值、方差、最可能尺寸）是什么？它们与模拟和理论值相比如何？

主要发现

• 所提分布 $ f_d(y) \propto y^{(3d-1)/2} \exp\left(-\frac{3d+1}{2}y\right) $，其中 $ a = (3d+1)/2 $，在二维和三维中对归一化细胞大小分布提供了良好近似，对应 $ a = 7/2 $ 和 $ a = 5 $。
• 在三维中，最可能的归一化体积为 $ y_{\text{prob}} = 4/5 = 0.8 $，二阶矩为 $ \langle y^2 \rangle = 6/5 = 1.2 $，与理论值 1.18 的相对误差为 1.7%。
• 由模拟数据计算的二阶矩为 $ \langle y^2 \rangle_{\text{sim}} = 1.19 $，相对误差为 0.85%，表明与模型高度一致。
• 由所提公式导出的累积分布函数在图示尺度下与模拟数据无显著差异，表明整体拟合效果极佳。
• 该单参数形式在简洁性上优于简单伽马拟合，在实际应用中（尤其在数据有限时）精度也具有竞争力。
• 该模型并非精确解，其精度低于三参数伽马拟合，但其简洁性使其对需要第一性原理近似的实验工作者极具价值。