[论文解读] Digital Nets and Sequences for Quasi-Monte Carlo Methods

Hee Sun Hong|arXiv (Cornell University)|Jul 27, 2022

Mathematical Approximation and Integration被引用 6

一句话总结

本论文通过引入随机化数字网、分析混淆网的不均匀性分布以及利用进化计算优化数字网质量，推进了准蒙特卡罗方法的发展。研究结果表明，混淆网可提升误差估计与性能表现，而经进化计算优化的网在高维积分中优于 Sobol’ 和 Niederreiter-Xing 序列，尤其在复杂测试函数下表现更优。

ABSTRACT

Quasi-Monte Carlo methods are a way of improving the efficiency of Monte Carlo methods. Digital nets and sequences are one of the low discrepancy point sets used in quasi-Monte Carlo methods. This thesis presents the three new results pertaining to digital nets and sequences: implementing randomized digital nets, finding the distribution of the discrepancy of scrambled digital nets, and obtaining better quality of digital nets through evolutionary computation. Finally, applications of scrambled and non-scrambled digital nets are provided.

研究动机与目标

• 开发随机化数字 $(t,s)$-序列的高效实现方法，以提升准蒙特卡罗积分的性能。
• 分析混淆数字网中平方不均匀性的统计分布，以实现严谨的误差估计。
• 利用进化计算优化数字网的 $t$-参数，以在高维空间中实现更优的均匀分布。
• 在基准多维积分问题上评估新旧数字网的性能表现。
• 在多种测试函数下对比混淆与非混淆数字网，以评估其鲁棒性与收敛速率。

提出的方法

• 通过将生成矩阵与随机矩阵相乘并应用数字平移，实现 Owen 的混淆方法与 Faure-Tezuka 的翻滚方法。
• 利用二进制算术与格雷码生成技术，加速基于 2 的序列计算，通过利用硬件效率降低时间复杂度。
• 将平方不均匀性经验分布拟合为卡方分布的混合模型，依据中心极限定理与 Loh 的理论框架。
• 应用进化计算搜索最优生成矩阵，以最小化数字 $(t,m,s)$-网的 $t$-参数。
• 基于多投影下的 $t$-值定义目标函数，并采用带有自适应环境设置的随机搜索策略。
• 采用拉丁超立方采样与 Genz 的测试函数套件，评估在振荡、尖峰、高斯及不连续被积函数上的性能表现。

实验结果

研究问题

• RQ1 如何高效实现随机化数字网，以在保持低不均匀性的同时实现方差估计？
• RQ2 混淆数字网中平方不均匀性的真实分布为何？其能否被理论分布良好拟合？
• RQ3 进化计算能否有效优化数字网的 $t$-参数，使其在性能上超越 Sobol’ 与 Niederreiter-Xing 等经典构造？
• RQ4 在多种高维积分问题中，混淆与非混淆数字网在收敛性与精度方面有何差异？
• RQ5 随着维度增加，新构造数字网的均匀分布特性（尤其是低维投影中的表现）改善程度如何？

主要发现

• 混淆数字网，特别是 Owen 混淆与 Faure-Tezuka 翻滚变体，在平均平方不均匀性上的方差显著低于非混淆版本。
• 混淆网中平方不均匀性的经验分布可良好拟合为卡方分布的混合模型，验证了理论近似，并支持统计误差估计。
• 在高维投影中，经进化的数字网（EC）的 $t$-参数低于 Sobol’ 与 Niederreiter-Xing 序列，尤其在 $s = 18$ 与 $s = 32$ 时表现更优，表明其具有更优的均匀分布特性。
• 在多维积分中，混淆 EC 网与混淆 Niederreiter-Xing 网持续优于非混淆版本及其他序列，且在振荡与不连续函数上表现出最佳的数字精度。
• Niederreiter-Xing 序列在低维下表现良好，但在高维中性能下降；而嵌入式 Niederreiter-Xing 变体表现较差，因其未针对嵌入设计。
• 生成混淆序列的开销约为非混淆序列的两倍，但其在收敛性与误差控制方面的性能提升足以证明该开销的合理性。