151 research outputs found
Recommended from our members
Structured Sub-Nyquist Sampling with Applications in Compressive Toeplitz Covariance Estimation, Super-Resolution and Phase Retrieval
Sub-Nyquist sampling has received a huge amount of interest in the past decade. In classical compressed sensing theory, if the measurement procedure satisfies a particular condition known as Restricted Isometry Property (RIP), we can achieve stable recovery of signals of low-dimensional intrinsic structures with an order-wise optimal sample size. Such low-dimensional structures include sparse and low rank for both vector and matrix cases. The main drawback of conventional compressed sensing theory is that random measurements are required to ensure the RIP property. However, in many applications such as imaging and array signal processing, applying independent random measurements may not be practical as the systems are deterministic. Moreover, random measurements based compressed sensing always exploits convex programs for signal recovery even in the noiseless case, and solving those programs is computationally intensive if the ambient dimension is large, especially in the matrix case. The main contribution of this dissertation is that we propose a deterministic sub-Nyquist sampling framework for compressing the structured signal and come up with computationally efficient algorithms. Besides widely studied sparse and low-rank structures, we particularly focus on the cases that the signals of interest are stationary or the measurements are of Fourier type. The key difference between our work from classical compressed sensing theory is that we explicitly exploit the second-order statistics of the signals, and study the equivalent quadratic measurement model in the correlation domain. The essential observation made in this dissertation is that a difference/sum coarray structure will arise from the quadratic model if the measurements are of Fourier type. With these observations, we are able to achieve a better compression rate for covariance estimation, identify more sources in array signal processing or recover the signals of larger sparsity. In this dissertation, we will first study the problem of Toeplitz covariance estimation. In particular, we will show how to achieve an order-wise optimal compression rate using the idea of sparse arrays in both general and low-rank cases. Then, an analysis framework of super-resolution with positivity constraint is established. We will present fundamental robustness guarantees, efficient algorithms and applications in practices. Next, we will study the problem of phase-retrieval for which we successfully apply the sparse array ideas by fully exploiting the quadratic measurement model. We achieve near-optimal sample complexity for both sparse and general cases with practical Fourier measurements and provide efficient and deterministic recovery algorithms. In the end, we will further elaborate on the essential role of non-negative constraint in underdetermined inverse problems. In particular, we will analyze the nonlinear co-array interpolation problem and develop a universal upper bound of the interpolation error. Bilinear problem with non-negative constraint will be considered next and the exact characterization of the ambiguous solutions will be established for the first time in literature. At last, we will show how to apply the nested array idea to solve real problems such as Kriging. Using spatial correlation information, we are able to have a stable estimate of the field of interest with fewer sensors than classic methodologies. Extensive numerical experiments are implemented to demonstrate our theoretical claims
GLRT-based threshold detection-estimation performance improvement and application to uniform circular antenna arrays
©2006 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE."This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder."The problem of estimating the number of independent Gaussian sources and their parameters impinging upon an antenna array is addressed for scenarios that are problematic for standard techniques, namely, under "threshold conditions" (where subspace techniques such as MUSIC experience an abrupt and dramatic performance breakdown). We propose an antenna geometry-invariant method that adopts the generalized-likelihood-ratio test (GLRT) methodology, supported by a maximum-likelihood-ratio lower-bound analysis that allows erroneous solutions ("outliers") to be found and rectified. Detection-estimation performance in both uniform circular and linear antenna arrays is shown to be significantly improved compared with conventional techniques but limited by the performance-breakdown phenomenon that is intrinsic to all such maximum-likelihood (ML) techniques.Yuri I. Abramovich, Nicholas K. Spencer, and Alexei Y. Gorokho
Toeplitz Inverse Eigenvalue Problem: Application to the Uniform Linear Antenna Array Calibration
The inverse Toeplitz eigenvalue problem (ToIEP) concerns finding a vector
that specifies the real-valued symmetric Toeplitz matrix with the prescribed
set of eigenvalues. Since phase "calibration" errors in uniform linear antenna
arrays (ULAs) do not change the covariance matrix eigenvalues and the moduli of
the covariance matrix elements, we formulate a number of the new ToIEP problems
of the Hermitian Toeplitz matrix reconstruction, given the moduli of the matrix
elements and the matrix eigenvalues. We demonstrate that for the real-valued
case, only two solutions to this problem exist, with the "non-physical" one
that in most practical cases could be easily disregarded. The computational
algorithm for the real-valued case is quite simple. For the complex-valued
case, we demonstrate that the family of solutions is broader and includes
solutions inappropriate for calibration. For this reason, we modified this
ToIEP problem to match the covariance matrix of the uncalibrated ULA. We
investigate the statistical convergence of the ad-hoc algorithm with the sample
matrices instead of the true ones. The proposed ad-hoc algorithms require the
so-called "strong" or "argumental" convergence, which means a large enough
required sample volume that reduces the errors in the estimated covariance
matrix elements. Along with the ULA arrays, we also considered the fully
augmentable minimum redundancy arrays that generate the same (full) set of
covariance lags as the uniform linear arrays, and we specified the conditions
when the ULA Toeplitz covariance matrix may be reconstructed given the
M-variate MRA covariance matrix.Comment: 27 pages, 40 figure
Subspace-based order estimation techniques in massive MIMO
Order estimation, also known as source enumeration, is a classical problem in array signal processing which consists in estimating the number of signals received by an array of sensors. In the last decades, numerous approaches to this problem have been proposed. However, the need of working with large-scale arrays (like in massive MIMO systems), low signal-to-noise- ratios, and poor sample regime scenarios, introduce new challenges to order estimation problems. For instance, most of the classical approaches are based on information theoretic criteria, which usually require a large sample size, typically several times larger than the number of sensors. Obtaining a number of samples several times larger than the number of sensors is not always possible with large-scale arrays. In addition, most of the methods found in literature assume that the noise is spatially white, which is very restrictive for many practical scenarios.
This dissertation deals with the problem of source enumeration for large-scale arrays, and proposes solutions that work robustly in the small sample regime under various noise models. The first part of the dissertation solves the problem by applying the idea of subspace averaging. The input data are modelled as subspaces, and an average or central subspace is computed. The source enumeration problem can be seen as an estimation of the dimension of the central subspace. A key element of the proposed method is to construct a bootstrap procedure, based on a newly proposed discrete distribution on the manifold of projection matrices, for stochastically generating subspaces from a function of experimentally determined eigenvalues. In this way, the proposed subspace averaging (SA) technique determines the order based on the eigenvalues of an average projection matrix, rather than on the likelihood of a covariance model, penalized by functions of the model order. The proposed SA criterion is especially effective in high-dimensional scenarios with low sample support for uniform linear arrays in the presence of white noise. Further, the proposed SA method is extended for: i) non-white noises, and ii) non-uniform linear arrays. The SA criterion is sensitive with the chosen dimension of extracted subspaces. To solve this problem, we combine the SA technique with a majority vote approach. The number of sources is detected for increasing dimensions of the SA technique and then a majority vote is applied to determine the final estimate. Further, to extend SA for arrays with arbitrary geometries, the SA is combined with a sparse reconstruction (SR) step. In the first step, each received snapshot is approximated by a sparse linear combination of the rest of snapshots. The SR problem is regularized by the logarithm-based surrogate of the l-0 norm and solved using a majorization-minimization approach. Based on the SR solution, a sampling mechanism is proposed in the second step to generate a collection of subspaces, all of which approximately span the same signal subspace. Finally, the dimension of the average of this collection of subspaces provides a robust estimate for the number of sources.
The second half of the dissertation introduces a completely different approach to the order estimation for uniform linear arrays, which is based on matrix completion algorithms. This part first discusses the problem of order estimation in the presence of noise whose spatial covariance structure is a diagonal matrix with possibly different variances. The diagonal terms of the sample covariance matrix are removed and, after applying Toeplitz rectification as a denoising step, the signal covariance matrix is reconstructed by using a low-rank matrix completion method adapted to enforce the Toeplitz structure of the sought solution. The proposed source enumeration criterion is based on the Frobenius norm of the reconstructed signal covariance matrix obtained for increasing rank values. The proposed method performs robustly for both small and large-scale arrays with few snapshots.
Finally, an approach to work with a reduced number of radio–frequency (RF) chains is proposed. The receiving array relies on antenna switching so that at every time instant only the signals received by a randomly selected subset of antennas are downconverted to baseband and sampled. Low-rank matrix completion (MC) techniques are then used to reconstruct the missing entries of the signal data matrix to keep the angular resolution of the original large-scale array. The proposed MC algorithm exploits not only the low- rank structure of the signal subspace, but also the shift-invariance property of uniform linear arrays, which results in a better estimation of the signal subspace. In addition, the effect of MC on DOA estimation is discussed under the perturbation theory framework. Further, this approach is extended to devise a novel order estimation criterion for missing data scenario. The proposed source enumeration criterion is based on the chordal subspace distance between two sub-matrices extracted from the reconstructed matrix after using MC for increasing rank values. We show that the proposed order estimation criterion performs consistently with a very few available entries in the data matrix.This work was supported by the Ministerio de Ciencia e Innovación (MICINN) of Spain, under grants TEC2016-75067-C4-4-R (CARMEN) and BES-2017-080542
Nested Sampling and its Applications in Stable Compressive Covariance Estimation and Phase Retrieval with Near-Minimal Measurements
Compressed covariance sensing using quadratic samplers is gaining increasing interest in recent literature. Covariance matrix often plays the role of a sufficient statistic in many signal and information processing tasks. However, owing to the large dimension of the data, it may become necessary to obtain a compressed sketch of the high dimensional covariance matrix to reduce the associated storage and communication costs. Nested sampling has been proposed in the past as an efficient sub-Nyquist sampling strategy that enables perfect reconstruction of the autocorrelation sequence of Wide-Sense Stationary (WSS) signals, as though it was sampled at the Nyquist rate. The key idea behind nested sampling is to exploit properties of the difference set that naturally arises in quadratic measurement model associated with covariance compression. In this thesis, we will focus on developing novel versions of nested sampling for low rank Toeplitz covariance estimation, and phase retrieval, where the latter problem finds many applications in high resolution optical imaging, X-ray crystallography and molecular imaging.
The problem of low rank compressive Toeplitz covariance estimation is first shown to be fundamentally related to that of line spectrum recovery. In absence if noise, this connection can be exploited to develop a particular kind of sampler called the Generalized Nested Sampler (GNS), that can achieve optimal compression rates. In presence of bounded noise, we develop a regularization-free algorithm that provably leads to stable recovery of the high dimensional Toeplitz matrix from its order-wise minimal sketch acquired using a GNS. Contrary to existing TV-norm and nuclear norm based reconstruction algorithms, our technique does not use any tuning parameters, which can be of great practical value.
The idea of nested sampling idea also finds a surprising use in the problem of phase retrieval, which has been of great interest in recent times for its convex formulation via PhaseLift, By using another modified version of nested sampling, namely the Partial Nested Fourier Sampler (PNFS), we show that with probability one, it is possible to achieve a certain conjectured lower bound on the necessary measurement size. Moreover, for sparse data, an l1 minimization based algorithm is proposed that can lead to stable phase retrieval using order-wise minimal number of measurements
Nested Arrays: A Novel Approach to Array Processing With Enhanced Degrees of Freedom
A new array geometry, which is capable of significantly
increasing the degrees of freedom of linear arrays, is
proposed. This structure is obtained by systematically nesting two
or more uniform linear arrays and can provide O(N^2) degrees
of freedom using only physical sensors when the second-order
statistics of the received data is used. The concept of nesting is
shown to be easily extensible to multiple stages and the structure
of the optimally nested array is found analytically. It is possible to
provide closed form expressions for the sensor locations and the
exact degrees of freedom obtainable from the proposed array as a
function of the total number of sensors. This cannot be done for
existing classes of arrays like minimum redundancy arrays which
have been used earlier for detecting more sources than the number
of physical sensors. In minimum-input–minimum-output (MIMO)
radar, the degrees of freedom are increased by constructing a
longer virtual array through active sensing. The method proposed
here, however, does not require active sensing and is capable of
providing increased degrees of freedom in a completely passive
setting. To utilize the degrees of freedom of the nested co-array, a
novel spatial smoothing based approach to DOA estimation is also
proposed, which does not require the inherent assumptions of the
traditional techniques based on fourth-order cumulants or quasi
stationary signals. As another potential application of the nested
array, a new approach to beamforming based on a nonlinear
preprocessing is also introduced, which can effectively utilize the
degrees of freedom offered by the nested arrays. The usefulness of
all the proposed methods is verified through extensive computer
simulations
- …