[論文レビュー] Computing a k-sparse n-length Discrete Fourier Transform using at most 4k samples and O(k log k) complexity

Sameer Pawar, Kannan Ramchandran|arXiv (Cornell University)|May 4, 2013

Sparse and Compressive Sensing Techniques参考文献 21被引用数 7

ひとこと要約

本稿では、中国剰余定理（CRT）を用いた均等なサブサンプリングによりアリアイジングを活用することで、O(k)のサンプルとO(k log k)の演算でkスパースな離散フーリエ変換（DFT）を計算する、FFASTアルゴリズムを提案する。この手法は、スパースグラフ符号の復号問題に変換することで、実用的なスパース度範囲において定数が4k未満の近似的に最適なサンプルおよび計算複雑度を達成する。

ABSTRACT

Given an $n$-length input signal $\mbf{x}$, it is well known that its Discrete Fourier Transform (DFT), $\mbf{X}$, can be computed in $O(n \log n)$ complexity using a Fast Fourier Transform (FFT). If the spectrum $\mbf{X}$ is exactly $k$-sparse (where $k<

研究の動機と目的

• k << n のとき、kスパース信号のDFTを計算するためのサンプル数と計算複雑度を低減すること。
• サブサンプリングによるアリアイジングを歪みではなく特徴として活用し、優れた消去訂正符号のパリティ制約に類似させる。
• kスパースDFTに対してO(k log k)の計算複雑度とO(k)のサンプル複雑度を達成し、定数を小さくすること。
• スペクトルがスパースでサポートが一様ランダムな場合に、従来のFFTを上回る実用的で低複雑度のアルゴリズムを設計すること。
• 密度推定を用いた符号理論の技術を用いて、CRTに基づくアリアイジングフレームワークとランダムなスパースグラフ符号を理論的に結びつけ、性能保証を示すこと。

提案手法

• FFASTアルゴリズムは、中国剰余定理（CRT）を各段階で用いて構造化されたアリアイジングパターンを生成する、均等なサブサンプリングを複数段階で行う。
• 非常にスパースな領域（k = O(n^δ), δ ≤ 1/3）では、互いに素な整数を用いてサブサンプリングパターンを設計し、ややスパースでない領域（δ > 1/3）では、循環的にシフトされた互いに素な整数を用いる。
• 得られるDFTドメインにおけるアリアイジングは、低密度パリティーチェック（LDPC）符号およびフォンテーヌン符号のパリティーチェック制約に類似しており、スパースDFTの復元を復号問題として扱える。
• FFASTのバックエンドでは、O(k)のサブサンプリング測定値からスパースDFTを回復するために、ピーリング型の復号アルゴリズムが適用される。
• 理論的分析では、符号理論における密度推定技術を用いて復号プロセスをモデル化し、一様ランダムなサポートの仮定下での性能を予測する。
• シミュレーションにより、理論的予測と非常に一致する結果が得られ、特にδ < 0.99の場合に顕著である。

実験結果

リサーチクエスチョン

• RQ1 nに比べてkが部分線形である場合、n未満のサンプルと線形時間未満の計算複雑度でkスパースDFTを計算可能か？
• RQ2 サブサンプリングによるアリアイジングを、DFTドメインに符号に類似した構造を生み出すために体系的に活用可能か？
• RQ3 構造化されたサブサンプリングを用いて、高い確率でkスパースDFTを復元するために必要な最小サンプル数は何か？
• RQ4 中国剰余定理を用いて、スパースDFT復元に有利なアリアイジングを誘発するサブサンプリングパターンをどのように設計できるか？
• RQ5 このようなフレームワークの理論的性能限界は何か？ また、情報理論的下界にどれほど近づけるか？

主な発見

• δ < 0.99 のとき、k = O(n^δ) に対して、FFASTは4k未満のサンプルでkスパースDFTを計算し、O(k)のサンプル複雑度の定数が小さい。
• FFASTの計算複雑度はO(k log k)であり、スパース信号に対して標準的なFFTのO(n log n)と比べて顕著に低い。
• 非ゼロ係数のサポートが一様ランダムである場合、kスパースDFTの高確率復元が達成される。
• 密度推定を用いた理論的分析により、ピーリング復号器が高確率で収束し、復号失敗確率がkとともに指数関数的に減少することが示された。
• シミュレーションにより理論的予測が確認され、さまざまなスパース度において解析的性能と実測性能の間にきわめて高い一致が得られた。
• このフレームワークは、スパースDFT計算とスパースグラフ符号の復号の間に強い関係を確立し、低複雑度のメッセージパッシングによる効率的回復を可能にした。