Reference sheet

Linear Algebra II reference sheet

The main definitions, theorems and factorizations of Linear Algebra II with their hypotheses, unit by unit, with links to the lessons that explain each topic in full.

For Linear Algebra II

This sheet follows the ten units of the course in order. Each unit starts with the big idea in a few sentences, then gives the definitions, theorems and algorithms you will look up most often, each with the hypotheses it needs. When you want a proof or more examples, follow the lesson links under each topic.

Every result names its field F\mathbb{F}, which is R\mathbb{R} or C\mathbb{C}, and all vector spaces are finite-dimensional unless a statement says otherwise. Vectors such as v,w,xv, w, x are written in plain italic, coordinates are columns, and a matrix A∈Fm×nA \in \mathbb{F}^{m \times n} represents a linear map T:V→WT: V \to W from an nn-dimensional space to an mm-dimensional one. The linear maps from VV to WW form the space Hom⁡(V,W)\operatorname{Hom}(V, W). Coordinates change by [v]C=PC←B[v]B[v]_{\mathcal{C}} = P_{\mathcal{C} \leftarrow \mathcal{B}}[v]_{\mathcal{B}}, and the matrix [T]D←B[T]_{\mathcal{D} \leftarrow \mathcal{B}} has as its jjth column the D\mathcal{D} coordinates of T(bj)T(b_j).

The transpose is ATA^{\mathsf{T}} and the conjugate transpose is A∗=A‾TA^{*} = \overline{A}^{\mathsf{T}}. Inner products are conjugate-linear in the first argument and linear in the second, with the standard inner product ⟨x,y⟩=x∗y\langle x, y \rangle = x^{*}y on Fn\mathbb{F}^n and ∥x∥=⟨x,x⟩\lVert x \rVert = \sqrt{\langle x, x \rangle}. The algebraic dual is V∨V^{\vee} and the dual map is T∨T^{\vee}, while T∗T^{*} always means the adjoint, which needs an inner product. The space Pd\mathcal{P}_d holds the polynomials of degree at most dd, including the zero polynomial. The characteristic polynomial is χA(t)=det⁡(tI−A)\chi_A(t) = \det(tI - A) and the minimal polynomial μA\mu_A is monic. Singular values are listed in decreasing order σ1≥⋯≥σp≥0\sigma_1 \ge \cdots \ge \sigma_p \ge 0 with p=min⁡(m,n)p = \min(m, n), while the eigenvalues of a Hermitian matrix are listed in increasing order λ1≤⋯≤λn\lambda_1 \le \cdots \le \lambda_n whenever the order matters. A real symmetric matrix counts as Hermitian.

Vector Spaces and Linear Maps

The first unit proves the facts that every later argument relies on. It sets out how to read and write a proof, fixes the field of scalars, shows that all bases of a space have the same size, and builds the space of linear maps, the rank-nullity theorem and direct sums with their projections.

Proofs and quantifiers

A universal claim ∀x∈S, P(x)\forall x \in S,\ P(x) is proved by an argument about an arbitrary member of SS and disproved by one counterexample. An existence claim ∃x∈S, P(x)\exists x \in S,\ P(x) is proved by producing one member with the property, or by an argument that one must exist. An implication P⇒QP \Rightarrow Q is false only when PP holds and QQ fails. Negation switches each quantifier.

StatementNegation
∀x∈S, P(x)\forall x \in S,\ P(x)∃x∈S, ¬P(x)\exists x \in S,\ \neg P(x)
∃x∈S, P(x)\exists x \in S,\ P(x)∀x∈S, ¬P(x)\forall x \in S,\ \neg P(x)
P⇒QP \Rightarrow QPP and ¬Q\neg Q
∀x∈S, (P(x)⇒Q(x))\forall x \in S,\ (P(x) \Rightarrow Q(x))some x∈Sx \in S has P(x)P(x) and ¬Q(x)\neg Q(x)
  • Converse and contrapositive. The converse of P⇒QP \Rightarrow Q is Q⇒PQ \Rightarrow P, a different statement. The contrapositive ¬Q⇒¬P\neg Q \Rightarrow \neg P is equivalent to the original. A statement P⇔QP \Leftrightarrow Q needs a proof of both directions.
  • Induction. If the base case P(k0)P(k_0) holds and P(k)P(k) implies P(k+1)P(k + 1) for every integer k≥k0k \ge k_0, then P(k)P(k) holds for every integer k≥k0k \ge k_0. Checking several values of kk proves nothing about the others.

For example, if N2=0N^2 = 0, then (I+N)k=I+kN(I + N)^k = I + kN for every integer k≥0k \ge 0. The base case is (I+N)0=I(I + N)^0 = I, and the step is

(I+kN)(I+N)=I+(k+1)N+kN2=I+(k+1)N.\begin{aligned} (I + kN)(I + N) &= I + (k + 1)N + kN^2 \\[0.4em] &= I + (k + 1)N. \end{aligned}

Lessons Proof Strategies in Linear Algebra

Fields, vector spaces and subspaces

A complex number z=a+biz = a + bi with a,b∈Ra, b \in \mathbb{R} has conjugate z‾=a−bi\overline{z} = a - bi and modulus ∣z∣=a2+b2=zz‾\lvert z \rvert = \sqrt{a^2 + b^2} = \sqrt{z\overline{z}}. Conjugation respects sums and products, ∣zw∣=∣z∣∣w∣\lvert zw \rvert = \lvert z \rvert \lvert w \rvert, and for z≠0z \ne 0

1z=z‾∣z∣2.\frac{1}{z} = \frac{\overline{z}}{\lvert z \rvert^2}.

A vector space over F\mathbb{F} allows scalars only from F\mathbb{F}. Its axioms ask that addition is commutative and associative with a zero vector and negatives, that a(bv)=(ab)va(bv) = (ab)v and 1v=v1v = v, and that both distributive laws hold. Everything proved from the axioms alone, such as 0v=00v = 0 and (−1)v=−v(-1)v = -v, holds in every vector space.

  • Subspace test. A subset UU of VV is a subspace exactly when it is nonempty and au+bv∈Uau + bv \in U for all u,v∈Uu, v \in U and a,b∈Fa, b \in \mathbb{F}. The line x+y=1x + y = 1 in R2\mathbb{R}^2 fails because it misses 00, and the union of the two axes fails because (1,0)+(0,1)(1, 0) + (0, 1) lies on neither.
  • Real and complex dimension. If v1,…,vnv_1, \dots, v_n is a basis of a complex space VV, then v1,iv1,…,vn,ivnv_1, iv_1, \dots, v_n, iv_n is a basis of VV as a real space, so dim⁡RV=2dim⁡CV\dim_{\mathbb{R}} V = 2\dim_{\mathbb{C}} V. A statement proved for real spaces does not transfer to complex scalars without a check.
  • Polynomials. The space Pd\mathcal{P}_d has basis 1,t,…,td1, t, \dots, t^d, so dim⁡Pd=d+1\dim \mathcal{P}_d = d + 1. The polynomials of degree exactly dd are not a subspace, since the zero polynomial is missing.

For example, the vector (2+i,3−2i)∈C2(2 + i, 3 - 2i) \in \mathbb{C}^2 has two complex coordinates in the standard basis, and over R\mathbb{R} it has the four coordinates (2,1,3,−2)(2, 1, 3, -2) in the basis (1,0),(i,0),(0,1),(0,i)(1, 0), (i, 0), (0, 1), (0, i).

Lessons Real and Complex Vector Spaces

Bases and dimension

The exchange lemma says that if v1,…,vmv_1, \dots, v_m are linearly independent in VV and w1,…,wnw_1, \dots, w_n span VV, then m≤nm \le n. Its key step replaces some wkw_k by a vector v=∑jajwjv = \sum_j a_j w_j with ak≠0a_k \ne 0, which keeps the span. When v1,v2,…v_1, v_2, \dots are inserted one at a time, independence always leaves such a nonzero coefficient on a remaining wkw_k. Applied in both directions to two bases it shows that all bases of VV have the same length, the dimension dim⁡V\dim V, with dim⁡{0}=0\dim \{0\} = 0.

A list of vectors in VV, with dim⁡V=n\dim V = nConclusion
more than nn vectorslinearly dependent
fewer than nn vectorsdoes not span VV
exactly nn independent vectorsa basis
exactly nn spanning vectorsa basis
  • Extracting. Every finite spanning list contains a basis. Remove, one at a time, a vector that is a combination of the others.
  • Extending. Every independent list in a finite-dimensional space extends to a basis. Append, one at a time, a vector outside the current span.
  • Subspaces. If UU is a subspace of VV, then dim⁡U≤dim⁡V\dim U \le \dim V, with equality only when U=VU = V.

For example, u1=(1,0,1)u_1 = (1, 0, 1) and u2=(0,1,1)u_2 = (0, 1, 1) extend to a basis of R3\mathbb{R}^3 by e3e_3, and (2,3,4)=2u1+3u2−e3(2, 3, 4) = 2u_1 + 3u_2 - e_3. The three vectors u1,u2,u1+u2u_1, u_2, u_1 + u_2 span only a plane, so three vectors in R3\mathbb{R}^3 need not form a basis.

Lessons Basis Exchange and Dimension

Linear maps and rank-nullity

A linear map is determined by its values on a basis. If b1,…,bnb_1, \dots, b_n is a basis of VV and w1,…,wn∈Ww_1, \dots, w_n \in W are arbitrary, exactly one linear T:V→WT: V \to W has T(bj)=wjT(b_j) = w_j for every jj. With pointwise addition and scaling, (S+T)(x)=S(x)+T(x)(S + T)(x) = S(x) + T(x) and (aT)(x)=aT(x)(aT)(x) = aT(x), the set Hom⁡(V,W)\operatorname{Hom}(V, W) is a vector space over F\mathbb{F}. For bases B=(b1,…,bn)\mathcal{B} = (b_1, \dots, b_n) of VV and C=(c1,…,cm)\mathcal{C} = (c_1, \dots, c_m) of WW, the matrix-unit maps with Eij(bk)=δjkciE_{ij}(b_k) = \delta_{jk}c_i form a basis, and the coordinates of TT in it are the entries of [T]C←B[T]_{\mathcal{C} \leftarrow \mathcal{B}}. Hence

dim⁡Hom⁡(V,W)=mn.\dim \operatorname{Hom}(V, W) = mn.

A map has rank one exactly when T(x)=φ(x) wT(x) = \varphi(x)\,w for a nonzero w∈Ww \in W and a nonzero linear φ:V→F\varphi: V \to \mathbb{F}. In matrices that is a nonzero column times a nonzero row. The kernel ker⁡T\ker T and the image im⁡T\operatorname{im} T are subspaces, and when VV is finite-dimensional the rank-nullity theorem says

dim⁡V=dim⁡ker⁡T+rank⁡T.\dim V = \dim \ker T + \operatorname{rank} T.
Property of T:V→WT: V \to WTestConsequence in finite dimensions
Injectiveker⁡T={0}\ker T = \{0\}dim⁡V≤dim⁡W\dim V \le \dim W
Surjectiverank⁡T=dim⁡W\operatorname{rank} T = \dim Wdim⁡W≤dim⁡V\dim W \le \dim V
Isomorphisminjective and surjectivedim⁡V=dim⁡W\dim V = \dim W, and T−1T^{-1} is linear

When dim⁡V=dim⁡W\dim V = \dim W is finite, an injective or a surjective linear map is already an isomorphism. With unequal dimensions this fails, as the inclusion of R2\mathbb{R}^2 in R3\mathbb{R}^3 shows. For composable maps, rank⁡(ST)≤min⁡(rank⁡S,rank⁡T)\operatorname{rank}(ST) \le \min(\operatorname{rank} S, \operatorname{rank} T).

For example, T(p)=(p(0),p(1))T(p) = (p(0), p(1)) from P2\mathcal{P}_2 to R2\mathbb{R}^2 sends a+bt+ct2a + bt + ct^2 to (a,a+b+c)(a, a + b + c). Its kernel is spanned by t2−tt^2 - t and it is onto, so 3=1+23 = 1 + 2. All polynomials with measurements (2,5)(2, 5) are 2+3t+s(t2−t)2 + 3t + s(t^2 - t) with s∈Rs \in \mathbb{R}.

Lessons The Vector Space of Linear Maps, Rank-Nullity and Isomorphisms

Sums, direct sums and projections

For subspaces UU and WW of VV, the sum U+W={u+w∣u∈U, w∈W}U + W = \{u + w \mid u \in U,\ w \in W\} and the intersection U∩WU \cap W are subspaces, and for finite-dimensional UU and WW

dim⁡(U+W)=dim⁡U+dim⁡W−dim⁡(U∩W).\dim(U + W) = \dim U + \dim W - \dim(U \cap W).

The union U∪WU \cup W is a subspace only when one of them contains the other. The product U×WU \times W of two spaces, with componentwise operations, has dimension dim⁡U+dim⁡W\dim U + \dim W.

  • Direct sum. The space VV is the direct sum V=U⊕WV = U \oplus W when U+W=VU + W = V and U∩W={0}U \cap W = \{0\}, which holds exactly when every v∈Vv \in V is u+wu + w with u∈Uu \in U and w∈Ww \in W in only one way. Then dim⁡V=dim⁡U+dim⁡W\dim V = \dim U + \dim W.
  • Complements. In a finite-dimensional space every subspace UU has a complement WW with V=U⊕WV = U \oplus W, found by extending a basis of UU. It is unique only when U={0}U = \{0\} or U=VU = V.
  • Projection along a complement. If V=U⊕WV = U \oplus W, the map P(u+w)=uP(u + w) = u is linear with P2=PP^2 = P, im⁡P=U\operatorname{im} P = U and ker⁡P=W\ker P = W, and I−PI - P projects onto WW along UU.
  • Idempotents. Conversely, every linear PP with P2=PP^2 = P gives V=im⁡P⊕ker⁡PV = \operatorname{im} P \oplus \ker P and is the projection onto im⁡P\operatorname{im} P along ker⁡P\ker P. The splitting V=im⁡P⊕ker⁡PV = \operatorname{im} P \oplus \ker P alone does not force P2=PP^2 = P, as P=2IP = 2I shows.

A projection needs no inner product. With the dot product on Rn\mathbb{R}^n, an idempotent real matrix is the orthogonal projection onto its image exactly when it is symmetric, and otherwise it is an oblique projection. For example, the projection onto U=span⁡(1,0)U = \operatorname{span}(1, 0) along W=span⁡(1,1)W = \operatorname{span}(1, 1) is P(x,y)=(x−y,0)P(x, y) = (x - y, 0), with matrix

P=[1−100],P = \begin{bmatrix} 1 & -1 \\ 0 & 0 \end{bmatrix},

and it splits (3,2)(3, 2) as (1,0)+(2,2)(1, 0) + (2, 2). The matrix is not symmetric, so the projection is oblique.

Lessons Sums and Intersections of Subspaces, Direct Sums and Projections

Quotient Spaces and Duality

A quotient space treats vectors that differ by an element of a fixed subspace as one object, and the dual space collects all linear measurements of vectors. Neither construction needs an inner product. Together they describe what a map cannot detect, explain why measurements change coordinates in the opposite way to vectors, and show why row rank equals column rank.

Quotient spaces

Let WW be a subspace of VV. Vectors vv and uu of VV are equivalent modulo WW when v−u∈Wv - u \in W. The coset of vv is

v+W={v+w∣w∈W},v + W = \{v + w \mid w \in W\},

and v+W=u+Wv + W = u + W exactly when v−u∈Wv - u \in W. Two cosets are equal or disjoint, and a coset is a subspace only when it is WW itself. The quotient space V/WV/W is the set of cosets with the operations

(v+W)+(u+W)=(v+u)+W,a(v+W)=av+W,\begin{aligned} (v + W) + (u + W) &= (v + u) + W, \\[0.4em] a(v + W) &= av + W, \end{aligned}

which are well defined because changing representatives changes the results by vectors of WW. Its zero vector is the coset WW.

  • Well-defined rules. A rule v+W↦f(v)v + W \mapsto f(v) defines a function on V/WV/W exactly when v−u∈Wv - u \in W implies f(v)=f(u)f(v) = f(u). For a linear ff this means W⊆ker⁡fW \subseteq \ker f.
  • Quotient basis. If w1,…,wrw_1, \dots, w_r is a basis of WW extended to a basis w1,…,wr,v1,…,vsw_1, \dots, w_r, v_1, \dots, v_s of VV, then the cosets v1+W,…,vs+Wv_1 + W, \dots, v_s + W form a basis of V/WV/W.

When VV is finite-dimensional, counting that basis gives

dim⁡(V/W)=dim⁡V−dim⁡W.\dim(V/W) = \dim V - \dim W.

For example, in R3\mathbb{R}^3 with W=span⁡(e3)W = \operatorname{span}(e_3), each coset is a vertical line (a,b,c)+W(a, b, c) + W. The function a+ba + b is constant on each line and descends to the quotient, while cc is not and does not.

Lessons Quotient Spaces

The first isomorphism theorem

