v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
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
- Page 511 and 512: 7.1. FIRST PREVALENT PROBLEM: 511 w
- Page 513 and 514: 7.1. FIRST PREVALENT PROBLEM: 513 T
- Page 515 and 516: 7.2. SECOND PREVALENT PROBLEM: 515
- Page 517 and 518: 7.2. SECOND PREVALENT PROBLEM: 517
- Page 519 and 520: 7.2. SECOND PREVALENT PROBLEM: 519
- Page 521 and 522: 7.2. SECOND PREVALENT PROBLEM: 521
- Page 523 and 524: 7.2. SECOND PREVALENT PROBLEM: 523
- Page 525 and 526: 7.3. THIRD PREVALENT PROBLEM: 525 g
- Page 527 and 528: 7.3. THIRD PREVALENT PROBLEM: 527 w
- Page 529 and 530: 7.3. THIRD PREVALENT PROBLEM: 529 7
- Page 531 and 532: 7.3. THIRD PREVALENT PROBLEM: 531 7
- Page 533 and 534: 7.3. THIRD PREVALENT PROBLEM: 533 O
- Page 535 and 536: 7.4. CONCLUSION 535 The rank constr
- Page 537 and 538: Appendix A Linear algebra A.1 Main-
- Page 539 and 540: A.1. MAIN-DIAGONAL δ OPERATOR, λ
- Page 541 and 542: A.1. MAIN-DIAGONAL δ OPERATOR, λ
- Page 543 and 544: A.2. SEMIDEFINITENESS: DOMAIN OF TE
- Page 545 and 546: A.3. PROPER STATEMENTS 545 (AB) T
- Page 547 and 548: A.3. PROPER STATEMENTS 547 A.3.1 Se
- Page 549 and 550: A.3. PROPER STATEMENTS 549 For A di
- Page 551 and 552: A.3. PROPER STATEMENTS 551 Diagonal
- Page 553 and 554: A.3. PROPER STATEMENTS 553 For A,B
- Page 555 and 556: A.3. PROPER STATEMENTS 555 A.3.1.0.
- Page 557 and 558: A.4. SCHUR COMPLEMENT 557 A.4 Schur
- Page 559 and 560: A.4. SCHUR COMPLEMENT 559 A.4.0.0.2
- Page 561: A.5. EIGEN DECOMPOSITION 561 When B
- Page 565 and 566: A.6. SINGULAR VALUE DECOMPOSITION,
- Page 567 and 568: A.6. SINGULAR VALUE DECOMPOSITION,
- Page 569 and 570: A.7. ZEROS 569 Given symmetric matr
- Page 571 and 572: A.7. ZEROS 571 (TRANSPOSE.) Likewis
- Page 573 and 574: A.7. ZEROS 573 For X,A∈ S M + [31
- Page 575 and 576: A.7. ZEROS 575 A.7.5.0.1 Propositio
- Page 577 and 578: Appendix B Simple matrices Mathemat
- Page 579 and 580: B.1. RANK-ONE MATRIX (DYAD) 579 R(v
- Page 581 and 582: B.1. RANK-ONE MATRIX (DYAD) 581 ran
- Page 583 and 584: B.2. DOUBLET 583 R([u v ]) R(Π)= R
- Page 585 and 586: B.3. ELEMENTARY MATRIX 585 If λ
- Page 587 and 588: B.4. AUXILIARY V -MATRICES 587 the
- Page 589 and 590: B.4. AUXILIARY V -MATRICES 589 18.
- Page 591 and 592: B.5. ORTHOGONAL MATRIX 591 B.5 Orth
- Page 593 and 594: B.5. ORTHOGONAL MATRIX 593 Figure 1
- Page 595 and 596: Appendix C Some analytical optimal
- Page 597 and 598: C.2. TRACE, SINGULAR AND EIGEN VALU
- Page 599 and 600: C.2. TRACE, SINGULAR AND EIGEN VALU
- Page 601 and 602: C.2. TRACE, SINGULAR AND EIGEN VALU
- Page 603 and 604: C.3. ORTHOGONAL PROCRUSTES PROBLEM
- Page 605 and 606: C.4. TWO-SIDED ORTHOGONAL PROCRUSTE
- Page 607 and 608: C.4. TWO-SIDED ORTHOGONAL PROCRUSTE
- Page 609 and 610: Appendix D Matrix calculus From too
- Page 611 and 612: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
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