Ludium
Sign In
Linear Algebra Done Right
Vector Spaces
01Complex Numbers as Pairs of Real Numbers02Lists and Fⁿ: Coordinate Addition and a First Proof03Vectors in R², Scalar Multiplication, and FieldsProblem set0/10Problem set 20/10
01Definition of a Vector Space: Fⁿ, F^∞ and Function Spaces02Unique Zero and Inverses, 0v = 0, and the Subspace Test03Sums of Subspaces: The Smallest Subspace Containing Them04Direct Sums: Only One Way to Write ZeroProblem set0/10Problem set 20/10
Finite-Dimensional Spaces
01Linear Combinations, Span, and Finite-Dimensional Spaces02Linear Independence and the Linear Dependence Lemma03Independent Lists Are No Longer Than Spanning ListsProblem set0/10Problem set 20/10
01Bases: Reduce a Spanning List, Extend an Independent List02Dimension: Why Every Basis Has the Same Length03The Dimension of a Sum of Two SubspacesProblem set0/10Problem set 20/10
Linear Maps
01Linear Maps and the Linear Map Lemma: Values on a Basis02Sums, Scalar Multiples and Products of Linear Maps03Null Space, Range, Injectivity and Surjectivity04The Fundamental Theorem of Linear Maps and Linear SystemsProblem set0/10Problem set 20/10
01The Matrix of a Linear Map: Each Column Is Tvₖ02Matrix Operations and Why Multiplication Is Row Times Column03Rows, Columns and Linear Combinations in Matrix Products04Column Rank Equals Row Rank: Column-Row FactorizationProblem set0/10Problem set 20/10
01Invertible Linear Maps: Injective and Surjective02Equal Dimensions: Injective or Surjective Means Invertible03Isomorphic Spaces and the Dimension of L(V, W)04The Matrix of a Vector: M(Tv) = M(T)M(v)05Change of Basis for an Operator: A = C⁻¹BCProblem set0/10Problem set 20/10
01Products of Vector Spaces: When Dimensions Add02Quotient Spaces V/U: Parallel Lines as Vectors03The Quotient Map and V/null T ≅ range TProblem set0/10Problem set 20/10
01The Dual Basis Reads Coefficients and T′ Reverses Arrows02The Annihilator U⁰ and dim U⁰ = dim V − dim U03Null T′ and Range T′ as Annihilators; M(T′) Is the TransposeProblem set0/10Problem set 20/10
Polynomials and Eigenvalues
01Conjugates, |z| and the Complex Triangle Inequality02Zeros, Factors and the Division Algorithm via a Basis03Fundamental Theorem of Algebra and Complete Factorization04Real Polynomials: Linear and Irreducible Quadratic FactorsProblem set0/10Problem set 20/10
01Invariant Subspaces and Eigenvalues: Tv = λv02Eigenvectors, T − λI and At Most dim V Eigenvalues03p(T) and Why Every Complex Operator Has an Eigenvalue04The Minimal Polynomial: Existence, Uniqueness and Computation05Eigenvalues Are the Zeros of the Minimal Polynomial06Why Odd Dimension Forces a Real EigenvalueProblem set0/10Problem set 20/10
01Upper-Triangular Matrices and Invariant Spans02Eigenvalues Sit on the Upper-Triangular Diagonal03The Minimal Polynomial Decides Upper-TriangularityProblem set0/10Problem set 20/10
01Diagonalizable Operators and Direct Sums of Eigenspaces02Four Equivalent Tests for Diagonalizability03Distinct Linear Factors Decide Diagonalizability04The Gershgorin Disk Theorem: Where Eigenvalues Live05Commuting Operators Share a Diagonalizing Basis06Commuting Operators Over C: One Upper-Triangular BasisProblem set0/10Problem set 20/10
Inner Product Spaces
01Inner Products: Why Complex Spaces Need a Conjugate02Norms, Orthogonality and the Pythagorean Theorem03Cauchy-Schwarz and the Triangle InequalityProblem set0/10Problem set 20/10
01Orthonormal Bases: Bessel's Inequality and Parseval02Gram-Schmidt, Orthonormal Bases and Schur's Theorem03Linear Functionals and the Riesz Representation TheoremProblem set0/10Problem set 20/10
01Orthogonal Complements and the Direct Sum V = U ⊕ U⊥02Orthogonal Projection onto U and Why U⊥⊥ = U03Riesz Representation Without a Basis: v ⊥ null φProblem set0/10Problem set 20/10
01The Closest Point Is the Projection: Fitting Sine on [−π, π]02The Pseudoinverse T† and Its Three Inverse-Like Properties03T†w: The Best Approximate Solution of Smallest NormProblem set0/10Problem set 20/10
Operators on Inner Product Spaces
01The Adjoint T*: Moving T into the Second Slot of ⟨Tv, w⟩02Null T* = (range T)⊥ and M(T*) as the Conjugate Transpose03Self-Adjoint Operators: Real Eigenvalues and Real ⟨Tv, v⟩04Normal Operators: ‖Tv‖ = ‖T*v‖ and Orthogonal EigenvectorsProblem set0/10Problem set 20/10
01Real Spectral Theorem: Self-Adjoint ⇔ Orthonormal Eigenbasis02The Complex Spectral Theorem: Normal Is Enough03Positive Operators: Six Equivalent Characterizations04The Unique Positive Square Root √TProblem set0/10Problem set 20/10
01Isometries: Five Equivalent Ways to Preserve Length02Unitary Operators: S⁻¹ = S*, Eigenvalues on |λ| = 103QR Factorization: Gram–Schmidt Recorded as A = QR04Cholesky Factorization: Positive Definite B = R*RProblem set0/10Problem set 20/10
01Singular Values: Square Roots of the Eigenvalues of T*T02Singular Value Decomposition: Rank, Isometries, Proof03Diagonal Matrices, Adjoints and Pseudoinverses via SVD04Matrix SVD: Computing It and Factoring A = BDC*Problem set0/10Problem set 20/10

