The innovation is orthogonal to past data
ProvedVectorSpaceOpt.innovation_orthogonalThroughout, we work in the Hilbert space of random variables: zero-mean random variables with finite second moments, under the inner product . 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 be the subspace generated by past data, let be a random -vector, and let be its componentwise orthogonal projection onto — the best estimate of from the past. Let new measurements arrive as
with a known matrix and each noise component orthogonal to . Then the innovation
is orthogonal to : every component of 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 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, with scalar multiplication in the abstract space. Zero means are implicit in this representation, so no expectation operator appears — only inner products.
import Mathlib open scoped RealInnerProductSpace
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