v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
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)
- Page 167 and 168: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 169 and 170: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 171 and 172: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 173 and 174: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 175 and 176: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 177 and 178: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 179 and 180: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 181 and 182: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 183 and 184: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 185 and 186: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 187 and 188: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 189 and 190: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 191 and 192: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 193 and 194: 2.13. DUAL CONE & GENERALIZED INEQU
- Page 195 and 196: Chapter 3 Geometry of convex functi
- Page 197 and 198: 3.1. CONVEX FUNCTION 197 f 1 (x) f
- Page 199 and 200: 3.1. CONVEX FUNCTION 199 3.1.3 norm
- Page 201 and 202: 3.1. CONVEX FUNCTION 201 A B 1 Figu
- Page 203 and 204: 3.1. CONVEX FUNCTION 203 k/m 1 0.9
- Page 205 and 206: k∑ i=1 3.1. CONVEX FUNCTION 205 S
- Page 207 and 208: 3.1. CONVEX FUNCTION 207 rather x >
- Page 209 and 210: 3.1. CONVEX FUNCTION 209 rather ] x
- Page 211 and 212: 3.1. CONVEX FUNCTION 211 3.1.6.0.2
- Page 213 and 214: 3.1. CONVEX FUNCTION 213 q(x) f(x)
- Page 215 and 216: 3.1. CONVEX FUNCTION 215 3.1.7.0.2
- Page 217: 3.1. CONVEX FUNCTION 217 We learned
- Page 221 and 222: 3.1. CONVEX FUNCTION 221 2 1.5 1 0.
- Page 223 and 224: 3.1. CONVEX FUNCTION 223 Setting th
- Page 225 and 226: 3.1. CONVEX FUNCTION 225 Similarly,
- Page 227 and 228: 3.1. CONVEX FUNCTION 227 For vector
- Page 229 and 230: 3.1. CONVEX FUNCTION 229 This means
- Page 231 and 232: 3.1. CONVEX FUNCTION 231 f(Y ) −
- Page 233 and 234: 3.2. MATRIX-VALUED CONVEX FUNCTION
- Page 235 and 236: 3.2. MATRIX-VALUED CONVEX FUNCTION
- Page 237 and 238: 3.2. MATRIX-VALUED CONVEX FUNCTION
- Page 239 and 240: 3.3. QUASICONVEX 239 exponential al
- Page 241 and 242: 3.3. QUASICONVEX 241 Unlike convex
- Page 243 and 244: 3.4. SALIENT PROPERTIES 243 6. (af
- Page 245 and 246: Chapter 4 Semidefinite programming
- Page 247 and 248: 4.1. CONIC PROBLEM 247 where K is a
- Page 249 and 250: 4.1. CONIC PROBLEM 249 4.1.1.2 Redu
- Page 251 and 252: 4.1. CONIC PROBLEM 251 In any SDP f
- Page 253 and 254: 4.1. CONIC PROBLEM 253 Proposition
- Page 255 and 256: 4.2. FRAMEWORK 255 sets are closed
- Page 257 and 258: 4.2. FRAMEWORK 257 4.2.1.1.3 Exampl
- Page 259 and 260: 4.2. FRAMEWORK 259 4.2.2 Duals The
- Page 261 and 262: 4.2. FRAMEWORK 261 When equality is
- Page 263 and 264: 4.2. FRAMEWORK 263 The pseudoinvers
- Page 265 and 266: 4.2. FRAMEWORK 265 For the data giv
- Page 267 and 268: 4.2. FRAMEWORK 267 minimizes an aff
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)