10.03.2015 Views

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

4.6. CARDINALITY AND RANK CONSTRAINT EXAMPLES 323<br />

4.6.0.0.9 Example. Cardinality/rank problem.<br />

d’Aspremont et alii [83] propose approximating a positive semidefinite matrix<br />

A∈ S N + by a rank-one matrix having a constraint on cardinality c : for<br />

0 < c < N<br />

minimize ‖A − zz T ‖ F<br />

z<br />

(739)<br />

subject to cardz ≤ c<br />

which, they explain, is a hard problem equivalent to<br />

maximize x T Ax<br />

x<br />

subject to ‖x‖ = 1<br />

cardx ≤ c<br />

(740)<br />

where z ∆ = √ λx and where optimal solution x ⋆ is a principal eigenvector<br />

(1575) (A.5) of A and λ = x ⋆T Ax ⋆ is the principal eigenvalue [134, p.331]<br />

when c is true cardinality of that eigenvector. This is principal component<br />

analysis with a cardinality constraint which controls solution sparsity. Define<br />

the matrix variable<br />

X ∆ = xx T ∈ S N (741)<br />

whose desired rank is 1, and whose desired diagonal cardinality<br />

cardδ(X) ≡ cardx (742)<br />

is equivalent to cardinality c of vector x . Then we can transform cardinality<br />

problem (740) to an equivalent problem in new variable X : 4.38<br />

maximize<br />

X∈S N 〈X , A〉<br />

subject to 〈X , I〉 = 1<br />

(X ≽ 0)<br />

rankX = 1<br />

cardδ(X) ≤ c<br />

(743)<br />

We transform problem (743) to an equivalent convex problem by<br />

introducing two direction matrices into regularization terms: W to achieve<br />

4.38 A semidefiniteness constraint X ≽ 0 is not required, theoretically, because positive<br />

semidefiniteness of a rank-1 matrix is enforced by symmetry. (Theorem A.3.1.0.7)

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

Saved successfully!

Ooh no, something went wrong!