DE-Sinc approximation for unilateral rapidly decreasing functions and its computational error bound

技術要約：片側急速減少関数に対するDE-Sinc近似

問題提起
Sinc近似法は、実軸の両端で急速に減少する関数に対して非常に効果的である。しかし、片側急速減少関数（$-\infty \to \infty$ の方向において、一方が代数的に減少し、他方が指数関数的に減少する関数）に対しては、標準的な変数変換を用いると、$O(\exp(-c\sqrt{n}))$ の根指数収束（root-exponential convergence）しか得られない。このクラスの関数のために改良された単一指数（SE）変換も存在するが、この漸近的な限界を克服するには至っていない。さらに、二重指数（DE）変換は他のケース（有限区間、対称な減衰など）には成功裏に適用され、ほぼ指数的な収束を実現しているが、片側急速減少関数（Stengerの分類におけるケース5）に対するDE-Sinc近似の厳密かつ計算可能な誤差界は、これまで確立されていなかった。

手法
本論文では、特定の二重指数変数変換に基づくDE-Sinc近似を提案する：
$$t = \phi_5(x) = 2 \sinh(\log(\log(1 + e^{\pi \sinh x})))$$
この変換は以前、数値積分のために導入されたものであるが、ここでは関数近似のために適応されている。手法は以下の通りである：

1. 変換： 元の関数 $f(t)$ を、実軸上で定義された変換関数 $F(x) = f(\phi_5(x))$ へ写像する。
2. 近似： $F(x)$ に対して標準的なSinc級数展開を適用し、その結果を $t$ ドメインへ写像する。
3. 誤差解析： 全誤差を離散化誤差と打ち切り誤差に分解する。
   • 離散化誤差は、ストリップ領域 $D_d$ 上の Hardy 空間 $H_1(D_d)$ における $F$ のノルムを用いて抑えられる。
   • 打ち切り誤差は、$F(x)$ の $\pm \infty$ における減衰特性を分析することで抑えられる。
4. 技術的課題： 証明には、$\phi_5$ の入れ子になった対数構造、複素平面におけるその分岐構造、および実軸の両端における非対称な減衰挙動による、新しい推定が必要となる。著者らは、定数を厳密に抑えるために、ストリップ領域上の変換関数の挙動に関する明示的な不等式を導出している。

主な貢献

• 新定理（定理 2.3）： 本論文は、片側急速減少関数のDE-Sinc近似に関する厳密な誤差界を確立した。この界は以下の形式をとる：
  $$\sup_{t \in \mathbb{R}} |f(t) - \text{Approx}(t)| \leq C e^{-\pi d n / \log(2dn/\mu)}$$
  ここで、$n$ はサンプリング点数、$d$ はストリップ幅、$\mu$ は減衰率に関連する。
• 明示的な定数： 非計算的な定数に依存する場合がある従来の理論的結果とは異なり、本論文は定数 $C$ を問題のパラメータ（$K_+, K_-, \alpha, \beta, d$）を用いて明示的な公式として提供している。これにより、この誤差界は計算可能であり、保証された精度を必要とするアルゴリズムに適したものとなっている。
• 収束速度： 本手法は、$O(\exp(-cn/\log n))$ のオーダーの**ほぼ指数的な収束（almost exponential convergence）**を達成する。これは、このクラスの関数に対する既存のSE-Sinc法の根指数収束 $O(\exp(-c\sqrt{n}))$ よりも大幅な改善である。
• サンプル複雑性： 精度 $\epsilon$ を達成するために、提案されたDE-Sinc法は $n = O(\log(1/\epsilon) \log \log(1/\epsilon))$ のサンプリング点を必要とする。これに対し、SE-Sinc近似では $n = O((\log(1/\epsilon))^2)$ が必要となる。

結果
本論文は、3つの特定の片側急速減少関数（$f_1, f_2, f_3$）を用いた数値例を通じて、理論的知見を検証している：

• 新定理の解析性の仮定を満たす関数 $f_1$ および $f_2$ について、観察された誤差は、SE-Sinc近似（定理 2.1 および 2.2）よりも、DE-Sinc近似（定理 2.3）において著しく速く収束した。
• これらの関数における観察された誤差は、導出された理論的誤差界の範囲内に収まっており、定数 $C$ と収束次数の妥当性を裏付けている。
• 新定理の厳密な解析性の仮定を満たさない関数 $f_3$ についても、DE-Sinc近似は競争力のある性能を示したが、特徴的なほぼ指数的な収束は示さず、理論的誤差界は適用できなかった。

意義
本論文は、片側急速減少関数のためのDE-Sinc近似に関する初の計算可能な誤差界を提供することにより、Sinc数値計算法における重要な空白を埋めたと主張している。主な意義は、この特定のクラスの問題に対して保証された精度での計算を可能にした点にある。ほぼ指数的な収束を確立することで、本手法は、既存の単一指数アプローチと比較して、高い精度を達成するために必要な計算コスト（サンプリング点数）を大幅に削減する。著者らは、Sinc近似は求積法やコロケーション法の基礎的な構成要素であるため、これらの結果が、片側減衰を伴う微分方程式や積分方程式に関する将来の保証精度アルゴリズムの解析的基盤となる可能性があると述べている。