v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization

convexoptimization.com
from convexoptimization.com More from this publisher
10.03.2015 Views

218 CHAPTER 3. GEOMETRY OF CONVEX FUNCTIONS variables which are later combined to solve the original problem (510). What makes this approach sound is that the constraints are separable, the partitioned feasible sets are not interdependent, and the fact that the original problem (though nonlinear) is convex simultaneously in both variables. 3.11 But partitioning alone does not guarantee a projector. To make orthogonal projector W a certainty, we must invoke a known analytical optimal solution to problem (512): Diagonalize optimal solution from problem (511) x ⋆ x ⋆T = ∆ QΛQ T (A.5.2) and set U ⋆ = Q(:, 1:k)∈ R n×k per (1581c); W = U ⋆ U ⋆T = x⋆ x ⋆T ‖x ⋆ ‖ 2 + Q(:, 2:k)Q(:, 2:k)T (513) Then optimal solution (x ⋆ , U ⋆ ) to problem (510) is found, for small ǫ , by iterating solution to problem (511) with optimal (projector) solution (513) to convex problem (512). Proof. Optimal vector x ⋆ is orthogonal to the last n −1 columns of orthogonal matrix Q , so f ⋆ (511) = ‖x ⋆ ‖ 2 (1 − (1 + ǫ) −1 ) (514) after each iteration. Convergence of f(511) ⋆ is proven with the observation that iteration (511) (512a) is a nonincreasing sequence that is bounded below by 0. Any bounded monotonic sequence in R is convergent. [222,1.2] [37,1.1] Expression (513) for optimal projector W holds at each iteration, therefore ‖x ⋆ ‖ 2 (1 − (1 + ǫ) −1 ) must also represent the optimal objective value f(511) ⋆ at convergence. Because the objective f (510) from problem (510) is also bounded below by 0 on the same domain, this convergent optimal objective value f(511) ⋆ (for positive ǫ arbitrarily close to 0) is necessarily optimal for (510); id est, by (1563), and f ⋆ (511) ≥ f ⋆ (510) ≥ 0 (515) lim ǫ→0 +f⋆ (511) = 0 (516) 3.11 A convex problem has convex feasible set, and the objective surface has one and only one global minimum.

3.1. CONVEX FUNCTION 219 Since optimal (x ⋆ , U ⋆ ) from problem (511) is feasible to problem (510), and because their objectives are equivalent for projectors by (507), then converged (x ⋆ , U ⋆ ) must also be optimal to (510) in the limit. Because problem (510) is convex, this represents a globally optimal solution. 3.1.7.2 Semidefinite program via Schur Schur complement (1410) can be used to convert a projection problem to an optimization problem in epigraph form. Suppose, for example, we are presented with the constrained projection problem studied by Hayden & Wells in [158] (who provide analytical solution): Given A∈ R M×M and some full-rank matrix S ∈ R M×L with L < M minimize ‖A − X‖ 2 X∈ S M F subject to S T XS ≽ 0 (517) Variable X is constrained to be positive semidefinite, but only on a subspace determined by S . First we write the epigraph form: minimize t X∈ S M , t∈R subject to ‖A − X‖ 2 F ≤ t S T XS ≽ 0 (518) Next we use Schur complement [239,6.4.3] [214] and matrix vectorization (2.2): minimize t X∈ S M , t∈R [ ] tI vec(A − X) subject to vec(A − X) T ≽ 0 (519) 1 S T XS ≽ 0 This semidefinite program (4) is an epigraph form in disguise, equivalent to (517); it demonstrates how a quadratic objective or constraint can be converted to a semidefinite constraint. Were problem (517) instead equivalently expressed without the square minimize ‖A − X‖ F X∈ S M subject to S T XS ≽ 0 (520)

3.1. CONVEX FUNCTION 219<br />

Since optimal (x ⋆ , U ⋆ ) from problem (511) is feasible to problem (510), and<br />

because their objectives are equivalent for projectors by (507), then converged<br />

(x ⋆ , U ⋆ ) must also be optimal to (510) in the limit. Because problem (510)<br />

is convex, this represents a globally optimal solution.<br />

<br />

3.1.7.2 Semidefinite program via Schur<br />

Schur complement (1410) can be used to convert a projection problem<br />

to an optimization problem in epigraph form. Suppose, for example,<br />

we are presented with the constrained projection problem studied by<br />

Hayden & Wells in [158] (who provide analytical solution): Given A∈ R M×M<br />

and some full-rank matrix S ∈ R M×L with L < M<br />

minimize ‖A − X‖ 2<br />

X∈ S M<br />

F<br />

subject to S T XS ≽ 0<br />

(517)<br />

Variable X is constrained to be positive semidefinite, but only on a subspace<br />

determined by S . First we write the epigraph form:<br />

minimize t<br />

X∈ S M , t∈R<br />

subject to ‖A − X‖ 2 F ≤ t<br />

S T XS ≽ 0<br />

(518)<br />

Next we use Schur complement [239,6.4.3] [214] and matrix vectorization<br />

(2.2):<br />

minimize t<br />

X∈ S M , t∈R<br />

[<br />

]<br />

tI vec(A − X)<br />

subject to<br />

vec(A − X) T ≽ 0 (519)<br />

1<br />

S T XS ≽ 0<br />

This semidefinite program (4) is an epigraph form in disguise, equivalent<br />

to (517); it demonstrates how a quadratic objective or constraint can be<br />

converted to a semidefinite constraint.<br />

Were problem (517) instead equivalently expressed without the square<br />

minimize ‖A − X‖ F<br />

X∈ S M<br />

subject to S T XS ≽ 0<br />

(520)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!