Positive solutions to linear inverse problems

被引:28
作者
de Villiers, GD
McNally, B
Pike, ER
机构
[1] DERA Malvern, Malvern WR14 3PS, Worcs, England
[2] Univ London Kings Coll, Dept Phys, London WC2R 2LS, England
关键词
D O I
10.1088/0266-5611/15/2/018
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We discuss two methods for incorporating the prior knowledge that the solution is positive into the truncated singular value decomposition method for solving linear inverse problems. The methods are based on mathematical programming techniques. One method can be viewed as a primal method and the other as its dual. Provided the singular functions are analytic these methods both deliver the same solution-namely the positive solution of minimum 2-norm which agrees with the truncated singular function expansion in its known terms-and this solution also appears to possess higher resolution. In the presence of noise both methods can sometimes fail to converge and in these situations we give simple remedies which yield approximate solutions. While there may be no reason to suppose the unknown object should be that of minimum 2-norm, our method has the advantage over other nonlinear methods that for a noiseless band-limited object the exact solution is delivered.
引用
收藏
页码:615 / 635
页数:21
相关论文
共 22 条
[1]   MAXIMUM-ENTROPY REGULARIZATION OF FREDHOLM INTEGRAL-EQUATIONS OF THE 1ST KIND [J].
AMATO, U ;
HUGHES, W .
INVERSE PROBLEMS, 1991, 7 (06) :793-808
[2]  
[Anonymous], 1976, P NONL PROGR SIAM AM
[3]  
[Anonymous], INTRO FUNCTIONAL ANA
[4]  
Bazaraa M.S., 2013, Nonlinear Programming-Theory and Algorithms, V3rd
[5]   RESOLUTION IN DIFFRACTION-LIMITED IMAGING, A SINGULAR VALUE ANALYSIS .1. THE CASE OF COHERENT ILLUMINATION [J].
BERTERO, M ;
PIKE, ER .
OPTICA ACTA, 1982, 29 (06) :727-746
[6]   LINEAR REGULARIZING ALGORITHMS FOR POSITIVE SOLUTIONS OF LINEAR INVERSE PROBLEMS [J].
BERTERO, M ;
BRIANZI, P ;
PIKE, ER ;
REBOLIA, L .
PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES, 1988, 415 (1849) :257-275
[7]   COMMUTING DIFFERENTIAL-OPERATORS FOR THE FINITE LAPLACE TRANSFORM [J].
BERTERO, M ;
GRUNBAUM, FA .
INVERSE PROBLEMS, 1985, 1 (03) :181-192
[8]   RESOLUTION IN DIFFRACTION-LIMITED IMAGING, A SINGULAR VALUE ANALYSIS .3. THE EFFECT OF SAMPLING AND TRUNCATION OF THE DATA [J].
BERTERO, M ;
BRIANZI, P ;
PARKER, P ;
PIKE, ER .
OPTICA ACTA, 1984, 31 (02) :181-201
[9]   ON THE RECOVERY AND RESOLUTION OF EXPONENTIAL RELAXATION RATES FROM EXPERIMENTAL-DATA - A SINGULAR-VALUE ANALYSIS OF THE LAPLACE TRANSFORM INVERSION IN THE PRESENCE OF NOISE [J].
BERTERO, M ;
BOCCACCI, P ;
PIKE, ER .
PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES, 1982, 383 (1784) :15-29
[10]   A SIMPLE CONSTRAINT QUALIFICATION IN INFINITE DIMENSIONAL PROGRAMMING [J].
BORWEIN, JM ;
WOLKOWICZ, H .
MATHEMATICAL PROGRAMMING, 1986, 35 (01) :83-96