Abstract
Given a linear system Ax as b over the real or complex field, where both A and b are subject to noise, the total least squares (TLS) problem seeks to find a correction matrix and a correction right-hand side vector of minimal norm which makes the linear system feasible. To avoid ill posedness, a regularization term is added to the objective function; this leads to the so-called regularized TLS problem. A further complication arises when the matrix A and correspondingly the correction matrix must have a specific structure. This is modeled by the regularized structured TLS (RSTLS) problem. In general this problem is nonconvex and hence difficult to solve. However, the RSTLS problem arising from image deblurring applications under reflexive or periodic boundary conditions possesses a special structure where all relevant matrices are simultaneously diagonalizable (SD). In this paper we introduce an algorithm for finding the global optimum of the RSTLS problem with this SD structure. The devised method is based on decomposing the problem into single variable problems and then transforming them into one-dimensional unimodal real-valued minimization problems which can be solved globally. Based on the uniqueness and attainment properties of the RSTLS solution we show that a constrained version of the problem possesses a strong duality result and can thus be solved via a sequence of RSTLS problems.
Original language | English |
---|---|
Pages (from-to) | 419-443 |
Number of pages | 25 |
Journal | SIAM Journal on Matrix Analysis and Applications |
Volume | 30 |
Issue number | 1 |
DOIs | |
State | Published - 2008 |
Externally published | Yes |
Keywords
- Image deblurring
- Nonconvex optimization
- Simultaneously diagonalizable matrices
- Structured total least squares
- Unimodal functions