Let T:V→WT: V \to W be linear over F\mathbb{F}. The fiber over yy is T−1(y)={v∣T(v)=y}T^{-1}(y) = \{v \mid T(v) = y\}, which is empty when y∉im⁡Ty \notin \operatorname{im} T and equals the coset v0+ker⁡Tv_0 + \ker T when T(v0)=yT(v_0) = y. The notation does not assume an inverse.

  • Quotient map. For a subspace UU of VV, the map q(v)=v+Uq(v) = v + U from VV to V/UV/U is linear and onto, with kernel UU.
  • Induced map. The rule T~(v+U)=T(v)\widetilde{T}(v + U) = T(v) defines a map on V/UV/U exactly when U⊆ker⁡TU \subseteq \ker T. It is then linear, satisfies T=T~∘qT = \widetilde{T} \circ q, and is the only map with that property.

The first isomorphism theorem says that the induced map

T~:V/ker⁡T→im⁡T,v+ker⁡T↦T(v),\begin{gathered} \widetilde{T}: V/\ker T \to \operatorname{im} T, \\[0.4em] v + \ker T \mapsto T(v), \end{gathered}

is a linear bijection, so V/ker⁡T≅im⁡TV/\ker T \cong \operatorname{im} T. The target is the image, not the codomain, and no finite-dimensional hypothesis is needed. In finite dimensions, counting dimensions recovers rank-nullity.

For example, T(x,y,z)=(x+y,z)T(x, y, z) = (x + y, z) on R3\mathbb{R}^3 has ker⁡T=span⁡(1,−1,0)\ker T = \operatorname{span}(1, -1, 0), and the coset of (x,y,z)(x, y, z) corresponds to (x+y,z)(x + y, z). The cosets of (1,0,0)(1, 0, 0) and (0,0,1)(0, 0, 1) form a quotient basis. Evaluation ev⁡s(p)=p(s)\operatorname{ev}_s(p) = p(s) maps Pd\mathcal{P}_d onto F\mathbb{F}, so Pd/ker⁡ev⁡s≅F\mathcal{P}_d / \ker \operatorname{ev}_s \cong \mathbb{F}, and for d≥1d \ge 1 the kernel consists of the multiples (t−s)r(t)(t - s)r(t) with r∈Pd−1r \in \mathcal{P}_{d-1}.

Lessons The First Isomorphism Theorem

Dual spaces and dual bases

A linear functional on VV is a linear map ℓ:V→F\ell: V \to \mathbb{F}. Over C\mathbb{C} scalars come out with no conjugate. The dual space is

V∨=Hom⁡(V,F),V^{\vee} = \operatorname{Hom}(V, \mathbb{F}),

with pointwise operations. Many books write V∗V^{*} for the dual space. This course keeps V∨V^{\vee} so that the dual is not confused with the adjoint T∗T^{*}. On Fn\mathbb{F}^n every functional is ℓ(x)=rx\ell(x) = rx for exactly one row rr, whose entries are ℓ(e1),…,ℓ(en)\ell(e_1), \dots, \ell(e_n), with no conjugation.

  • Dual basis. For a basis b1,…,bnb_1, \dots, b_n of VV, the coordinate functionals β1,…,βn\beta_1, \dots, \beta_n with βi(bj)=δij\beta_i(b_j) = \delta_{ij} form a basis of V∨V^{\vee}, and every functional expands as ℓ=∑iℓ(bi) βi\ell = \sum_i \ell(b_i)\,\beta_i. So dim⁡V∨=dim⁡V\dim V^{\vee} = \dim V. On Fn\mathbb{F}^n, if the bjb_j are the columns of BB, the βi\beta_i are the rows of B−1B^{-1}.
  • Evaluation functionals. On Pd\mathcal{P}_d with basis 1,t,…,td1, t, \dots, t^d and dual basis β1,…,βd+1\beta_1, \dots, \beta_{d+1}, evaluation at a fixed ss is ev⁡s=∑k=1d+1sk−1βk\operatorname{ev}_s = \sum_{k=1}^{d+1} s^{k-1}\beta_k.
  • Not functionals. On a complex space, a conjugate-linear map ff, with f(av)=a‾f(v)f(av) = \overline{a}f(v), is not in V∨V^{\vee} unless it is zero. On a nonzero space the norm x↦∥x∥x \mapsto \lVert x \rVert is not linear either.

The equality dim⁡V∨=dim⁡V\dim V^{\vee} = \dim V does not single out an isomorphism between the two spaces. Matching bib_i with βi\beta_i depends on the basis. For example, the basis (1,1),(1,−1)(1, 1), (1, -1) of R2\mathbb{R}^2 has dual basis β1(x,y)=(x+y)/2\beta_1(x, y) = (x + y)/2 and β2(x,y)=(x−y)/2\beta_2(x, y) = (x - y)/2, and the functional φ(x,y)=2x−y\varphi(x, y) = 2x - y equals β1+3β2\beta_1 + 3\beta_2. The vector (3,1)(3, 1) has coordinates (2,1)(2, 1) in this basis, and both descriptions give φ(3,1)=5\varphi(3, 1) = 5.

Lessons Dual Spaces and Dual Bases

Dual maps and covector coordinates

A covector is another name for a linear functional. The dual map of a linear T:V→WT: V \to W pulls measurements back,

T∨:W∨→V∨,T∨(ℓ)=ℓ∘T.T^{\vee}: W^{\vee} \to V^{\vee}, \qquad T^{\vee}(\ell) = \ell \circ T.

It is linear over R\mathbb{R} and over C\mathbb{C}, it uses no inner product, it reverses compositions as (ST)∨=T∨S∨(ST)^{\vee} = T^{\vee}S^{\vee}, and it satisfies (T−1)∨=(T∨)−1(T^{-1})^{\vee} = (T^{\vee})^{-1} for an isomorphism TT. If A=[T]C←BA = [T]_{\mathcal{C} \leftarrow \mathcal{B}} and B∨\mathcal{B}^{\vee}, C∨\mathcal{C}^{\vee} are the dual bases, then

[T∨]B∨←C∨=AT,[T^{\vee}]_{\mathcal{B}^{\vee} \leftarrow \mathcal{C}^{\vee}} = A^{\mathsf{T}},

the plain transpose even over C\mathbb{C}. In rows, a functional with evaluation row rr on WW pulls back to the row rArA on VV.

Object, with P=PC←BP = P_{\mathcal{C} \leftarrow \mathcal{B}}Coordinates in B\mathcal{B}Coordinates in C\mathcal{C}
Vector vvcolumn [v]B[v]_{\mathcal{B}}[v]C=P[v]B[v]_{\mathcal{C}} = P[v]_{\mathcal{B}}
Covector ℓ\ell as an evaluation rowrow rBr_{\mathcal{B}} with entries ℓ(bj)\ell(b_j)rC=rBP−1r_{\mathcal{C}} = r_{\mathcal{B}}P^{-1}
Covector ℓ\ell in the dual basiscolumn [ℓ]B∨=rBT[\ell]_{\mathcal{B}^{\vee}} = r_{\mathcal{B}}^{\mathsf{T}}[ℓ]C∨=(P−1)T[ℓ]B∨[\ell]_{\mathcal{C}^{\vee}} = (P^{-1})^{\mathsf{T}}[\ell]_{\mathcal{B}^{\vee}}

The reading ℓ(v)=rC[v]C=rB[v]B\ell(v) = r_{\mathcal{C}}[v]_{\mathcal{C}} = r_{\mathcal{B}}[v]_{\mathcal{B}} does not depend on the basis, which is why covector coordinates use the inverse. Multiplying both the coefficient column of ℓ\ell and the coordinate column of vv by PP would keep every reading only when PTP=IP^{\mathsf{T}}P = I.

For example, with AA having rows (1,2)(1, 2) and (0,1)(0, 1) and the output sensor r=(3,−1)r = (3, -1), the pullback is rA=(3,5)rA = (3, 5). At v=(1,2)v = (1, 2) both routes give 1313, since Av=(5,2)Av = (5, 2). For φ(x,y)=2x−y\varphi(x, y) = 2x - y, standard B\mathcal{B} and C=((1,1),(1,−1))\mathcal{C} = ((1, 1), (1, -1)), the row rB=(2,−1)r_{\mathcal{B}} = (2, -1) becomes rC=(1,3)r_{\mathcal{C}} = (1, 3), the coefficients of φ\varphi in the dual basis above.

Lessons Dual Maps and Covectors

Annihilators and the four subspaces

The annihilator of a subspace UU of VV is

U0={ℓ∈V∨∣ℓ(u)=0 for all u∈U},U^{0} = \{\ell \in V^{\vee} \mid \ell(u) = 0 \text{ for all } u \in U\},

a subspace of V∨V^{\vee}, not of VV. For finite-dimensional VV, dim⁡U0=dim⁡V−dim⁡U\dim U^{0} = \dim V - \dim U. For a linear T:V→WT: V \to W,

ker⁡T∨=(im⁡T)0,im⁡T∨=(ker⁡T)0,\begin{aligned} \ker T^{\vee} &= (\operatorname{im} T)^{0}, \\[0.4em] \operatorname{im} T^{\vee} &= (\ker T)^{0}, \end{aligned}

where the first identity holds for all vector spaces and the second follows by counting dimensions when VV and WW are finite-dimensional. Hence rank⁡T∨=rank⁡T\operatorname{rank} T^{\vee} = \operatorname{rank} T, and T∨T^{\vee} is injective exactly when TT is surjective. For a matrix A∈Fm×nA \in \mathbb{F}^{m \times n} of rank rr, the dual map sends an output row yy to yAyA, which gives the row space as (Null⁡A)0(\operatorname{Null} A)^{0} and the left nullspace as (Col⁡A)0(\operatorname{Col} A)^{0}.

SubspaceLives inDimension
Column space Col⁡A\operatorname{Col} AFm\mathbb{F}^mrr
Nullspace Null⁡A\operatorname{Null} AFn\mathbb{F}^nn−rn - r
Row space Row⁡A\operatorname{Row} Arows of length nnrr
Left nullspace {y∣yA=0}\{y \mid yA = 0\}rows of length mmm−rm - r

An annihilator needs no extra structure, while an orthogonal complement U⊥U^{\perp} lies in VV and depends on an inner product. In a finite-dimensional inner product space, the functional ⟨v,⋅⟩\langle v, \cdot \rangle lies in U0U^{0} exactly when v∈U⊥v \in U^{\perp}, and the correspondence is conjugate-linear over C\mathbb{C}. For example, AA with rows (1,1,0)(1, 1, 0) and (0,0,1)(0, 0, 1) has kernel spanned by (1,−1,0)(1, -1, 0). The rows that vanish on it are (a,a,b)(a, a, b), exactly the row space, and since AA maps onto R2\mathbb{R}^2 the left nullspace is zero.

Lessons Annihilators and Fundamental Subspaces

Determinants and Bilinear Forms

The determinant is the only function of the columns of a square matrix that is multilinear, alternating and equal to 11 at the identity, and its product rule follows from that uniqueness. Trace, determinant and characteristic polynomial stay the same under similarity, the rule by which the matrix of an operator changes, while the matrix of a bilinear form changes by congruence. The same square array can therefore follow two different coordinate laws.

The determinant from alternating forms

A function D(v1,…,vk)D(v_1, \dots, v_k) of kk vectors is multilinear if it is linear in each slot with the others fixed, and alternating if it is 00 whenever two arguments are equal. An alternating multilinear form changes sign when two arguments are swapped, is unchanged when a multiple of one argument is added to another, and vanishes on linearly dependent arguments. The sign of a permutation σ\sigma of {1,…,n}\{1, \dots, n\} is sgn⁡(σ)=(−1)N(σ)\operatorname{sgn}(\sigma) = (-1)^{N(\sigma)}, where N(σ)N(\sigma) counts the pairs i<ji < j with σ(i)>σ(j)\sigma(i) > \sigma(j), and every exchange of two entries flips it. For A=[aij]∈Fn×nA = [a_{ij}] \in \mathbb{F}^{n \times n},

det⁡A=∑σ∈Snsgn⁡(σ)∏j=1naσ(j)j.\det A = \sum_{\sigma \in S_n} \operatorname{sgn}(\sigma) \prod_{j=1}^{n} a_{\sigma(j) j}.

The uniqueness theorem says that every function DD of the columns that is multilinear and alternating satisfies

D(A)=D(I) det⁡A,D(A) = D(I)\,\det A,

so det⁡\det is the only one with D(I)=1D(I) = 1. Applying it to D(B)=det⁡(AB)D(B) = \det(AB) gives det⁡(AB)=det⁡Adet⁡B\det(AB) = \det A \det B for all square AA and BB of the same size, singular or not. Also det⁡AT=det⁡A\det A^{\mathsf{T}} = \det A, so every column rule holds for rows, and AA is invertible exactly when det⁡A≠0\det A \ne 0, with det⁡(A−1)=1/det⁡A\det(A^{-1}) = 1/\det A.

Operation on the columns or the rowsEffect on det⁡A\det A
swap two of themchanges the sign
multiply one by ccmultiplies it by cc
add a multiple of one to anothernone
make them linearly dependentgives 00

Lessons Determinants and Multilinear Forms

Trace, similarity and the characteristic polynomial

The trace tr⁡A\operatorname{tr} A is the sum of the diagonal entries. For A∈Fm×nA \in \mathbb{F}^{m \times n} and B∈Fn×mB \in \mathbb{F}^{n \times m},

tr⁡(AB)=tr⁡(BA),\operatorname{tr}(AB) = \operatorname{tr}(BA),

even when the two products have different sizes. With three factors only the cyclic rotations tr⁡(ABC)=tr⁡(BCA)=tr⁡(CAB)\operatorname{tr}(ABC) = \operatorname{tr}(BCA) = \operatorname{tr}(CAB) are allowed. Matrices A,B∈Fn×nA, B \in \mathbb{F}^{n \times n} are similar when B=S−1ASB = S^{-1}AS for an invertible S∈Fn×nS \in \mathbb{F}^{n \times n}. Similar matrices are exactly the matrices of one operator in two bases. If S=PB←CS = P_{\mathcal{B} \leftarrow \mathcal{C}}, then [T]C←C=S−1[T]B←BS[T]_{\mathcal{C} \leftarrow \mathcal{C}} = S^{-1}[T]_{\mathcal{B} \leftarrow \mathcal{B}}S.

The characteristic polynomial of A∈Fn×nA \in \mathbb{F}^{n \times n} is monic of degree nn. For n≥2n \ge 2 it expands as

χA(t)=det⁡(tI−A)=tn−(tr⁡A) tn−1+⋯+(−1)ndet⁡A,\begin{aligned} \chi_A(t) &= \det(tI - A) \\[0.4em] &= t^n - (\operatorname{tr} A)\,t^{n-1} \\[0.4em] &\quad + \cdots + (-1)^n \det A, \end{aligned}

and for n=2n = 2 it is t2−(tr⁡A) t+det⁡At^2 - (\operatorname{tr} A)\,t + \det A. A scalar λ∈F\lambda \in \mathbb{F} is an eigenvalue of AA exactly when χA(λ)=0\chi_A(\lambda) = 0. Trace, determinant, characteristic polynomial, rank and the dimension of every eigenspace are similarity invariants, so they are defined for an operator on a finite-dimensional space. Equal invariants do not imply similarity. The identity and the shear with rows (1,1)(1, 1) and (0,1)(0, 1) share trace 22, determinant 11 and characteristic polynomial (t−1)2(t - 1)^2, but the shear is not similar to II, since S−1IS=IS^{-1}IS = I for every invertible SS.

Over C\mathbb{C} a determinant is not a signed volume. Multiplication by A∈Cn×nA \in \mathbb{C}^{n \times n}, viewed as a real map on R2n\mathbb{R}^{2n}, scales real volume by ∣det⁡A∣2\lvert \det A \rvert^2. For example, D=diag⁡(2,3)D = \operatorname{diag}(2, 3) and SS with rows (1,1)(1, 1) and (0,1)(0, 1) give

S−1DS=[2−103],S^{-1}DS = \begin{bmatrix} 2 & -1 \\ 0 & 3 \end{bmatrix},

with the same determinant 66, trace 55 and characteristic polynomial t2−5t+6t^2 - 5t + 6 as DD.

Lessons Trace, Determinant and Similarity

Bilinear forms and congruence

A bilinear form on a real space VV is a function b:V×V→Rb: V \times V \to \mathbb{R} that is linear in each argument. In an ordered basis B=(v1,…,vn)\mathcal{B} = (v_1, \dots, v_n) its matrix MBM_{\mathcal{B}} has entries mij=b(vi,vj)m_{ij} = b(v_i, v_j), and with ξ=[x]B\xi = [x]_{\mathcal{B}} and η=[y]B\eta = [y]_{\mathcal{B}}

b(x,y)=ξTMB η.b(x, y) = \xi^{\mathsf{T}} M_{\mathcal{B}}\, \eta.

The form is symmetric exactly when MBM_{\mathcal{B}} is, and nondegenerate exactly when MBM_{\mathcal{B}} is invertible. If the columns of S=PB←CS = P_{\mathcal{B} \leftarrow \mathcal{C}} hold the new basis vectors in old coordinates, then

