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.

A.6. SINGULAR VALUE DECOMPOSITION, SVD 565<br />

Square matrix Σ is diagonal (A.1.1)<br />

δ 2 (Σ) = Σ ∈ R η×η (1452)<br />

holding the singular values {σ i ∈ R} of A which are always arranged in<br />

nonincreasing order by convention and are related to eigenvalues by A.14<br />

⎧<br />

⎨<br />

σ(A) i = σ(A T ) i =<br />

⎩<br />

√<br />

λ(AT A) i = √ (√<br />

λ(AA T ) i = λ AT A)<br />

i<br />

)<br />

= λ(√<br />

AA<br />

T<br />

> 0, 1 ≤ i ≤ ρ<br />

i<br />

0, ρ < i ≤ η<br />

(1453)<br />

of which the last η −ρ are 0 , A.15 where<br />

ρ ∆ = rankA = rank Σ (1454)<br />

A point sometimes lost: Any real matrix may be decomposed in terms of<br />

its real singular values σ(A) ∈ R η and real matrices U and Q as in (1450),<br />

where [134,2.5.3]<br />

R{u i |σ i ≠0} = R(A)<br />

R{u i |σ i =0} ⊆ N(A T )<br />

R{q i |σ i ≠0} = R(A T (1455)<br />

)<br />

R{q i |σ i =0} ⊆ N(A)<br />

A.6.2<br />

Subcompact SVD<br />

Some authors allow only nonzero singular values. In that case the compact<br />

decomposition can be made smaller; it can be redimensioned in terms of rank<br />

ρ because, for any A∈ R m×n<br />

ρ = rankA = rank Σ = max {i∈{1... η} | σ i ≠ 0} ≤ η (1456)<br />

There are η singular values. For any flavor SVD, rank is equivalent to<br />

the number of nonzero singular values on the main diagonal of Σ.<br />

A.14 When A is normal,<br />

√<br />

σ(A) = |λ(A)|. [344,8.1]<br />

)<br />

A.15 For η = n , σ(A) = λ(A T A) = λ(√<br />

AT A where λ denotes eigenvalues.<br />

√<br />

)<br />

For η = m , σ(A) = λ(AA T ) = λ(√<br />

AA<br />

T<br />

.

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

Saved successfully!

Ooh no, something went wrong!