Documentation

LeanPool.DavisKahan.ForTauCeti.Analysis.InnerProductSpace.SchurHorn

Schur–Horn majorization (forward direction, Karamata form) #

Let T be a symmetric operator on a finite-dimensional inner product space over π•œ = ℝ, β„‚, with sorted eigenvalues Ξ» (hT.eigenvalues hn) and orthonormal eigenbasis v (hT.eigenvectorBasis hn). Fix any orthonormal basis e. The "diagonal" of T in e is the tuple d k = re βŸͺT (e k), e k⟫.

The forward direction of the Schur–Horn theorem (due to Schur, 1923) says the diagonal is majorized by the spectrum, diag T β‰Ί spec T. We prove the equivalent Hardy–Littlewood–PΓ³lya / Karamata characterisation: βˆ‘ Ο† (d k) ≀ βˆ‘ Ο† (Ξ» i) for every convex Ο† defined on a set containing the eigenvalues.

The mechanism is the doubly-stochastic weight matrix w i k = β€–βŸͺvα΅’, e kβŸ«β€–Β² (schurWeight): its rows and columns sum to 1 by Parseval, and the diagonal is its image of the spectrum, d k = βˆ‘ i, Ξ» i * w i k. Row-wise Jensen followed by a sum swap over the column sums gives the inequality.

Mathlib has the spectral theorem and Birkhoff's theorem but no majorization predicate and no Schur–Horn theorem (only a comment in Mathlib/Analysis/InnerProductSpace/Spectrum.lean); this file supplies the forward direction in the self-contained convex-function form.

Main results #

References #

Provenance #

noncomputable def TauCeti.schurWeight {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) (i k : Fin n) :

The doubly-stochastic weight w i k = β€–βŸͺvα΅’, e kβŸ«β€–Β² of the i-th eigenvector vα΅’ of T against the k-th vector of a chosen orthonormal basis e.

Equations
Instances For
    theorem TauCeti.schurWeight_nonneg {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) (i k : Fin n) :
    0 ≀ schurWeight hT hn e i k

    Schur--Horn weights are nonnegative, being squared moduli of basis coefficients.

    theorem TauCeti.schurWeight_row_sum {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) (k : Fin n) :
    βˆ‘ i : Fin n, schurWeight hT hn e i k = 1

    Rows sum to one. By Parseval for the eigenbasis v, βˆ‘ i, β€–βŸͺvα΅’, e kβŸ«β€–Β² = β€–e kβ€–Β² = 1.

    theorem TauCeti.schurWeight_col_sum {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) (i : Fin n) :
    βˆ‘ k : Fin n, schurWeight hT hn e i k = 1

    Columns sum to one. By Parseval for the basis e, βˆ‘ k, β€–βŸͺvα΅’, e kβŸ«β€–Β² = β€–vα΅’β€–Β² = 1.

    theorem TauCeti.re_inner_orthonormalBasis_self_eq_sum_eigenvalues_mul {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) (k : Fin n) :
    RCLike.re (inner π•œ (T (e k)) (e k)) = βˆ‘ i : Fin n, hT.eigenvalues hn i * schurWeight hT hn e i k

    Diagonal = doubly-stochastic image of the spectrum. The diagonal entry re βŸͺT (e k), e k⟫ of T in the basis e is the convex combination βˆ‘ i, Ξ» i * w i k of the eigenvalues. Immediate from the diagonalisation of the quadratic form together with vβ±Ό.repr (e k) i = βŸͺvα΅’, e k⟫.

    theorem TauCeti.convexOn_sum_re_inner_orthonormalBasis_self_le {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) {Ο† : ℝ β†’ ℝ} {s : Set ℝ} (hΟ† : ConvexOn ℝ s Ο†) (hmem : βˆ€ (i : Fin n), hT.eigenvalues hn i ∈ s) :
    βˆ‘ k : Fin n, Ο† (RCLike.re (inner π•œ (T (e k)) (e k))) ≀ βˆ‘ i : Fin n, Ο† (hT.eigenvalues hn i)

    Forward Schur–Horn theorem (convex / Karamata form). For any convex function Ο† on a set s containing all eigenvalues of the symmetric operator T, the diagonal of T in any orthonormal basis e satisfies βˆ‘ k, Ο† (re βŸͺT (e k), e k⟫) ≀ βˆ‘ i, Ο† (Ξ» i), i.e. diag T β‰Ί spec T.

    Row-by-row Jensen against the doubly-stochastic weight matrix schurWeight, followed by a sum swap collapsing the column sums.

    theorem TauCeti.sum_re_inner_orthonormalBasis_self_eq_sum_eigenvalues {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) :
    βˆ‘ k : Fin n, RCLike.re (inner π•œ (T (e k)) (e k)) = βˆ‘ i : Fin n, hT.eigenvalues hn i

    Basis independence of the trace (the equality case of Schur majorization). The sum of the diagonal entries of T in any orthonormal basis equals the sum of its eigenvalues. No convexity needed.

    theorem TauCeti.sum_sq_re_inner_orthonormalBasis_self_le_sum_sq_eigenvalues {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [InnerProductSpace π•œ E] [FiniteDimensional π•œ E] {n : β„•} {T : E β†’β‚—[π•œ] E} (hT : T.IsSymmetric) (hn : Module.finrank π•œ E = n) (e : OrthonormalBasis (Fin n) π•œ E) :
    βˆ‘ k : Fin n, RCLike.re (inner π•œ (T (e k)) (e k)) ^ 2 ≀ βˆ‘ i : Fin n, hT.eigenvalues hn i ^ 2

    Sum-of-squares bound (the Ο† = (Β·)Β² instance of Schur majorization). The diagonal of T in any orthonormal basis has Euclidean length no larger than the spectrum: βˆ‘ k, (re βŸͺT (e k), e k⟫)Β² ≀ βˆ‘ i, (Ξ» i)Β².