Documentation

Mathlib.Analysis.InnerProductSpace.Spectrum

Spectral theory of self-adjoint operators #

This file covers the spectral theory of self-adjoint operators on an inner product space.

The first part of the file covers general properties, true without any condition on boundedness or compactness of the operator or finite-dimensionality of the underlying space, notably:

The second part of the file covers properties of self-adjoint operators in finite dimension. Letting T be a self-adjoint operator on a finite-dimensional inner product space T,

These are forms of the diagonalization theorem for self-adjoint operators on finite-dimensional inner product spaces.

TODO #

Spectral theory for compact self-adjoint operators, bounded self-adjoint operators.

Tags #

self-adjoint operator, spectral theorem, diagonalization theorem

theorem LinearMap.IsSymmetric.invariant_orthogonalComplement_eigenspace {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) (ฮผ : ๐•œ) (v : E) (hv : v โˆˆ (Module.End.eigenspace T ฮผ)แ—ฎ) :

A self-adjoint operator preserves orthogonal complements of its eigenspaces.

theorem LinearMap.IsSymmetric.conj_eigenvalue_eq_self {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) {ฮผ : ๐•œ} (hฮผ : Module.End.HasEigenvalue T ฮผ) :
(starRingEnd ๐•œ) ฮผ = ฮผ

The eigenvalues of a self-adjoint operator are real.

theorem LinearMap.IsSymmetric.orthogonalFamily_eigenspaces {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) :
OrthogonalFamily ๐•œ (fun (ฮผ : ๐•œ) => โ†ฅ(Module.End.eigenspace T ฮผ)) fun (ฮผ : ๐•œ) => Submodule.subtypeโ‚—แตข (Module.End.eigenspace T ฮผ)

The eigenspaces of a self-adjoint operator are mutually orthogonal.

theorem LinearMap.IsSymmetric.orthogonalFamily_eigenspaces' {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) :
OrthogonalFamily ๐•œ (fun (ฮผ : Module.End.Eigenvalues T) => โ†ฅ(Module.End.eigenspace T (โ†‘T ฮผ))) fun (ฮผ : Module.End.Eigenvalues T) => Submodule.subtypeโ‚—แตข (Module.End.eigenspace T (โ†‘T ฮผ))
theorem LinearMap.IsSymmetric.orthogonalComplement_iSup_eigenspaces_invariant {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) โฆƒv : Eโฆ„ (hv : v โˆˆ (โจ† (ฮผ : ๐•œ), Module.End.eigenspace T ฮผ)แ—ฎ) :
T v โˆˆ (โจ† (ฮผ : ๐•œ), Module.End.eigenspace T ฮผ)แ—ฎ

The mutual orthogonal complement of the eigenspaces of a self-adjoint operator on an inner product space is an invariant subspace of the operator.

theorem LinearMap.IsSymmetric.orthogonalComplement_iSup_eigenspaces {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) (ฮผ : ๐•œ) :

The mutual orthogonal complement of the eigenspaces of a self-adjoint operator on an inner product space has no eigenvalues.

Finite-dimensional theory #

theorem LinearMap.IsSymmetric.orthogonalComplement_iSup_eigenspaces_eq_bot {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] :
(โจ† (ฮผ : ๐•œ), Module.End.eigenspace T ฮผ)แ—ฎ = โŠฅ

The mutual orthogonal complement of the eigenspaces of a self-adjoint operator on a finite-dimensional inner product space is trivial.

theorem LinearMap.IsSymmetric.orthogonalComplement_iSup_eigenspaces_eq_bot' {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] :
(โจ† (ฮผ : Module.End.Eigenvalues T), Module.End.eigenspace T (โ†‘T ฮผ))แ—ฎ = โŠฅ
noncomputable instance LinearMap.IsSymmetric.directSumDecomposition {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} [FiniteDimensional ๐•œ E] [hT : Fact (LinearMap.IsSymmetric T)] :

The eigenspaces of a self-adjoint operator on a finite-dimensional inner product space E gives an internal direct sum decomposition of E.

Note this takes hT as a Fact to allow it to be an instance.

Equations
theorem LinearMap.IsSymmetric.directSum_decompose_apply {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} [FiniteDimensional ๐•œ E] [_hT : Fact (LinearMap.IsSymmetric T)] (x : E) (ฮผ : Module.End.Eigenvalues T) :
((DirectSum.decompose fun (ฮผ : Module.End.Eigenvalues T) => Module.End.eigenspace T (โ†‘T ฮผ)) x) ฮผ = (orthogonalProjection (Module.End.eigenspace T (โ†‘T ฮผ))) x
theorem LinearMap.IsSymmetric.direct_sum_isInternal {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] :
DirectSum.IsInternal fun (ฮผ : Module.End.Eigenvalues T) => Module.End.eigenspace T (โ†‘T ฮผ)

The eigenspaces of a self-adjoint operator on a finite-dimensional inner product space E gives an internal direct sum decomposition of E.