MC=STMB S,M_{\mathcal{C}} = S^{\mathsf{T}} M_{\mathcal{B}}\, S,

which is the congruence rule. With P=PC←B=S−1P = P_{\mathcal{C} \leftarrow \mathcal{B}} = S^{-1} it reads MC=(P−1)TMB P−1M_{\mathcal{C}} = (P^{-1})^{\mathsf{T}} M_{\mathcal{B}}\, P^{-1}.

Object, new basis given by the columns of SSNew matrixUnchanged
Operator T:V→VT: V \to VS−1ASS^{-1}AS, similarityeigenvalues, trace, determinant, χA\chi_A
Bilinear form on a real spaceSTMSS^{\mathsf{T}}MS, congruencesymmetry, rank, the sign of det⁡M\det M, the inertia of a symmetric form
Hermitian form on a complex spaceS∗HSS^{*}HS, conjugate congruencethe Hermitian property, rank, inertia
  • Quadratic form. The function q(v)=b(v,v)q(v) = b(v, v) satisfies q(cv)=c2q(v)q(cv) = c^2 q(v).

  • Polarization. For a symmetric bb, the polarization identity

    b(x,y)=12(q(x+y)−q(x)−q(y))b(x, y) = \frac{1}{2}\bigl(q(x + y) - q(x) - q(y)\bigr)

    recovers bb from qq, so symmetric forms and quadratic forms correspond one to one.

  • Symmetric and skew parts. Every real square MM is H+KH + K with H=12(M+MT)H = \tfrac{1}{2}(M + M^{\mathsf{T}}) symmetric and K=12(M−MT)K = \tfrac{1}{2}(M - M^{\mathsf{T}}) skew-symmetric. Since xTKx=0x^{\mathsf{T}}Kx = 0, the quadratic form sees only HH.

For example, MB=diag⁡(1,4)M_{\mathcal{B}} = \operatorname{diag}(1, 4) and the new basis (2,0),(0,1)(2, 0), (0, 1) give S=diag⁡(2,1)S = \operatorname{diag}(2, 1) and MC=diag⁡(4,4)M_{\mathcal{C}} = \operatorname{diag}(4, 4). The eigenvalues changed, which similarity never does. The vector (2,1)(2, 1) has C\mathcal{C} coordinates (1,1)(1, 1), and both matrices give q=8q = 8. The matrix with rows (0,1)(0, 1) and (−1,0)(-1, 0) is nonzero but has q=0q = 0 everywhere.

Lessons Bilinear Forms and Congruence

Inner Product Spaces

An inner product measures lengths and angles. Over the complex numbers it needs a conjugate, which keeps the squared length of every nonzero vector positive. This unit proves the Cauchy-Schwarz inequality, builds orthonormal bases and expansions, matches every linear functional with a vector, defines the adjoint, and finds best approximations, including weighted least squares, where changing the inner product can move the best fit.

Inner products, norms and Gram matrices

An inner product on a vector space VV over F\mathbb{F} assigns a scalar ⟨x,y⟩\langle x, y \rangle to each pair so that, for all vectors and scalars,

⟨x,ay+bz⟩=a⟨x,y⟩+b⟨x,z⟩,⟨y,x⟩=⟨x,y⟩‾,⟨x,x⟩>0  for x≠0.\begin{aligned} &\langle x, ay + bz \rangle = a\langle x, y \rangle + b\langle x, z \rangle, \\[0.4em] &\langle y, x \rangle = \overline{\langle x, y \rangle}, \\[0.4em] &\langle x, x \rangle > 0 \ \text{ for } x \ne 0. \end{aligned}

These are linearity in the second argument, conjugate symmetry and positive definiteness. They force ⟨ax,y⟩=a‾⟨x,y⟩\langle ax, y \rangle = \overline{a}\langle x, y \rangle, so the first argument is conjugate-linear. The standard inner product on Fn\mathbb{F}^n is ⟨x,y⟩=x∗y=∑jxj‾ yj\langle x, y \rangle = x^{*}y = \sum_j \overline{x_j}\,y_j.

  • Norm. The norm ∥x∥=⟨x,x⟩\lVert x \rVert = \sqrt{\langle x, x \rangle} has ∥ax∥=∣a∣ ∥x∥\lVert ax \rVert = \lvert a \rvert\,\lVert x \rVert, and ∥x+y∥2=∥x∥2+∥y∥2+2Re⁡⟨x,y⟩\lVert x + y \rVert^2 = \lVert x \rVert^2 + \lVert y \rVert^2 + 2\operatorname{Re}\langle x, y \rangle. Orthogonal vectors, with ⟨x,y⟩=0\langle x, y \rangle = 0, satisfy the Pythagorean identity.
  • Sesquilinear forms. A map that is linear in the second argument and conjugate-linear in the first is sesquilinear, and Hermitian if s(y,x)=s(x,y)‾s(y, x) = \overline{s(x, y)}. A complex inner product is exactly a positive definite Hermitian sesquilinear form. On Cn\mathbb{C}^n these forms are s(x,y)=x∗Hys(x, y) = x^{*}Hy, Hermitian exactly when H∗=HH^{*} = H.
  • Gram matrix. For a basis b1,…,bnb_1, \dots, b_n the Gram matrix has entries Gij=⟨bi,bj⟩G_{ij} = \langle b_i, b_j \rangle, and with coordinate columns ξ\xi and η\eta we have ⟨x,y⟩=ξ∗Gη\langle x, y \rangle = \xi^{*}G\eta. It is Hermitian and positive definite, it equals II exactly for an orthonormal basis, and a change of basis replaces it by S∗GSS^{*}GS, where the columns of SS hold the new basis vectors in old coordinates.

For example, x=(1,i)x = (1, i) has xTx=1+i2=0x^{\mathsf{T}}x = 1 + i^2 = 0 although x≠0x \ne 0, while x∗x=2x^{*}x = 2. Scalars leave the first slot conjugated, so ⟨ix,y⟩=−i⟨x,y⟩\langle ix, y \rangle = -i\langle x, y \rangle and ⟨x,iy⟩=i⟨x,y⟩\langle x, iy \rangle = i\langle x, y \rangle.

Lessons Complex Inner Products

Cauchy-Schwarz and orthonormal expansions

For all x,yx, y in an inner product space over F\mathbb{F}, in any dimension,

∣⟨x,y⟩∣≤∥x∥ ∥y∥,\lvert \langle x, y \rangle \rvert \le \lVert x \rVert\,\lVert y \rVert,

with equality exactly when xx and yy are linearly dependent. For x≠0x \ne 0 the proof expands ∥y−cx∥2≥0\lVert y - cx \rVert^2 \ge 0 with c=⟨x,y⟩/∥x∥2c = \langle x, y \rangle / \lVert x \rVert^2, which makes y−cxy - cx orthogonal to xx. It gives the triangle inequality ∥x+y∥≤∥x∥+∥y∥\lVert x + y \rVert \le \lVert x \rVert + \lVert y \rVert, with equality exactly when one vector is a nonnegative real multiple of the other.

The complex Gram-Schmidt process turns linearly independent a1,…,ana_1, \dots, a_n into an orthonormal list with the same spans span⁡(q1,…,qk)=span⁡(a1,…,ak)\operatorname{span}(q_1, \dots, q_k) = \operatorname{span}(a_1, \dots, a_k) for every kk. Set q1=a1/∥a1∥q_1 = a_1/\lVert a_1 \rVert and, for k≥2k \ge 2,

rk=ak−∑j=1k−1⟨qj,ak⟩ qj,qk=rk∥rk∥.\begin{aligned} r_k &= a_k - \sum_{j=1}^{k-1} \langle q_j, a_k \rangle\, q_j, \\[0.4em] q_k &= \frac{r_k}{\lVert r_k \rVert}. \end{aligned}

The basis vector sits in the first slot of each coefficient. With dependent inputs, the first aka_k in the span of the earlier ones gives rk=0r_k = 0, or is itself zero when k=1k = 1, and that zero signals the dependence. For an orthonormal basis q1,…,qnq_1, \dots, q_n of VV and every v∈Vv \in V,

v=∑j=1n⟨qj,v⟩ qj,∥v∥2=∑j=1n∣⟨qj,v⟩∣2.\begin{gathered} v = \sum_{j=1}^{n} \langle q_j, v \rangle\, q_j, \\[0.4em] \lVert v \rVert^2 = \sum_{j=1}^{n} \lvert \langle q_j, v \rangle \rvert^2. \end{gathered}

The second line is Parseval’s identity. For a finite orthonormal list q1,…,qkq_1, \dots, q_k that need not be a basis, Bessel’s inequality says

∑j=1k∣⟨qj,v⟩∣2≤∥v∥2,\sum_{j=1}^{k} \lvert \langle q_j, v \rangle \rvert^2 \le \lVert v \rVert^2,

with equality exactly when vv lies in the span of the list.

For example, x=(1,i)x = (1, i) and y=(1,1)y = (1, 1) have ∥x∥2=∥y∥2=2\lVert x \rVert^2 = \lVert y \rVert^2 = 2 and ⟨x,y⟩=1−i\langle x, y \rangle = 1 - i, of modulus 2<2\sqrt{2} < 2. The remainder y−1−i2x=(1+i2,1−i2)y - \tfrac{1 - i}{2}x = \bigl(\tfrac{1 + i}{2}, \tfrac{1 - i}{2}\bigr) is orthogonal to xx.

Lessons Cauchy-Schwarz and Orthonormal Expansions

Riesz representation and adjoints

The Riesz representation theorem says that if VV is a finite-dimensional inner product space over F\mathbb{F}, then for every ℓ∈V∨\ell \in V^{\vee} there is exactly one z∈Vz \in V with

ℓ(x)=⟨z,x⟩for all x∈V.\ell(x) = \langle z, x \rangle \quad \text{for all } x \in V.

In an orthonormal basis e1,…,ene_1, \dots, e_n the vector is z=∑jℓ(ej)‾ ejz = \sum_j \overline{\ell(e_j)}\,e_j. On Fn\mathbb{F}^n with the inner product x∗Gyx^{*}Gy for a Hermitian positive definite GG, the functional ℓ(x)=rx\ell(x) = rx is represented by z=G−1r∗z = G^{-1}r^{*}, which is r∗r^{*} for the standard inner product. The vector depends on the inner product, and the identification z↦⟨z,⋅⟩z \mapsto \langle z, \cdot \rangle is conjugate-linear over C\mathbb{C}.

For finite-dimensional inner product spaces VV and WW over F\mathbb{F} and a linear T:V→WT: V \to W, the adjoint is the unique linear map T∗:W→VT^{*}: W \to V with

⟨Tv,w⟩W=⟨v,T∗w⟩Vfor all v∈V and w∈W.\begin{gathered} \langle Tv, w \rangle_W = \langle v, T^{*}w \rangle_V \\[0.4em] \text{for all } v \in V \text{ and } w \in W. \end{gathered}

It obeys (ST)∗=T∗S∗(ST)^{*} = T^{*}S^{*}, (T∗)∗=T(T^{*})^{*} = T, (T1+T2)∗=T1∗+T2∗(T_1 + T_2)^{*} = T_1^{*} + T_2^{*} and (aT)∗=a‾ T∗(aT)^{*} = \overline{a}\,T^{*}. In orthonormal bases its matrix is A∗A^{*}. In bases with Gram matrices GVG_V and GWG_W, where A=[T]C←BA = [T]_{\mathcal{C} \leftarrow \mathcal{B}},

[T∗]B←C=GV−1A∗GW.[T^{*}]_{\mathcal{B} \leftarrow \mathcal{C}} = G_V^{-1}A^{*}G_W.

The orthogonal complement U⊥={v∣⟨u,v⟩=0 for all u∈U}U^{\perp} = \{v \mid \langle u, v \rangle = 0 \text{ for all } u \in U\} has dimension dim⁡V−dim⁡U\dim V - \dim U, and (U⊥)⊥=U(U^{\perp})^{\perp} = U in finite dimensions. For linear T:V→WT: V \to W,

(im⁡T)⊥=ker⁡T∗,(ker⁡T)⊥=im⁡T∗,\begin{aligned} (\operatorname{im} T)^{\perp} &= \ker T^{*}, \\[0.4em] (\ker T)^{\perp} &= \operatorname{im} T^{*}, \end{aligned}

so W=im⁡T⊕ker⁡T∗W = \operatorname{im} T \oplus \ker T^{*}, V=ker⁡T⊕im⁡T∗V = \ker T \oplus \operatorname{im} T^{*} and rank⁡T∗=rank⁡T\operatorname{rank} T^{*} = \operatorname{rank} T.

Dual map T∨T^{\vee}Adjoint T∗T^{*}
Actsfrom W∨W^{\vee} to V∨V^{\vee}from WW to VV
Needsnothing beyond TTinner products on VV and WW
MatrixATA^{\mathsf{T}} in dual basesA∗A^{*} in orthonormal bases, GV−1A∗GWG_V^{-1}A^{*}G_W in general

For example, AA with rows (1,i)(1, i) and (0,2)(0, 2) has A∗A^{*} with rows (1,0)(1, 0) and (−i,2)(-i, 2). With x=(1,i)x = (1, i) and y=(1,1)y = (1, 1), Ax=(0,2i)Ax = (0, 2i), and both ⟨Ax,y⟩\langle Ax, y \rangle and ⟨x,A∗y⟩\langle x, A^{*}y \rangle equal −2i-2i. With the domain Gram matrix GV=diag⁡(2,1)G_V = \operatorname{diag}(2, 1), the Euclidean target and A=(1  1)A = (1 \ \ 1), the adjoint has coordinate column (1/2,1)(1/2, 1), not A∗=(1,1)A^{*} = (1, 1).

Lessons Riesz Representation in Finite Dimensions, Adjoints and Gram Matrices

Orthogonal projection and best approximation

Let UU be a subspace of a finite-dimensional inner product space VV over F\mathbb{F}. Then V=U⊕U⊥V = U \oplus U^{\perp}, and the orthogonal projection PUvP_U v is the unique vector with PUv∈UP_U v \in U and v−PUv∈U⊥v - P_U v \in U^{\perp}. For an orthonormal basis q1,…,qrq_1, \dots, q_r of UU,

PUv=∑j=1r⟨qj,v⟩ qj.P_U v = \sum_{j=1}^{r} \langle q_j, v \rangle\, q_j.

For a basis u1,…,uru_1, \dots, u_r that is not orthonormal, PUv=∑jcjujP_U v = \sum_j c_j u_j where Gc=hGc = h, with Gij=⟨ui,uj⟩G_{ij} = \langle u_i, u_j \rangle and hi=⟨ui,v⟩h_i = \langle u_i, v \rangle.

The projection theorem says that PUvP_U v is the only best approximation to vv from UU, and that u0∈Uu_0 \in U is that best approximation exactly when v−u0∈U⊥v - u_0 \in U^{\perp}. For every u∈Uu \in U,

∥v−u∥2=∥v−PUv∥2+∥PUv−u∥2,\begin{aligned} \lVert v - u \rVert^2 &= \lVert v - P_U v \rVert^2 \\[0.4em] &\quad + \lVert P_U v - u \rVert^2, \end{aligned}

and the squared distance from vv to UU is ∥v∥2−∥PUv∥2\lVert v \rVert^2 - \lVert P_U v \rVert^2. A linear P:V→VP: V \to V is an orthogonal projector exactly when P2=PP^2 = P and P∗=PP^{*} = P. If the orthonormal columns of QQ span UU in orthonormal coordinates, its matrix is QQ∗QQ^{*}.

On Pd\mathcal{P}_d, distinct real nodes t1,…,tNt_1, \dots, t_N with N≥d+1N \ge d + 1 give the inner product ⟨p,q⟩=∑kp(tk)‾ q(tk)\langle p, q \rangle = \sum_k \overline{p(t_k)}\,q(t_k), and the best polynomial approximation from a subspace depends on that choice. For example, with nodes −1,0,1-1, 0, 1 the projection of 1+t1 + t onto the constants is 11, and the residual tt has samples (−1,0,1)(-1, 0, 1), orthogonal to the constant samples (1,1,1)(1, 1, 1).

Lessons Orthogonal Projection and Best Approximation

Weighted least squares

A weight matrix W∈Fm×mW \in \mathbb{F}^{m \times m} is Hermitian and positive definite, and a diagonal weight needs positive diagonal entries. It defines the inner product ⟨x,y⟩W=x∗Wy\langle x, y \rangle_W = x^{*}Wy with norm ∥x∥W\lVert x \rVert_W, and every inner product on Fm\mathbb{F}^m has this form. For A∈Fm×nA \in \mathbb{F}^{m \times n} and b∈Fmb \in \mathbb{F}^m, a vector x^\hat{x} minimizes ∥b−Ax∥W2=(b−Ax)∗W(b−Ax)\lVert b - Ax \rVert_W^2 = (b - Ax)^{*}W(b - Ax) exactly when it solves the weighted normal equations

A∗WA x^=A∗Wb,A^{*}WA\,\hat{x} = A^{*}Wb,

