Direction of Arrival Estimation and Tracking with Sparse Arrays

Direction of Arrival (DOA) estimation and tracking of a plane wave or multiple plane waves impinging on an array of sensors from noisy data are two of the most important tasks in array signal processing, which have attracted tremendous research interest over the past several decades. It is well-known that the estimation accuracy, angular resolution, tracking capacity, computational complexity, and hardware implementation cost of a DOA estimation and/or tracking technique depend largely on the array geometry. Large arrays with many sensors provide accurate DOA estimation and perfect target tracking, but they usually suffer from a high cost for hardware implementation. Sparse arrays can yield similar DOA estimates and tracking performance with fewer elements for the same-size array aperture as compared to the traditional uniform arrays. In addition, the signals of interest may have rich temporal information that can be exploited to effectively eliminate background noise and significantly improve the performance and capacity of DOA estimation and tracking, and/or even dramatically reduce the computational burden of estimation and tracking algorithms. Therefore, this thesis aims to provide some solutions to improving the DOA estimation and tracking performance by designing sparse arrays and exploiting prior knowledge of the incident signals such as AR modeled sources and known waveforms. First, we design two sparse linear arrays to efficiently extend the array aperture and improve the DOA estimation performance. One scheme is called minimum redundancy sparse subarrays (MRSSA), where the subarrays are used to obtain an extended correlation matrix according to the principle of minimum redundancy linear array (MRLA). The other linear array is constructed using two sparse ULAs, where the inter-sensor spacing within the same ULA is much larger than half wavelength. Moreover, we propose a 2-D DOA estimation method based on sparse L-shaped arrays, where the signal subspace is selected from the noise-free correlation matrix without requiring the eigen-decomposition to estimate the elevation angle, while the azimuth angles are estimated based on the modified total least squares (TLS) technique. Second, we develop two DOA estimation and tracking methods for autoregressive (AR) modeled signal source using sparse linear arrays together with Kalman filter and LS-based techniques. The proposed methods consist of two common stages: in the first stage, the sources modeled by AR processes are estimated by the celebrated Kalman filter and in the second stage, the efficient LS or TLS techniques are employed to estimate the DOAs and AR coefficients simultaneously. The AR-modeled sources can provide useful temporal information to handle cases such as the ones, where the number of sources is larger than the number of antennas. In the first method, we exploit the symmetric array to transfer a complex-valued nonlinear problem to a real-valued linear one, which can reduce the computational complexity, while in the second method, we use the ordinary sparse arrays to provide a more accurate DOA estimation. Finally, we study the problem of estimating and tracking the direction of arrivals (DOAs) of multiple moving targets with known signal source waveforms and unknown gains in the presence of Gaussian noise using a sparse sensor array. The core idea is to consider the output of each sensor as a linear regression model, each of whose coefficients contains a pair of DOAs and gain information corresponding to one target. These coefficients are determined by solving a linear least squares problem and then updating recursively, based on a block QR decomposition recursive least squares (QRD-RLS) technique or a block regularized LS technique. It is shown that the coefficients from different sensors have the same amplitude, but variable phase information for the same signal. Then, simple algebraic manipulations and the well-known generalized least squares (GLS) are used to obtain an asymptotically-optimal DOA estimate without requiring a search over a large region of the parameter space

Broadband DOA estimation using Convolutional neural networks trained with noise signals

A convolution neural network (CNN) based classification method for broadband DOA estimation is proposed, where the phase component of the short-time Fourier transform coefficients of the received microphone signals are directly fed into the CNN and the features required for DOA estimation are learnt during training. Since only the phase component of the input is used, the CNN can be trained with synthesized noise signals, thereby making the preparation of the training data set easier compared to using speech signals. Through experimental evaluation, the ability of the proposed noise trained CNN framework to generalize to speech sources is demonstrated. In addition, the robustness of the system to noise, small perturbations in microphone positions, as well as its ability to adapt to different acoustic conditions is investigated using experiments with simulated and real data.Comment: Published in Proceedings of IEEE Workshop on Applications of Signal Processing to Audio and Acoustics (WASPAA) 201

