# Convex Iteration

### From Wikimization

(→constraining cardinality of signed variable) |
(→constraining cardinality of signed variable) |
||

Line 15: | Line 15: | ||

\mbox{subject to}&A_{}x=b\\ | \mbox{subject to}&A_{}x=b\\ | ||

&\|x\|_0\leq_{}k | &\|x\|_0\leq_{}k | ||

- | \end{array} | + | \end{array}</math> (1) |

where <math>\|x\|_{0\!}\leq_{_{}\!}k</math> means vector <math>x</math> has at most <math>k</math> nonzero entries; | where <math>\|x\|_{0\!}\leq_{_{}\!}k</math> means vector <math>x</math> has at most <math>k</math> nonzero entries; | ||

Line 27: | Line 27: | ||

\mbox{subject to}&A_{}x=b\\ | \mbox{subject to}&A_{}x=b\\ | ||

&-t\preceq x\preceq_{_{}}t | &-t\preceq x\preceq_{_{}}t | ||

- | \end{array} | + | \end{array}</math> (2) |

is iterated with | is iterated with | ||

Line 34: | Line 34: | ||

\mbox{subject to}&0\preceq y\preceq\mathbf{1}\\ | \mbox{subject to}&0\preceq y\preceq\mathbf{1}\\ | ||

&y^T\mathbf{1}=n-_{}k | &y^T\mathbf{1}=n-_{}k | ||

- | \end{array} | + | \end{array}</math> (3) |

to find ''direction vector'' <math>y</math>. | to find ''direction vector'' <math>y</math>. | ||

Line 53: | Line 53: | ||

\mbox{subject to}&A_{}x=b | \mbox{subject to}&A_{}x=b | ||

\end{array} | \end{array} | ||

- | \end{array} | + | \end{array}</math> (4) |

Subsequent iterations of problem (2) engaging cardinality term <math>\langle t_{}\,,\,y\rangle</math> | Subsequent iterations of problem (2) engaging cardinality term <math>\langle t_{}\,,\,y\rangle</math> | ||

Line 67: | Line 67: | ||

===equivalent=== | ===equivalent=== | ||

Several simple equivalents to linear programs (2) (3) are easily devised, | Several simple equivalents to linear programs (2) (3) are easily devised, | ||

- | but their geometrical interpretation is not as apparent: ''e.g.'', equivalent in the limit <math>\,\varepsilon\!\ | + | but their geometrical interpretation is not as apparent: ''e.g.'', equivalent in the limit <math>\,\varepsilon\!\rightarrow0^+\,</math> |

<math>\begin{array}{cl}\text{minimize}_{x_{}\in_{_{}}\mathbb{R}^n,~t_{}\in_{_{}}\mathbb{R}^n}&\langle t\,,\,y\rangle\\ | <math>\begin{array}{cl}\text{minimize}_{x_{}\in_{_{}}\mathbb{R}^n,~t_{}\in_{_{}}\mathbb{R}^n}&\langle t\,,\,y\rangle\\ |

## Revision as of 23:08, 4 December 2011

Convex iteration is method for constraining rank or cardinality in an otherwise convex optimization problem.

A rank or cardinality constraint is replaced by a linear regularization term added to the objective.

Then two convex problems are iterated until convergence where, ideally, solution to the original problem is found.

## constraining cardinality of signed variable

Excert from Chapter 4 Semidefinite Programming:

Consider a feasibility problem equivalent to the classical problem from linear algebra
,
but with an upper bound on cardinality :
for vector

(1)

where means vector has at most nonzero entries;
such a vector is presumed existent in the feasible set.

Convex iteration works with a nonnegative variable; absolute value is therefore needed.

We propose that nonconvex problem (1) can be equivalently written
as a convex iteration:
for a relatively small positive constant,

(2)

is iterated with

(3)

to find *direction vector* .
The cardinality constraint has been moved to the objective as a linear regularization.

The term in (2) is necessary to determine absolute value
because vector can take zero values in its entries.

By initializing direction vector to ,
the first iteration of problem (2) is a 1-norm problem; *i.e.*,

(4)

Subsequent iterations of problem (2) engaging cardinality term
can be interpreted as corrections to this 1-norm problem leading to a -norm solution.

Iteration (2) (3) always converges to a locally optimal solution by virtue of
a monotonically nonincreasing objective sequence.

There can be no proof of global optimality, defined by an optimal objective equal to 0.

Local solutions are therefore detected by nonzero optimal objective.

Heuristics for reinitializing direction vector can lead to a globally optimal solution.

### equivalent

Several simple equivalents to linear programs (2) (3) are easily devised,
but their geometrical interpretation is not as apparent: *e.g.*, equivalent in the limit

For a coded numerical example, see Candes.m