which say that the residual is WW-orthogonal to every column of AA. They always have a solution.

  • Fitted vector. The fitted vector Ax^A\hat{x} is the same for every minimizer. It is the projection of bb onto the column space in the inner product of WW.
  • Coefficients. The minimizers form the set x^0+ker⁡A\hat{x}_0 + \ker A. They are unique exactly when the columns of AA are independent, and then x^=(A∗WA)−1A∗Wb\hat{x} = (A^{*}WA)^{-1}A^{*}Wb.
  • Weights. Multiplying WW by a constant c>0c > 0 does not change the fit. A singular positive semidefinite WW or an indefinite WW does not define a distance.

With positive weights wkw_k and distinct real nodes, ⟨p,q⟩=∑k=1Nwk p(tk)‾ q(tk)\langle p, q \rangle = \sum_{k=1}^{N} w_k\,\overline{p(t_k)}\,q(t_k) is an inner product on Pd\mathcal{P}_d exactly when N≥d+1N \ge d + 1, and projecting in it is weighted least squares on the sample values. For example, approximating b=(1,0)b = (1, 0) by t(1,1)t(1, 1) gives t=1/2t = 1/2 with W=IW = I, but t=4/5t = 4/5 with W=diag⁡(4,1)W = \operatorname{diag}(4, 1), where the weighted squared error 4(1−t)2+t24(1 - t)^2 + t^2 is 4/54/5.

Lessons Weighted Least Squares

Spectral Theory

Spectral theory asks which operators become simple in an orthonormal basis and how simple they can be made. Invariant subspaces give block triangular matrices, and Schur’s theorem makes every complex square matrix upper triangular in an orthonormal basis, with a real version that uses diagonal blocks of size one and two. Over the complex numbers, normal operators are exactly the ones with an orthonormal eigenbasis.

Invariant subspaces and block matrices

A subspace UU is invariant under an operator T:V→VT: V \to V if T(U)⊆UT(U) \subseteq U, which can be checked on a spanning list of UU. The subspaces {0}\{0\}, VV, ker⁡T\ker T, im⁡T\operatorname{im} T and every eigenspace are invariant, and a line span⁡(v)\operatorname{span}(v) is invariant exactly when vv is an eigenvector. The restriction T∣UT|_U is then an operator on UU. Over C\mathbb{C} every nonzero finite-dimensional invariant subspace contains an eigenvector, which can fail over R\mathbb{R}.

In a basis (u1,…,ur,w1,…,wn−r)(u_1, \dots, u_r, w_1, \dots, w_{n-r}) whose first rr vectors span UU, the subspace UU is invariant exactly when

[T]B=[BC0D],[T]_{\mathcal{B}} = \begin{bmatrix} B & C \\ 0 & D \end{bmatrix},

with BB the matrix of T∣UT|_U. Then χT=χB χD\chi_T = \chi_B\,\chi_D. The matrix is block diagonal, with C=0C = 0, exactly when the complement spanned by the wjw_j is invariant too, and an invariant complement need not exist.

  • Adjoints. In an inner product space, if UU is invariant under TT, then U⊥U^{\perp} is invariant under T∗T^{*}.
  • Reducing subspaces. The subspace UU reduces TT when both UU and U⊥U^{\perp} are invariant under TT, which holds exactly when UU is invariant under both TT and T∗T^{*}. Then TT is block diagonal in an orthonormal basis adapted to V=U⊕U⊥V = U \oplus U^{\perp}. For a self-adjoint TT every invariant subspace reduces it.

For example, AA with rows (2,1,0)(2, 1, 0), (0,3,0)(0, 3, 0) and (0,0,4)(0, 0, 4) leaves span⁡(e1)\operatorname{span}(e_1) invariant, but not its orthogonal complement, since Ae2=e1+3e2Ae_2 = e_1 + 3e_2. The invariant subspaces span⁡(e1,e2)\operatorname{span}(e_1, e_2) and span⁡(e3)\operatorname{span}(e_3) give a block diagonal description with blocks of sizes 22 and 11.

Lessons Invariant Subspaces and Block Matrices

Schur decomposition and real Schur form

A matrix Q∈Cn×nQ \in \mathbb{C}^{n \times n} is unitary if Q∗Q=IQ^{*}Q = I, which says that its columns are orthonormal. Then Q−1=Q∗Q^{-1} = Q^{*}, ∣det⁡Q∣=1\lvert \det Q \rvert = 1, and a real unitary matrix is orthogonal. Schur’s theorem says that every A∈Cn×nA \in \mathbb{C}^{n \times n} has a factorization

A=QRQ∗,Q∗AQ=R,A = QRQ^{*}, \qquad Q^{*}AQ = R,

with QQ unitary and RR upper triangular. The diagonal of RR lists the eigenvalues of AA, each as often as its algebraic multiplicity, in any chosen order, and the factorization is not unique. The first kk columns of QQ span an invariant subspace for each kk, but only the first column is guaranteed to be an eigenvector. The proof takes an eigenvalue from the fundamental theorem of algebra, completes a unit eigenvector to an orthonormal basis and applies induction to the lower right block, the compression of AA to the orthogonal complement of that eigenvector. If AA is real with only real eigenvalues, QQ and RR can be taken real.

  • Multiplicities. The algebraic multiplicity mjm_j of λj\lambda_j is its exponent in χA\chi_A, the geometric multiplicity is gj=dim⁡ker⁡(A−λjI)g_j = \dim \ker(A - \lambda_j I), and 1≤gj≤mj1 \le g_j \le m_j.
  • Trace and determinant. With the eigenvalues λ1,…,λn\lambda_1, \dots, \lambda_n listed with algebraic multiplicity, tr⁡A=λ1+⋯+λn\operatorname{tr} A = \lambda_1 + \cdots + \lambda_n and det⁡A=λ1⋯λn\det A = \lambda_1 \cdots \lambda_n, with no diagonalizability needed.

The real Schur theorem says that every A∈Rn×nA \in \mathbb{R}^{n \times n} has QTAQ=TQ^{\mathsf{T}}AQ = T for a real orthogonal QQ, where TT is real and block upper triangular with diagonal blocks of size 11, holding the real eigenvalues, and of size 22, each holding a pair of nonreal conjugate eigenvalues. The 2×22 \times 2 blocks come from real invariant planes. If Az=λzAz = \lambda z with λ=a+ib\lambda = a + ib, b≠0b \ne 0 and a nonzero z=x+iyz = x + iy for real xx and yy, then xx and yy are independent and

Ax=ax−by,Ay=bx+ay.\begin{aligned} Ax &= ax - by, \\[0.4em] Ay &= bx + ay. \end{aligned}

A real 2×22 \times 2 matrix with nonreal eigenvalues has the rotation-scaling form with rows (a,−b)(a, -b) and (b,a)(b, a) exactly when it is normal. So a general real Schur block is not a rotation followed by a uniform stretch. For example, CC with rows (1,−2)(1, -2) and (1/2,1)(1/2, 1) has eigenvalues 1±i1 \pm i, but CTCC^{\mathsf{T}}C has rows (5/4,−3/2)(5/4, -3/2) and (−3/2,5)(-3/2, 5), which differs from CCTCC^{\mathsf{T}}, so CC is not normal.

Lessons Schur Decomposition, Real Schur Form

Normal operators and the spectral theorem

An operator TT on a finite-dimensional inner product space is normal if TT∗=T∗TTT^{*} = T^{*}T, and a square matrix is normal if AA∗=A∗AAA^{*} = A^{*}A. The matrix test describes the operator only in an orthonormal basis. The complex spectral theorem says that for a finite-dimensional complex inner product space VV and a linear T:V→VT: V \to V, the space has an orthonormal basis of eigenvectors of TT if and only if TT is normal. In matrices, A∈Cn×nA \in \mathbb{C}^{n \times n} is normal if and only if

A=QΛQ∗=∑j=1nλj qjqj∗A = Q\Lambda Q^{*} = \sum_{j=1}^{n} \lambda_j\, q_j q_j^{*}

with QQ unitary and Λ\Lambda diagonal. The proof applies Schur’s theorem and shows that an upper triangular normal matrix is diagonal. Eigenvectors of a normal operator for distinct eigenvalues are orthogonal. The matrix QQ is not unique, since each column can be multiplied by a scalar of modulus 11 and a repeated eigenvalue allows any orthonormal basis of its eigenspace.

Class of A∈Cn×nA \in \mathbb{C}^{n \times n}DefinitionEigenvaluesOrthonormal eigenbasis
NormalAA∗=A∗AAA^{*} = A^{*}Aany complex numbersyes, over C\mathbb{C}
HermitianA∗=AA^{*} = Arealyes
UnitaryA∗A=IA^{*}A = Imodulus 11yes
Real symmetricAT=AA^{\mathsf{T}} = A, real entriesrealyes, with QQ real orthogonal
Not normalAA∗≠A∗AAA^{*} \ne A^{*}Aanyno

On a complex inner product space an operator is Hermitian exactly when it has an orthonormal eigenbasis with real eigenvalues, and unitary exactly when it has an orthonormal eigenbasis with eigenvalues on the unit circle. For a real symmetric matrix this gives A=QΛQTA = Q\Lambda Q^{\mathsf{T}} with QQ real orthogonal and Λ\Lambda real. A real normal matrix that is not symmetric always has nonreal eigenvalues, since real eigenvalues alone would make it real symmetric, and over R\mathbb{R} it keeps 2×22 \times 2 rotation-scaling blocks.

For example, HH with rows (2,i)(2, i) and (−i,2)(-i, 2) is Hermitian with eigenvalues 11 and 33 and orthonormal eigenvectors (−i,1)/2(-i, 1)/\sqrt{2} and (i,1)/2(i, 1)/\sqrt{2}. The shear AA with rows (1,1)(1, 1) and (0,1)(0, 1) is not normal, since AA∗AA^{*} has rows (2,1)(2, 1) and (1,1)(1, 1) while A∗AA^{*}A has rows (1,1)(1, 1) and (1,2)(1, 2).

Lessons The Spectral Theorem for Normal Operators

Simultaneous diagonalization

Operators SS and TT on the same space commute if ST=TSST = TS. A family is simultaneously diagonalizable when one basis consists of eigenvectors of every member, possibly with different eigenvalues. Commuting is necessary, because diagonal matrices commute. If ST=TSST = TS and Tv=λvTv = \lambda v, then T(Sv)=λSvT(Sv) = \lambda Sv, so SS maps every eigenspace of TT into itself.

The theorem on common orthonormal eigenbases says that for a finite family T1,…,TkT_1, \dots, T_k of operators on a finite-dimensional complex inner product space, there is an orthonormal basis of common eigenvectors if and only if every TjT_j is normal and the TjT_j commute in pairs. In matrices, one unitary QQ then makes every Q∗AjQQ^{*}A_jQ diagonal. The proof splits the space into the orthogonal eigenspaces of T1T_1, which reduce the other operators, and uses that the restriction of a normal operator to a reducing subspace is normal.

Commuting alone is not enough for an orthonormal basis. For example, AA with rows (1,1)(1, 1) and (0,2)(0, 2) commutes with A2A^2, and the two share the eigenbasis (1,0),(1,1)(1, 0), (1, 1), but AA is not normal, so no orthonormal common eigenbasis exists. A matrix that is not diagonalizable commutes with II and has no eigenbasis at all.

Lessons Simultaneous Diagonalization

Quadratic Forms

A Hermitian matrix AA defines a quadratic form x↦x∗Axx \mapsto x^{*}Ax with real values. This unit classifies such forms by the signs they take, factors positive definite matrices by Cholesky’s method, and describes every eigenvalue of a Hermitian matrix as the answer to a minimum or maximum problem.

Congruence and Sylvester’s law of inertia

Let A∈Fn×nA \in \mathbb{F}^{n \times n} be Hermitian, which for F=R\mathbb{F} = \mathbb{R} means real symmetric. Under an invertible change of coordinates x=Syx = Sy, the form becomes y∗(S∗AS) yy^{*}(S^{*}AS)\,y, so the new matrix is S∗ASS^{*}AS, conjugate congruence, which over R\mathbb{R} is the congruence STASS^{\mathsf{T}}AS. Congruence is not similarity and can change the eigenvalues.

  • Diagonalization by congruence. There is an invertible SS with S∗AS=diag⁡(Ip,−Iq,0z)S^{*}AS = \operatorname{diag}(I_p, -I_q, 0_z). One choice is S=QDS = QD, from A=QΛQ∗A = Q\Lambda Q^{*} with the positive eigenvalues listed first, then the negative ones, then the zeros, and a diagonal DD with entries 1/∣λi∣1/\sqrt{\lvert \lambda_i \rvert} for λi≠0\lambda_i \ne 0 and 11 otherwise.
  • Inertia. The triple (p,q,z)(p, q, z) is the inertia of AA. The numbers pp and qq count the positive and the negative eigenvalues with multiplicity, and z=dim⁡ker⁡A=n−rank⁡Az = \dim \ker A = n - \operatorname{rank} A.
  • Positive subspaces. The number pp is the largest dimension of a subspace on which x∗Ax>0x^{*}Ax > 0 for every nonzero xx, and qq is the same for negative values.

Sylvester’s law of inertia says that if A,B∈Fn×nA, B \in \mathbb{F}^{n \times n} are both real symmetric or both complex Hermitian, then B=S∗ASB = S^{*}AS for some invertible S∈Fn×nS \in \mathbb{F}^{n \times n} if and only if AA and BB have the same inertia. Over C\mathbb{C} this is conjugate congruence, since STASS^{\mathsf{T}}AS with a complex SS need not even be Hermitian. Every diagonal matrix S∗ASS^{*}AS with SS invertible has pp positive, qq negative and zz zero diagonal entries.

Form x∗Axx^{*}AxConditionInertia
Positive definitex∗Ax>0x^{*}Ax > 0 for x≠0x \ne 0(n,0,0)(n, 0, 0)
Positive semidefinitex∗Ax≥0x^{*}Ax \ge 0 for all xxq=0q = 0
Negative definitex∗Ax<0x^{*}Ax < 0 for x≠0x \ne 0(0,n,0)(0, n, 0)
Negative semidefinitex∗Ax≤0x^{*}Ax \le 0 for all xxp=0p = 0
Indefinitetakes both signsp≥1p \ge 1 and q≥1q \ge 1

For example, the real form x12+2x1x2+x22−x32=(x1+x2)2−x32x_1^2 + 2x_1x_2 + x_2^2 - x_3^2 = (x_1 + x_2)^2 - x_3^2 has inertia (1,1,1)(1, 1, 1), and its matrix has kernel spanned by (1,−1,0)(1, -1, 0). The invertible matrix with rows (0,1)(0, 1) and (1,0)(1, 0) has inertia (1,1,0)(1, 1, 0), and its form 2x1x22x_1x_2 still vanishes on nonzero vectors such as (1,0)(1, 0).

Lessons Quadratic Forms and Sylvester’s Law of Inertia

Positive definite matrices and Cholesky factorization

A matrix A∈Fn×nA \in \mathbb{F}^{n \times n} is positive definite if it is Hermitian and x∗Ax>0x^{*}Ax > 0 for every nonzero xx, and positive semidefinite if it is Hermitian and x∗Ax≥0x^{*}Ax \ge 0 for every xx. Over R\mathbb{R} symmetry is part of the definition, because xTAxx^{\mathsf{T}}Ax sees only the symmetric part.

For a Hermitian AA, positive definite is equivalent toPositive semidefinite is equivalent to
every eigenvalue positiveevery eigenvalue nonnegative
inertia (n,0,0)(n, 0, 0)inertia (r,0,n−r)(r, 0, n - r) with r=rank⁡Ar = \operatorname{rank} A
A=B∗BA = B^{*}B for an invertible BBA=B∗BA = B^{*}B for some BB with nn columns
every leading principal minor positivenot decided by the leading minors

Completing the square splits off the first coordinate. If AA is Hermitian with first entry a≠0a \ne 0, first column (a,b)(a, b), first row (a,b∗)(a, b^{*}) and lower right block CC, then for x=(x1,y)x = (x_1, y)

x∗Ax=a ∣x1+b∗ya∣2+y∗Sy,S=C−1a bb∗.\begin{aligned} x^{*}Ax &= a\,\Bigl\lvert x_1 + \frac{b^{*}y}{a} \Bigr\rvert^2 + y^{*}Sy, \\[0.4em] S &= C - \frac{1}{a}\,bb^{*}. \end{aligned}

The matrix SS is the Schur complement of aa, and det⁡A=adet⁡S\det A = a \det S. If AA is positive definite, then a>0a > 0 and SS is positive definite.

The Cholesky factorization theorem says that every positive definite A∈Fn×nA \in \mathbb{F}^{n \times n} has exactly one factorization

A=R∗RA = R^{*}R

