[論文レビュー] Numerical Analysis for Iterative Filtering with New Efficient Implementations Based on FFT

Antonio Cicone, Haomin Zhou|arXiv (Cornell University)|Feb 5, 2018

Image and Signal Denoising Methods参考文献 5被引用数 13

ひとこと要約

本稿では、非定常的かつ非線形信号を固有モード関数（IMF）に分解するための反復フィルタリング（IF）法について、包括的な数値解析を提示する。高速フーリエ変換（FFT）に基づく実装を活用することで、収束を加速し、計算複雑度を低減する直接計算手法を導出し、信号の即時的分解を可能にした。実際の応用例では、従来のIF法が数日を要する信号（数千万点規模）を、1時間未塔で処理可能であることが示された。

ABSTRACT

Real life signals are in general non--stationary and non--linear. The development of methods able to extract their hidden features in a fast and reliable way is of high importance in many research fields. In this work we tackle the problem of further analyzing the convergence of the Iterative Filtering method both in a continuous and a discrete setting in order to provide a comprehensive analysis of its behavior. Based on these results we provide new ideas for efficient implementations of Iterative Filtering algorithm which are based on Fast Fourier Transform (FFT), and the reduction of the original iterative algorithm to a direct method.

研究の動機と目的

• 連続的および離散的設定における反復フィルタリング（IF）法の収束性を厳密に分析すること。
• IFにおける未解決の問題、特に停止基準の影響と、全IMF分解における外側ループの収束性を解消すること。
• 計算コストを著しく低減しながらも精度を保持する、効率的なFFTベースの実装を開発すること。
• 巡回行列の性質とDFT/IDFT演算を用いて、IMFの直接的・1ステップ計算を確立すること。
• 反復的IFを高速で安定的かつスケーラブルなアルゴリズムに変換することで、実用的かつ大規模な信号解析を可能にすること。

提案手法

• 内側ループおよび外側ループに焦点を当てた、連続的および離散的設定におけるIFアルゴリズムの形式的収束解析。
• 離散逆フーリエ変換（DIF）を用いたIMFの明示的公式の導出：IMF = IDFT((I - D)^N₀ DFT(s))。
• O(n log n) 時間でDFTおよびIDFT演算をFFTを用いて計算することで、反復的畳み込みの置き換えを実現。
• さまざまなフィルタスケーリングにおける固有値λₖの事前計算により、アルゴリズムのさらなる高速化を実現。
• スペクトル分解を用いて1ステップでIMFを計算する、高速反復フィルタリング（FIF）法の導入。
• フィルタおよび信号の性質に基づいた反復回数N₀の事前推定を実装し、不要な反復を回避。

実験結果

リサーチクエスチョン

• RQ1 反復フィルタリングにおける内側ループの収束性および精度に、停止基準の選択がどのように影響するか？
• RQ2 IFアルゴリズムにおける外側ループの有限収束を保証する条件は何か？ また、有限個のIMFの後で分解は終了するか？
• RQ3 巡回行列のスペクトル的性質を用いて、反復的IFプロセスを直接的かつ反復のない方法に変換可能か？
• RQ4 FFTベースの技術をどのように活用すれば、大規模信号におけるIFの計算複雑度を著しく低減できるか？
• RQ5 異なるフィルタ関数がスペクトル分解および得られるIMFの品質に与える影響は何か？

主な発見

• IFアルゴリズムは連続的および離散的設定の両方で、偽の振動を発生させることなく、有限個のIMFに収束する。
• DIF法により、DFT係数の線形結合として各IMFを明示的に直接計算する公式が得られ、1ステップ計算が可能となる。
• 高速反復フィルタリング（FIF）アルゴリズムにより、数千万点規模の信号処理において、従来の数日間を要する処理を1時間未塔に短縮した。
• フィルタ固有値λₖの事前計算により、繰り返し処理やバッチ処理においてFIFアルゴリズムのさらなる高速化が可能となった。
• 信号の摂動に対しても安定しており、連続的および離散的定式化の両方で偽の振動が生じないことが証明された。
• オンラインで利用可能なFIF実装は、FFTを用いて大規模な非定常信号のほぼ即時の分解を達成している。