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.

3.1. CONVEX FUNCTION 215<br />

3.1.7.0.2 Example. Matrix pseudofractional function.<br />

Consider a real function of two variables<br />

f(A, x) : S n × R n → R = x T A † x (498)<br />

on domf = S+× n R(A). This function is convex simultaneously in both<br />

variables when variable matrix A belongs to the entire positive semidefinite<br />

cone S n + and variable vector x is confined to range R(A) of matrix A .<br />

To explain this, we need only demonstrate that the function epigraph is<br />

convex. Consider Schur-form (1410) fromA.4: for t ∈ R<br />

[ ] A z<br />

G(A, z , t) =<br />

z T ≽ 0<br />

t<br />

⇔<br />

z T (I − AA † ) = 0<br />

t − z T A † z ≥ 0<br />

A ≽ 0<br />

(499)<br />

Inverse image of the positive semidefinite cone S n+1<br />

+ under affine mapping<br />

G(A, z , t) is convex by Theorem 2.1.9.0.1. Of the equivalent conditions for<br />

positive semidefiniteness of G , the first is an equality demanding vector z<br />

belong to R(A). Function f(A, z)=z T A † z is convex on convex domain<br />

S+× n R(A) because the Cartesian product constituting its epigraph<br />

epi f(A, z) = { (A, z , t) | A ≽ 0, z ∈ R(A), z T A † z ≤ t } = G −1( S n+1<br />

+<br />

is convex.<br />

)<br />

(500)<br />

3.1.7.0.3 Exercise. Matrix product function.<br />

Continue Example 3.1.7.0.2 by introducing vector variable x and making<br />

the substitution z ←Ax . Because of matrix symmetry (E), for all x∈ R n<br />

f(A, z(x)) = z T A † z = x T A T A † Ax = x T Ax = f(A, x) (501)<br />

whose epigraph is<br />

epif(A, x) = { (A, x, t) | A ≽ 0, x T Ax ≤ t } (502)<br />

Provide two simple explanations why f(A, x) = x T Ax is not a function<br />

convex simultaneously in positive semidefinite matrix A and vector x on<br />

domf = S+× n R n .

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

Saved successfully!

Ooh no, something went wrong!