[論文レビュー] Adaptive sub-linear Fourier algorithms

Debbie A. Lawlor, Yan Wang+1 more|arXiv (Cornell University)|Jul 26, 2012

Sparse and Compressive Sensing Techniques参考文献 25被引用数 15

ひとこと要約

本稿では、$k$ 個の顕著な周波数の数に対して線形スケーリングを達成する決定的で適応的な非線形フーリエ法を提示する。従来の決定的手法の $O(k^2)$ の複雑さに比べて優れた性能を発揮する。時間遅延付DFTにおける位相関係を活用することでアリエイジョンを検出し、アリエイジョンのない周波数を回復可能となる。平均ケースにおいて、既存手法に比べて桁違いの高速化を実現する。

ABSTRACT

We present a new deterministic algorithm for the sparse Fourier transform problem, in which we seek to identify k << N significant Fourier coefficients from a signal of bandwidth N. Previous deterministic algorithms exhibit quadratic runtime scaling, while our algorithm scales linearly with k in the average case. Underlying our algorithm are a few simple observations relating the Fourier coefficients of time-shifted samples to unshifted samples of the input function. This allows us to detect when aliasing between two or more frequencies has occurred, as well as to determine the value of unaliased frequencies. We show that empirically our algorithm is orders of magnitude faster than competing algorithms.

研究の動機と目的

• 標準的なFFTの $\Omega(N)$ の壁を打ち破ることを目的とした、$k \ll N$ の周波数係数のみが顕著なスパースフーリエ変換問題に対処する。
• 既存の非線形手法に比べ、平均ケースの実行時間とサンプリング複雑度を向上させた決定的アルゴリズムを開発する。
• 増加するシフトサイズ $\varepsilon_j$ を用いたマルチスケール手法を導入することで、ノイジーな環境下でも頑健な周波数推定を可能にする。
• 非整数周波数への適用可能性を拡張し、本研究の文脈において画期的な貢献を果たす。
• 周波数局在化と係数忠実度の両面を考慮した、カスタムコスト関数を用いた地球移動距離（EMD）メトリックを用いて近似品質を評価する。

提案手法

• 時間遅延付きサンプルを活用し、元の信号のフーリエ係数と遅延信号の係数との関係を確立することで、アリエイジョンの検出とアリエイジョンのない周波数の回復を可能にする。
• 互いに素な長さの短いDFTに対して中国剰余定理を適用し、異なる長さにおける剰余から主要周波数を再構成する。
• 複数のシフトにおける位相の一貫性に基づく丸め手順を適用し、周波数成分を段階的に高精度で推定する。
• シフトサイズ $\varepsilon_j$ を段階的に増加させるマルチスケール戦略を導入し、上位ビット（MSB）から下位ビット（LSB）へと周波数成分を1ビットずつ抽出する。
• 周波数と振幅の両方の精度をバランスさせるために、コスト関数 $d_1(\omega, \tilde{\omega}) = \frac{|\omega - \tilde{\omega}|}{N} + |c_\omega - c_{\tilde{\omega}}|$ を用いた地球移動距離（EMD）を用いて近似品質を評価する。
• 正しさを保証するがランダムな実行時間を持つラスベガス方式のアルゴリズムを設計し、適応的サンプリングと位相ベースのフィルタリングにより性能を最適化する。

実験結果

リサーチクエスチョン

• RQ1 スパースフーリエ変換問題において、$k$ に関して $O(k^2)$ の既存の境界を上回る、決定的アルゴリズムが $k$ に関して線形スケーリングの実行時間で達成可能かどうか。
• RQ2 時間遅延付DFTからの位相情報は、どのようにしてアリエイジョンを検出し、アリエイジョンのない周波数を回復するのに活用できるか。
• RQ3 増加するシフトサイズ $\varepsilon_j$ を用いたマルチスケール手法は、ノイジーな状況下での周波数推定精度をどの程度向上できるか。
• RQ4 周波数と振幅のハイブリッドコスト関数を用いた地球移動距離（EMD）メトリックは、$\ell_2$ ベースの誤差メトリックと比較して、スパースフーリエ近似品質の評価においてどのように異なるか。
• RQ5 非整数周波数を処理できるように、このアルゴリズムを拡張可能か。

主な発見

• 提案手法は平均ケースにおいて $k$ に関して線形スケーリングを達成し、従来の決定的手法の $O(k^2)$ 複雑度に比べて顕著な改善を示す。
• 実験的結果では、$k \ll N$ のスパarsityレベルにおいて、FFTと同様の手法（FFTW）に比べて、本手法は桁違いに高速である。
• ノイズなしの状況では、周波数回復に十分なパラメータ $c_1 = 5$ で十分であったが、ノイズがある状況では低EMD(1)誤差を維持するためには $c_1 \approx 256$ 必要であった。
• ノイズレベル $\sigma$ が上昇するにつれて、EMD(1)誤差は線形に増加する傾向を示し、入力ノイズに対して高い耐性を示す。
• 増加するシフトサイズ $\varepsilon_j$ を用いたマルチスケール手法により、1ビットずつ周波数を推定可能となり、ノイジーな環境下での精度が向上する。
• 時間遅延による位相の一貫性を用いたアリエイジョン検出により、従来手法よりも少ないサンプル数でより多くの情報を抽出可能となる。