Prove2Me
Navigate
DiscoverFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

The innovation is orthogonal to past data

Proved
VectorSpaceOpt.innovation_orthogonal

by Shuze Chen · Aug 23, 2026 · Mathlib 0df444a (Lean v4.33.1)

estimationhilbert-spacekalman-filter

Throughout, we work in the Hilbert space of random variables: zero-mean random variables with finite second moments, under the inner product ⟨a,b⟩=E[ab]\langle a, b\rangle = E[ab]⟨a,b⟩=E[ab]. In this space, orthogonality is uncorrelatedness, and the best linear estimate of a variable given some data is its orthogonal projection onto the subspace the data generates.

Let SSS be the subspace generated by past data, let β\betaβ be a random nnn-vector, and let β^\hat\betaβ^​ be its componentwise orthogonal projection onto SSS — the best estimate of β\betaβ from the past. Let new measurements arrive as

y=Wβ+ε,y = W\beta + \varepsilon,y=Wβ+ε,

with WWW a known m×nm \times nm×n matrix and each noise component εi\varepsilon_iεi​ orthogonal to SSS. Then the innovation

y~=y−Wβ^\tilde y = y - W\hat\betay~​=y−Wβ^​

is orthogonal to SSS: every component of y~\tilde yy~​ is uncorrelated with every element of the past-data subspace.

The innovation is the part of the new measurement that could not have been anticipated from the old data. Its orthogonality to SSS is what makes recursive estimation possible: updating an estimate requires only the innovation, never a recomputation over the whole data history.

Formalization Note. The action of a matrix on a random vector is componentwise, (Wβ)i=∑jWij βj(W\beta)_i = \sum_j W_{ij}\,\beta_j(Wβ)i​=∑j​Wij​βj​ with scalar multiplication in the abstract space. Zero means are implicit in this representation, so no expectation operator appears — only inner products.

Preamble
import Mathlib
open scoped RealInnerProductSpace
Formal statement
namespace VectorSpaceOpt

theorem innovation_orthogonal {H : Type} [NormedAddCommGroup H]
    [InnerProductSpace ℝ H] {n m : ℕ} (S : Submodule ℝ H)
    (β βh : Fin n → H)
    (hproj : ∀ i, ∀ s ∈ S, ⟪β i - βh i, s⟫ = 0)
    (W : Matrix (Fin m) (Fin n) ℝ) (ε : Fin m → H)
    (hεS : ∀ i, ∀ s ∈ S, ⟪ε i, s⟫ = 0)
    (y : Fin m → H) (hy : ∀ i, y i = ∑ j, W i j • β j + ε i) :
    ∀ i, ∀ s ∈ S, ⟪y i - ∑ j, W i j • βh j, s⟫ = 0 := by sorry

end VectorSpaceOpt
Source
David G. Luenberger, Optimization by Vector Space Methods, John Wiley & Sons, 1969, §4.6 (updating) and §4.7, pp. 92–93

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me