v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization

convexoptimization.com
from convexoptimization.com More from this publisher
10.03.2015 Views

562 APPENDIX A. LINEAR ALGEBRA A.5.0.0.1 Theorem. Real eigenvector. Eigenvectors of a real matrix corresponding to real eigenvalues must be real. ⋄ Proof. Ax = λx . Given λ=λ ∗ , x H Ax = λx H x = λ‖x‖ 2 = x T Ax ∗ x = x ∗ , where x H =x ∗T . The converse is equally simple. ⇒ A.5.0.1 Uniqueness From the fundamental theorem of algebra it follows: eigenvalues, including their multiplicity, for a given square matrix are unique; meaning, there is no other set of eigenvalues for that matrix. (Conversely, many different matrices may share the same unique set of eigenvalues; e.g., for any X , λ(X) = λ(X T ).) Uniqueness of eigenvectors, in contrast, disallows multiplicity of the same direction: A.5.0.1.1 Definition. Unique eigenvectors. When eigenvectors are unique, we mean: unique to within a real nonzero scaling, and their directions are distinct. △ If S is a matrix of eigenvectors of X as in (1438), for example, then −S is certainly another matrix of eigenvectors decomposing X with the same eigenvalues. For any square matrix, the eigenvector corresponding to a distinct eigenvalue is unique; [285, p.220] distinct eigenvalues ⇒ eigenvectors unique (1442) Eigenvectors corresponding to a repeated eigenvalue are not unique for a diagonalizable matrix; repeated eigenvalue ⇒ eigenvectors not unique (1443) Proof follows from the observation: any linear combination of distinct eigenvectors of diagonalizable X , corresponding to a particular eigenvalue, produces another eigenvector. For eigenvalue λ whose multiplicity A.13 A.13 A matrix is diagonalizable iff algebraic multiplicity (number of occurrences of same eigenvalue) equals geometric multiplicity dim N(X −λI) = m − rank(X −λI) [285, p.15] (number of Jordan blocks w.r.t λ or number of corresponding l.i. eigenvectors).

A.5. EIGEN DECOMPOSITION 563 dim N(X −λI) exceeds 1, in other words, any choice of independent vectors from N(X −λI) (of the same multiplicity) constitutes eigenvectors corresponding to λ . Caveat diagonalizability insures linear independence which implies existence of distinct eigenvectors. We may conclude, for diagonalizable matrices, distinct eigenvalues ⇔ eigenvectors unique (1444) A.5.1 eigenmatrix The (right-)eigenvectors {s i } (1438) are naturally orthogonal w T i s j =0 to left-eigenvectors {w i } except, for i=1... m , w T i s i =1 ; called a biorthogonality condition [318,2.2.4] [176] because neither set of left or right eigenvectors is necessarily an orthogonal set. Consequently, each dyad from a diagonalization is an independent (B.1.1) nonorthogonal projector because s i w T i s i w T i = s i w T i (1445) (whereas the dyads of singular value decomposition are not inherently projectors (confer (1450))). The dyads of eigen decomposition can be termed eigenmatrices because Sum of the eigenmatrices is the identity; X s i w T i = λ i s i w T i (1446) m∑ s i wi T = I (1447) i=1 A.5.2 Symmetric matrix diagonalization The set of normal matrices is, precisely, that set of all real matrices having a complete orthonormal set of eigenvectors; [344,8.1] [289, prob.10.2.31] id est, any matrix X for which XX T = X T X ; [134,7.1.3] [285, p.3] e.g., orthogonal and circulant matrices [142]. All normal matrices are diagonalizable. A symmetric matrix is a special normal matrix whose

A.5. EIGEN DECOMPOSITION 563<br />

dim N(X −λI) exceeds 1, in other words, any choice of independent<br />

vectors from N(X −λI) (of the same multiplicity) constitutes eigenvectors<br />

corresponding to λ .<br />

<br />

Caveat diagonalizability insures linear independence which implies<br />

existence of distinct eigenvectors. We may conclude, for diagonalizable<br />

matrices,<br />

distinct eigenvalues ⇔ eigenvectors unique (1444)<br />

A.5.1<br />

eigenmatrix<br />

The (right-)eigenvectors {s i } (1438) are naturally orthogonal w T i s j =0<br />

to left-eigenvectors {w i } except, for i=1... m , w T i s i =1 ; called a<br />

biorthogonality condition [318,2.2.4] [176] because neither set of left or right<br />

eigenvectors is necessarily an orthogonal set. Consequently, each dyad from a<br />

diagonalization is an independent (B.1.1) nonorthogonal projector because<br />

s i w T i s i w T i = s i w T i (1445)<br />

(whereas the dyads of singular value decomposition are not inherently<br />

projectors (confer (1450))).<br />

The dyads of eigen decomposition can be termed eigenmatrices because<br />

Sum of the eigenmatrices is the identity;<br />

X s i w T i = λ i s i w T i (1446)<br />

m∑<br />

s i wi T = I (1447)<br />

i=1<br />

A.5.2<br />

Symmetric matrix diagonalization<br />

The set of normal matrices is, precisely, that set of all real matrices having<br />

a complete orthonormal set of eigenvectors; [344,8.1] [289, prob.10.2.31]<br />

id est, any matrix X for which XX T = X T X ; [134,7.1.3] [285, p.3]<br />

e.g., orthogonal and circulant matrices [142]. All normal matrices are<br />

diagonalizable. A symmetric matrix is a special normal matrix whose

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

Saved successfully!

Ooh no, something went wrong!