A new algorithm in spectral analysis and band-limited extrapolation

If only a segment of a function is given, then its Fourier spectrum is estimated either as the transform of the product of with a time-limited window , or by certain techniques based on various a priori assumptions. In the following, a new algorithm is proposed for computing the transform of a band-limited function. The algorithm is a simple iteration involving only the fast Fourier transform (FFT). The effect of noise and the error due to aliasing are determined and it is shown that they can be controlled by early termination of the iteration. The proposed method can also be used to extrapolate bandlimited functions.