noncomputable def LinearMap.IsSymmetric.diagonalization {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] :
E โ‰ƒโ‚—แตข[๐•œ] PiLp 2 fun (ฮผ : Module.End.Eigenvalues T) => โ†ฅ(Module.End.eigenspace T (โ†‘T ฮผ))

Isometry from an inner product space E to the direct sum of the eigenspaces of some self-adjoint operator T on E.

Equations
Instances For
    @[simp]
    theorem LinearMap.IsSymmetric.diagonalization_symm_apply {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] (w : PiLp 2 fun (ฮผ : Module.End.Eigenvalues T) => โ†ฅ(Module.End.eigenspace T (โ†‘T ฮผ))) :
    (LinearIsometryEquiv.symm (LinearMap.IsSymmetric.diagonalization hT)) w = Finset.sum Finset.univ fun (ฮผ : Module.End.Eigenvalues T) => โ†‘(w ฮผ)
    theorem LinearMap.IsSymmetric.diagonalization_apply_self_apply {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] (v : E) (ฮผ : Module.End.Eigenvalues T) :

    Diagonalization theorem, spectral theorem; version 1: A self-adjoint operator T on a finite-dimensional inner product space E acts diagonally on the decomposition of E into the direct sum of the eigenspaces of T.

    @[irreducible]
    noncomputable def LinearMap.IsSymmetric.eigenvectorBasis {๐•œ : Type u_3} [RCLike ๐•œ] {E : Type u_4} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] {n : โ„•} (hn : FiniteDimensional.finrank ๐•œ E = n) :
    OrthonormalBasis (Fin n) ๐•œ E

    A choice of orthonormal basis of eigenvectors for self-adjoint operator T on a finite-dimensional inner product space E.

    TODO Postcompose with a permutation so that these eigenvectors are listed in increasing order of eigenvalue.

    Equations
    Instances For
      theorem LinearMap.IsSymmetric.eigenvalues_def {๐•œ : Type u_3} [RCLike ๐•œ] {E : Type u_4} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] {n : โ„•} (hn : FiniteDimensional.finrank ๐•œ E = n) (i : Fin n) :
      @[irreducible]
      noncomputable def LinearMap.IsSymmetric.eigenvalues {๐•œ : Type u_3} [RCLike ๐•œ] {E : Type u_4} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] {n : โ„•} (hn : FiniteDimensional.finrank ๐•œ E = n) (i : Fin n) :

      The sequence of real eigenvalues associated to the standard orthonormal basis of eigenvectors for a self-adjoint operator T on E.

      TODO Postcompose with a permutation so that these eigenvalues are listed in increasing order.

      Equations
      Instances For
        theorem LinearMap.IsSymmetric.hasEigenvalue_eigenvalues {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] {n : โ„•} (hn : FiniteDimensional.finrank ๐•œ E = n) (i : Fin n) :
        @[simp]
        theorem LinearMap.IsSymmetric.apply_eigenvectorBasis {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] {n : โ„•} (hn : FiniteDimensional.finrank ๐•œ E = n) (i : Fin n) :
        theorem LinearMap.IsSymmetric.eigenvectorBasis_apply_self_apply {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {T : E โ†’โ‚—[๐•œ] E} (hT : LinearMap.IsSymmetric T) [FiniteDimensional ๐•œ E] {n : โ„•} (hn : FiniteDimensional.finrank ๐•œ E = n) (v : E) (i : Fin n) :

        Diagonalization theorem, spectral theorem; version 2: A self-adjoint operator T on a finite-dimensional inner product space E acts diagonally on the identification of E with Euclidean space induced by an orthonormal basis of eigenvectors of T.

        @[simp]
        theorem inner_product_apply_eigenvector {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {ฮผ : ๐•œ} {v : E} {T : E โ†’โ‚—[๐•œ] E} (h : v โˆˆ Module.End.eigenspace T ฮผ) :
        โŸชv, T vโŸซ_๐•œ = ฮผ * โ†‘โ€–vโ€– ^ 2
        theorem eigenvalue_nonneg_of_nonneg {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {ฮผ : โ„} {T : E โ†’โ‚—[๐•œ] E} (hฮผ : Module.End.HasEigenvalue T โ†‘ฮผ) (hnn : โˆ€ (x : E), 0 โ‰ค RCLike.re โŸชx, T xโŸซ_๐•œ) :
        0 โ‰ค ฮผ
        theorem eigenvalue_pos_of_pos {๐•œ : Type u_1} [RCLike ๐•œ] {E : Type u_2} [NormedAddCommGroup E] [InnerProductSpace ๐•œ E] {ฮผ : โ„} {T : E โ†’โ‚—[๐•œ] E} (hฮผ : Module.End.HasEigenvalue T โ†‘ฮผ) (hnn : โˆ€ (x : E), 0 < RCLike.re โŸชx, T xโŸซ_๐•œ) :
        0 < ฮผ