A practical approach for the design of nonuniform lapped transforms

We propose a simple method for the design of lapped transforms with nonuniform frequency resolution and good time localization. The method is a generalization of an approach previously proposed by Princen, where the nonuniform filter bank is obtained by joining uniform cosine-modulated filter banks (CMFBs) using a transition filter. We use several transition filters to obtain a near perfect-reconstruction (PR) nonuniform lapped transform with significantly reduced overall distortion. The main advantage of the proposed method is in reducing the length of the transition filters, which leads to a reduction in processing delay that can be useful for applications such as real-time audio coding

Optimized Nonuniform FFTs and Their Application to Array Factor Computation

We deal with developing an optimized approach for implementing nonuniform fast Fourier transform (NUFFT) algorithms under a general and new perspective for 1-D transformations. The computations of nonequispaced results, nonequispaced data, and Type-3 nonuniform discrete Fourier transforms are tackled in a unified way. They exploit “uniformly sampled” exponentials to interpolate the “nonuniformly sampled” ones involved in the nonuniform discrete Fourier transforms (NUFDTs), so as to enable the use of standard fast Fourier transforms, and an optimized window. The computational costs and the memory requirements are analyzed, and their convenient performance is assessed also by comparing them with other approaches in the literature. Numerical results demonstrate that the method is more accurate and does not introduce any additional computational or memory burden. The computation of the window functions amounts to that of a Legendre polynomial expansion, i.e., a simple polynomial evaluation. This is convenient in terms of computational burden and of the proper arrangement of the calculations. A case study of electromagnetic interest has been carried out by applying the developed NUFFTs to the radiation of linear regular or irregular arrays onto a set of regular or irregular spectral points. Guidelines for multidimensional extension of the proposed approach are also presented

Nonuniform Fast Fourier Transforms Using Min-Max Interpolation

The fast Fourier transform (FFT) is used widely in signal processing for efficient computation of the FT of finite-length signals over a set of uniformly spaced frequency locations. However, in many applications, one requires nonuniform sampling in the frequency domain, i.e., a nonuniform FT. Several papers have described fast approximations for the nonuniform FT based on interpolating an oversampled FFT. This paper presents an interpolation method for the nonuniform FT that is optimal in the min-max sense of minimizing the worst-case approximation error over all signals of unit norm. The proposed method easily generalizes to multidimensional signals. Numerical results show that the min-max approach provides substantially lower approximation errors than conventional interpolation methods. The min-max criterion is also useful for optimizing the parameters of interpolation kernels such as the Kaiser-Bessel function.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85840/1/Fessler70.pd

Iterative Tomographic Image Reconstruction Using Fourier-Based Forward and Back-Projectors

Iterative image reconstruction algorithms play an increasingly important role in modern tomographic systems, especially in emission tomography. With the fast increase of the sizes of the tomographic data, reduction of the computation demands of the reconstruction algorithms is of great importance. Fourier-based forward and back-projection methods have the potential to considerably reduce the computation time in iterative reconstruction. Additional substantial speed-up of those approaches can be obtained utilizing powerful and cheap off-the-shelf fast Fourier transform (FFT) processing hardware. The Fourier reconstruction approaches are based on the relationship between the Fourier transform of the image and Fourier transformation of the parallel-ray projections. The critical two steps are the estimations of the samples of the projection transform, on the central section through the origin of Fourier space, from the samples of the transform of the image, and vice versa for back-projection. Interpolation errors are a limitation of Fourier-based reconstruction methods. We have applied min-max optimized Kaiser-Bessel interpolation within the nonuniform FFT (NUFFT) framework and devised ways of incorporation of resolution models into the Fourier-based iterative approaches. Numerical and computer simulation results show that the min-max NUFFT approach provides substantially lower approximation errors in tomographic forward and back-projection than conventional interpolation methods. Our studies have further confirmed that Fourier-based projectors using the NUFFT approach provide accurate approximations to their space-based counterparts but with about ten times faster computation, and that they are viable candidates for fast iterative image reconstruction.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85804/1/Fessler62.pd

Fine structures of solar radio type III bursts and their possible relationship with coronal density turbulence

Solar radio type III bursts are believed to be the most sensitive signatures of near-relativistic electron beam propagation in the corona. A solar radio type IIIb-III pair burst with fine frequency structures, observed by the Low Frequency Array (LOFAR) with high temporal (~10 ms) and spectral (12.5 kHz) resolutions at 30–80 MHz, is presented. The observations show that the type III burst consists of many striae, which have a frequency scale of about 0.1 MHz in both the fundamental (plasma) and the harmonic (double plasma) emission. We investigate the effects of background density fluctuations based on the observation of striae structure to estimate the density perturbation in the solar corona. It is found that the spectral index of the density fluctuation spectrum is about −1.7, and the characteristic spatial scale of the density perturbation is around 700 km. This spectral index is very close to a Kolmogorov turbulence spectral index of −5/3, consistent with a turbulent cascade. This fact indicates that the coronal turbulence may play the important role of modulating the time structures of solar radio type III bursts, and the fine structure of radio type III bursts could provide a useful and unique tool to diagnose the turbulence in the solar corona