10.03.2015 Views

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

218 CHAPTER 3. GEOMETRY OF CONVEX FUNCTIONS<br />

variables which are later combined to solve the original problem (510).<br />

What makes this approach sound is that the constraints are separable, the<br />

partitioned feasible sets are not interdependent, and the fact that the original<br />

problem (though nonlinear) is convex simultaneously in both variables. 3.11<br />

But partitioning alone does not guarantee a projector. To make<br />

orthogonal projector W a certainty, we must invoke a known analytical<br />

optimal solution to problem (512): Diagonalize optimal solution from<br />

problem (511) x ⋆ x ⋆T = ∆ QΛQ T (A.5.2) and set U ⋆ = Q(:, 1:k)∈ R n×k<br />

per (1581c);<br />

W = U ⋆ U ⋆T = x⋆ x ⋆T<br />

‖x ⋆ ‖ 2 + Q(:, 2:k)Q(:, 2:k)T (513)<br />

Then optimal solution (x ⋆ , U ⋆ ) to problem (510) is found, for small ǫ , by<br />

iterating solution to problem (511) with optimal (projector) solution (513)<br />

to convex problem (512).<br />

Proof. Optimal vector x ⋆ is orthogonal to the last n −1 columns of<br />

orthogonal matrix Q , so<br />

f ⋆ (511) = ‖x ⋆ ‖ 2 (1 − (1 + ǫ) −1 ) (514)<br />

after each iteration. Convergence of f(511) ⋆ is proven with the observation that<br />

iteration (511) (512a) is a nonincreasing sequence that is bounded below by 0.<br />

Any bounded monotonic sequence in R is convergent. [222,1.2] [37,1.1]<br />

Expression (513) for optimal projector W holds at each iteration, therefore<br />

‖x ⋆ ‖ 2 (1 − (1 + ǫ) −1 ) must also represent the optimal objective value f(511)<br />

⋆<br />

at convergence.<br />

Because the objective f (510) from problem (510) is also bounded below<br />

by 0 on the same domain, this convergent optimal objective value f(511) ⋆ (for<br />

positive ǫ arbitrarily close to 0) is necessarily optimal for (510); id est,<br />

by (1563), and<br />

f ⋆ (511) ≥ f ⋆ (510) ≥ 0 (515)<br />

lim<br />

ǫ→0 +f⋆ (511) = 0 (516)<br />

3.11 A convex problem has convex feasible set, and the objective surface has one and only<br />

one global minimum.

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

Saved successfully!

Ooh no, something went wrong!