Quadratic form of a symmetric operator in a Hilbert basis of eigenvectors
ProvedAsymptoticOperator.hasSum_inner_self_of_hilbertBasis_eigenvectorsLet be a real inner product space and let be a linear operator defined on a subspace that is symmetric: for all . Let be a Hilbert basis of consisting of eigenvectors of in , with for real numbers . Then for every the series converges unconditionally and
Role: expresses the quadratic form of a symmetric operator with an orthonormal eigenbasis through its eigenvalues. Applies to any such operator, for example the asymptotic operators on .
Formalization note: E →ₗ.[ℝ] E is Mathlib's LinearPMap (a linear map on the submodule A.domain), A.IsFormalAdjoint A is the symmetry condition, HilbertBasis ι ℝ E is Mathlib's Hilbert basis, and convergence is HasSum (unconditional summation over ι). The index type and E live in Type; no completeness of is assumed.
import Mathlib.Analysis.InnerProductSpace.l2Space import Mathlib.Analysis.InnerProductSpace.LinearPMap
namespace AsymptoticOperator
/-- Let `A` be a symmetric partially defined operator on a real inner product space, and let `b`
be a Hilbert basis of eigenvectors of `A` lying in its domain, `A bᵢ = evᵢ bᵢ`. Then for every
`x` in the domain the quadratic form of `A` is `⟨A x, x⟩ = ∑ᵢ evᵢ ⟨x, bᵢ⟩²`. -/
theorem hasSum_inner_self_of_hilbertBasis_eigenvectors {E ι : Type} [NormedAddCommGroup E]
[InnerProductSpace ℝ E] (A : E →ₗ.[ℝ] E) (hA : A.IsFormalAdjoint A)
(b : HilbertBasis ι ℝ E) (ev : ι → ℝ)
(hb : ∀ i, ∃ h : b i ∈ A.domain, A ⟨b i, h⟩ = ev i • b i) (x : A.domain) :
HasSum (fun i => ev i * inner ℝ (x : E) (b i) ^ 2) (inner ℝ (A x) (x : E)) := by sorry
end AsymptoticOperator