with RR upper triangular and every diagonal entry real and positive. Over R\mathbb{R} it reads A=RTRA = R^{\mathsf{T}}R, and many books write A=LL∗A = LL^{*} with L=R∗L = R^{*}. Sylvester’s criterion says that a Hermitian AA is positive definite if and only if every leading principal minor Δk\Delta_k, the determinant of the upper left k×kk \times k block, is positive, and then rkk2=Δk/Δk−1r_{kk}^2 = \Delta_k / \Delta_{k-1} with Δ0=1\Delta_0 = 1. The Hermitian hypothesis matters, since the matrix with rows (1,−3)(1, -3) and (0,1)(0, 1) has Δ1=Δ2=1\Delta_1 = \Delta_2 = 1 but gives −1-1 at x=(1,1)x = (1, 1). Nonnegative leading minors do not prove semidefiniteness, as diag⁡(0,−1)\operatorname{diag}(0, -1) shows.

For example, AA with rows (2,1)(2, 1) and (1,2)(1, 2) has xTAx=2(x1+x2/2)2+32x22x^{\mathsf{T}}Ax = 2(x_1 + x_2/2)^2 + \tfrac{3}{2}x_2^2, leading minors 22 and 33, and

R=[21/203/2].R = \begin{bmatrix} \sqrt{2} & 1/\sqrt{2} \\ 0 & \sqrt{3/2} \end{bmatrix}.

Lessons Positive Definite Matrices and Cholesky Factorization

Rayleigh quotients and the min-max principle

For a Hermitian A∈Fn×nA \in \mathbb{F}^{n \times n} and nonzero x∈Fnx \in \mathbb{F}^n, the Rayleigh quotient is

RA(x)=x∗Axx∗x.R_A(x) = \frac{x^{*}Ax}{x^{*}x}.

It is real, and RA(cx)=RA(x)R_A(cx) = R_A(x) for every nonzero scalar cc. With eigenvalues λ1≤⋯≤λn\lambda_1 \le \cdots \le \lambda_n and an orthonormal eigenbasis q1,…,qnq_1, \dots, q_n in the same order, it is a weighted average of the eigenvalues,

RA(x)=∑j=1nwjλj,wj=∣qj∗x∣2∥x∥2,R_A(x) = \sum_{j=1}^{n} w_j \lambda_j, \qquad w_j = \frac{\lvert q_j^{*}x \rvert^2}{\lVert x \rVert^2},

with wj≥0w_j \ge 0 and ∑jwj=1\sum_j w_j = 1. Hence λ1≤RA(x)≤λn\lambda_1 \le R_A(x) \le \lambda_n, with RA(x)=λ1R_A(x) = \lambda_1 exactly when Ax=λ1xAx = \lambda_1 x and RA(x)=λnR_A(x) = \lambda_n exactly when Ax=λnxAx = \lambda_n x. When an extreme eigenvalue is repeated, its extremizers are not one line but all nonzero vectors of an eigenspace of dimension at least 22. The Hermitian hypothesis is needed, since for a real nonsymmetric matrix the quotient sees only the symmetric part.

Two subspaces of an nn-dimensional space satisfy dim⁡(U∩W)≥dim⁡U+dim⁡W−n\dim(U \cap W) \ge \dim U + \dim W - n, so they share a nonzero vector when their dimensions add up to more than nn. This dimension count proves the min-max principle, also called the Courant-Fischer theorem. For a Hermitian AA and k=1,…,nk = 1, \dots, n,

λk=min⁡dim⁡S=k max⁡0≠x∈SRA(x)=max⁡dim⁡L=n−k+1 min⁡0≠x∈LRA(x),\begin{aligned} \lambda_k &= \min_{\dim S = k} \ \max_{0 \ne x \in S} R_A(x) \\[0.4em] &= \max_{\dim L = n - k + 1} \ \min_{0 \ne x \in L} R_A(x), \end{aligned}

where SS and LL run over subspaces of Fn\mathbb{F}^n. The subspaces span⁡(q1,…,qk)\operatorname{span}(q_1, \dots, q_k) and span⁡(qk,…,qn)\operatorname{span}(q_k, \dots, q_n) attain the two values, but attaining subspaces are not unique in general.

For example, A=diag⁡(1,4,9)A = \operatorname{diag}(1, 4, 9) has RA(x)=5/2R_A(x) = 5/2 at x=(1,1,0)/2x = (1, 1, 0)/\sqrt{2}. On span⁡(e1,e2)\operatorname{span}(e_1, e_2) the largest value is 44, and every two-dimensional subspace contains a nonzero vector of span⁡(e2,e3)\operatorname{span}(e_2, e_3), where the quotient is at least 44, so λ2=4\lambda_2 = 4.

Lessons Rayleigh Quotients, The Min-Max Principle

Singular Value Decomposition

The singular value decomposition writes every rectangular matrix as an orthonormal change of input coordinates, a scaling by nonnegative singular values and an orthonormal change of output coordinates. On it rest the pseudoinverse with its four orthogonal projectors, the best approximations of lower rank, principal component analysis of centered data, positive square roots and the polar decomposition.

The singular value decomposition

Let A∈Fm×nA \in \mathbb{F}^{m \times n} have rank rr, and let p=min⁡(m,n)p = \min(m, n). The Gram operator A∗AA^{*}A is Hermitian and positive semidefinite, since x∗A∗Ax=∥Ax∥2x^{*}A^{*}Ax = \lVert Ax \rVert^2, with ker⁡A∗A=ker⁡A\ker A^{*}A = \ker A and rank⁡A∗A=r\operatorname{rank} A^{*}A = r. The singular values σ1≥⋯≥σp≥0\sigma_1 \ge \cdots \ge \sigma_p \ge 0 are the nonnegative square roots of the pp largest eigenvalues of A∗AA^{*}A, and exactly rr of them are positive.

The theorem says that every A∈Fm×nA \in \mathbb{F}^{m \times n} has a singular value decomposition

A=UΣV∗A = U\Sigma V^{*}

with U∈Fm×mU \in \mathbb{F}^{m \times m} and V∈Fn×nV \in \mathbb{F}^{n \times n} unitary, which means real orthogonal when F=R\mathbb{F} = \mathbb{R}, and Σ∈Rm×n\Sigma \in \mathbb{R}^{m \times n} zero except for Σii=σi\Sigma_{ii} = \sigma_i, i≤pi \le p. To build it, take an orthonormal eigenbasis v1,…,vnv_1, \dots, v_n of A∗AA^{*}A ordered by decreasing eigenvalue, so that A∗Avi=σi2viA^{*}Av_i = \sigma_i^2 v_i for i≤pi \le p, set ui=Avi/σiu_i = Av_i/\sigma_i for i≤ri \le r, and complete the uiu_i to an orthonormal basis of Fm\mathbb{F}^m. Then Avi=σiuiAv_i = \sigma_i u_i and A∗ui=σiviA^{*}u_i = \sigma_i v_i for i≤pi \le p, and Avi=0Av_i = 0 for i>ri > r.

FormLeft factorMiddle factorRight factor
FullUU, size m×mm \times mΣ\Sigma, size m×nm \times nVV, size n×nn \times n
ThinUpU_p, size m×pm \times pΣp\Sigma_p, size p×pp \times pVpV_p, size n×pn \times p
CompactUrU_r, size m×rm \times rΣr\Sigma_r, size r×rr \times rVrV_r, size n×rn \times r

Here UkU_k and VkV_k hold the first kk columns and Σk=diag⁡(σ1,…,σk)\Sigma_k = \operatorname{diag}(\sigma_1, \dots, \sigma_k), and in each form AA is the left factor times the middle factor times the adjoint of the right factor. The compact form is A=∑i=1rσi uivi∗A = \sum_{i=1}^{r} \sigma_i\, u_i v_i^{*}, and the thin form may still contain zero singular values. The singular vectors give orthonormal bases of the four subspaces.

SubspaceOrthonormal basis
im⁡A\operatorname{im} Au1,…,uru_1, \dots, u_r
ker⁡A∗\ker A^{*}ur+1,…,umu_{r+1}, \dots, u_m
im⁡A∗=(ker⁡A)⊥\operatorname{im} A^{*} = (\ker A)^{\perp}v1,…,vrv_1, \dots, v_r
ker⁡A\ker Avr+1,…,vnv_{r+1}, \dots, v_n

The singular values are determined by AA, but the singular vectors are not. For a positive singular value σ\sigma repeated kk times, the kk matching columns UσU_\sigma of UU and VσV_\sigma of VV can be replaced by UσQU_\sigma Q and VσQV_\sigma Q for any unitary Q∈Fk×kQ \in \mathbb{F}^{k \times k}, and for k=1k = 1 this multiplies uiu_i and viv_i by the same scalar of modulus 11. The remaining columns can be any orthonormal bases of ker⁡A∗\ker A^{*} and ker⁡A\ker A. Singular values are not the moduli of eigenvalues in general. They agree for normal matrices, but the shear with rows (1,1)(1, 1) and (0,1)(0, 1) has eigenvalues 1,11, 1 and singular values (5±1)/2(\sqrt{5} \pm 1)/2.

For example, AA with rows (3,0)(3, 0), (0,0)(0, 0) and (0,2)(0, 2) has singular values 33 and 22, with V=I2V = I_2, UU with columns e1,e3,e2e_1, e_3, e_2, and Σ\Sigma with rows (3,0)(3, 0), (0,2)(0, 2) and (0,0)(0, 0). The left null direction e2e_2 has no positive singular value.

Lessons The Singular Value Decomposition

The pseudoinverse and the four projectors

For A∈Fm×nA \in \mathbb{F}^{m \times n} of rank rr with A=UΣV∗A = U\Sigma V^{*}, let Σ+∈Rn×m\Sigma^{+} \in \mathbb{R}^{n \times m} hold 1/σi1/\sigma_i for i≤ri \le r and zeros elsewhere. The pseudoinverse is

A+=VΣ+U∗=∑i=1r1σi viui∗∈Fn×m,\begin{aligned} A^{+} &= V\Sigma^{+}U^{*} \\[0.4em] &= \sum_{i=1}^{r} \frac{1}{\sigma_i}\, v_i u_i^{*} \in \mathbb{F}^{n \times m}, \end{aligned}

and it does not depend on which singular value decomposition is used.

CasePseudoinverse
AA square and invertibleA−1A^{-1}
full column rank, r=nr = n(A∗A)−1A∗(A^{*}A)^{-1}A^{*}
full row rank, r=mr = mA∗(AA∗)−1A^{*}(AA^{*})^{-1}
A=0A = 0the n×mn \times m zero matrix

A vector x^\hat{x} minimizes ∥Ax−b∥\lVert Ax - b \rVert exactly when Ax^=AA+bA\hat{x} = AA^{+}b, or equivalently when A∗Ax^=A∗bA^{*}A\hat{x} = A^{*}b. Least-squares solutions always exist and form the set A+b+ker⁡AA^{+}b + \ker A. Among them A+bA^{+}b is the only one of smallest norm, because it lies in (ker⁡A)⊥(\ker A)^{\perp} and ∥A+b+z∥2=∥A+b∥2+∥z∥2\lVert A^{+}b + z \rVert^2 = \lVert A^{+}b \rVert^2 + \lVert z \rVert^2 for z∈ker⁡Az \in \ker A.

ProjectorSizeProjects orthogonally onto
AA+=UrUr∗AA^{+} = U_rU_r^{*}m×mm \times mim⁡A\operatorname{im} A, the column space
Im−AA+I_m - AA^{+}m×mm \times mker⁡A∗\ker A^{*}, where the residuals lie
A+A=VrVr∗A^{+}A = V_rV_r^{*}n×nn \times nim⁡A∗=(ker⁡A)⊥\operatorname{im} A^{*} = (\ker A)^{\perp}, spanned by the conjugated rows of AA
In−A+AI_n - A^{+}An×nn \times nker⁡A\ker A

The Moore-Penrose theorem says that X=A+X = A^{+} is the only matrix X∈Fn×mX \in \mathbb{F}^{n \times m} with

AXA=A,XAX=X,(AX)∗=AX,(XA)∗=XA.\begin{aligned} &AXA = A, \qquad XAX = X, \\[0.4em] &(AX)^{*} = AX, \qquad (XA)^{*} = XA. \end{aligned}

The first condition makes AXAX and XAXA idempotent, the last two make them Hermitian and hence orthogonal projectors, and all four are needed. A matrix with only AXA=AAXA = A is a generalized inverse. It solves every consistent system, but it is unique only when AA is square and invertible. For A=diag⁡(1,0)A = \operatorname{diag}(1, 0), every XX with rows (1,t)(1, t) and (0,0)(0, 0) satisfies the first two conditions, but AXAX is not symmetric when t≠0t \ne 0.

For example, AA with rows (1,1)(1, 1), (0,0)(0, 0) and (0,0)(0, 0) has A+A^{+} with rows (1/2,0,0)(1/2, 0, 0) and (1/2,0,0)(1/2, 0, 0). For b=(2,3,0)b = (2, 3, 0) the least-squares solutions are (1+s,1−s)(1 + s, 1 - s), the minimum-norm one is A+b=(1,1)A^{+}b = (1, 1), the fitted vector is (2,0,0)(2, 0, 0) and the residual is (0,3,0)(0, 3, 0).

Lessons The Pseudoinverse and Four Projectors, The Moore-Penrose Conditions

Best low-rank approximation

Two matrix norms measure an error E∈Fm×nE \in \mathbb{F}^{m \times n} with singular values σi(E)\sigma_i(E). Both are unchanged by unitary factors on either side.

∥E∥2=max⁡∥x∥=1∥Ex∥=σ1(E),∥E∥F2=∑i,j∣eij∣2=∑iσi(E)2.\begin{aligned} \lVert E \rVert_2 &= \max_{\lVert x \rVert = 1} \lVert Ex \rVert = \sigma_1(E), \\[0.4em] \lVert E \rVert_F^2 &= \sum_{i,j} \lvert e_{ij} \rvert^2 = \sum_{i} \sigma_i(E)^2. \end{aligned}

The spectral norm is the error in the worst input direction and the squared Frobenius norm is the total squared error, with ∥E∥2≤∥E∥F≤rank⁡E ∥E∥2\lVert E \rVert_2 \le \lVert E \rVert_F \le \sqrt{\operatorname{rank} E}\,\lVert E \rVert_2. The truncated singular value decomposition keeps the first kk terms,

Ak=∑i=1kσi uivi∗,0≤k≤p.A_k = \sum_{i=1}^{k} \sigma_i\, u_i v_i^{*}, \qquad 0 \le k \le p.

The Eckart-Young theorem says that for every B∈Fm×nB \in \mathbb{F}^{m \times n} with rank⁡B≤k\operatorname{rank} B \le k,

∥A−B∥2≥σk+1,∥A−B∥F2≥σk+12+⋯+σp2,\begin{aligned} \lVert A - B \rVert_2 &\ge \sigma_{k+1}, \\[0.4em] \lVert A - B \rVert_F^2 &\ge \sigma_{k+1}^2 + \cdots + \sigma_p^2, \end{aligned}

with σp+1=0\sigma_{p+1} = 0, and that B=AkB = A_k attains both bounds. For 1≤k<p1 \le k < p the matrices of rank at most kk do not form a subspace, so the projection theorem does not apply directly.

  • Frobenius uniqueness. For 1≤k<p1 \le k < p, the truncation AkA_k is the only best approximation in the Frobenius norm when σk>σk+1\sigma_k > \sigma_{k+1}, and there are infinitely many when σk=σk+1>0\sigma_k = \sigma_{k+1} > 0.
  • Spectral nonuniqueness. When k≥1k \ge 1 and σk+1>0\sigma_{k+1} > 0 the spectral norm always has infinitely many best approximations, since changing each σi\sigma_i with i≤ki \le k by at most σk+1\sigma_{k+1} keeps the error σk+1\sigma_{k+1}.

For example, A=diag⁡(3,1)A = \operatorname{diag}(3, 1) and k=1k = 1 give spectral error 11 for every Bt=diag⁡(t,0)B_t = \operatorname{diag}(t, 0) with 2≤t≤42 \le t \le 4, while the Frobenius error (3−t)2+1\sqrt{(3 - t)^2 + 1} is smallest only at t=3t = 3. For I2I_2 every rank-one orthogonal projector ww∗ww^{*} is a best approximation in both norms.

Lessons The Eckart-Young Theorem

Principal component analysis

Let x1,…,xN∈Rdx_1, \dots, x_N \in \mathbb{R}^d be N>1N > 1 observations with mean xˉ\bar{x}. The centered data matrix X∈RN×dX \in \mathbb{R}^{N \times d} has rows (xi−xˉ)T(x_i - \bar{x})^{\mathsf{T}}, so rows are observations, columns are measured coordinates and every column sums to 00. The sample covariance matrix is

C=1N−1XTX,C = \frac{1}{N - 1}X^{\mathsf{T}}X,

