[论文解读] A reduced fast component-by-component construction of lattice point sets with small weighted star discrepancy

Ralph Kritzinger, Helene Laimer|arXiv (Cornell University)|Jan 28, 2015

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

一句话总结

本文提出了一种用于构建具有小加权星不等性偏差的格点集的简化快速分量-分量（CBC）算法，当权重迅速衰减时，显著降低了构造成本。通过利用结构化的生成向量并优化分量选择，该方法在权重和移位参数满足可 summability 条件时实现了强多项式可 tractability，从而在高维情况下实现了高效的准蒙特卡洛积分。

ABSTRACT

The weighted star discrepancy of point sets appears in the weighted Koksma-Hlawka inequality and thus is a measure for the quality of point sets with respect to their performance in quasi-Monte Carlo algorithms. A special choice of point sets are lattice point sets whose generating vector can be obtained one component at a time such that the resulting lattice point set has a small weighted star discrepancy. In this paper we consider a reduced fast component-by-component algorithm which significantly reduces the construction cost for such generating vectors provided that the weights decrease fast enough.

研究动机与目标

• 降低构造具有小加权星不等性偏差的格点集的计算成本，以实现准蒙特卡洛积分。
• 通过利用结构化的生成向量，开发一种在维度和样本大小上均能高效扩展的分量-分量算法。
• 建立在何种条件下，所得的格点规则在高维积分问题中实现强多项式可 tractability。
• 通过减少所需运算次数，在保持最优偏差界的同时，改进现有CBC构造方法。

提出的方法

• 该方法使用一种简化快速CBC算法，按分量逐个构造生成向量，仅基于权重序列和移位参数关注相关分量。
• 它引入了基于 $ b $ 进制数字和移位参数 $ w_j $ 的生成向量变换，从而更高效地计算加权星不等性偏差。
• 该算法利用恒等式 $ R_{N,\boldsymbol{\nu}}^s(\textbf{z}) = \frac{1}{N} \sum_{k=0}^{N-1} \eta_s(k) $ 高效计算与偏差相关的项。
• 它应用调和和的界 $ S_N \leq 4\log N $ 来控制偏差估计中的误差项。
• 构造成本降低至 $ O(N\log N + \min\{s,t\}N + N\sum_{d=1}^{\min\{s,t\}} (m - w_d)b^{-w_d}) $，远低于标准CBC方法的 $ O(sN\log N) $ 成本。
• 它使用乘积权重 $ \gamma_{\mathfrak{u}} = \prod_{j\in\mathfrak{u}} \gamma_j $，并假设 $ \gamma_j \leq 1 $，非增，并且 $ \sum_j \gamma_j b^{w_j} < \infty $。

实验结果

研究问题

• RQ1 能否在保持低偏差的同时降低构造具有小加权星不等性偏差的格点集的计算成本？
• RQ2 在何种权重和移位参数条件下，简化CBC算法能实现强多项式可 tractability？
• RQ3 与标准CBC算法相比，简化CBC算法在计算复杂度方面的效率如何？
• RQ4 移位参数 $ w_j $ 对偏差和构造成本有何影响？
• RQ5 能否通过变换后的生成向量和调和和估计，紧密控制偏差界？

主要发现

• 简化快速CBC算法将构造成本降低至 $ O(N\log N + \min\{s,t\}N + N\sum_{d=1}^{\min\{s,t\}} (m - w_d)b^{-w_d}) $，远低于标准CBC方法的 $ O(sN\log N) $。
• 在条件 $ \sum_{j=1}^\infty \gamma_j b^{w_j} < \infty $ 下，该方法实现了强多项式可 tractability，且对任意 $ \delta > 0 $，有 $ R_{N,\boldsymbol{\gamma}}^s(\mathbf{z}) \leq c_{\boldsymbol{\gamma},\delta} N^{\delta - 1} $。
• 加权星不等性偏差满足 $ D_{N,\boldsymbol{\gamma}}^*(\mathbf{z}) \leq O(1/N) + \frac{1}{2} R_{N,\boldsymbol{\gamma}}^s(\mathbf{z}) $，在可求和条件下确保收敛速度与 $ 1/N $ 成正比。
• 当 $ N \geq (c_{\boldsymbol{\gamma},\delta} \varepsilon^{-1})^{1/(1 - \delta)} $ 时，该方法保证 $ R_{N,\boldsymbol{\gamma}}^s(\mathbf{z}) \leq \varepsilon $，从而确认了强多项式可 tractability。
• 不等式 $ \prod_{j=1}^s (1 + 13\gamma_j b^{w_j} \log N) \leq (1 + \sigma_d^{-1})^d N^{(\sigma_0 + 1)\sigma_d} $ 使得在可求和权重下推导出可 tractability 结果成为可能。
• 理论框架证实，简化CBC算法在大幅降低计算量的同时，仍能保持最优偏差界，尤其在权重快速衰减时优势显著。