Multiple and single snapshot compressive beamforming

For a sound field observed on a sensor array, compressive sensing (CS) reconstructs the direction-of-arrival (DOA) of multiple sources using a sparsity constraint. The DOA estimation is posed as an underdetermined problem by expressing the acoustic pressure at each sensor as a phase-lagged superposition of source amplitudes at all hypothetical DOAs. Regularizing with an $\ell_1$-norm constraint renders the problem solvable with convex optimization, and promoting sparsity gives high-resolution DOA maps. Here, the sparse source distribution is derived using maximum a posteriori (MAP) estimates for both single and multiple snapshots. CS does not require inversion of the data covariance matrix and thus works well even for a single snapshot where it gives higher resolution than conventional beamforming. For multiple snapshots, CS outperforms conventional high-resolution methods, even with coherent arrivals and at low signal-to-noise ratio. The superior resolution of CS is demonstrated with vertical array data from the SWellEx96 experiment for coherent multi-paths.Comment: In press Journal of Acoustical Society of Americ

Biomimetic direction of arrival estimation for resolving front-back confusions in hearing aids

Sound sources at the same angle in front or behind a two-microphone array (e.g., bilateral hearing aids) produce the same time delay and two estimates for the direction of arrival: A front-back confusion. The auditory system can resolve this issue using head movements. To resolve front-back confusion for hearing-aid algorithms, head movement was measured using an inertial sensor. Successive time-delay estimates between the microphones are shifted clockwise and counterclockwise by the head movement between estimates and aggregated in two histograms. The histogram with the largest peak after multiple estimates predicted the correct hemifield for the source, eliminating the front-back confusions

Parametric high resolution techniques for radio astronomical imaging

The increased sensitivity of future radio telescopes will result in requirements for higher dynamic range within the image as well as better resolution and immunity to interference. In this paper we propose a new matrix formulation of the imaging equation in the cases of non co-planar arrays and polarimetric measurements. Then we improve our parametric imaging techniques in terms of resolution and estimation accuracy. This is done by enhancing both the MVDR parametric imaging, introducing alternative dirty images and by introducing better power estimates based on least squares, with positive semi-definite constraints. We also discuss the use of robust Capon beamforming and semi-definite programming for solving the self-calibration problem. Additionally we provide statistical analysis of the bias of the MVDR beamformer for the case of moving array, which serves as a first step in analyzing iterative approaches such as CLEAN and the techniques proposed in this paper. Finally we demonstrate a full deconvolution process based on the parametric imaging techniques and show its improved resolution and sensitivity compared to the CLEAN method.Comment: To appear in IEEE Journal of Selected Topics in Signal Processing, Special issue on Signal Processing for Astronomy and space research. 30 page

Clutter rejection for MTI radar using a single antenna and a long integration time

Moving Target Indicators (MTI) are airborne radar systems designed to detect and track moving vehicles or aircrafts. In this paper, we address the problem of detecting hazardous collision targets to avoid them. One of the best known solutions to solve this problem is given by the so-called Space-Time Adaptive Processing (STAP) algorithms which optimally filter the target signal from interference and noise exploiting the specific relationship between Direction Of Arrival (DOA) and Doppler for the ground clutter. However, these algorithms require an antenna array and multiple reception channels that increase cost and complexity. The authors propose an alternative solution using a single antenna only. In addition to the standard Doppler shift related to the radial speed, the orthoradial speed of any target can be estimated if using a long integration time. Dangerous targets and ground clutter have different signatures in the radial-orthoradial velocity plane. An optimal detector is then proposed based on the oblique projection onto the signal subspace orthogonal to the clutter subspace. The theoretical performances of this detector are derived and a realistic radar scene simulation shows the benefits of this new MTI detector