which is symmetric and positive semidefinite. The divisor N−1N - 1 is a convention, and dividing by NN changes every variance by the same factor and none of the directions. Along a unit vector vv the scores XvXv have sample variance ∥Xv∥2/(N−1)=vTCv\lVert Xv \rVert^2/(N - 1) = v^{\mathsf{T}}Cv, a Rayleigh quotient of CC.

  • Principal directions. With a full real singular value decomposition X=UΣVTX = U\Sigma V^{\mathsf{T}}, the right singular vectors v1,v2,…v_1, v_2, \dots are the principal directions. For i≤p=min⁡(N,d)i \le p = \min(N, d) the variances along them are σi2/(N−1)\sigma_i^2/(N - 1) and the score vectors are Xvi=σiuiXv_i = \sigma_i u_i.
  • Reconstruction. For a subspace with orthonormal basis given by the columns of VkV_k and projector P=VkVkTP = V_kV_k^{\mathsf{T}}, the identity ∥X∥F2=∥XP∥F2+∥X−XP∥F2\lVert X \rVert_F^2 = \lVert XP \rVert_F^2 + \lVert X - XP \rVert_F^2 shows that maximizing retained variance minimizes reconstruction error. The top kk directions give the smallest error, σk+12+⋯+σp2\sigma_{k+1}^2 + \cdots + \sigma_p^2, by the Eckart-Young theorem.
  • Choices. Without centering, XrawTXraw=XTX+NxˉxˉTX_{\text{raw}}^{\mathsf{T}}X_{\text{raw}} = X^{\mathsf{T}}X + N\bar{x}\bar{x}^{\mathsf{T}}, and the leading direction can describe the mean instead of the spread. Changing units changes the inner product and can change the answer. Each viv_i can be replaced by −vi-v_i, and when singular values tie, any orthonormal basis of the tied eigenspace of CC serves equally well.

For example, the centered rows (2,0)(2, 0), (−2,0)(-2, 0), (0,1)(0, 1) and (0,−1)(0, -1) give XTX=diag⁡(8,2)X^{\mathsf{T}}X = \operatorname{diag}(8, 2), variances 8/38/3 and 2/32/3, and the horizontal line as the best one-dimensional fit, with squared error 22. Shifting every point up by 1010 leaves the centered answer unchanged but makes the uncentered Gram matrix diag⁡(8,402)\operatorname{diag}(8, 402), whose leading direction is vertical.

Lessons Principal Component Analysis

Positive square roots and polar decomposition

The positive square root theorem says that a Hermitian positive semidefinite H∈Fn×nH \in \mathbb{F}^{n \times n} has exactly one Hermitian positive semidefinite SS with S2=HS^2 = H. If H=QΛQ∗H = Q\Lambda Q^{*}, it is

H1/2=Q diag⁡(λ1,…,λn) Q∗,H^{1/2} = Q\,\operatorname{diag}\bigl(\sqrt{\lambda_1}, \dots, \sqrt{\lambda_n}\bigr)\,Q^{*},

and ker⁡H1/2=ker⁡H\ker H^{1/2} = \ker H. Without the positivity requirement square roots need not be unique, and a Hermitian matrix with a negative eigenvalue has no Hermitian square root at all. A Cholesky factor solves a different equation, H=R∗RH = R^{*}R with RR triangular, and is usually not Hermitian.

The absolute value of A∈Fm×nA \in \mathbb{F}^{m \times n} is ∣A∣=(A∗A)1/2∈Fn×n\lvert A \rvert = (A^{*}A)^{1/2} \in \mathbb{F}^{n \times n}. It satisfies ∥∣A∣x∥=∥Ax∥\lVert \lvert A \rvert x \rVert = \lVert Ax \rVert for every xx, so ker⁡∣A∣=ker⁡A\ker \lvert A \rvert = \ker A. The polar decomposition theorem says that exactly one U∈Fm×nU \in \mathbb{F}^{m \times n} satisfies

A=U∣A∣,Ux=0  for every x∈ker⁡A.\begin{gathered} A = U\lvert A \rvert, \\[0.4em] Ux = 0 \ \text{ for every } x \in \ker A. \end{gathered}

In a compact singular value decomposition, ∣A∣=∑i≤rσivivi∗\lvert A \rvert = \sum_{i \le r} \sigma_i v_iv_i^{*} and U=∑i≤ruivi∗U = \sum_{i \le r} u_iv_i^{*}.

FactorFormulaUniqueness
∣A∣\lvert A \rvert, size n×nn \times n(A∗A)1/2(A^{*}A)^{1/2}determined by AA alone
Canonical UU, size m×nm \times n∑i≤ruivi∗\sum_{i \le r} u_iv_i^{*}unique with U=0U = 0 on ker⁡A\ker A
Unitary WW with A=W∣A∣A = W\lvert A \rvertexists when m=nm = nunique only when AA is invertible
KK in A=KUA = KU, size m×mm \times m(AA∗)1/2=U∣A∣U∗(AA^{*})^{1/2} = U\lvert A \rvert U^{*}determined by AA alone

The canonical UU is a partial isometry, with U∗U=A+AU^{*}U = A^{+}A and UU∗=AA+UU^{*} = AA^{+}. It is an isometry, U∗U=InU^{*}U = I_n, exactly when AA has full column rank, and unitary exactly when AA is square and invertible, in which case U=A∣A∣−1U = A\lvert A \rvert^{-1}.

For example, HH with rows (5,4)(5, 4) and (4,5)(4, 5) has eigenvalues 99 and 11, and its positive root has rows (2,1)(2, 1) and (1,2)(1, 2). The matrix with rows (1,2)(1, 2) and (2,1)(2, 1) also squares to HH but is indefinite. The matrix AA with rows (0,−1)(0, -1) and (2,0)(2, 0) has A∗A=diag⁡(4,1)A^{*}A = \operatorname{diag}(4, 1), ∣A∣=diag⁡(2,1)\lvert A \rvert = \operatorname{diag}(2, 1) and the quarter-turn UU with rows (0,−1)(0, -1) and (1,0)(1, 0). For the singular A=diag⁡(2,0)A = \operatorname{diag}(2, 0), the canonical factor is diag⁡(1,0)\operatorname{diag}(1, 0), and both diag⁡(1,1)\operatorname{diag}(1, 1) and diag⁡(1,−1)\operatorname{diag}(1, -1) are unitary extensions.

Lessons Positive Square Roots of Matrices, Polar Decomposition

Canonical Forms

Some operators have too few eigenvectors to be diagonalized, and polynomials in the operator measure exactly what is missing. The minimal polynomial decides diagonalizability, polynomial projectors split the space into generalized eigenspaces, and chains of generalized eigenvectors give the Jordan form. It classifies complex matrices up to similarity, but it is unreliable when the entries are known only approximately.

The minimal polynomial and the Cayley-Hamilton theorem

For f(t)=a0+a1t+⋯+adtdf(t) = a_0 + a_1 t + \cdots + a_d t^d with coefficients in F\mathbb{F} and A∈Fn×nA \in \mathbb{F}^{n \times n}, the matrix f(A)=a0I+a1A+⋯+adAdf(A) = a_0 I + a_1 A + \cdots + a_d A^d satisfies (fg)(A)=f(A)g(A)(fg)(A) = f(A)g(A), so polynomials in the same AA commute. A nonzero ff with f(A)=0f(A) = 0 is an annihilating polynomial. One exists because the n2+1n^2 + 1 matrices I,A,…,An2I, A, \dots, A^{n^2} are dependent.

  • Minimal polynomial. The minimal polynomial μA\mu_A is the unique monic annihilating polynomial of least degree, and f(A)=0f(A) = 0 exactly when μA\mu_A divides ff.
  • Roots. A scalar λ∈F\lambda \in \mathbb{F} is a root of μA\mu_A exactly when it is an eigenvalue of AA.
  • Cayley-Hamilton. Every A∈Fn×nA \in \mathbb{F}^{n \times n} satisfies χA(A)=0\chi_A(A) = 0, so μA\mu_A divides χA\chi_A and deg⁡μA≤n\deg \mu_A \le n. The proof compares coefficients in (tI−A)adj⁡(tI−A)=χA(t) I(tI - A)\operatorname{adj}(tI - A) = \chi_A(t)\,I. Substituting t=At = A into det⁡(tI−A)\det(tI - A) is not a proof.
MatrixχA(t)\chi_A(t)μA(t)\mu_A(t)
I2I_2(t−1)2(t - 1)^2t−1t - 1
shear with rows (1,1)(1, 1) and (0,1)(0, 1)(t−1)2(t - 1)^2(t−1)2(t - 1)^2
diag⁡(2,2,3)\operatorname{diag}(2, 2, 3)(t−2)2(t−3)(t - 2)^2(t - 3)(t−2)(t−3)(t - 2)(t - 3)
quarter-turn with rows (0,−1)(0, -1) and (1,0)(1, 0)t2+1t^2 + 1t2+1t^2 + 1

A repeated root in some annihilating polynomial does not show that a matrix is defective. Only the minimal polynomial decides. High powers reduce by division. For AA with rows (2,1)(2, 1) and (0,3)(0, 3), Cayley-Hamilton gives A2=5A−6IA^2 = 5A - 6I, then A3=19A−30IA^3 = 19A - 30I and

A4=65A−114I=[1665081].A^4 = 65A - 114I = \begin{bmatrix} 16 & 65 \\ 0 & 81 \end{bmatrix}.

Lessons Matrix Polynomials and the Minimal Polynomial, The Cayley-Hamilton Theorem

Polynomial projectors and diagonalizability

Repeated division gives the greatest common divisor of two polynomials as the last nonzero remainder, made monic, and back substitution gives Bézout’s identity. If aa and bb are coprime, there are polynomials uu and vv with

u(t)a(t)+v(t)b(t)=1.u(t)a(t) + v(t)b(t) = 1.

If moreover a(A) b(A)=0a(A)\,b(A) = 0 for A∈Fn×nA \in \mathbb{F}^{n \times n}, then

Fn=ker⁡a(A)⊕ker⁡b(A),\mathbb{F}^n = \ker a(A) \oplus \ker b(A),

and P=v(A)b(A)P = v(A)b(A) is the projection onto ker⁡a(A)\ker a(A) along ker⁡b(A)\ker b(A), with I−P=u(A)a(A)I - P = u(A)a(A). These projectors are polynomials in AA and need not be orthogonal.

The diagonalizability criterion says that A∈Fn×nA \in \mathbb{F}^{n \times n} is diagonalizable over F\mathbb{F} if and only if μA\mu_A splits over F\mathbb{F} into distinct linear factors, μA(t)=(t−λ1)⋯(t−λs)\mu_A(t) = (t - \lambda_1) \cdots (t - \lambda_s). The eigenspace projectors are then the Lagrange polynomials

Pj=∏k≠jA−λkIλj−λk,P_j = \prod_{k \ne j} \frac{A - \lambda_k I}{\lambda_j - \lambda_k},

with ∑jPj=I\sum_j P_j = I, PjPk=0P_jP_k = 0 for j≠kj \ne k and A=∑jλjPjA = \sum_j \lambda_j P_j. Splitting is part of the hypothesis. The real quarter-turn has μ=t2+1\mu = t^2 + 1 with distinct complex roots and is diagonalizable over C\mathbb{C} but not over R\mathbb{R}. A finite family of commuting operators, each diagonalizable over F\mathbb{F}, has a common eigenbasis, which need not be orthonormal.

For example, AA with rows (1,1)(1, 1) and (0,2)(0, 2) has P1=2I−AP_1 = 2I - A and P2=A−IP_2 = A - I, with rows (1,−1),(0,0)(1, -1), (0, 0) and (0,1),(0,1)(0, 1), (0, 1). They split x=(2,3)x = (2, 3) as (−1,0)+(3,3)(-1, 0) + (3, 3), an eigenvector for 11 plus an eigenvector for 22, along lines that are not perpendicular.

Lessons Polynomial Projectors and Diagonalizability

Cyclic subspaces and companion matrices

The orbit of a nonzero vv under TT is v,Tv,T2v,…v, Tv, T^2v, \dots. If dd is the smallest number for which v,Tv,…,Tdvv, Tv, \dots, T^dv are dependent, the first dependence

Tdv=−a0v−a1Tv−⋯−ad−1Td−1vT^d v = -a_0 v - a_1 Tv - \cdots - a_{d-1}T^{d-1}v

defines the monic polynomial p(t)=td+ad−1td−1+⋯+a0p(t) = t^d + a_{d-1}t^{d-1} + \cdots + a_0. The cyclic subspace CT(v)=span⁡(v,Tv,T2v,… )\mathcal{C}_T(v) = \operatorname{span}(v, Tv, T^2v, \dots) has basis v,Tv,…,Td−1vv, Tv, \dots, T^{d-1}v, is invariant under TT, and is the smallest invariant subspace containing vv. A polynomial qq satisfies q(T)v=0q(T)v = 0 exactly when pp divides qq, so pp divides μT\mu_T. In the orbit basis the restriction of TT has the companion matrix

Cp=[00⋯0−a010⋯0−a101⋯0−a2⋮⋮⋱⋮⋮00⋯1−ad−1],C_p = \begin{bmatrix} 0 & 0 & \cdots & 0 & -a_0 \\ 1 & 0 & \cdots & 0 & -a_1 \\ 0 & 1 & \cdots & 0 & -a_2 \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ 0 & 0 & \cdots & 1 & -a_{d-1} \end{bmatrix},

whose characteristic and minimal polynomials both equal pp.

  • Cyclic vectors. If CT(v)=V\mathcal{C}_T(v) = V, then vv is a cyclic vector and χT=μT=p\chi_T = \mu_T = p. A scalar matrix cIcI on a space of dimension at least 22 has no cyclic vector.
  • Recurrences. The orbit vectors xk=Tkvx_k = T^kv satisfy xk+d=−ad−1xk+d−1−⋯−a0xkx_{k+d} = -a_{d-1}x_{k+d-1} - \cdots - a_0x_k, and so does every scalar sequence ℓ(Tkv)\ell(T^kv) with ℓ∈V∨\ell \in V^{\vee}.

For example, AA with rows (0,−2)(0, -2) and (1,3)(1, 3) sends e1e_1 to e2e_2 and e2e_2 to −2e1+3e2-2e_1 + 3e_2, so A2−3A+2I=0A^2 - 3A + 2I = 0 and AA is the companion matrix of t2−3t+2t^2 - 3t + 2. The seed (1,−1)(1, -1) is an eigenvector for 22 and generates only a line.

Lessons Cyclic Subspaces and Companion Matrices

Generalized eigenspaces

Fix an eigenvalue λ∈F\lambda \in \mathbb{F} of an operator TT on an nn-dimensional space and let N=T−λIN = T - \lambda I. A generalized eigenvector is a nonzero vv with Nkv=0N^kv = 0 for some k≥1k \ge 1, and its rank is the smallest such kk. Rank one means an ordinary eigenvector. A vector of rank kk generates a Jordan chain, listed bottom first,

v1=Nk−1v, …, vk=v,Nv1=0,Nvj=vj−1,\begin{gathered} v_1 = N^{k-1}v, \ \dots, \ v_k = v, \\[0.4em] Nv_1 = 0, \quad Nv_j = v_{j-1}, \end{gathered}

so Tv1=λv1Tv_1 = \lambda v_1 and Tvj=λvj+vj−1Tv_j = \lambda v_j + v_{j-1} for j≥2j \ge 2. Every chain is linearly independent. The kernels ker⁡N⊆ker⁡N2⊆⋯\ker N \subseteq \ker N^2 \subseteq \cdots stop growing after the first equality, and by step nn at the latest, so the generalized eigenspace is

Gλ=ker⁡(T−λI)n.G_\lambda = \ker(T - \lambda I)^n.

It is invariant under TT and contains the eigenspace.

The primary decomposition theorem says that if μT(t)=(t−λ1)s1⋯(t−λr)sr\mu_T(t) = (t - \lambda_1)^{s_1} \cdots (t - \lambda_r)^{s_r} splits over F\mathbb{F}, with distinct λj\lambda_j, which always holds over C\mathbb{C}, then

V=Gλ1⊕⋯⊕Gλr,V = G_{\lambda_1} \oplus \cdots \oplus G_{\lambda_r},

with Gλj=ker⁡(T−λjI)sjG_{\lambda_j} = \ker(T - \lambda_j I)^{s_j}. On GλjG_{\lambda_j} the operator T−λjIT - \lambda_j I is nilpotent with index exactly sjs_j, and dim⁡Gλj\dim G_{\lambda_j} is the algebraic multiplicity of λj\lambda_j. The projectors come from Bézout’s identity for the coprime factors (t−λj)sj(t - \lambda_j)^{s_j}. Over R\mathbb{R} the splitting hypothesis matters, since a real quarter-turn has no real eigenvalue and no nonzero real generalized eigenspace.

For example, AA with rows (2,1,0)(2, 1, 0), (0,2,1)(0, 2, 1) and (0,0,2)(0, 0, 2) has the chain e3↦e2↦e1↦0e_3 \mapsto e_2 \mapsto e_1 \mapsto 0 under N=A−2IN = A - 2I, kernel dimensions 1,2,31, 2, 3, and A2e3=(4I+4N+N2)e3=(1,4,4)A^2e_3 = (4I + 4N + N^2)e_3 = (1, 4, 4). The matrix with rows (1,1,0)(1, 1, 0), (0,1,0)(0, 1, 0) and (0,0,2)(0, 0, 2) has μA=(t−1)2(t−2)\mu_A = (t - 1)^2(t - 2), G1=span⁡(e1,e2)G_1 = \operatorname{span}(e_1, e_2) and G2=span⁡(e3)G_2 = \operatorname{span}(e_3), and (A−I)2=diag⁡(0,0,1)(A - I)^2 = \operatorname{diag}(0, 0, 1) projects onto G2G_2.

