Unified commutation-pruning technique for efficient computation of composite DFTs

An efficient computation of a composite length discrete Fourier transform (DFT), as well as a fast Fourier transform (FFT) of both time and space data sequences in uncertain (non-sparse or sparse) computational scenarios, requires specific processing algorithms. Traditional algorithms typically employ some pruning methods without any commutations, which prevents them from attaining the potential computational efficiency. In this paper, we propose an alternative unified approach with automatic commutations between three computational modalities aimed at efficient computations of the pruned DFTs adapted for variable composite lengths of the non-sparse input-output data. The first modality is an implementation of the direct computation of a composite length DFT, the second one employs the second-order recursive filtering method, and the third one performs the new pruned decomposed transform. The pruned decomposed transform algorithm performs the decimation in time or space (DIT) data acquisition domain and, then, decimation in frequency (DIF). The unified combination of these three algorithms is addressed as the DFTCOMM technique. Based on the treatment of the combinational-type hypotheses testing optimization problem of preferable allocations between all feasible commuting-pruning modalities, we have found the global optimal solution to the pruning problem that always requires a fewer or, at most, the same number of arithmetic operations than other feasible modalities. The DFTCOMM method outperforms the existing competing pruning techniques in the sense of attainable savings in the number of required arithmetic operations. It requires fewer or at most the same number of arithmetic operations for its execution than any other of the competing pruning methods reported in the literature. Finally, we provide the comparison of the DFTCOMM with the recently developed sparse fast Fourier transform (SFFT) algorithmic family. We feature that, in the sensing scenarios with sparse/non-sparse data Fourier spectrum, the DFTCOMM technique manifests robustness against such model uncertainties in the sense of insensitivity for sparsity/non-sparsity restrictions and the variability of the operating parameters.

[1]  Franz Franchetti,et al.  Polar format synthetic aperture radar in energy efficient application-specific logic-in-memory , 2012, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[2]  Jörn Schumacher High performance Sparse Fast Fourier Transform , 2013 .

[3]  Christian Roche A split-radix partial input/output fast Fourier transform algorithm , 1992, IEEE Trans. Signal Process..

[4]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

[5]  Markus Püschel,et al.  High-performance sparse fast Fourier transforms , 2014 .

[6]  T. V. Sreenivas,et al.  FFT algorithm for both input and output pruning , 1979 .

[7]  Mark A. Iwen,et al.  Empirical evaluation of a sub-linear time sparse DFT algorithm , 2007 .

[8]  Modesto Guadalupe Medina Melendrez,et al.  Input and/or output pruning of composite length FFTs using a DIF-DIT transform decomposition , 2009 .

[9]  Xiaofang Zhou,et al.  Generic Mixed-Radix FFT Pruning , 2012, IEEE Signal Processing Letters.

[10]  慧 廣瀬 A Mathematical Introduction to Compressive Sensing , 2015 .

[11]  Piotr Indyk,et al.  Nearly optimal sparse fourier transform , 2012, STOC '12.

[12]  P. Townsend Principles and Applications of Imaging Radar: Manual of Remote Sensing , 2000 .

[13]  Yihu Xu,et al.  Split-radix FFT pruning for the reduction of computational complexity in OFDM based Cognitive Radio system , 2010, Proceedings of 2010 IEEE International Symposium on Circuits and Systems.

[14]  Chih-Peng Fan,et al.  Pruning fast Fourier transform algorithm design using group-based method , 2007, Signal Process..

[15]  D. Skinner Pruning the decimation in-time FFT algorithm , 1976 .

[16]  M. Omair Ahmad,et al.  Efficient pruning algorithms for the DFT computation for a subset of output samples , 2003, Proceedings of the 2003 International Symposium on Circuits and Systems, 2003. ISCAS '03..

[17]  Richard G. Baraniuk,et al.  Sparsity and Structure in Hyperspectral Imaging , 2014 .

[18]  Yuriy Shkvarko,et al.  ℓ2-ℓ1 Structured descriptive experiment design regularization based enhancement of fractional SAR imagery , 2013, Signal Process..

[19]  J. Markel,et al.  FFT pruning , 1971 .

[20]  Piotr Indyk,et al.  Simple and practical algorithm for sparse Fourier transform , 2012, SODA.

[21]  S. Rickard,et al.  Source Separation on Seismic Data , 2022 .

[22]  C. Sidney Burrus,et al.  On computing the split-radix FFT , 1986, IEEE Trans. Acoust. Speech Signal Process..

[23]  Holger Rauhut,et al.  A Mathematical Introduction to Compressive Sensing , 2013, Applied and Numerical Harmonic Analysis.

[24]  Scott T. Rickard,et al.  Source Separation on Seismic Data: Application in a Geophysical Setting , 2012, IEEE Signal Processing Magazine.

[25]  Yaowu Chen,et al.  Pruning split-radix FFT with time shift , 2011, 2011 International Conference on Electronics, Communications and Control (ICECC).

[26]  Ken'iti Kido,et al.  A new FFT algorithm of radix 3, 6, and 12 , 1986, IEEE Trans. Acoust. Speech Signal Process..

[27]  Richard G. Baraniuk,et al.  Sparsity and Structure in Hyperspectral Imaging : Sensing, Reconstruction, and Target Detection , 2014, IEEE Signal Processing Magazine.

[28]  Yuriy Shkvarko Unifying Experiment Design and Convex Regularization Techniques for Enhanced Imaging With Uncertain Remote Sensing Data—Part I: Theory , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[29]  Vasanthan Raghavan,et al.  Complexity of pruning strategies for the frequency domain LMS algorithm , 2006, Signal Process..

[30]  Yuriy Shkvarko,et al.  High-resolution imaging with uncertain radar measurement data: A doubly regularized compressive sensing experiment design approach , 2012, 2012 IEEE International Geoscience and Remote Sensing Symposium.

[31]  Alan V. Oppenheim,et al.  Discrete-Time Signal Pro-cessing , 1989 .

[32]  C. Sidney Burrus,et al.  Efficient computation of the DFT with only a subset of input or output points , 1993, IEEE Trans. Signal Process..

[33]  Yuriy Shkvarko,et al.  Unifying Experiment Design and Convex Regularization Techniques for Enhanced Imaging With Uncertain Remote Sensing Data—Part II: Adaptive Implementation and Performance Issues , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[34]  Harrison H. Barrett,et al.  Foundations of Image Science , 2003, J. Electronic Imaging.