Robust uncertainty principles:: Exact signal reconstruction from highly incomplete frequency information

被引:11197
作者
Candès, EJ
Romberg, J
Tao, T
机构
[1] CALTECH, Dept Appl & Computat Math, Pasadena, CA 91125 USA
[2] Univ Calif Los Angeles, Dept Math, Los Angeles, CA 90095 USA
基金
美国国家科学基金会;
关键词
convex optimization; duality in optimization; free probability; image reconstruction; linear programming; random matrices; sparsity; total-variation minimization; trigonometric expansions; uncertainty principle;
D O I
10.1109/TIT.2005.862083
中图分类号
TP [自动化技术、计算机技术];
学科分类号
0812 ;
摘要
This paper considers the model problem of reconstructing an object from incomplete frequency samples. Consider a discrete-time signal f is an element of C-N and a randomly chosen set of frequencies Q. Is it possible to reconstruct f from the partial knowledge of its Fourier coefficients on the set Q? A typical result of this paper is as follows. Suppose that f is a superposition of vertical bar T vertical bar spikes f(t) = E-tau is an element of T f(tau)delta(t - tau) obeying vertical bar T vertical bar <= C-M (.) (logN)(-1 .) vertical bar ohm vertical bar for some constant C-M > 0. We do not know the locations of the spikes nor their amplitudes. Then with probability at least 1 - O(N-m), f can be reconstructed exactly as the solution to the l(1) minimization problem min/g Sigma(N-1)/t=0 vertical bar g(t)vertical bar, s.t. (g) over cap(omega) = (f) over cap(omega) for al omega is an element of ohm In short, exact recovery may be obtained by solving a convex optimization problem. We give numerical values for Cm which depend on the desired probability of success. Our result may be interpreted as a novel kind of nonlinear sampling theorem. In effect, it says that any signal made out of vertical bar T vertical bar spikes may be recovered by convex programming from almost every set of frequencies of size O(vertical bar T vertical bar (.) log N). Moreover, this is nearly optimal in the sense that any method succeeding with probability 1 - O(N-M) would in general require a number of frequency samples at least proportional to vertical bar T vertical bar (.) log N. The methodology extends to a variety of other situations and higher dimensions. For example, we show how one can reconstruct a piecewise constant (one- or two-dimensional) object from incomplete frequency samples-provided that the number of jumps (discontinuities) obeys the condition above-by minimizing other convex functionals such as the total variation of f.
引用
收藏
页码:489 / 509
页数:21
相关论文
共 32 条
[1]  
Boucheron S, 2000, RANDOM STRUCT ALGOR, V16, P277, DOI 10.1002/(SICI)1098-2418(200005)16:3<277::AID-RSA4>3.0.CO
[2]  
2-1
[3]  
CANDES EJ, 2002, IMAGE RECONSTRUCTION
[4]  
Chen Y, 1998, NONCON OPTIM ITS APP, V20, P1
[5]   A fast and accurate Fourier algorithm for iterative parallel-beam tomography [J].
Delaney, AH ;
Bresler, Y .
IEEE TRANSACTIONS ON IMAGE PROCESSING, 1996, 5 (05) :740-753
[6]   Recovery of blocky images from noisy and blurred data [J].
Dobson, DC ;
Santosa, F .
SIAM JOURNAL ON APPLIED MATHEMATICS, 1996, 56 (04) :1181-1198
[7]   UNCERTAINTY PRINCIPLES AND SIGNAL RECOVERY [J].
DONOHO, DL ;
STARK, PB .
SIAM JOURNAL ON APPLIED MATHEMATICS, 1989, 49 (03) :906-931
[8]   Optimally sparse representation in general (nonorthogonal) dictionaries via l1 minimization [J].
Donoho, DL ;
Elad, M .
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA, 2003, 100 (05) :2197-2202
[9]   SIGNAL RECOVERY AND THE LARGE SIEVE [J].
DONOHO, DL ;
LOGAN, BF .
SIAM JOURNAL ON APPLIED MATHEMATICS, 1992, 52 (02) :577-591
[10]   Uncertainty principles and ideal atomic decomposition [J].
Donoho, DL ;
Huo, XM .
IEEE TRANSACTIONS ON INFORMATION THEORY, 2001, 47 (07) :2845-2862