Abstract

Computing the Sparse Fast Fourier Transform(sFFT) has emerged as a critical topic for a long time. The sFFT algorithms decrease the runtime and sampling complexity by taking advantage of the signal's inherent characteristics that a large number of signals are sparse in the frequency domain. More than ten sFFT algorithms have been proposed, which can be classified into many types according to filter, framework, method of location, method of estimation. In this paper, the technology of these algorithms is completely analyzed in theory. The performance of them is thoroughly tested and verified in practice. The techniques involved in different sFFT algorithms include the following contents: five operations of signal, three methods of frequency bucketization, five methods of location, four methods of estimation, two problems caused by bucketization, three methods to solve these two problems, four algorithmic frameworks. All the above technologies and methods are introduced in detail and examples are given to illustrate the above research. After theoretical research, we make experiments for computing the signals of different SNR, length, sparsity by a standard testing platform and record the run time, percentage of the signal sampled and L <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">0</sub> ,L <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">1</sub> ,L <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sub> error with eight different typical sFFT algorithms. The result of experiments satisfies the inferences obtained in theory. It helps us to optimize these algorithms and use them selectively.

Highlights

  • The Discrete Fourier Transform(DFT) is one of the most important and widely used techniques in mathematical computing

  • The filtered signal is equal to the product of the original signal and the spike train filter in the time domain, it can be represented by x′ = DLx and it is equivalent to aliasing in the frequency domain

  • The filtered signal is equal to the convolution of the original signal and the spike train filter in the time domain, it can be represented by x′ = ULx and it is equivalent to subsampled in the frequency domain

Read more

Summary

Introduction

The Discrete Fourier Transform(DFT) is one of the most important and widely used techniques in mathematical computing. The most popular algorithm to compute DFT is the fast Fourier transform(FFT) invented by Cooley and Tukey, which can compute a signal of size N in O(N logN ) time and use O(N ) samples. The performance of the Ann Arbor fast Fourier transform(AAFFT0.5 [1]) algorithm was later improved in the AAFFT0.9 [2], [3] algorithm through the use of unequally-spaced FFTs and binary search technique for spectrum reconstruction. We present an empirical analysis of the performance of the algorithms in theory and in practice

Problem Statement of sFFT
Example and character of the above three operations
Method description
The second stage of sFFT
The first method of location
The second method of location
The third method of location
The fourth method of location
The first method of estimation
The second method of estimation
The third method of estimation
The fourth method of estimation
The first problem caused by bucketization
The second problem caused by bucketization
The first method to solve these problems
The second method to solve these problems
The third method to solve these problems
Experiment results and analysis
Conclusion

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.