[論文レビュー] A stability barrier for reconstructions from Fourier samples

Ben Adcock, Anders C. Hansen+1 more|arXiv (Cornell University)|Oct 29, 2012

Mathematical functions and polynomials参考文献 44被引用数 4

ひとこと要約

本稿は、解析的かつ非周期的関数を最初の $ m $ 個のフーリエ係数から再構成する際の根本的な安定性の限界を確立する。任意の安定な手法は、$ m $ に関してせいぜいルート指数的収束しか達成できず、指数的収束は指数的悪条件性を意味する。著者らは、フーリエ級数が有界な多項式の最大挙動に関する新しい推定を用いてこの結果を証明し、フーリエ拡張が機械精度まで到達する安定的かつ超代数的収束を実現する代替手法であることを示している。

ABSTRACT

We prove that any stable method for resolving the Gibbs phenomenon - that is, recovering high-order accuracy from the first $m$ Fourier coefficients of an analytic and nonperiodic function - can converge at best root-exponentially fast in $m$. Any method with faster convergence must also be unstable, and in particular, exponential convergence implies exponential ill-conditioning. This result is analogous to a recent theorem of Platte, Trefethen & Kuijlaars concerning recovery from pointwise function values on an equispaced $m$-grid. The main step in our proof is an estimate for the maximal behaviour of a polynomial of degree $n$ with bounded $m$-term Fourier series, which is related to a conjecture of Hrycak & Groechenig. In the second part of the paper we discuss the implications of our main theorem to polynomial-based interpolation and least-squares approaches for overcoming the Gibbs phenomenon. Finally, we consider the use of so-called Fourier extensions as an attractive alternative for this problem. We present numerical results demonstrating rapid convergence in a stable manner.

研究の動機と目的

• フーリエ係数から生じるギブス現象を解消する安定な手法の収束速度に根本的な限界を確立すること。
• 既存の幾何的収束を示す手法がなぜ悪条件性を示すのか、および安定な手法が通常はゆっくり収束する理由を説明すること。
• 多項式に基づく手法（例：逆多項式再構成や最小二乗法）の理論的限界を分析すること。
• フーリエ拡張を提案し、解析的関数をフーリエデータから再構成するための安定的かつ高次精度の代替手法として評価すること。

提案手法

• 次数 $ n $ の多項式が有界な $ m $ 項フーリエ級数を持つ場合の最大挙動に関する新しい推定を導出し、コッパースミスとリヴリンの結果をフーリエ設定に拡張する。
• この推定を用いて、収束速度が $ \rho^{-m^\tau} $（$ \rho > 1 $、$ \tau \in (\frac{1}{2}, 1] $）である任意の手法は、$ \rho^{m^{2\tau - 1}} $ のオーダーの悪条件性を持つ必要があることを証明し、ルート指数的障壁を確立する。
• 逆多項式再構成法（IPRM）を分析し、指数的悪条件性であることを示し、最小二乗法による安定化により収束速度がルート指数的へ低下することを示す。
• フーリエ拡張を安定な代替手法として提案する。ここで近似は、関数をより大きな区間に拡張し、そのフーリエ級数を計算することで構築される。
• 数値実験を用いて、さまざまな特異性を有するテスト関数に対して、多項式最小二乗法とフーリエ拡張を比較する。
• 拡張区間長を制御するパラメータ $ T $ を用い、フーリエ拡張法における条件数と収束速度のバランスをとる。

実験結果

リサーチクエスチョン

• RQ1 解析的かつ非周期的関数を最初の $ m $ 個のフーリエ係数から再構成する安定な手法が達成可能な最良の収束速度は何か？
• RQ2 実際の応用では幾何的収束を示す手法がなぜ重度の悪条件性を示すのか？
• RQ3 多項式に基づく最小二乗法は、安定性を保ちつつルート指数的収束を上回る収束を達成できるか？
• RQ4 フーリエ拡張は、機械精度まで収束することで、安定性の障壁をどれほど回避できるか？
• RQ5 境界層、実または複素の特異性を有する関数などの異なる関数タイプは、再構成手法の性能にどのように影響を与えるか？

主な発見

• 解析的かつ非周期的関数を最初の $ m $ 個のフーリエ係数から再構成する安定な手法は、$ m $ に関してせいぜいルート指数的収束しか達成できず、すなわちある $ \rho > 1 $ に対して $ \rho^{-m^{1/2}} $ のオーダーとなる。
• 指数的収束は必然的に指数的悪条件性を意味するため、このような手法は数値的に不安定である。
• 逆多項式再構成法（IPRM）は指数的悪条件性であることが示され、高精度な用途には不適切である。
• IPRM を最小二乗法で安定化すると、収束速度はルート指数的へ低下し、理論的な安定性の限界と一致する。
• フーリエ拡張は、超代数的（しばしば幾何的）収束を機械精度まで実現し、複雑な特異性や高周波数振動を有する関数に対して、多項式最小二乗法を上回る性能を示す。
• 実際には、拡張パラメータ $ T $ の選択が誤差に与える影響は最小限に抑えられる。なぜなら、自由度の増加が各自由度あたりの収束速度の低下と相殺されるからである。