Greedy solution of ill-posed problems: error bounds and exact inversion

被引:25
作者
Denis, L. [1 ,2 ]
Lorenz, D. A. [3 ]
Trede, D. [4 ]
机构
[1] Ecole Super Chim Phys Elect Lyon, F-69616 Lyon, France
[2] Univ Lyon, CNRS, UMR 5516, Lab Hubert Curien, F-42000 St Etienne, France
[3] TU Braunschweig, D-38092 Braunschweig, Germany
[4] Univ Bremen, Zentrum Technomath, D-28334 Bremen, Germany
关键词
PARTICLE DIGITAL HOLOGRAPHY; CONVERGENCE-RATES; TIKHONOV REGULARIZATION; SIGNAL RECOVERY; BANACH-SPACES; SPARSITY; REPRESENTATIONS; MINIMIZATION; SHRINKAGE;
D O I
10.1088/0266-5611/25/11/115017
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The orthogonal matching pursuit (OMP) is a greedy algorithm to solve sparse approximation problems. Sufficient conditions for exact recovery are known with and without noise. In this paper we investigate the applicability of the OMP for the solution of ill-posed inverse problems in general, and in particular for two deconvolution examples from mass spectrometry and digital holography, respectively. In sparse approximation problems one often has to deal with the problem of redundancy of a dictionary, i.e. the atoms are not linearly independent. However, one expects them to be approximatively orthogonal and this is quantified by the so-called incoherence. This idea cannot be transferred to ill-posed inverse problems since here the atoms are typically far from orthogonal. The ill-posedness of the operator probably causes the correlation of two distinct atoms to become huge, i.e. that two atoms look much alike. Therefore, one needs conditions which take the structure of the problem into account and work without the concept of coherence. In this paper we develop results for the exact recovery of the support of noisy signals. In the two examples, mass spectrometry and digital holography, we show that our results lead to practically relevant estimates such that one may check a priori if the experimental setup guarantees exact deconvolution with OMP. Especially in the example from digital holography, our analysis may be regarded as a first step to calculate the resolution power of droplet holography.
引用
收藏
页数:24
相关论文
共 45 条
[1]  
[Anonymous], 2008, Introduction to Fourier optics
[2]  
[Anonymous], 2009, WAVELET TOUR SIGNAL
[3]   A generalized conditional gradient method for nonlinear operator equations with sparsity constraints [J].
Bonesky, Thomas ;
Bredies, Kristian ;
Lorenz, Dirk A. ;
Maass, Peter .
INVERSE PROBLEMS, 2007, 23 (05) :2041-2058
[4]   Morozov's discrepancy principle and Tikhonov-type functionals [J].
Bonesky, Thomas .
INVERSE PROBLEMS, 2009, 25 (01)
[5]  
BREDIES K, 2008, J PHYS C SER, V124
[6]   Regularization with non-convex separable constraints [J].
Bredies, Kristian ;
Lorenz, Dirk A. .
INVERSE PROBLEMS, 2009, 25 (08)
[7]   Iterated hard shrinkage for minimization problems with sparsity constraints [J].
Bredies, Kristian ;
Lorenz, Dirk A. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2008, 30 (02) :657-683
[9]   On the Uniqueness of Nonnegative Sparse Solutions to Underdetermined Systems of Equations [J].
Bruckstein, Alfred M. ;
Elad, Michael ;
Zibulevsky, Michael .
IEEE TRANSACTIONS ON INFORMATION THEORY, 2008, 54 (11) :4813-4820
[10]   Convergence rates of convex variational regularization [J].
Burger, M ;
Osher, S .
INVERSE PROBLEMS, 2004, 20 (05) :1411-1421