[论文解读] Low discrepancy constructions in the triangle

Kinjal Basu, Art B. Owen|arXiv (Cornell University)|Mar 11, 2014

Mathematical Approximation and Integration参考文献 16被引用 6

一句话总结

本文提出了两种用于单位三角形上数值积分的新型准蒙特卡罗构造，实现了消失的不均匀性。第一种是可扩展的数字范德科普特型序列，其不均匀性被限制在 $12/\sqrt{N}$ 以内；第二种是将二维整数格点以二次无理数切线的角度旋转后，仅保留位于三角形内的点，实现了最优的 $O(\log N / N)$ 不均匀性。两种构造均能使黎曼可积函数的积分收敛，且无需有界变差条件。

ABSTRACT

Most quasi-Monte Carlo research focuses on sampling from the unit cube. Many problems, especially in computer graphics, are defined via quadrature over the unit triangle. Quasi-Monte Carlo methods for the triangle have been developed by Pillands and Cools (2005) and by Brandolini et al. (2013). This paper presents two QMC constructions in the triangle with a vanishing discrepancy. The first is a version of the van der Corput sequence customized to the unit triangle. It is an extensible digital construction that attains a discrepancy below 12/sqrt(N). The second construction rotates an integer lattice through an angle whose tangent is a quadratic irrational number. It attains a discrepancy of O(log(N)/N) which is the best possible rate. Previous work strongly indicated that such a discrepancy was possible, but no constructions were available. Scrambling the digits of the first construction improves its accuracy for integration of smooth functions. Both constructions also yield convergent estimates for integrands that are Riemann integrable on the triangle without requiring bounded variation.

研究动机与目标

• 在三角形中构造具有消失不均匀性的准蒙特卡罗点集，以实现黎曼可积函数的收敛。
• 将科克斯马-克拉夫卡不等式框架扩展至单纯形（三角形），为三角形积分中的误差界提供理论基础。
• 构造显式、确定性的序列，其不均匀性可被严格证明为低值，从而弥补先前工作中此类构造的缺失。
• 确保构造在计算机图形学和数值积分应用中具有实用性，具备可扩展性与混淆支持。

提出的方法

• 通过将定义域递归地二等分划分为子三角形，将范德科普特序列适配到三角形中。
• 使用以4为基的数字映射到三角形子区域，构建数字序列，确保可扩展性并支持数字混淆。
• 将缩放后的整数格点按具有二次无理数切线的角进行旋转，仅保留位于三角形内的点。
• 利用Chen和Travaglini（2007）以及Brandolini等（2013）的定理，证明旋转格点实现了 $O(\log N / N)$ 的不均匀性。
• 建立Pillards与Cools的不均匀性与Brandolini等的平行四边形不均匀性之间的等价性，表明两者在 $N \to \infty$ 时均趋于消失。
• 采用带符号的不均匀性函数，并通过将区域分解为平行四边形来界定子三角形的不均匀性。

实验结果

研究问题

• RQ1 能否直接在三角形中构造低不均匀性序列，而非通过从单位立方体映射得到？
• RQ2 在三角形中，确定性点集可实现的最优不均匀性速率是多少？
• RQ3 准蒙特卡罗积分能否在不依赖有界变差条件的情况下，对三角形上的黎曼可积函数实现收敛？
• RQ4 能否对基于三角形的QMC序列应用数字混淆，以提升光滑函数的收敛性？
• RQ5 是否存在一种直接构造的类似格点的序列，其在三角形中具有最优的 $O(\log N / N)$ 不均匀性？

主要发现

• 三角形中的数字范德科普特构造实现了至多 $12/\sqrt{N}$ 的平行四边形不均匀性，当 $N = 4^k$ 时，其精确值为 $2/(3\sqrt{N}) - 1/(9N)$。
• 旋转格点构造实现了最优的 $O(\log N / N)$ 不均匀性速率，与三角形中可能达到的最佳速率一致。
• 对数字构造的数字进行混淆，可提升光滑函数的积分精度，当 $N = 4^k$ 时，实现 $O(N^{-1})$ 的均方根误差。
• 两种构造均确保QMC估计对任意三角形上的黎曼可积函数收敛到真实积分，即使在无有界变差条件下亦成立。
• Pillards与Cools度量的不均匀性至多为Brandolini等的平行四边形不均匀性的两倍，因此两者在 $N \to \infty$ 时均趋于消失。
• 理论框架支持对任意非退化三角形的收敛性，而不仅限于标准单位三角形。