Matrix SVD: Computing It and Factoring A = BDC*

How do orthonormal eigenvectors of T*T produce an SVD, and why does every rank-m matrix factor as A = BDC*?


Loading…

←Previous Diagonal Matrices, Adjoints and Pseudoinverses via SVD

Your summary note

    1. 1

      Finding an SVD from eigenvectors of T∗TT^*TT∗T

      State the recipe, then carry it through for T(x1,x2,x3,x4)=(−5x4,0,x1+x2)T(x_1,x_2,x_3,x_4) = (-5x_4, 0, x_1+x_2)T(x1​,x2​,x3​,x4​)=(−5x4​,0,x1​+x2​): unit eigenvectors e1,e2e_1, e_2e1​,e2​ for eigenvalues 25,225, 225,2, then f1,f2f_1, f_2f1​,f2​ from fk=Tek/skf_k = Te_k/s_kfk​=Tek​/sk​.

    2. 2

      Matrix SVD: statement and the three factors

      State the theorem for a ppp-by-nnn matrix AAA of rank mmm, giving the sizes and properties of BBB, DDD and CCC, and build them from f1,…,fmf_1,\dots,f_mf1​,…,fm​, s1,…,sms_1,\dots,s_ms1​,…,sm​ and e1,…,eme_1,\dots,e_me1​,…,em​.

    3. 3

      Proof that A=BDC∗A = BDC^*A=BDC∗

      Show AC=BDAC = BDAC=BD column by column, multiply by C∗C^*C∗, then prove ACC∗=AACC^* = AACC∗=A by checking agreement on span(e1,…,em)\text{span}(e_1,\dots,e_m)span(e1​,…,em​) and on its orthogonal complement.

    4. 4

      Storage cost of a low-rank factorization

      Count the m(p+m+n)m(p+m+n)m(p+m+n) entries of BBB, DDD and CCC against the pnpnpn entries of AAA, and record when the factorization saves storage.

    Attempt 1 of 2