Suppose that we have @math sensors and each one intends to send a function @math (e.g. a signal or an image) to a receiver common to all @math sensors. During transmission, each @math gets convolved with a function @math . The receiver records the function @math , given by the sum of all these convolved signals. When and under which conditions is it possible to recover the individual signals @math and the blurring functions @math from just one received signal @math ? This challenging problem, which intertwines blind deconvolution with blind demixing, appears in a variety of applications, such as audio processing, image processing, neuroscience, spectroscopy, and astronomy. It is also expected to play a central role in connection with the future Internet-of-Things. We will prove that under reasonable and practical assumptions, it is possible to solve this otherwise highly ill-posed problem and recover the @math transmitted functions @math and the impulse responses @math in a robust, reliable, and efficient manner from just one single received function @math by solving a semidefinite program. We derive explicit bounds on the number of measurements needed for successful recovery and prove that our method is robust in the presence of noise. Our theory is actually sub-optimal, since numerical experiments demonstrate that, quite remarkably, recovery is still possible if the number of measurements is close to the number of degrees of freedom.
Problems of the type or are ubiquitous in many applied scientific disciplines and in applications, see e.g @cite_19 @cite_24 @cite_25 @cite_37 @cite_16 @cite_23 @cite_29 @cite_45 @cite_22 @cite_11 @cite_43 . Thus, there is a large body of works to solve different versions of these problems. Most of the existing works however require the availability of multiple received signals @math . And indeed, it is not hard to imagine that for instance an SVD-based approach will succeed if @math (and must fail if @math ). A sparsity-based approach can be found in @cite_8 . However, in this paper we are interested in the case where we have only one single received signal @math -- a single snapshot, in the jargon of array processing. Hence, there is little overlap between these methods heavily relying on multiple snapshots (many of which do not come with any theory) and the work presented here.
@cite_37 @cite_22 @cite_8 @cite_29 @cite_24 @cite_19 @cite_43 @cite_45 @cite_23 @cite_16 @cite_25 @cite_11
1976709621 2962959718 2568623927 2275123856
A full-rank matrix @math with @math generates an underdetermined system of linear equations @math having infinitely many solutions. Suppose we seek the sparsest solution, i.e., the one with the fewest nonzero entries. Can it ever be unique? If so, when? As optimization of sparsity is combinatorial in nature, are there efficient methods for finding the sparsest solution? These questions have been answered positively and constructively in recent years, exposing a wide variety of surprising phenomena, in particular the existence of easily verifiable conditions under which optimally sparse solutions can be found by concrete, effective computational methods. Such theoretical results inspire a bold perspective on some important practical problems in signal and image processing. Several well-known signal and image processing problems can be cast as demanding solutions of undetermined systems of equations. Such problems have previously seemed, to many, intractable, but there is considerable evidence that these problems often have sparse solutions. Hence, advances in finding sparse solutions to underdetermined systems have energized research on such signal and image processing problems—to striking effect. In this paper we review the theoretical results on sparse solutions of linear systems, empirical results on sparse modeling of signals and images, and recent applications in inverse problems and compression in image processing. This work lies at the intersection of signal processing and applied mathematics, and arose initially from the wavelets and harmonic analysis research communities. The aim of this paper is to introduce a few key notions and applications connected to sparsity, targeting newcomers interested in either the mathematical aspects of this area or its applications. The idea that many important classes of signals can be well-represented by linear combinations of a small set of atoms selected from a given dictionary has had dramatic impact on the theory and practice of signal processing. For practical problems in which an appropriate sparsifying dictionary is not known ahead of time, a very popular and successful heuristic is to search for a dictionary that minimizes an appropriate sparsity surrogate over a given set of sample data. While this idea is appealing, the behavior of these algorithms is largely a mystery; although there is a body of empirical evidence suggesting they do learn very eective representations, there is little theory to guarantee when they will behave correctly, or when the learned dictionary can be expected to generalize. In this paper, we take a step towards such a theory. We show that under mild hypotheses, the dictionary learning problem is locally well-posed: the desired solution is indeed a local minimum of the ‘ 1 norm. Namely, if A 2 R m n is an incoherent (and possibly overcomplete) dictionary, and the coecients X 2 R n p follow a random sparse model, then with high probability (A;X) is a local minimum of the ‘ 1 norm over the manifold of factorizations (A 0 ;X 0 ) satisfying A 0 X 0 = Y , provided the number of samples p = Consider the problem of recovering an unknown signal from undersampled measurements, given the knowledge that the signal has a sparse representation in a specified dictionary @math . This problem is now understood to be well-posed and efficiently solvable under suitable assumptions on the measurements and dictionary if the number of measurements scales roughly with the sparsity level. One sufficient condition for such is the @math -restricted isometry property ( @math -RIP), which asks that the sampling matrix approximately preserve the norm of all signals which are sufficiently sparse in @math . While many classes of random matrices are known to satisfy such conditions, such matrices are not representative of the structural constraints imposed by practical sensing systems. We close this gap in the theory by demonstrating that one can subsample a fixed orthogonal matrix in such a way that the @math -RIP will hold, provided this basis is sufficiently incoherent with the sparsifying dictionary $ ... Recent development of compressive sensing has greatly benefited radar imaging problems. In this paper, we investigate the problem of obtaining enhanced targets such as ships and airplanes, where targets often exhibit structured sparsity. A novel structured sparsity-driven autofocus algorithm is proposed based on sparse Bayesian framework.The structured sparse prior is imposed on the target scene in a statistical manner. Based on a statistical framework, the proposed algorithm can simultaneously cope with structured sparse recovery and phase error correction problem. The focused high-resolution radar image can be obtained by iteratively estimating scattering coefficients and phase. Due to the structured sparse constraint, the proposed algorithm can desirably preserve the target region and alleviate over-shrinkage problem, compared to previous sparsity-driven auto-focus approaches.Moreover, to accelerate convergence rate of the algorithm, we propose to adaptively eliminate noise-only range cells in estimating phase errors. The selection is conveniently conducted based on the parameters controlling sparsity degree of the signal in the proposed hierarchical model.The simulated and real data experimental results demonstrate that the proposed algorithm can obtain more concentrated images with much smaller number of iterations, particularly in low SNR and highly under-sampling scenarios. HighlightsA novel structured sparsity-driven autofocus algorithm is proposed based on sparse Bayesian framework.To accelerate convergence rate of the algorithm, we propose to adaptively eliminate noise-only range cells in phase error estimation stage.The proposed algorithm can obtain more concentrated images with much smaller number of iterations, particularly in low SNR and highly under-sampling scenarios.