Lessons Generalized Eigenvectors, Generalized Eigenspace Decomposition

Jordan form

A Jordan block of size kk has λ\lambda on the diagonal and 11 directly above it,

Jk(λ)=[λ1λ⋱⋱1λ],J_k(\lambda) = \begin{bmatrix} \lambda & 1 & & \\ & \lambda & \ddots & \\ & & \ddots & 1 \\ & & & \lambda \end{bmatrix},

which is the matrix of the operator on the span of one Jordan chain, ordered bottom first. The Jordan form theorem says that if χA\chi_A splits over F\mathbb{F}, which always holds for F=C\mathbb{F} = \mathbb{C}, then P−1AP=JP^{-1}AP = J for an invertible P∈Fn×nP \in \mathbb{F}^{n \times n} and a block diagonal JJ of Jordan blocks. The number of blocks of each size for each eigenvalue is unique, so JJ is unique up to the order of its blocks, while PP is not unique. The columns of PP are Jordan chains, and checking AP=PJAP = PJ with PP invertible verifies a proposed answer.

  • Construction. On GλG_\lambda with N=A−λIN = A - \lambda I, the spaces Ek=ker⁡N∩im⁡Nk−1E_k = \ker N \cap \operatorname{im} N^{k-1} decrease in kk. Choose a basis of the highest nonzero EkE_k and extend it downward. Each vector bb first added at level kk is the bottom of one chain of length kk, lifted through a ww with Nk−1w=bN^{k-1}w = b, and dim⁡Ek=rank⁡Nk−1−rank⁡Nk\dim E_k = \operatorname{rank} N^{k-1} - \operatorname{rank} N^k shows that the chains span.

  • Counting. With dk=dim⁡ker⁡(A−λI)kd_k = \dim \ker(A - \lambda I)^k and d0=0d_0 = 0, the number ckc_k of blocks of size at least kk and the number bk=ck−ck+1b_k = c_k - c_{k+1} of blocks of size exactly kk are

    ck=dk−dk−1,bk=2dk−dk−1−dk+1.\begin{aligned} c_k &= d_k - d_{k-1}, \\[0.4em] b_k &= 2d_k - d_{k-1} - d_{k+1}. \end{aligned}
Quantity for an eigenvalue λ\lambdaMeaning in the Jordan form
dim⁡ker⁡(A−λI)\dim \ker(A - \lambda I), the geometric multiplicitynumber of blocks for λ\lambda
dim⁡Gλ\dim G_\lambda, the algebraic multiplicitytotal size of those blocks
exponent of t−λt - \lambda in μA\mu_Asize of the largest block

For example, a nilpotent operator on a four-dimensional space with d1,d2,d3=2,3,4d_1, d_2, d_3 = 2, 3, 4 has c1,c2,c3=2,1,1c_1, c_2, c_3 = 2, 1, 1 and blocks of sizes 33 and 11. Kernel dimensions 2,4,5,52, 4, 5, 5 on a five-dimensional space give blocks of sizes 33 and 22, minimal polynomial t3t^3 and rank⁡N=3\operatorname{rank} N = 3.

Exact Jordan structure is fragile. The matrix AεA_\varepsilon with rows (1,1)(1, 1) and (ε,1)(\varepsilon, 1) is one block of size 22 at ε=0\varepsilon = 0, but for every ε>0\varepsilon > 0 it has the distinct eigenvalues 1±ε1 \pm \sqrt{\varepsilon}, with eigenvectors (1,±ε)(1, \pm\sqrt{\varepsilon}) that are nearly parallel when ε\varepsilon is small. A tiny computed singular value cannot establish an exact block without an error model, and numerical work prefers the unitary Schur form.

Lessons Constructing Jordan Form, Computing Jordan Form and Its Limitations

Numerical Linear Algebra

Exact algebra and floating-point computation answer different questions. This unit separates the sensitivity of a problem, measured by a condition number, from the error a method adds, measured by backward error. That distinction explains why a small residual need not mean an accurate answer, guides the choice between QR, the normal equations and Cholesky, and shows when power iteration finds a dominant eigenvector.

Matrix norms and condition numbers

Given vector norms on the domain and the codomain, the induced matrix norm is

∥A∥=max⁡x≠0∥Ax∥∥x∥.\lVert A \rVert = \max_{x \ne 0} \frac{\lVert Ax \rVert}{\lVert x \rVert}.

It satisfies ∥Ax∥≤∥A∥ ∥x∥\lVert Ax \rVert \le \lVert A \rVert\,\lVert x \rVert, ∥AB∥≤∥A∥ ∥B∥\lVert AB \rVert \le \lVert A \rVert\,\lVert B \rVert when the middle norms match, and ∥I∥=1\lVert I \rVert = 1 when one norm is used on both sides.

NormValueComment
∥A∥1\lVert A \rVert_1largest column sum of the ∣aij∣\lvert a_{ij} \rvertinduced by ∥x∥1\lVert x \rVert_1
∥A∥∞\lVert A \rVert_\inftylargest row sum of the ∣aij∣\lvert a_{ij} \rvertinduced by ∥x∥∞\lVert x \rVert_\infty
∥A∥2\lVert A \rVert_2σ1\sigma_1induced by the Euclidean norm, equals ∥A∗∥2\lVert A^{*} \rVert_2
∥A∥F\lVert A \rVert_Ftr⁡(A∗A)=∑iσi2\sqrt{\operatorname{tr}(A^{*}A)} = \sqrt{\sum_i \sigma_i^2}for square sizes n≥2n \ge 2 not induced by any one vector norm, since ∥In∥F=n\lVert I_n \rVert_F = \sqrt{n}

The Frobenius norm satisfies ∥AB∥F≤∥A∥2 ∥B∥F≤∥A∥F ∥B∥F\lVert AB \rVert_F \le \lVert A \rVert_2\,\lVert B \rVert_F \le \lVert A \rVert_F\,\lVert B \rVert_F. For a square invertible AA and an induced norm, the condition number is

κ(A)=∥A∥ ∥A−1∥≥1,κ2(A)=σ1σn.\begin{gathered} \kappa(A) = \lVert A \rVert\,\lVert A^{-1} \rVert \ge 1, \\[0.4em] \kappa_2(A) = \frac{\sigma_1}{\sigma_n}. \end{gathered}

If Ax=bAx = b with b≠0b \ne 0 and the data change to b+δbb + \delta b, the solution changes by δx=A−1δb\delta x = A^{-1}\delta b, and

∥δx∥∥x∥≤κ(A) ∥δb∥∥b∥.\frac{\lVert \delta x \rVert}{\lVert x \rVert} \le \kappa(A)\,\frac{\lVert \delta b \rVert}{\lVert b \rVert}.

This is a worst case, and many perturbations are amplified far less. The number κ2\kappa_2 is unchanged by unitary factors and by scaling AA, but rescaling coordinates can change it. A small determinant is no test, since εI\varepsilon I has κ2=1\kappa_2 = 1 for every ε≠0\varepsilon \ne 0. For a singular matrix the inverse problem has no finite condition number.

For example, A=diag⁡(1,10−3)A = \operatorname{diag}(1, 10^{-3}) has κ2(A)=1000\kappa_2(A) = 1000. With b=(1,0)b = (1, 0), the perturbation δb=(0,10−6)\delta b = (0, 10^{-6}) gives δx=(0,10−3)\delta x = (0, 10^{-3}), a thousand times larger in relative size, while the same perturbation along e1e_1 is not amplified. A matrix with singular values 55, 22 and 1/21/2 has ∥A∥2=5\lVert A \rVert_2 = 5, ∥A−1∥2=2\lVert A^{-1} \rVert_2 = 2, κ2=10\kappa_2 = 10 and ∥A∥F=117/2\lVert A \rVert_F = \sqrt{117}/2.

Lessons Matrix Norms and Condition Numbers

Residuals, forward error and backward error

Let AA be square and invertible, b≠0b \ne 0, xx the exact solution of Ax=bAx = b and x^\hat{x} a computed one, with all norms consistent and induced.

QuantityDefinitionWhat it measures
Residualr=b−Ax^r = b - A\hat{x}mismatch in the equation, computable
Relative forward error∥x^−x∥/∥x∥\lVert \hat{x} - x \rVert / \lVert x \rVertdistance to the true answer
Backward error, data bb only∥r∥/∥b∥\lVert r \rVert / \lVert b \rVertx^\hat{x} solves exactly Ax^=b−rA\hat{x} = b - r
Backward error, matrix only∥r∥2/∥x^∥2\lVert r \rVert_2 / \lVert \hat{x} \rVert_2 for x^≠0\hat{x} \ne 0smallest ∥δA∥2\lVert \delta A \rVert_2 with (A+δA)x^=b(A + \delta A)\hat{x} = b

The matrix-only bound is attained by the rank-one change δA=rx^∗/∥x^∥22\delta A = r\hat{x}^{*}/\lVert \hat{x} \rVert_2^2. In the model that changes only bb, subtracting the equations gives x^−x=−A−1r\hat{x} - x = -A^{-1}r, and

1κ(A) ∥r∥∥b∥≤∥x^−x∥∥x∥≤κ(A) ∥r∥∥b∥.\begin{aligned} \frac{1}{\kappa(A)}\,\frac{\lVert r \rVert}{\lVert b \rVert} &\le \frac{\lVert \hat{x} - x \rVert}{\lVert x \rVert} \\[0.4em] &\le \kappa(A)\,\frac{\lVert r \rVert}{\lVert b \rVert}. \end{aligned}
  • Backward stability. An algorithm is backward stable when its computed answer is the exact answer for data changed by a modest multiple of the unit roundoff uu. In the usual model each rounded operation gives the exact result times 1+δ1 + \delta with ∣δ∣≤u\lvert \delta \rvert \le u, ignoring overflow and underflow, and in IEEE double precision u=2−53≈1.1×10−16u = 2^{-53} \approx 1.1 \times 10^{-16}. A backward stable method cannot remove sensitivity that the problem already has.
  • Numerical rank. Exact rank counts nonzero singular values. Numerical rank counts singular values above a tolerance that must be tied to the data uncertainty or the arithmetic precision.

For example, A=diag⁡(1,10−3)A = \operatorname{diag}(1, 10^{-3}), b=(1,0)b = (1, 0) and x^=(1,1)\hat{x} = (1, 1) give r=(0,−10−3)r = (0, -10^{-3}). The relative residual is 10−310^{-3}, yet the relative forward error is 11, exactly κ2(A)\kappa_2(A) times the relative residual.

Lessons Residuals, Forward Error and Backward Error

Householder QR factorization

For a nonzero real vector vv, the Householder reflector

H=I−2 vvTvTvH = I - \frac{2\,vv^{\mathsf{T}}}{v^{\mathsf{T}}v}

is symmetric and orthogonal with H2=IH^2 = I. It reverses the component along vv and fixes everything perpendicular to vv. To send a real column aa with a nonzero tail to a multiple of e1e_1, let α=∥a∥\alpha = \lVert a \rVert and s=sign⁡(a1)s = \operatorname{sign}(a_1) with sign⁡(0)=1\operatorname{sign}(0) = 1. Then

v=a+sαe1givesHa=−sαe1.v = a + s\alpha e_1 \quad \text{gives} \quad Ha = -s\alpha e_1.

This sign adds magnitudes in the first entry and avoids cancellation. If the tail of aa is already zero, use the identity, and never divide by vTv=0v^{\mathsf{T}}v = 0. Complex data need unitary reflectors with a phase, not this real sign rule.

Applying reflectors column by column, each acting only on the rows not yet finished, gives A=QRA = QR for A∈Rm×nA \in \mathbb{R}^{m \times n} with m≥nm \ge n. In the full form Q=[Q1 Q2]Q = [Q_1 \ Q_2] is m×mm \times m orthogonal and RR is m×nm \times n with an upper triangular top block R1R_1 and zero rows below, and the reduced form is A=Q1R1A = Q_1R_1. With full column rank R1R_1 is invertible, and because orthogonal matrices keep lengths,

∥Ax−b∥2=∥R1x−Q1Tb∥2+∥Q2Tb∥2.\begin{aligned} \lVert Ax - b \rVert^2 &= \lVert R_1x - Q_1^{\mathsf{T}}b \rVert^2 \\[0.4em] &\quad + \lVert Q_2^{\mathsf{T}}b \rVert^2. \end{aligned}

So the least-squares solution solves R1x^=Q1TbR_1\hat{x} = Q_1^{\mathsf{T}}b by back substitution, and the smallest residual norm is ∥Q2Tb∥=∥b−Q1Q1Tb∥\lVert Q_2^{\mathsf{T}}b \rVert = \lVert b - Q_1Q_1^{\mathsf{T}}b \rVert. A reduced calculation must keep that last term in mind.

For example, a=(3,4,0)a = (3, 4, 0) has α=5\alpha = 5 and v=(8,4,0)v = (8, 4, 0), and the reflector sends aa to (−5,0,0)(-5, 0, 0). For AA with rows (1,0)(1, 0), (1,1)(1, 1) and (1,2)(1, 2) and b=(1,1,2)b = (1, 1, 2), the least-squares solution is x^=(5/6,1/2)\hat{x} = (5/6, 1/2) with residual (1/6,−1/3,1/6)(1/6, -1/3, 1/6), which is orthogonal to both columns.

Lessons Householder QR Factorization

QR, normal equations and Cholesky

For a real A∈Rm×nA \in \mathbb{R}^{m \times n} with full column rank, the condition number is κ2(A)=∥A∥2 ∥A+∥2=σ1/σn\kappa_2(A) = \lVert A \rVert_2\,\lVert A^{+} \rVert_2 = \sigma_1/\sigma_n, and

κ2(ATA)=κ2(A)2.\kappa_2(A^{\mathsf{T}}A) = \kappa_2(A)^2.

The normal equations ATAx^=ATbA^{\mathsf{T}}A\hat{x} = A^{\mathsf{T}}b are exactly correct, but they solve a system with the squared condition number, and forming ATAA^{\mathsf{T}}A rounds away information. For the columns (1,1,1)(1, 1, 1) and (1,1+ε,1−ε)(1, 1 + \varepsilon, 1 - \varepsilon) with ε=10−4\varepsilon = 10^{-4}, the exact Gram matrix has rows (3,3)(3, 3) and (3,3+2⋅10−8)(3, 3 + 2 \cdot 10^{-8}) and determinant 6⋅10−86 \cdot 10^{-8}, but rounding each entry to six significant digits makes it singular, although the columns are independent.

  • Cholesky solve. For a given positive definite system Ax=bAx = b, factor A=RTRA = R^{\mathsf{T}}R and solve RTy=bR^{\mathsf{T}}y = b and then Rx=yRx = y by substitution.
  • Partial pivoting. For a general square system, elimination chooses as the pivot the entry of largest modulus in the current column, on or below the diagonal, so every multiplier has modulus at most 11, and the same row swaps are applied to bb. For the matrix with rows (ε,1)(\varepsilon, 1) and (1,1)(1, 1), eliminating without a swap uses the multiplier 1/ε1/\varepsilon.
ProblemMethodReason
Least squares with full column rankHouseholder QRworks with κ2(A)\kappa_2(A), not κ2(A)2\kappa_2(A)^2
A given positive definite systemCholeskyuses the structure that is already there
A general square systemelimination with partial pivotingkeeps the multipliers at most 11 in modulus
Rank in doubt or noisy dataSVD or pivoted QR with a stated toleranceshows small singular values or nearly dependent columns

For example, with AA with rows (1,0)(1, 0), (1,1)(1, 1), (1,2)(1, 2) and b=(1,1,2)b = (1, 1, 2), the normal equations have matrix rows (3,3)(3, 3), (3,5)(3, 5) and right side (4,5)(4, 5), and give the same x^=(5/6,1/2)\hat{x} = (5/6, 1/2) as QR. In the Cholesky factorization RTRR^{\mathsf{T}}R of that Gram matrix, the lower triangular RTR^{\mathsf{T}} has rows (3,0)(\sqrt{3}, 0) and (3,2)(\sqrt{3}, \sqrt{2}).

Lessons QR, Normal Equations and Cholesky

Power iteration

Power iteration starts from a nonzero x0x_0 and repeats yk+1=Axky_{k+1} = Ax_k and xk+1=yk+1/∥yk+1∥x_{k+1} = y_{k+1}/\lVert y_{k+1} \rVert. If some yk+1=0y_{k+1} = 0, stop and report that xkx_k lies in ker⁡A\ker A.

