[論文レビュー] Bounds on Walsh coefficients by dyadic difference and a new Koksma-Hlawka type inequality for Quasi-Monte Carlo integration

Takehito Yoshiki|arXiv (Cornell University)|Apr 10, 2015

Mathematical Approximation and Integration被引用数 6

ひとこと要約

本稿は、二進差分作用素を用いてウォルシュ係数の tighter な上限を導出し、Quasi-Monte Carlo 統合における Koksma-Hlawka 型不等式を新たに導入する。関数の混合偏導関数が順序 α まで連続である場合、従来の推定値と比較して定数を著しく削減でき、特に不確実性評価やデジタルネットの CBC 構成において顕著な利点を示す。

ABSTRACT

In this paper we give a new Koksma-Hlawka type inequality for Quasi-Monte Carlo (QMC) integration. QMC integration of a function $f\colon[0,1)^s ightarrow \mathbb{R}$ by a finite point set $\mathcal{P}\subset [0,1)^s$ is the approximation of the integral $I(f):=\int_{[0,1)^s}f(\mathbf{x})\,d\mathbf{x}$ by the average $I_{\mathcal{P}}(f):=\frac{1}{|\mathcal{P}|}\sum_{\mathbf{x} \in \mathcal{P}}f(\mathbf{x})$. We treat a certain class of point sets $\mathcal{P}$ called digital nets. A Koksma-Hlawka type inequality is an inequality bounding the integration error $ ext{Err}(f;\mathcal{P}):=I(f)-I_{\mathcal{P}}(f)$ by a bound of the form $| ext{Err}(f;\mathcal{P})|\le C\cdot \|f\|\cdot D(\mathcal{P})$. We can obtain a Koksma-Hlawka type inequality by estimating bounds on $|\hat{f}(\mathbf{k})|$, where $\hat{f}(\mathbf{k})$ is a generalized Fourier coefficient with respect to the Walsh system. In this paper we prove bounds on Walsh coefficients $\hat{f}(\mathbf{k})$ by introducing an operator called `dyadic difference' $\partial_{i,n}$. By converting dyadic differences $\partial_{i,n}$ to derivatives $\frac{\partial }{\partial x_i}$, we get a new bound on $|\hat{f}(\mathbf{k})|$ for a function $f$ whose mixed partial derivatives up to order $α$ in each variable are continuous. This new bound is smaller than the known bound on $|\hat{f}(\mathbf{k})|$ under some condition. The new Koksma-Hlawka inequality is derived using this new bound on the Walsh coefficients.

研究の動機と目的

• Quasi-Monte Carlo 統合における既存の Koksma-Hlawka 型不等式の大きな定数を低減すること、特に不確実性評価の文脈において。
• 相互にねじれ型多項式格子ルールの CBC 構成における現実的でないほど大きな定数の問題を解決すること。
• 離散的な二進構造と連続的微分の間の関係を確立する新たな枠組みを構築し、ウォルシュ係数を二進差分作用素を用いて評価すること。
• 一般化されたフーリエ係数の推定を改善することで、基数 b=2 のデジタルネットに対してよりタイトな Koksma-Hlawka 不等式を導出すること。

提案手法

• 関数の離散的二進シフトにおける挙動を分析することで、ウォルシュ係数 |f^(k)| を二進差分作用素 ∂i,n で上界付ける。
• 極限過程を通じて、二進差分 ∂i,n と古典的偏微分 ∂/∂xi 間の関係を確立する。
• 二進差分作用素を用いて、各変数について混合偏導関数が順序 α まで連続である関数 f に対する |f^(k)| の新しい上限を導出する。
• ホルダーの不等式と Lp-ノルム推定を用いて、二進差分カーネルの L1 および L∞ ノルムを制御し、L∞ ノルムが常に 2 である一様な上限を得る。
• 改善されたウォルシュ係数の上限とデジタルネットの双対ネット構造を組み合わせることで、新しい Koksma-Hlawka 型不等式を導出する。
• 多変数インデックス k の二進構造に関する帰納法を用いて、境界を検証し、すべての k ∈ ℕ に対してカーネルの L1 ノルムが 1、L∞ ノルムが 2 であることを証明する。

実験結果

リサーチクエスチョン

• RQ1 古典的微分に基づく手法と比較して、二進差分を用いることでウォルシュ係数のよりタイトな上限を導出可能か？
• RQ2 デジタルネットの既存の Koksma-Hlawka 不等式における定数と、二進差分解析によって達成可能な定数はどのように比較されるか？
• RQ3 二進差分作用素は、滑らかで周期的でない関数に対する QMC 統合の最悪誤差上限をどの程度改善できるか？
• RQ4 ウォルシュ係数に対する新しい上限は、デジタルネット統合における収束速度や定数因子に意味的な改善をもたらすか？
• RQ5 改善された上限は、相互にねじれ型多項式格子ルールの CBC 構成における計算コスト低減に実用的に応用可能か？

主な発見

• 本稿は、任意の k ∈ ℕ に対して、二進差分カーネル W(k) の L1 ノルムが正確に 1、L∞ ノルムが正確に 2 であることを確立しており、最終的な境界において極めて重要である。
• 特定の条件下では、混合偏導関数が順序 α まで連続である場合、既知の境界よりも |f^(k)| の新しい境界が厳密に小さくなる。
• 導出された Koksma-Hlawka 不等式は、特に高次元デジタルネットの文脈において、従来の結果と比較して著しく定数が小さくなっている。
• 古典的微分に基づく推定を、二進差分に基づくアプローチに置き換えることで、二進ネットの構造をよりよく捉えることができた。
• CBC 構成における相互にねじれ型多項式格子ルールに、より小さい定数を使用する理論的根拠を提供し、不確実性評価の応用における問題を解決した。
• 境界は k に依存せず一様に成立し、カーネル W(k) の L∞ ノルムが常に 2 で抑えられ、k とは無関係である。