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 , which is or , and all vector spaces are finite-dimensional unless a statement says otherwise. Vectors such as are written in plain italic, coordinates are columns, and a matrix represents a linear map from an -dimensional space to an -dimensional one. The linear maps from to form the space . Coordinates change by , and the matrix has as its th column the coordinates of .
The transpose is and the conjugate transpose is . Inner products are conjugate-linear in the first argument and linear in the second, with the standard inner product on and . The algebraic dual is and the dual map is , while always means the adjoint, which needs an inner product. The space holds the polynomials of degree at most , including the zero polynomial. The characteristic polynomial is and the minimal polynomial is monic. Singular values are listed in decreasing order with , while the eigenvalues of a Hermitian matrix are listed in increasing order 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 is proved by an argument about an arbitrary member of and disproved by one counterexample. An existence claim is proved by producing one member with the property, or by an argument that one must exist. An implication is false only when holds and fails. Negation switches each quantifier.
| Statement | Negation |
|---|---|
| and | |
| some has and |
- Converse and contrapositive. The converse of is , a different statement. The contrapositive is equivalent to the original. A statement needs a proof of both directions.
- Induction. If the base case holds and implies for every integer , then holds for every integer . Checking several values of proves nothing about the others.
For example, if , then for every integer . The base case is , and the step is
Lessons Proof Strategies in Linear Algebra
Fields, vector spaces and subspaces
A complex number with has conjugate and modulus . Conjugation respects sums and products, , and for
A vector space over allows scalars only from . Its axioms ask that addition is commutative and associative with a zero vector and negatives, that and , and that both distributive laws hold. Everything proved from the axioms alone, such as and , holds in every vector space.
- Subspace test. A subset of is a subspace exactly when it is nonempty and for all and . The line in fails because it misses , and the union of the two axes fails because lies on neither.
- Real and complex dimension. If is a basis of a complex space , then is a basis of as a real space, so . A statement proved for real spaces does not transfer to complex scalars without a check.
- Polynomials. The space has basis , so . The polynomials of degree exactly are not a subspace, since the zero polynomial is missing.
For example, the vector has two complex coordinates in the standard basis, and over it has the four coordinates in the basis .
Lessons Real and Complex Vector Spaces
Bases and dimension
The exchange lemma says that if are linearly independent in and span , then . Its key step replaces some by a vector with , which keeps the span. When are inserted one at a time, independence always leaves such a nonzero coefficient on a remaining . Applied in both directions to two bases it shows that all bases of have the same length, the dimension , with .
| A list of vectors in , with | Conclusion |
|---|---|
| more than vectors | linearly dependent |
| fewer than vectors | does not span |
| exactly independent vectors | a basis |
| exactly spanning vectors | a 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 is a subspace of , then , with equality only when .
For example, and extend to a basis of by , and . The three vectors span only a plane, so three vectors in 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 is a basis of and are arbitrary, exactly one linear has for every . With pointwise addition and scaling, and , the set is a vector space over . For bases of and of , the matrix-unit maps with form a basis, and the coordinates of in it are the entries of . Hence
A map has rank one exactly when for a nonzero and a nonzero linear . In matrices that is a nonzero column times a nonzero row. The kernel and the image are subspaces, and when is finite-dimensional the rank-nullity theorem says
| Property of | Test | Consequence in finite dimensions |
|---|---|---|
| Injective | ||
| Surjective | ||
| Isomorphism | injective and surjective | , and is linear |
When is finite, an injective or a surjective linear map is already an isomorphism. With unequal dimensions this fails, as the inclusion of in shows. For composable maps, .
For example, from to sends to . Its kernel is spanned by and it is onto, so . All polynomials with measurements are with .
Lessons The Vector Space of Linear Maps, Rank-Nullity and Isomorphisms
Sums, direct sums and projections
For subspaces and of , the sum and the intersection are subspaces, and for finite-dimensional and
The union is a subspace only when one of them contains the other. The product of two spaces, with componentwise operations, has dimension .
- Direct sum. The space is the direct sum when and , which holds exactly when every is with and in only one way. Then .
- Complements. In a finite-dimensional space every subspace has a complement with , found by extending a basis of . It is unique only when or .
- Projection along a complement. If , the map is linear with , and , and projects onto along .
- Idempotents. Conversely, every linear with gives and is the projection onto along . The splitting alone does not force , as shows.
A projection needs no inner product. With the dot product on , 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 along is , with matrix
and it splits as . 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 be a subspace of . Vectors and of are equivalent modulo when . The coset of is
and exactly when . Two cosets are equal or disjoint, and a coset is a subspace only when it is itself. The quotient space is the set of cosets with the operations
which are well defined because changing representatives changes the results by vectors of . Its zero vector is the coset .
- Well-defined rules. A rule defines a function on exactly when implies . For a linear this means .
- Quotient basis. If is a basis of extended to a basis of , then the cosets form a basis of .
When is finite-dimensional, counting that basis gives
For example, in with , each coset is a vertical line . The function is constant on each line and descends to the quotient, while is not and does not.
Lessons Quotient Spaces
The first isomorphism theorem
Let be linear over . The fiber over is , which is empty when and equals the coset when . The notation does not assume an inverse.
- Quotient map. For a subspace of , the map from to is linear and onto, with kernel .
- Induced map. The rule defines a map on exactly when . It is then linear, satisfies , and is the only map with that property.
The first isomorphism theorem says that the induced map
is a linear bijection, so . 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, on has , and the coset of corresponds to . The cosets of and form a quotient basis. Evaluation maps onto , so , and for the kernel consists of the multiples with .
Lessons The First Isomorphism Theorem
Dual spaces and dual bases
A linear functional on is a linear map . Over scalars come out with no conjugate. The dual space is
with pointwise operations. Many books write for the dual space. This course keeps so that the dual is not confused with the adjoint . On every functional is for exactly one row , whose entries are , with no conjugation.
- Dual basis. For a basis of , the coordinate functionals with form a basis of , and every functional expands as . So . On , if the are the columns of , the are the rows of .
- Evaluation functionals. On with basis and dual basis , evaluation at a fixed is .
- Not functionals. On a complex space, a conjugate-linear map , with , is not in unless it is zero. On a nonzero space the norm is not linear either.
The equality does not single out an isomorphism between the two spaces. Matching with depends on the basis. For example, the basis of has dual basis and , and the functional equals . The vector has coordinates in this basis, and both descriptions give .
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 pulls measurements back,
It is linear over and over , it uses no inner product, it reverses compositions as , and it satisfies for an isomorphism . If and , are the dual bases, then
the plain transpose even over . In rows, a functional with evaluation row on pulls back to the row on .
| Object, with | Coordinates in | Coordinates in |
|---|---|---|
| Vector | column | |
| Covector as an evaluation row | row with entries | |
| Covector in the dual basis | column |
The reading does not depend on the basis, which is why covector coordinates use the inverse. Multiplying both the coefficient column of and the coordinate column of by would keep every reading only when .
For example, with having rows and and the output sensor , the pullback is . At both routes give , since . For , standard and , the row becomes , the coefficients of in the dual basis above.
Lessons Dual Maps and Covectors
Annihilators and the four subspaces
The annihilator of a subspace of is
a subspace of , not of . For finite-dimensional , . For a linear ,
where the first identity holds for all vector spaces and the second follows by counting dimensions when and are finite-dimensional. Hence , and is injective exactly when is surjective. For a matrix of rank , the dual map sends an output row to , which gives the row space as and the left nullspace as .
| Subspace | Lives in | Dimension |
|---|---|---|
| Column space | ||
| Nullspace | ||
| Row space | rows of length | |
| Left nullspace | rows of length |
An annihilator needs no extra structure, while an orthogonal complement lies in and depends on an inner product. In a finite-dimensional inner product space, the functional lies in exactly when , and the correspondence is conjugate-linear over . For example, with rows and has kernel spanned by . The rows that vanish on it are , exactly the row space, and since maps onto 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 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 of vectors is multilinear if it is linear in each slot with the others fixed, and alternating if it is 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 of is , where counts the pairs with , and every exchange of two entries flips it. For ,
The uniqueness theorem says that every function of the columns that is multilinear and alternating satisfies
so is the only one with . Applying it to gives for all square and of the same size, singular or not. Also , so every column rule holds for rows, and is invertible exactly when , with .
| Operation on the columns or the rows | Effect on |
|---|---|
| swap two of them | changes the sign |
| multiply one by | multiplies it by |
| add a multiple of one to another | none |
| make them linearly dependent | gives |
Lessons Determinants and Multilinear Forms
Trace, similarity and the characteristic polynomial
The trace is the sum of the diagonal entries. For and ,
even when the two products have different sizes. With three factors only the cyclic rotations are allowed. Matrices are similar when for an invertible . Similar matrices are exactly the matrices of one operator in two bases. If , then .
The characteristic polynomial of is monic of degree . For it expands as
and for it is . A scalar is an eigenvalue of exactly when . 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 and share trace , determinant and characteristic polynomial , but the shear is not similar to , since for every invertible .
Over a determinant is not a signed volume. Multiplication by , viewed as a real map on , scales real volume by . For example, and with rows and give
with the same determinant , trace and characteristic polynomial as .
Lessons Trace, Determinant and Similarity
Bilinear forms and congruence
A bilinear form on a real space is a function that is linear in each argument. In an ordered basis its matrix has entries , and with and
The form is symmetric exactly when is, and nondegenerate exactly when is invertible. If the columns of hold the new basis vectors in old coordinates, then
which is the congruence rule. With it reads .
| Object, new basis given by the columns of | New matrix | Unchanged |
|---|---|---|
| Operator | , similarity | eigenvalues, trace, determinant, |
| Bilinear form on a real space | , congruence | symmetry, rank, the sign of , the inertia of a symmetric form |
| Hermitian form on a complex space | , conjugate congruence | the Hermitian property, rank, inertia |
-
Quadratic form. The function satisfies .
-
Polarization. For a symmetric , the polarization identity
recovers from , so symmetric forms and quadratic forms correspond one to one.
-
Symmetric and skew parts. Every real square is with symmetric and skew-symmetric. Since , the quadratic form sees only .
For example, and the new basis give and . The eigenvalues changed, which similarity never does. The vector has coordinates , and both matrices give . The matrix with rows and is nonzero but has 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 over assigns a scalar to each pair so that, for all vectors and scalars,
These are linearity in the second argument, conjugate symmetry and positive definiteness. They force , so the first argument is conjugate-linear. The standard inner product on is .
- Norm. The norm has , and . Orthogonal vectors, with , 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 . A complex inner product is exactly a positive definite Hermitian sesquilinear form. On these forms are , Hermitian exactly when .
- Gram matrix. For a basis the Gram matrix has entries , and with coordinate columns and we have . It is Hermitian and positive definite, it equals exactly for an orthonormal basis, and a change of basis replaces it by , where the columns of hold the new basis vectors in old coordinates.
For example, has although , while . Scalars leave the first slot conjugated, so and .
Lessons Complex Inner Products
Cauchy-Schwarz and orthonormal expansions
For all in an inner product space over , in any dimension,
with equality exactly when and are linearly dependent. For the proof expands with , which makes orthogonal to . It gives the triangle inequality , with equality exactly when one vector is a nonnegative real multiple of the other.
The complex Gram-Schmidt process turns linearly independent into an orthonormal list with the same spans for every . Set and, for ,
The basis vector sits in the first slot of each coefficient. With dependent inputs, the first in the span of the earlier ones gives , or is itself zero when , and that zero signals the dependence. For an orthonormal basis of and every ,
The second line is Parseval’s identity. For a finite orthonormal list that need not be a basis, Bessel’s inequality says
with equality exactly when lies in the span of the list.
For example, and have and , of modulus . The remainder is orthogonal to .
Lessons Cauchy-Schwarz and Orthonormal Expansions
Riesz representation and adjoints
The Riesz representation theorem says that if is a finite-dimensional inner product space over , then for every there is exactly one with
In an orthonormal basis the vector is . On with the inner product for a Hermitian positive definite , the functional is represented by , which is for the standard inner product. The vector depends on the inner product, and the identification is conjugate-linear over .
For finite-dimensional inner product spaces and over and a linear , the adjoint is the unique linear map with
It obeys , , and . In orthonormal bases its matrix is . In bases with Gram matrices and , where ,
The orthogonal complement has dimension , and in finite dimensions. For linear ,
so , and .
| Dual map | Adjoint | |
|---|---|---|
| Acts | from to | from to |
| Needs | nothing beyond | inner products on and |
| Matrix | in dual bases | in orthonormal bases, in general |
For example, with rows and has with rows and . With and , , and both and equal . With the domain Gram matrix , the Euclidean target and , the adjoint has coordinate column , not .
Lessons Riesz Representation in Finite Dimensions, Adjoints and Gram Matrices
Orthogonal projection and best approximation
Let be a subspace of a finite-dimensional inner product space over . Then , and the orthogonal projection is the unique vector with and . For an orthonormal basis of ,
For a basis that is not orthonormal, where , with and .
The projection theorem says that is the only best approximation to from , and that is that best approximation exactly when . For every ,
and the squared distance from to is . A linear is an orthogonal projector exactly when and . If the orthonormal columns of span in orthonormal coordinates, its matrix is .
On , distinct real nodes with give the inner product , and the best polynomial approximation from a subspace depends on that choice. For example, with nodes the projection of onto the constants is , and the residual has samples , orthogonal to the constant samples .
Lessons Orthogonal Projection and Best Approximation
Weighted least squares
A weight matrix is Hermitian and positive definite, and a diagonal weight needs positive diagonal entries. It defines the inner product with norm , and every inner product on has this form. For and , a vector minimizes exactly when it solves the weighted normal equations
which say that the residual is -orthogonal to every column of . They always have a solution.
- Fitted vector. The fitted vector is the same for every minimizer. It is the projection of onto the column space in the inner product of .
- Coefficients. The minimizers form the set . They are unique exactly when the columns of are independent, and then .
- Weights. Multiplying by a constant does not change the fit. A singular positive semidefinite or an indefinite does not define a distance.
With positive weights and distinct real nodes, is an inner product on exactly when , and projecting in it is weighted least squares on the sample values. For example, approximating by gives with , but with , where the weighted squared error is .
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 is invariant under an operator if , which can be checked on a spanning list of . The subspaces , , , and every eigenspace are invariant, and a line is invariant exactly when is an eigenvector. The restriction is then an operator on . Over every nonzero finite-dimensional invariant subspace contains an eigenvector, which can fail over .
In a basis whose first vectors span , the subspace is invariant exactly when
with the matrix of . Then . The matrix is block diagonal, with , exactly when the complement spanned by the is invariant too, and an invariant complement need not exist.
- Adjoints. In an inner product space, if is invariant under , then is invariant under .
- Reducing subspaces. The subspace reduces when both and are invariant under , which holds exactly when is invariant under both and . Then is block diagonal in an orthonormal basis adapted to . For a self-adjoint every invariant subspace reduces it.
For example, with rows , and leaves invariant, but not its orthogonal complement, since . The invariant subspaces and give a block diagonal description with blocks of sizes and .
Lessons Invariant Subspaces and Block Matrices
Schur decomposition and real Schur form
A matrix is unitary if , which says that its columns are orthonormal. Then , , and a real unitary matrix is orthogonal. Schur’s theorem says that every has a factorization
with unitary and upper triangular. The diagonal of lists the eigenvalues of , each as often as its algebraic multiplicity, in any chosen order, and the factorization is not unique. The first columns of span an invariant subspace for each , 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 to the orthogonal complement of that eigenvector. If is real with only real eigenvalues, and can be taken real.
- Multiplicities. The algebraic multiplicity of is its exponent in , the geometric multiplicity is , and .
- Trace and determinant. With the eigenvalues listed with algebraic multiplicity, and , with no diagonalizability needed.
The real Schur theorem says that every has for a real orthogonal , where is real and block upper triangular with diagonal blocks of size , holding the real eigenvalues, and of size , each holding a pair of nonreal conjugate eigenvalues. The blocks come from real invariant planes. If with , and a nonzero for real and , then and are independent and
A real matrix with nonreal eigenvalues has the rotation-scaling form with rows and exactly when it is normal. So a general real Schur block is not a rotation followed by a uniform stretch. For example, with rows and has eigenvalues , but has rows and , which differs from , so is not normal.
Lessons Schur Decomposition, Real Schur Form
Normal operators and the spectral theorem
An operator on a finite-dimensional inner product space is normal if , and a square matrix is normal if . 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 and a linear , the space has an orthonormal basis of eigenvectors of if and only if is normal. In matrices, is normal if and only if
with unitary and 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 is not unique, since each column can be multiplied by a scalar of modulus and a repeated eigenvalue allows any orthonormal basis of its eigenspace.
| Class of | Definition | Eigenvalues | Orthonormal eigenbasis |
|---|---|---|---|
| Normal | any complex numbers | yes, over | |
| Hermitian | real | yes | |
| Unitary | modulus | yes | |
| Real symmetric | , real entries | real | yes, with real orthogonal |
| Not normal | any | no |
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 with real orthogonal and real. A real normal matrix that is not symmetric always has nonreal eigenvalues, since real eigenvalues alone would make it real symmetric, and over it keeps rotation-scaling blocks.
For example, with rows and is Hermitian with eigenvalues and and orthonormal eigenvectors and . The shear with rows and is not normal, since has rows and while has rows and .
Lessons The Spectral Theorem for Normal Operators
Simultaneous diagonalization
Operators and on the same space commute if . 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 and , then , so maps every eigenspace of into itself.
The theorem on common orthonormal eigenbases says that for a finite family of operators on a finite-dimensional complex inner product space, there is an orthonormal basis of common eigenvectors if and only if every is normal and the commute in pairs. In matrices, one unitary then makes every diagonal. The proof splits the space into the orthogonal eigenspaces of , 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, with rows and commutes with , and the two share the eigenbasis , but is not normal, so no orthonormal common eigenbasis exists. A matrix that is not diagonalizable commutes with and has no eigenbasis at all.
Lessons Simultaneous Diagonalization
Quadratic Forms
A Hermitian matrix defines a quadratic form 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 be Hermitian, which for means real symmetric. Under an invertible change of coordinates , the form becomes , so the new matrix is , conjugate congruence, which over is the congruence . Congruence is not similarity and can change the eigenvalues.
- Diagonalization by congruence. There is an invertible with . One choice is , from with the positive eigenvalues listed first, then the negative ones, then the zeros, and a diagonal with entries for and otherwise.
- Inertia. The triple is the inertia of . The numbers and count the positive and the negative eigenvalues with multiplicity, and .
- Positive subspaces. The number is the largest dimension of a subspace on which for every nonzero , and is the same for negative values.
Sylvester’s law of inertia says that if are both real symmetric or both complex Hermitian, then for some invertible if and only if and have the same inertia. Over this is conjugate congruence, since with a complex need not even be Hermitian. Every diagonal matrix with invertible has positive, negative and zero diagonal entries.
| Form | Condition | Inertia |
|---|---|---|
| Positive definite | for | |
| Positive semidefinite | for all | |
| Negative definite | for | |
| Negative semidefinite | for all | |
| Indefinite | takes both signs | and |
For example, the real form has inertia , and its matrix has kernel spanned by . The invertible matrix with rows and has inertia , and its form still vanishes on nonzero vectors such as .
Lessons Quadratic Forms and Sylvester’s Law of Inertia
Positive definite matrices and Cholesky factorization
A matrix is positive definite if it is Hermitian and for every nonzero , and positive semidefinite if it is Hermitian and for every . Over symmetry is part of the definition, because sees only the symmetric part.
| For a Hermitian , positive definite is equivalent to | Positive semidefinite is equivalent to |
|---|---|
| every eigenvalue positive | every eigenvalue nonnegative |
| inertia | inertia with |
| for an invertible | for some with columns |
| every leading principal minor positive | not decided by the leading minors |
Completing the square splits off the first coordinate. If is Hermitian with first entry , first column , first row and lower right block , then for
The matrix is the Schur complement of , and . If is positive definite, then and is positive definite.
The Cholesky factorization theorem says that every positive definite has exactly one factorization
with upper triangular and every diagonal entry real and positive. Over it reads , and many books write with . Sylvester’s criterion says that a Hermitian is positive definite if and only if every leading principal minor , the determinant of the upper left block, is positive, and then with . The Hermitian hypothesis matters, since the matrix with rows and has but gives at . Nonnegative leading minors do not prove semidefiniteness, as shows.
For example, with rows and has , leading minors and , and
Lessons Positive Definite Matrices and Cholesky Factorization
Rayleigh quotients and the min-max principle
For a Hermitian and nonzero , the Rayleigh quotient is
It is real, and for every nonzero scalar . With eigenvalues and an orthonormal eigenbasis in the same order, it is a weighted average of the eigenvalues,
with and . Hence , with exactly when and exactly when . When an extreme eigenvalue is repeated, its extremizers are not one line but all nonzero vectors of an eigenspace of dimension at least . The Hermitian hypothesis is needed, since for a real nonsymmetric matrix the quotient sees only the symmetric part.
Two subspaces of an -dimensional space satisfy , so they share a nonzero vector when their dimensions add up to more than . This dimension count proves the min-max principle, also called the Courant-Fischer theorem. For a Hermitian and ,
where and run over subspaces of . The subspaces and attain the two values, but attaining subspaces are not unique in general.
For example, has at . On the largest value is , and every two-dimensional subspace contains a nonzero vector of , where the quotient is at least , so .
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 have rank , and let . The Gram operator is Hermitian and positive semidefinite, since , with and . The singular values are the nonnegative square roots of the largest eigenvalues of , and exactly of them are positive.
The theorem says that every has a singular value decomposition
with and unitary, which means real orthogonal when , and zero except for , . To build it, take an orthonormal eigenbasis of ordered by decreasing eigenvalue, so that for , set for , and complete the to an orthonormal basis of . Then and for , and for .
| Form | Left factor | Middle factor | Right factor |
|---|---|---|---|
| Full | , size | , size | , size |
| Thin | , size | , size | , size |
| Compact | , size | , size | , size |
Here and hold the first columns and , and in each form is the left factor times the middle factor times the adjoint of the right factor. The compact form is , and the thin form may still contain zero singular values. The singular vectors give orthonormal bases of the four subspaces.
| Subspace | Orthonormal basis |
|---|---|
The singular values are determined by , but the singular vectors are not. For a positive singular value repeated times, the matching columns of and of can be replaced by and for any unitary , and for this multiplies and by the same scalar of modulus . The remaining columns can be any orthonormal bases of and . Singular values are not the moduli of eigenvalues in general. They agree for normal matrices, but the shear with rows and has eigenvalues and singular values .
For example, with rows , and has singular values and , with , with columns , and with rows , and . The left null direction has no positive singular value.
Lessons The Singular Value Decomposition
The pseudoinverse and the four projectors
For of rank with , let hold for and zeros elsewhere. The pseudoinverse is
and it does not depend on which singular value decomposition is used.
| Case | Pseudoinverse |
|---|---|
| square and invertible | |
| full column rank, | |
| full row rank, | |
| the zero matrix |
A vector minimizes exactly when , or equivalently when . Least-squares solutions always exist and form the set . Among them is the only one of smallest norm, because it lies in and for .
| Projector | Size | Projects orthogonally onto |
|---|---|---|
| , the column space | ||
| , where the residuals lie | ||
| , spanned by the conjugated rows of | ||
The Moore-Penrose theorem says that is the only matrix with
The first condition makes and idempotent, the last two make them Hermitian and hence orthogonal projectors, and all four are needed. A matrix with only is a generalized inverse. It solves every consistent system, but it is unique only when is square and invertible. For , every with rows and satisfies the first two conditions, but is not symmetric when .
For example, with rows , and has with rows and . For the least-squares solutions are , the minimum-norm one is , the fitted vector is and the residual is .
Lessons The Pseudoinverse and Four Projectors, The Moore-Penrose Conditions
Best low-rank approximation
Two matrix norms measure an error with singular values . Both are unchanged by unitary factors on either side.
The spectral norm is the error in the worst input direction and the squared Frobenius norm is the total squared error, with . The truncated singular value decomposition keeps the first terms,
The Eckart-Young theorem says that for every with ,
with , and that attains both bounds. For the matrices of rank at most do not form a subspace, so the projection theorem does not apply directly.
- Frobenius uniqueness. For , the truncation is the only best approximation in the Frobenius norm when , and there are infinitely many when .
- Spectral nonuniqueness. When and the spectral norm always has infinitely many best approximations, since changing each with by at most keeps the error .
For example, and give spectral error for every with , while the Frobenius error is smallest only at . For every rank-one orthogonal projector is a best approximation in both norms.
Lessons The Eckart-Young Theorem
Principal component analysis
Let be observations with mean . The centered data matrix has rows , so rows are observations, columns are measured coordinates and every column sums to . The sample covariance matrix is
which is symmetric and positive semidefinite. The divisor is a convention, and dividing by changes every variance by the same factor and none of the directions. Along a unit vector the scores have sample variance , a Rayleigh quotient of .
- Principal directions. With a full real singular value decomposition , the right singular vectors are the principal directions. For the variances along them are and the score vectors are .
- Reconstruction. For a subspace with orthonormal basis given by the columns of and projector , the identity shows that maximizing retained variance minimizes reconstruction error. The top directions give the smallest error, , by the Eckart-Young theorem.
- Choices. Without centering, , and the leading direction can describe the mean instead of the spread. Changing units changes the inner product and can change the answer. Each can be replaced by , and when singular values tie, any orthonormal basis of the tied eigenspace of serves equally well.
For example, the centered rows , , and give , variances and , and the horizontal line as the best one-dimensional fit, with squared error . Shifting every point up by leaves the centered answer unchanged but makes the uncentered Gram matrix , 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 has exactly one Hermitian positive semidefinite with . If , it is
and . 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, with triangular, and is usually not Hermitian.
The absolute value of is . It satisfies for every , so . The polar decomposition theorem says that exactly one satisfies
In a compact singular value decomposition, and .
| Factor | Formula | Uniqueness |
|---|---|---|
| , size | determined by alone | |
| Canonical , size | unique with on | |
| Unitary with | exists when | unique only when is invertible |
| in , size | determined by alone |
The canonical is a partial isometry, with and . It is an isometry, , exactly when has full column rank, and unitary exactly when is square and invertible, in which case .
For example, with rows and has eigenvalues and , and its positive root has rows and . The matrix with rows and also squares to but is indefinite. The matrix with rows and has , and the quarter-turn with rows and . For the singular , the canonical factor is , and both and 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 with coefficients in and , the matrix satisfies , so polynomials in the same commute. A nonzero with is an annihilating polynomial. One exists because the matrices are dependent.
- Minimal polynomial. The minimal polynomial is the unique monic annihilating polynomial of least degree, and exactly when divides .
- Roots. A scalar is a root of exactly when it is an eigenvalue of .
- Cayley-Hamilton. Every satisfies , so divides and . The proof compares coefficients in . Substituting into is not a proof.
| Matrix | ||
|---|---|---|
| shear with rows and | ||
| quarter-turn with rows and |
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 with rows and , Cayley-Hamilton gives , then and
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 and are coprime, there are polynomials and with
If moreover for , then
and is the projection onto along , with . These projectors are polynomials in and need not be orthogonal.
The diagonalizability criterion says that is diagonalizable over if and only if splits over into distinct linear factors, . The eigenspace projectors are then the Lagrange polynomials
with , for and . Splitting is part of the hypothesis. The real quarter-turn has with distinct complex roots and is diagonalizable over but not over . A finite family of commuting operators, each diagonalizable over , has a common eigenbasis, which need not be orthonormal.
For example, with rows and has and , with rows and . They split as , an eigenvector for plus an eigenvector for , along lines that are not perpendicular.
Lessons Polynomial Projectors and Diagonalizability
Cyclic subspaces and companion matrices
The orbit of a nonzero under is . If is the smallest number for which are dependent, the first dependence
defines the monic polynomial . The cyclic subspace has basis , is invariant under , and is the smallest invariant subspace containing . A polynomial satisfies exactly when divides , so divides . In the orbit basis the restriction of has the companion matrix
whose characteristic and minimal polynomials both equal .
- Cyclic vectors. If , then is a cyclic vector and . A scalar matrix on a space of dimension at least has no cyclic vector.
- Recurrences. The orbit vectors satisfy , and so does every scalar sequence with .
For example, with rows and sends to and to , so and is the companion matrix of . The seed is an eigenvector for and generates only a line.
Lessons Cyclic Subspaces and Companion Matrices
Generalized eigenspaces
Fix an eigenvalue of an operator on an -dimensional space and let . A generalized eigenvector is a nonzero with for some , and its rank is the smallest such . Rank one means an ordinary eigenvector. A vector of rank generates a Jordan chain, listed bottom first,
so and for . Every chain is linearly independent. The kernels stop growing after the first equality, and by step at the latest, so the generalized eigenspace is
It is invariant under and contains the eigenspace.
The primary decomposition theorem says that if splits over , with distinct , which always holds over , then
with . On the operator is nilpotent with index exactly , and is the algebraic multiplicity of . The projectors come from Bézout’s identity for the coprime factors . Over the splitting hypothesis matters, since a real quarter-turn has no real eigenvalue and no nonzero real generalized eigenspace.
For example, with rows , and has the chain under , kernel dimensions , and . The matrix with rows , and has , and , and projects onto .
Lessons Generalized Eigenvectors, Generalized Eigenspace Decomposition
Jordan form
A Jordan block of size has on the diagonal and directly above it,
which is the matrix of the operator on the span of one Jordan chain, ordered bottom first. The Jordan form theorem says that if splits over , which always holds for , then for an invertible and a block diagonal of Jordan blocks. The number of blocks of each size for each eigenvalue is unique, so is unique up to the order of its blocks, while is not unique. The columns of are Jordan chains, and checking with invertible verifies a proposed answer.
-
Construction. On with , the spaces decrease in . Choose a basis of the highest nonzero and extend it downward. Each vector first added at level is the bottom of one chain of length , lifted through a with , and shows that the chains span.
-
Counting. With and , the number of blocks of size at least and the number of blocks of size exactly are
| Quantity for an eigenvalue | Meaning in the Jordan form |
|---|---|
| , the geometric multiplicity | number of blocks for |
| , the algebraic multiplicity | total size of those blocks |
| exponent of in | size of the largest block |
For example, a nilpotent operator on a four-dimensional space with has and blocks of sizes and . Kernel dimensions on a five-dimensional space give blocks of sizes and , minimal polynomial and .
Exact Jordan structure is fragile. The matrix with rows and is one block of size at , but for every it has the distinct eigenvalues , with eigenvectors that are nearly parallel when 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
It satisfies , when the middle norms match, and when one norm is used on both sides.
| Norm | Value | Comment |
|---|---|---|
| largest column sum of the | induced by | |
| largest row sum of the | induced by | |
| induced by the Euclidean norm, equals | ||
| for square sizes not induced by any one vector norm, since |
The Frobenius norm satisfies . For a square invertible and an induced norm, the condition number is
If with and the data change to , the solution changes by , and
This is a worst case, and many perturbations are amplified far less. The number is unchanged by unitary factors and by scaling , but rescaling coordinates can change it. A small determinant is no test, since has for every . For a singular matrix the inverse problem has no finite condition number.
For example, has . With , the perturbation gives , a thousand times larger in relative size, while the same perturbation along is not amplified. A matrix with singular values , and has , , and .
Lessons Matrix Norms and Condition Numbers
Residuals, forward error and backward error
Let be square and invertible, , the exact solution of and a computed one, with all norms consistent and induced.
| Quantity | Definition | What it measures |
|---|---|---|
| Residual | mismatch in the equation, computable | |
| Relative forward error | distance to the true answer | |
| Backward error, data only | solves exactly | |
| Backward error, matrix only | for | smallest with |
The matrix-only bound is attained by the rank-one change . In the model that changes only , subtracting the equations gives , and
- 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 . In the usual model each rounded operation gives the exact result times with , ignoring overflow and underflow, and in IEEE double precision . 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, , and give . The relative residual is , yet the relative forward error is , exactly times the relative residual.
Lessons Residuals, Forward Error and Backward Error
Householder QR factorization
For a nonzero real vector , the Householder reflector
is symmetric and orthogonal with . It reverses the component along and fixes everything perpendicular to . To send a real column with a nonzero tail to a multiple of , let and with . Then
This sign adds magnitudes in the first entry and avoids cancellation. If the tail of is already zero, use the identity, and never divide by . 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 for with . In the full form is orthogonal and is with an upper triangular top block and zero rows below, and the reduced form is . With full column rank is invertible, and because orthogonal matrices keep lengths,
So the least-squares solution solves by back substitution, and the smallest residual norm is . A reduced calculation must keep that last term in mind.
For example, has and , and the reflector sends to . For with rows , and and , the least-squares solution is with residual , which is orthogonal to both columns.
Lessons Householder QR Factorization
QR, normal equations and Cholesky
For a real with full column rank, the condition number is , and
The normal equations are exactly correct, but they solve a system with the squared condition number, and forming rounds away information. For the columns and with , the exact Gram matrix has rows and and determinant , but rounding each entry to six significant digits makes it singular, although the columns are independent.
- Cholesky solve. For a given positive definite system , factor and solve and then 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 , and the same row swaps are applied to . For the matrix with rows and , eliminating without a swap uses the multiplier .
| Problem | Method | Reason |
|---|---|---|
| Least squares with full column rank | Householder QR | works with , not |
| A given positive definite system | Cholesky | uses the structure that is already there |
| A general square system | elimination with partial pivoting | keeps the multipliers at most in modulus |
| Rank in doubt or noisy data | SVD or pivoted QR with a stated tolerance | shows small singular values or nearly dependent columns |
For example, with with rows , , and , the normal equations have matrix rows , and right side , and give the same as QR. In the Cholesky factorization of that Gram matrix, the lower triangular has rows and .
Lessons QR, Normal Equations and Cholesky
Power iteration
Power iteration starts from a nonzero and repeats and . If some , stop and report that lies in .
The convergence theorem assumes that , a real matrix being viewed as complex, is diagonalizable with eigenbasis of , where , and that . If with , then
so the line through approaches the line through , with error of order . A negative or nonreal makes change sign or phase while its line converges. The coefficient is a coordinate in the eigenbasis. It equals when the eigenbasis is orthonormal, but not in general.
| Situation | What happens |
|---|---|
| in exact arithmetic the iterates never approach | |
| with | directions can cycle, as for |
| not diagonalizable | the proof does not apply, and a defective dominant eigenvalue can slow convergence to a rate like |
| nearly dependent eigenbasis | large constants and slow visible progress |
For a unit vector and a scalar , the residual checks an approximate eigenpair. If is normal, then
but the bound fails for other matrices. The matrix with rows and , and proportional to give a residual , although the only eigenvalue is at distance .
For example, from gives , a component ratio of . For from a start with every component nonzero the error ratio is , and the start converges to instead, with ratio . For the matrix with rows and , the start is an eigenvector for and stays on its line, even though its dot product with the dominant eigenvector 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 , divide by the minimal polynomial, with . Then , and two polynomials give the same matrix exactly when divides their difference. Agreement at the eigenvalues is not enough for a defective matrix. The shear with rows and has the single eigenvalue , yet the polynomials and give .
Polynomial functions respect similarity, , so a diagonalizable has , with applied to the diagonal entries of . A matrix function is not computed entry by entry. For a Jordan block with , the two parts commute and the binomial theorem gives, for every ,
For the entries are times polynomials in of degree at most , so a long chain adds polynomial growth to the rate .
| Eigenvalues of | Powers as |
|---|---|
| every | tend to , possibly after transient growth |
| every , and blocks with of size | stay bounded |
| some in a Jordan block of size at least | unbounded |
| some | unbounded |
For example, has rows and , so drifts without bound although every eigenvalue is . The Jordan block of size with eigenvalue has , and with has . The matrix with rows and has powers with upper right entry , which tend to , but it sends to a vector of length greater than .
Lessons Matrix Powers and Polynomial Functions
The matrix exponential
For and a scalar , the matrix exponential is the series
which converges absolutely because in an induced norm, and .
| Property | Condition |
|---|---|
| always | |
| invertible | |
| always | |
| always | |
| real , using calculus |
Without commuting matrices the product rule can fail. For with rows and with rows , the product has rows and , while has both diagonal entries equal to . For a Jordan block with ,
and a diagonalizable has . For real , . If every eigenvalue has negative real part, then as , but can grow first. With calculus, is the solution of with .
For example, with rows and has
At it sends to , of length , although every trajectory tends to .
Lessons The Matrix Exponential
Coupled updates and eigenmodes
A linear update has . If is real symmetric with orthonormal eigenvectors and eigenvalues , each mode evolves on its own,
For a diagonalizable that is not symmetric, the coefficients come from , not from dot products, and the modes need not be orthogonal. A positive definite form can measure the size of the state. In an oblique basis with matrix the operator changes to by similarity while the form changes to by congruence, and the two rules agree only for orthogonal . 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 with rows and . The modes decay by the factors and , and from
Here with having rows and , and is positive for , with having eigenvalues and . The update has factors and , so it decays for every start exactly when . Changing the upper right coupling to gives eigenvalues with eigenvectors , which are not orthogonal.
Lessons Coupled Systems and Eigenmodes
Inverse problems and regularization
An inverse problem asks for from noisy data . With a compact singular value decomposition and , the usual reconstructions have the form
where the filter factors decide how much of each direction is kept. Small amplify noise, and components in are never measured.
| Method | Filter factor | Effect |
|---|---|---|
| Pseudoinverse | smallest residual, divides noise by small | |
| Truncated SVD at threshold | when , otherwise | drops the weak directions |
| Tikhonov with parameter | shrinks the weak directions gradually |
Tikhonov regularization minimizes , where is a penalty parameter, not an eigenvalue. The matrix is positive definite even when is rank deficient, so the minimizer is unique,
Increasing never increases the solution norm and never decreases the residual. Choosing or needs a noise model, a known tolerance or outside validation, and changing units changes what the penalty measures. All three methods set the component to zero by convention, which is not a measurement of it.
For example, and give the pseudoinverse solution with residual , the truncation at gives with residual norm , and Tikhonov with gives . A change in the second data entry reaches the second coefficient with gain , and in the three methods.
Lessons Inverse Problems and Regularization
Results that combine the course
- Dual of a quotient. For a subspace of a finite-dimensional , the pullback through the quotient map is an isomorphism from onto . A functional descends to exactly when it vanishes on . For in , the descending functionals are the combinations of and , such as .
- Normality test. Compare with before using the spectral theorem. The matrix with rows and has with rows and with rows , so it is not normal, and its one-dimensional eigenspace shows it is not diagonalizable either.
- Verifying a Jordan basis. A proposed and are right when is invertible and , checked column by column. For with rows , , , the columns give , and projects onto the eigenvector line for along the generalized eigenspace for , obliquely, since it is not symmetric.
- Orthogonal projectors. An idempotent is an orthogonal projector exactly when for the chosen inner product. The matrix with rows and is idempotent, but its kernel, spanned by , is not perpendicular to its image.
- Forward differences. The operator on lowers the degree of every nonconstant polynomial by one, so while , and its Jordan form is a single nilpotent block of size . On a chain is .
- 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 .
Lessons Proofs and Problem Solving