v2010.10.26 - Convex Optimization
v2010.10.26 - Convex Optimization v2010.10.26 - Convex Optimization
56 CHAPTER 2. CONVEX GEOMETRYConcurrently, consider injective linear operator Py=A † y : R m → R nwhere R(A † )= R(A T ). P(Ax)= PTx achieves projection of vector x onthe row space R(A T ). (E.3.1) This means vector Ax can be succinctlyinterpreted as coefficients of orthogonal projection.Pseudoinverse matrix A † is skinny and full-rank, so operator Py is a linearbijection with respect to its range R(A † ). By Definition 2.2.1.0.1, imageP(T B) of projection PT(B) on R(A T ) in R n must therefore be isomorphicwith the set of projection coefficients T B = {Ax |x∈ B} in R m and have thesame affine dimension by (49). To illustrate, we present a three-dimensionalEuclidean body B in Figure 19 where any point x in the nullspace N(A)maps to the origin.2.2.2 Symmetric matrices2.2.2.0.1 Definition. Symmetric matrix subspace.Define a subspace of R M×M : the convex set of all symmetric M×M matrices;S M { A∈ R M×M | A=A T} ⊆ R M×M (50)This subspace comprising symmetric matrices S M is isomorphic with thevector space R M(M+1)/2 whose dimension is the number of free variables in asymmetric M ×M matrix. The orthogonal complement [331] [250] of S M isS M⊥ { A∈ R M×M | A=−A T} ⊂ R M×M (51)the subspace of antisymmetric matrices in R M×M ; id est,S M ⊕ S M⊥ = R M×M (52)where unique vector sum ⊕ is defined on page 772.△All antisymmetric matrices are hollow by definition (have 0 maindiagonal). Any square matrix A∈ R M×M can be written as a sum of itssymmetric and antisymmetric parts: respectively,A = 1 2 (A +AT ) + 1 2 (A −AT ) (53)The symmetric part is orthogonal in R M2 to the antisymmetric part; videlicet,tr ( (A +A T )(A −A T ) ) = 0 (54)
2.2. VECTORIZED-MATRIX INNER PRODUCT 57In the ambient space of real matrices, the antisymmetric matrix subspacecan be described{ }1S M⊥ =2 (A −AT ) | A∈ R M×M ⊂ R M×M (55)because any matrix in S M is orthogonal to any matrix in S M⊥ . Furtherconfined to the ambient subspace of symmetric matrices, because ofantisymmetry, S M⊥ would become trivial.2.2.2.1 Isomorphism of symmetric matrix subspaceWhen a matrix is symmetric in S M , we may still employ the vectorizationtransformation (37) to R M2 ; vec , an isometric isomorphism. We mightinstead choose to realize in the lower-dimensional subspace R M(M+1)/2 byignoring redundant entries (below the main diagonal) during transformation.Such a realization would remain isomorphic but not isometric. Lack ofisometry is a spatial distortion due now to disparity in metric between R M2and R M(M+1)/2 . To realize isometrically in R M(M+1)/2 , we must make acorrection: For Y = [Y ij ]∈ S M we take symmetric vectorization [215,2.2.1]⎡ ⎤√2Y12Y 11Y 22√2Y13svec Y √2Y23∈ R M(M+1)/2 (56)⎢ Y 33 ⎥⎣ ⎦.Y MMwhere all entries off the main diagonal have been scaled. Now for Z ∈ S M〈Y , Z〉 tr(Y T Z) = vec(Y ) T vec Z = svec(Y ) T svec Z (57)Then because the metrics become equivalent, for X ∈ S M‖ svec X − svec Y ‖ 2 = ‖X − Y ‖ F (58)and because symmetric vectorization (56) is a linear bijective mapping, thensvec is an isometric isomorphism of the symmetric matrix subspace. In other
- Page 5 and 6: for Jennie Columba♦Antonio♦♦&
- Page 7 and 8: PreludeThe constant demands of my d
- Page 9 and 10: Convex Optimization&Euclidean Dista
- Page 11 and 12: CONVEX OPTIMIZATION & EUCLIDEAN DIS
- Page 13 and 14: List of Figures1 Overview 211 Orion
- Page 15 and 16: LIST OF FIGURES 1562 Shrouded polyh
- Page 17: LIST OF FIGURES 17130 Elliptope E 3
- Page 21 and 22: Chapter 1OverviewConvex Optimizatio
- Page 23 and 24: ˇx 4ˇx 3ˇx 2Figure 2: Applicatio
- Page 25 and 26: 25Figure 4: This coarsely discretiz
- Page 27 and 28: (biorthogonal expansion) is examine
- Page 29 and 30: 29cardinality Boolean solution to a
- Page 31 and 32: 31Figure 8: Robotic vehicles in con
- Page 33 and 34: an elaborate exposition offering in
- Page 35 and 36: Chapter 2Convex geometryConvexity h
- Page 37 and 38: 2.1. CONVEX SET 372.1.2 linear inde
- Page 39 and 40: 2.1. CONVEX SET 392.1.6 empty set v
- Page 41 and 42: 2.1. CONVEX SET 41(a)R(b)R 2(c)R 3F
- Page 43 and 44: 2.1. CONVEX SET 43where Q∈ R 3×3
- Page 45 and 46: 2.1. CONVEX SET 45Now let’s move
- Page 47 and 48: 2.1. CONVEX SET 47By additive inver
- Page 49 and 50: 2.1. CONVEX SET 49R nR mR(A T )x px
- Page 51 and 52: 2.2. VECTORIZED-MATRIX INNER PRODUC
- Page 53 and 54: 2.2. VECTORIZED-MATRIX INNER PRODUC
- Page 55: 2.2. VECTORIZED-MATRIX INNER PRODUC
- Page 59 and 60: 2.2. VECTORIZED-MATRIX INNER PRODUC
- Page 61 and 62: 2.3. HULLS 61Figure 20: Convex hull
- Page 63 and 64: 2.3. HULLS 63Aaffine hull (drawn tr
- Page 65 and 66: 2.3. HULLS 65subset of the affine h
- Page 67 and 68: 2.3. HULLS 672.3.2.0.2 Example. Nuc
- Page 69 and 70: 2.3. HULLS 692.3.2.0.3 Exercise. Co
- Page 71 and 72: 2.3. HULLS 71Figure 24: A simplicia
- Page 73 and 74: 2.4. HALFSPACE, HYPERPLANE 73H + =
- Page 75 and 76: 2.4. HALFSPACE, HYPERPLANE 7511−1
- Page 77 and 78: 2.4. HALFSPACE, HYPERPLANE 772.4.2.
- Page 79 and 80: 2.4. HALFSPACE, HYPERPLANE 792.4.2.
- Page 81 and 82: 2.4. HALFSPACE, HYPERPLANE 81tradit
- Page 83 and 84: 2.4. HALFSPACE, HYPERPLANE 832.4.2.
- Page 85 and 86: 2.5. SUBSPACE REPRESENTATIONS 852.5
- Page 87 and 88: 2.5. SUBSPACE REPRESENTATIONS 87(Ex
- Page 89 and 90: 2.5. SUBSPACE REPRESENTATIONS 89(a)
- Page 91 and 92: 2.5. SUBSPACE REPRESENTATIONS 91are
- Page 93 and 94: 2.6. EXTREME, EXPOSED 93A one-dimen
- Page 95 and 96: 2.6. EXTREME, EXPOSED 952.6.1.1 Den
- Page 97 and 98: 2.7. CONES 97X(a)00(b)XFigure 33: (
- Page 99 and 100: 2.7. CONES 99XXFigure 37: Truncated
- Page 101 and 102: 2.7. CONES 101Figure 39: Not a cone
- Page 103 and 104: 2.7. CONES 103cone that is a halfli
- Page 105 and 106: 2.7. CONES 105A pointed closed conv
2.2. VECTORIZED-MATRIX INNER PRODUCT 57In the ambient space of real matrices, the antisymmetric matrix subspacecan be described{ }1S M⊥ =2 (A −AT ) | A∈ R M×M ⊂ R M×M (55)because any matrix in S M is orthogonal to any matrix in S M⊥ . Furtherconfined to the ambient subspace of symmetric matrices, because ofantisymmetry, S M⊥ would become trivial.2.2.2.1 Isomorphism of symmetric matrix subspaceWhen a matrix is symmetric in S M , we may still employ the vectorizationtransformation (37) to R M2 ; vec , an isometric isomorphism. We mightinstead choose to realize in the lower-dimensional subspace R M(M+1)/2 byignoring redundant entries (below the main diagonal) during transformation.Such a realization would remain isomorphic but not isometric. Lack ofisometry is a spatial distortion due now to disparity in metric between R M2and R M(M+1)/2 . To realize isometrically in R M(M+1)/2 , we must make acorrection: For Y = [Y ij ]∈ S M we take symmetric vectorization [215,2.2.1]⎡ ⎤√2Y12Y 11Y 22√2Y13svec Y √2Y23∈ R M(M+1)/2 (56)⎢ Y 33 ⎥⎣ ⎦.Y MMwhere all entries off the main diagonal have been scaled. Now for Z ∈ S M〈Y , Z〉 tr(Y T Z) = vec(Y ) T vec Z = svec(Y ) T svec Z (57)Then because the metrics become equivalent, for X ∈ S M‖ svec X − svec Y ‖ 2 = ‖X − Y ‖ F (58)and because symmetric vectorization (56) is a linear bijective mapping, thensvec is an isometric isomorphism of the symmetric matrix subspace. In other