The convergence theorem assumes that A∈Cn×nA \in \mathbb{C}^{n \times n}, a real matrix being viewed as complex, is diagonalizable with eigenbasis v1,…,vnv_1, \dots, v_n of Cn\mathbb{C}^n, where Avj=λjvjAv_j = \lambda_j v_j, and that ∣λ1∣>∣λ2∣≥⋯≥∣λn∣\lvert \lambda_1 \rvert > \lvert \lambda_2 \rvert \ge \cdots \ge \lvert \lambda_n \rvert. If x0=∑jcjvjx_0 = \sum_j c_jv_j with c1≠0c_1 \ne 0, then

Akx0λ1k=c1v1+∑j≥2cj(λjλ1)kvj→c1v1,\begin{aligned} \frac{A^kx_0}{\lambda_1^k} &= c_1v_1 + \sum_{j \ge 2} c_j \Bigl(\frac{\lambda_j}{\lambda_1}\Bigr)^k v_j \\[0.4em] &\to c_1v_1, \end{aligned}

so the line through xkx_k approaches the line through v1v_1, with error of order ∣λ2/λ1∣k\lvert \lambda_2/\lambda_1 \rvert^k. A negative or nonreal λ1\lambda_1 makes xkx_k change sign or phase while its line converges. The coefficient c1c_1 is a coordinate in the eigenbasis. It equals ⟨v1,x0⟩\langle v_1, x_0 \rangle when the eigenbasis is orthonormal, but not in general.

SituationWhat happens
c1=0c_1 = 0in exact arithmetic the iterates never approach v1v_1
∣λ1∣=∣λ2∣\lvert \lambda_1 \rvert = \lvert \lambda_2 \rvert with λ1≠λ2\lambda_1 \ne \lambda_2directions can cycle, as for diag⁡(1,−1)\operatorname{diag}(1, -1)
AA not diagonalizablethe proof does not apply, and a defective dominant eigenvalue can slow convergence to a rate like 1/k1/k
nearly dependent eigenbasislarge constants and slow visible progress

For a unit vector vv and a scalar θ\theta, the residual ∥Av−θv∥\lVert Av - \theta v \rVert checks an approximate eigenpair. If AA is normal, then

min⁡j∣λj−θ∣≤∥Av−θv∥,\min_j \lvert \lambda_j - \theta \rvert \le \lVert Av - \theta v \rVert,

but the bound fails for other matrices. The matrix with rows (0,100)(0, 100) and (0,0)(0, 0), θ=1\theta = 1 and vv proportional to (1,1/100)(1, 1/100) give a residual 1/10001<1/1001/\sqrt{10001} < 1/100, although the only eigenvalue 00 is at distance 11.

For example, diag⁡(3,1)\operatorname{diag}(3, 1) from x0=(1,1)x_0 = (1, 1) gives Akx0=(3k,1)A^kx_0 = (3^k, 1), a component ratio of 3−k3^{-k}. For diag⁡(5,4,1)\operatorname{diag}(5, 4, 1) from a start with every component nonzero the error ratio is (4/5)k(4/5)^k, and the start (0,1,1)(0, 1, 1) converges to e2e_2 instead, with ratio (1/4)k(1/4)^k. For the matrix with rows (3,1)(3, 1) and (0,1)(0, 1), the start (−1/2,1)(-1/2, 1) is an eigenvector for 11 and stays on its line, even though its dot product with the dominant eigenvector (1,0)(1, 0) is nonzero.

Lessons Power Iteration

Matrix Functions and Applications

Polynomials and convergent power series apply to square matrices as well as to numbers. Matrix powers and the matrix exponential show how Jordan chains add polynomial factors to the rates set by the eigenvalues. The unit ends with two applications, independent modes of a coupled update and regularized reconstruction from noisy data, and with problems that combine the whole course.

Matrix powers and polynomial functions

For any polynomial pp, divide by the minimal polynomial, p=qμA+rp = q\mu_A + r with deg⁡r<deg⁡μA\deg r < \deg \mu_A. Then p(A)=r(A)p(A) = r(A), and two polynomials give the same matrix exactly when μA\mu_A divides their difference. Agreement at the eigenvalues is not enough for a defective matrix. The shear JJ with rows (1,1)(1, 1) and (0,1)(0, 1) has the single eigenvalue 11, yet the polynomials tt and 11 give J≠IJ \ne I.

Polynomial functions respect similarity, p(S−1AS)=S−1p(A)Sp(S^{-1}AS) = S^{-1}p(A)S, so a diagonalizable A=SΛS−1A = S\Lambda S^{-1} has p(A)=S p(Λ) S−1p(A) = S\,p(\Lambda)\,S^{-1}, with pp applied to the diagonal entries of Λ\Lambda. A matrix function is not computed entry by entry. For a Jordan block λI+N\lambda I + N with Ns=0N^s = 0, the two parts commute and the binomial theorem gives, for every k≥0k \ge 0,

(λI+N)k=∑j=0min⁡(k, s−1)(kj)λk−jNj.(\lambda I + N)^k = \sum_{j=0}^{\min(k,\, s-1)} \binom{k}{j}\lambda^{k-j}N^j.

For λ≠0\lambda \ne 0 the entries are λk\lambda^k times polynomials in kk of degree at most s−1s - 1, so a long chain adds polynomial growth to the rate ∣λ∣k\lvert \lambda \rvert^k.

Eigenvalues λ\lambda of AAPowers AkA^k as k→∞k \to \infty
every ∣λ∣<1\lvert \lambda \rvert < 1tend to 00, possibly after transient growth
every ∣λ∣≤1\lvert \lambda \rvert \le 1, and blocks with ∣λ∣=1\lvert \lambda \rvert = 1 of size 11stay bounded
some ∣λ∣=1\lvert \lambda \rvert = 1 in a Jordan block of size at least 22unbounded
some ∣λ∣>1\lvert \lambda \rvert > 1unbounded

For example, JkJ^k has rows (1,k)(1, k) and (0,1)(0, 1), so Jke2=(k,1)J^ke_2 = (k, 1) drifts without bound although every eigenvalue is 11. The Jordan block of size 33 with eigenvalue 11 has J4=I+4N+6N2J^4 = I + 4N + 6N^2, and A=2I+NA = 2I + N with N2=0N^2 = 0 has A4=16I+32N=32A−48IA^4 = 16I + 32N = 32A - 48I. The matrix with rows (1/2,10)(1/2, 10) and (0,1/2)(0, 1/2) has powers with upper right entry 20k/2k20k/2^k, which tend to 00, but it sends e2e_2 to a vector of length greater than 1010.

Lessons Matrix Powers and Polynomial Functions

The matrix exponential

For A∈Fn×nA \in \mathbb{F}^{n \times n} and a scalar tt, the matrix exponential is the series

etA=∑k=0∞tkAkk!,e^{tA} = \sum_{k=0}^{\infty} \frac{t^k A^k}{k!},

which converges absolutely because ∥Ak∥≤∥A∥k\lVert A^k \rVert \le \lVert A \rVert^k in an induced norm, and ∥etA∥≤e∣t∣ ∥A∥\lVert e^{tA} \rVert \le e^{\lvert t \rvert\,\lVert A \rVert}.

PropertyCondition
e0=Ie^{0} = Ialways
eS−1AS=S−1eASe^{S^{-1}AS} = S^{-1}e^{A}SSS invertible
eA+B=eAeBe^{A + B} = e^{A}e^{B}AB=BAAB = BA
(eA)−1=e−A(e^{A})^{-1} = e^{-A}always
e(s+t)A=esAetAe^{(s + t)A} = e^{sA}e^{tA}always
ddtetA=AetA=etAA\dfrac{d}{dt}e^{tA} = Ae^{tA} = e^{tA}Areal tt, using calculus

Without commuting matrices the product rule can fail. For AA with rows (0,1),(0,0)(0, 1), (0, 0) and CC with rows (0,0),(1,0)(0, 0), (1, 0), the product eAeC=(I+A)(I+C)e^{A}e^{C} = (I + A)(I + C) has rows (2,1)(2, 1) and (1,1)(1, 1), while eA+Ce^{A + C} has both diagonal entries equal to cosh⁡1\cosh 1. For a Jordan block with Ns=0N^s = 0,

et(λI+N)=eλt∑j=0s−1tjj!Nj,e^{t(\lambda I + N)} = e^{\lambda t} \sum_{j=0}^{s-1} \frac{t^j}{j!}N^j,

and a diagonalizable A=SΛS−1A = S\Lambda S^{-1} has etA=S diag⁡(etλ1,…,etλn) S−1e^{tA} = S\,\operatorname{diag}(e^{t\lambda_1}, \dots, e^{t\lambda_n})\,S^{-1}. For real tt, ∣etλ∣=etRe⁡λ\lvert e^{t\lambda} \rvert = e^{t\operatorname{Re}\lambda}. If every eigenvalue has negative real part, then etA→0e^{tA} \to 0 as t→∞t \to \infty, but ∥etAx∥\lVert e^{tA}x \rVert can grow first. With calculus, x(t)=etAx0x(t) = e^{tA}x_0 is the solution of x′=Axx' = Ax with x(0)=x0x(0) = x_0.

For example, BB with rows (−1,4)(-1, 4) and (0,−1)(0, -1) has

etB=e−t[14t01].e^{tB} = e^{-t}\begin{bmatrix} 1 & 4t \\ 0 & 1 \end{bmatrix}.

At t=1/2t = 1/2 it sends e2e_2 to e−1/2(2,1)e^{-1/2}(2, 1), of length 5/e≈1.36>1\sqrt{5/e} \approx 1.36 > 1, although every trajectory tends to 00.

Lessons The Matrix Exponential

Coupled updates and eigenmodes

A linear update xk+1=Bxkx_{k+1} = Bx_k has xk=Bkx0x_k = B^kx_0. If BB is real symmetric with orthonormal eigenvectors qjq_j and eigenvalues λj\lambda_j, each mode evolves on its own,

xk=∑jλjk (qjTx0) qj.x_k = \sum_j \lambda_j^k\,(q_j^{\mathsf{T}}x_0)\,q_j.

For a diagonalizable B=SΛS−1B = S\Lambda S^{-1} that is not symmetric, the coefficients come from S−1x0S^{-1}x_0, not from dot products, and the modes need not be orthogonal. A positive definite form xTHxx^{\mathsf{T}}Hx can measure the size of the state. In an oblique basis with matrix PP the operator changes to P−1BPP^{-1}BP by similarity while the form changes to PTHPP^{\mathsf{T}}HP by congruence, and the two rules agree only for orthogonal PP. Results that used symmetry must be checked again when the coupling becomes nonsymmetric or defective.

For example, two signals that keep half their value and receive a quarter of the other follow BB with rows (1/2,1/4)(1/2, 1/4) and (1/4,1/2)(1/4, 1/2). The modes q±=(1,±1)/2q_{\pm} = (1, \pm 1)/\sqrt{2} decay by the factors 3/43/4 and 1/41/4, and from x0=(1,0)x_0 = (1, 0)

xk=12[(3/4)k+(1/4)k(3/4)k−(1/4)k].x_k = \frac{1}{2}\begin{bmatrix} (3/4)^k + (1/4)^k \\ (3/4)^k - (1/4)^k \end{bmatrix}.

Here B=I−H/4B = I - H/4 with HH having rows (2,−1)(2, -1) and (−1,2)(-1, 2), and xTHx=(x1−x2)2+x12+x22x^{\mathsf{T}}Hx = (x_1 - x_2)^2 + x_1^2 + x_2^2 is positive for x≠0x \ne 0, with HH having eigenvalues 11 and 33. The update I−hHI - hH has factors 1−h1 - h and 1−3h1 - 3h, so it decays for every start exactly when 0<h<2/30 < h < 2/3. Changing the upper right coupling to 1/21/2 gives eigenvalues 1/2±1/(22)1/2 \pm 1/(2\sqrt{2}) with eigenvectors (±2,1)(\pm\sqrt{2}, 1), which are not orthogonal.

Lessons Coupled Systems and Eigenmodes

Inverse problems and regularization

An inverse problem asks for xx from noisy data b≈Axb \approx Ax. With a compact singular value decomposition A=∑i=1rσiuivi∗A = \sum_{i=1}^{r} \sigma_i u_iv_i^{*} and βi=ui∗b\beta_i = u_i^{*}b, the usual reconstructions have the form

x=∑i=1rφi βiσi vi,x = \sum_{i=1}^{r} \varphi_i\,\frac{\beta_i}{\sigma_i}\,v_i,

where the filter factors φi\varphi_i decide how much of each direction is kept. Small σi\sigma_i amplify noise, and components in ker⁡A\ker A are never measured.

MethodFilter factor φi\varphi_iEffect
Pseudoinverse11smallest residual, divides noise by small σi\sigma_i
Truncated SVD at threshold τ\tau11 when σi>τ\sigma_i > \tau, 00 otherwisedrops the weak directions
Tikhonov with parameter λ>0\lambda > 0σi2σi2+λ\dfrac{\sigma_i^2}{\sigma_i^2 + \lambda}shrinks the weak directions gradually

Tikhonov regularization minimizes ∥Ax−b∥2+λ∥x∥2\lVert Ax - b \rVert^2 + \lambda\lVert x \rVert^2, where λ>0\lambda > 0 is a penalty parameter, not an eigenvalue. The matrix A∗A+λIA^{*}A + \lambda I is positive definite even when AA is rank deficient, so the minimizer is unique,

(A∗A+λI) xλ=A∗b,xλ=∑i=1rσi βiσi2+λ vi.\begin{gathered} (A^{*}A + \lambda I)\,x_\lambda = A^{*}b, \\[0.4em] x_\lambda = \sum_{i=1}^{r} \frac{\sigma_i\,\beta_i}{\sigma_i^2 + \lambda}\,v_i. \end{gathered}

Increasing λ\lambda never increases the solution norm and never decreases the residual. Choosing τ\tau or λ\lambda needs a noise model, a known tolerance or outside validation, and changing units changes what the penalty measures. All three methods set the ker⁡A\ker A component to zero by convention, which is not a measurement of it.

For example, A=diag⁡(1,1/10,0)A = \operatorname{diag}(1, 1/10, 0) and b=(1,1/50,3/100)b = (1, 1/50, 3/100) give the pseudoinverse solution (1,1/5,0)(1, 1/5, 0) with residual (0,0,3/100)(0, 0, 3/100), the truncation at τ=3/20\tau = 3/20 gives (1,0,0)(1, 0, 0) with residual norm 13/100\sqrt{13}/100, and Tikhonov with λ=1/100\lambda = 1/100 gives (100/101,1/10,0)(100/101, 1/10, 0). A change in the second data entry reaches the second coefficient with gain 1010, 00 and 55 in the three methods.

Lessons Inverse Problems and Regularization

Results that combine the course

  • Dual of a quotient. For a subspace WW of a finite-dimensional VV, the pullback ℓ↦ℓ∘q\ell \mapsto \ell \circ q through the quotient map is an isomorphism from (V/W)∨(V/W)^{\vee} onto W0W^{0}. A functional descends to V/WV/W exactly when it vanishes on WW. For W=span⁡(1,1,0)W = \operatorname{span}(1, 1, 0) in R3\mathbb{R}^3, the descending functionals are the combinations of x−yx - y and zz, such as x−y+2zx - y + 2z.
  • Normality test. Compare AA∗AA^{*} with A∗AA^{*}A before using the spectral theorem. The matrix CC with rows (1,i)(1, i) and (0,1)(0, 1) has CC∗CC^{*} with rows (2,i),(−i,1)(2, i), (-i, 1) and C∗CC^{*}C with rows (1,i),(−i,2)(1, i), (-i, 2), so it is not normal, and its one-dimensional eigenspace shows it is not diagonalizable either.
  • Verifying a Jordan basis. A proposed PP and JJ are right when PP is invertible and AP=PJAP = PJ, checked column by column. For AA with rows (2,1,1)(2, 1, 1), (0,2,0)(0, 2, 0), (0,0,3)(0, 0, 3), the columns e1,e2,(1,0,1)e_1, e_2, (1, 0, 1) give J=diag⁡(J2(2),3)J = \operatorname{diag}(J_2(2), 3), and (A−2I)2(A - 2I)^2 projects onto the eigenvector line for 33 along the generalized eigenspace for 22, obliquely, since it is not symmetric.
  • Orthogonal projectors. An idempotent PP is an orthogonal projector exactly when P∗=PP^{*} = P for the chosen inner product. The matrix with rows (1,2)(1, 2) and (0,0)(0, 0) is idempotent, but its kernel, spanned by (−2,1)(-2, 1), is not perpendicular to its image.
  • Forward differences. The operator (Δp)(t)=p(t+1)−p(t)(\Delta p)(t) = p(t + 1) - p(t) on Pd\mathcal{P}_d lowers the degree of every nonconstant polynomial by one, so Δd+1=0\Delta^{d+1} = 0 while Δd≠0\Delta^d \ne 0, and its Jordan form is a single nilpotent block of size d+1d + 1. On P2\mathcal{P}_2 a chain is 1,t,(t2−t)/21, t, (t^2 - t)/2.
  • Method choice under noise. A reconstruction is justified by the noise level, the residual, the solution norm and the sensitivity of each direction together. No method recovers a component in ker⁡A\ker A.

Lessons Proofs and Problem Solving