Fix an admissible control α \alpha α ; write Y : = X α Y:=X^{\alpha} Y := X α , X ^ : = X ^ ( α ) \widehat X:=\widehat X(\alpha) X := X ( α ) , and let e t = X t − m t f e_t=X_t-m^{\mathrm f}_t e t = X t − m t f be the estimation error, so that by claim 3 of Conditional Expectation and Estimation Error of the Controlled State : Y t − X ^ t = e t Y_t-\widehat X_t=e_t Y t − X t = e t almost surely componentwise, E [ e t i ] = 0 \mathbb{E}[e^{i}_t]=0 E [ e t i ] = 0 , the covariance matrix of e t e_t e t is Π ( t ) \Pi(t) Π ( t ) , and σ ( e t 1 , … , e t l ) \sigma(e^{1}_t,\dots,e^{l}_t) σ ( e t 1 , … , e t l ) is independent of G t \mathcal{G}_t G t . Throughout, integrands are continuous in t t t and Riemann integrals of continuous functions are linear and monotone, via Agreement of the Riemann and Lebesgue Integrals for Continuous Functions on a Closed Interval and Linearity and Monotonicity of the Lebesgue Integral .
Claim 1. By Completion of Squares for the Linear-Quadratic-Gaussian Cost ,
J [ α ] = tr ( Z ( 0 ) P 0 ) + E [ ξ ] ⋅ ( Z ( 0 ) E [ ξ ] ) + ∫ 0 T tr ( Z ( t ) Θ ( t ) ) d t + ∫ 0 T E [ ( α t − Γ ( t ) Y t ) ⋅ ( R ( t ) ( α t − Γ ( t ) Y t ) ) ] d t . J[\alpha]=\operatorname{tr}\bigl(Z(0)P_0\bigr)+\mathbb{E}[\xi]\cdot\bigl(Z(0)\mathbb{E}[\xi]\bigr)+\int_0^T\operatorname{tr}\bigl(Z(t)\Theta(t)\bigr)dt+\int_0^T\mathbb{E}\bigl[(\alpha_t-\Gamma(t)Y_t)\cdot\bigl(R(t)(\alpha_t-\Gamma(t)Y_t)\bigr)\bigr]dt . J [ α ] = tr ( Z ( 0 ) P 0 ) + E [ ξ ] ⋅ ( Z ( 0 ) E [ ξ ] ) + ∫ 0 T tr ( Z ( t ) Θ ( t ) ) d t + ∫ 0 T E [ ( α t − Γ ( t ) Y t ) ⋅ ( R ( t ) ( α t − Γ ( t ) Y t ) ) ] d t .
Fix t t t and decompose, componentwise and almost surely,
α t − Γ ( t ) Y t = ( α t − Γ ( t ) X ^ t ) + ( − Γ ( t ) e t ) = : ζ + ρ . \alpha_t-\Gamma(t)Y_t=\bigl(\alpha_t-\Gamma(t)\widehat X_t\bigr)+\bigl(-\Gamma(t)e_t\bigr)=:\zeta+\rho . α t − Γ ( t ) Y t = ( α t − Γ ( t ) X t ) + ( − Γ ( t ) e t ) =: ζ + ρ .
We apply Orthogonal Decomposition of Expected Quadratic Forms under Independence with H : = G t \mathcal{H}:=\mathcal{G}_t H := G t and M : = R ( t ) M:=R(t) M := R ( t ) . Hypothesis (i): each component of ζ \zeta ζ is almost surely equal to a G t \mathcal{G}_t G t -measurable square-integrable random variable, by condition (ii) of Admissible Control for the Linear-Gaussian State-Observation Model , claim 1 of Conditional Expectation and Estimation Error of the Controlled State , and closure under finite linear combinations (claim 2 of The Closed Mean-Square Span of a Family of Random Variables ). Hypothesis (ii): each ρ κ = − ∑ i Γ κ i ( t ) e t i \rho^{\kappa}=-\sum_i\Gamma_{\kappa i}(t)e^{i}_t ρ κ = − ∑ i Γ κi ( t ) e t i has expectation 0 0 0 . Hypothesis (iii): each ρ κ \rho^{\kappa} ρ κ is a Borel (linear) function of the tuple ( e t 1 , … , e t l ) (e^{1}_t,\dots,e^{l}_t) ( e t 1 , … , e t l ) , so every preimage of a Borel set under a component of ρ \rho ρ lies in σ ( e t 1 , … , e t l ) \sigma(e^{1}_t,\dots,e^{l}_t) σ ( e t 1 , … , e t l ) ; hence σ ( ρ 1 , … , ρ k ) ⊆ σ ( e t 1 , … , e t l ) \sigma(\rho^{1},\dots,\rho^{k})\subseteq\sigma(e^{1}_t,\dots,e^{l}_t) σ ( ρ 1 , … , ρ k ) ⊆ σ ( e t 1 , … , e t l ) , which is independent of G t \mathcal{G}_t G t . (The expectation of the quadratic form depends on α t − Γ ( t ) Y t \alpha_t-\Gamma(t)Y_t α t − Γ ( t ) Y t only through its almost sure class, so the decomposition may be used after modification on a null set.) The lemma gives
E [ ( α t − Γ ( t ) Y t ) ⋅ ( R ( t ) ( α t − Γ ( t ) Y t ) ) ] = E [ ( α t − Γ ( t ) X ^ t ) ⋅ ( R ( t ) ( α t − Γ ( t ) X ^ t ) ) ] + tr ( R ( t ) C ρ ( t ) ) , \mathbb{E}\bigl[(\alpha_t-\Gamma(t)Y_t)\cdot\bigl(R(t)(\alpha_t-\Gamma(t)Y_t)\bigr)\bigr]=\mathbb{E}\bigl[(\alpha_t-\Gamma(t)\widehat X_t)\cdot\bigl(R(t)(\alpha_t-\Gamma(t)\widehat X_t)\bigr)\bigr]+\operatorname{tr}\bigl(R(t)C_{\rho}(t)\bigr), E [ ( α t − Γ ( t ) Y t ) ⋅ ( R ( t ) ( α t − Γ ( t ) Y t ) ) ] = E [ ( α t − Γ ( t ) X t ) ⋅ ( R ( t ) ( α t − Γ ( t ) X t ) ) ] + tr ( R ( t ) C ρ ( t ) ) ,
where C ρ ( t ) C_{\rho}(t) C ρ ( t ) is the covariance matrix of ρ = − Γ ( t ) e t \rho=-\Gamma(t)e_t ρ = − Γ ( t ) e t . By bilinearity of the covariance (linearity of the expectation applied to its defining formula),
C ρ ( t ) κ λ = ∑ i , j Γ κ i ( t ) Γ λ j ( t ) Cov ( e t i , e t j ) = ( Γ ( t ) Π ( t ) Γ ( t ) ⊤ ) κ λ , C_{\rho}(t)_{\kappa\lambda}=\sum_{i,j}\Gamma_{\kappa i}(t)\Gamma_{\lambda j}(t)\operatorname{Cov}(e^{i}_t,e^{j}_t)=\bigl(\Gamma(t)\Pi(t)\Gamma(t)^{\top}\bigr)_{\kappa\lambda}, C ρ ( t ) κλ = i , j ∑ Γ κi ( t ) Γ λj ( t ) Cov ( e t i , e t j ) = ( Γ ( t ) Π ( t ) Γ ( t ) ⊤ ) κλ ,
by the index formulas of the matrix product and transpose . Now apply the cyclic property (claim 3 of Basic Properties of the Trace ) to the pair of matrices R ( t ) Γ ( t ) Π ( t ) R(t)\Gamma(t)\Pi(t) R ( t ) Γ ( t ) Π ( t ) (size k × l k\times l k × l ) and Γ ( t ) ⊤ \Gamma(t)^{\top} Γ ( t ) ⊤ (size l × k l\times k l × k ), and compute, from the definition of the feedback gain in Completion of Squares for the Linear-Quadratic-Gaussian Cost : R ( t ) Γ ( t ) = − ( Z ( t ) B ( t ) + V ( t ) ) ⊤ R(t)\Gamma(t)=-(Z(t)B(t)+V(t))^{\top} R ( t ) Γ ( t ) = − ( Z ( t ) B ( t ) + V ( t ) ) ⊤ (since R ( t ) R ( t ) − 1 R(t)R(t)^{-1} R ( t ) R ( t ) − 1 is the identity, by associativity), and Γ ( t ) ⊤ = − ( Z ( t ) B ( t ) + V ( t ) ) R ( t ) − 1 \Gamma(t)^{\top}=-(Z(t)B(t)+V(t))R(t)^{-1} Γ ( t ) ⊤ = − ( Z ( t ) B ( t ) + V ( t )) R ( t ) − 1 (transpose rules of claim 3 of Componentwise Estimates, Transpose Identities, and Indefinite Riemann Integrals and symmetry of R ( t ) − 1 R(t)^{-1} R ( t ) − 1 , Invertibility of Symmetric Positive Definite Matrices ), so that
Γ ( t ) ⊤ R ( t ) Γ ( t ) = ( Z ( t ) B ( t ) + V ( t ) ) R ( t ) − 1 R ( t ) R ( t ) − 1 ( Z ( t ) B ( t ) + V ( t ) ) ⊤ = ( Z ( t ) B ( t ) + V ( t ) ) R ( t ) − 1 ( Z ( t ) B ( t ) + V ( t ) ) ⊤ . \Gamma(t)^{\top}R(t)\Gamma(t)=\bigl(Z(t)B(t)+V(t)\bigr)R(t)^{-1}R(t)R(t)^{-1}\bigl(Z(t)B(t)+V(t)\bigr)^{\top}=\bigl(Z(t)B(t)+V(t)\bigr)R(t)^{-1}\bigl(Z(t)B(t)+V(t)\bigr)^{\top}. Γ ( t ) ⊤ R ( t ) Γ ( t ) = ( Z ( t ) B ( t ) + V ( t ) ) R ( t ) − 1 R ( t ) R ( t ) − 1 ( Z ( t ) B ( t ) + V ( t ) ) ⊤ = ( Z ( t ) B ( t ) + V ( t ) ) R ( t ) − 1 ( Z ( t ) B ( t ) + V ( t ) ) ⊤ .
Hence
tr ( R ( t ) C ρ ( t ) ) = tr ( Γ ( t ) ⊤ R ( t ) Γ ( t ) Π ( t ) ) = tr ( ( Z ( t ) B ( t ) + V ( t ) ) R ( t ) − 1 ( Z ( t ) B ( t ) + V ( t ) ) ⊤ Π ( t ) ) . \operatorname{tr}\bigl(R(t)C_{\rho}(t)\bigr)=\operatorname{tr}\bigl(\Gamma(t)^{\top}R(t)\Gamma(t)\Pi(t)\bigr)=\operatorname{tr}\Bigl(\bigl(Z(t)B(t)+V(t)\bigr)R(t)^{-1}\bigl(Z(t)B(t)+V(t)\bigr)^{\top}\Pi(t)\Bigr). tr ( R ( t ) C ρ ( t ) ) = tr ( Γ ( t ) ⊤ R ( t ) Γ ( t ) Π ( t ) ) = tr ( ( Z ( t ) B ( t ) + V ( t ) ) R ( t ) − 1 ( Z ( t ) B ( t ) + V ( t ) ) ⊤ Π ( t ) ) .
Both t ↦ E [ ( α t − Γ ( t ) X ^ t ) ⋅ ( R ( t ) ( α t − Γ ( t ) X ^ t ) ) ] t\mapsto\mathbb{E}[(\alpha_t-\Gamma(t)\widehat X_t)\cdot(R(t)(\alpha_t-\Gamma(t)\widehat X_t))] t ↦ E [( α t − Γ ( t ) X t ) ⋅ ( R ( t ) ( α t − Γ ( t ) X t ))] (continuous by claim 2 of Expected Bilinear Forms: Trace Formula and Mean-Square Continuity , the components of α − Γ X ^ \alpha-\Gamma\widehat X α − Γ X being mean-square continuous by claim 1 of Conditional Expectation and Estimation Error of the Controlled State and claims 1-2 of Basic Properties of the Mean-Square Riemann Integral ) and the trace term (entries continuous) are continuous, and the first is nonnegative at each t t t : pointwise q ⋅ ( R ( t ) q ) ≥ 0 q\cdot(R(t)q)\ge0 q ⋅ ( R ( t ) q ) ≥ 0 since R ( t ) R(t) R ( t ) is positive definite , and expectations of nonnegative random variables are nonnegative (Linearity and Monotonicity of the Lebesgue Integral ). Integrating the pointwise identity over [ 0 , T ] [0,T] [ 0 , T ] and splitting the integral by linearity yields claim 1 with the stated V ∗ V^{*} V ∗ .
Claim 2. The inequality J [ α ] ≥ V ∗ J[\alpha]\ge V^{*} J [ α ] ≥ V ∗ is immediate from claim 1, the integrand being continuous and nonnegative and the Riemann integral monotone. If α t = Γ ( t ) X ^ t \alpha_t=\Gamma(t)\widehat X_t α t = Γ ( t ) X t almost surely for every t t t , the integrand vanishes identically and J [ α ] = V ∗ J[\alpha]=V^{*} J [ α ] = V ∗ . Conversely suppose J [ α ] = V ∗ J[\alpha]=V^{*} J [ α ] = V ∗ , i.e. ∫ 0 T g ( t ) d t = 0 \int_0^T g(t)\,dt=0 ∫ 0 T g ( t ) d t = 0 for the continuous nonnegative integrand g g g . If g ( t 0 ) > 0 g(t_0)>0 g ( t 0 ) > 0 for some t 0 t_0 t 0 , then by continuity g > g ( t 0 ) / 2 g>g(t_0)/2 g > g ( t 0 ) /2 on a nondegenerate subinterval, and splitting the integral by additivity (claim 4 of Componentwise Estimates, Transpose Identities, and Indefinite Riemann Integrals ) and monotonicity gives ∫ 0 T g > 0 \int_0^Tg>0 ∫ 0 T g > 0 , a contradiction; hence g ≡ 0 g\equiv0 g ≡ 0 . Fix t t t and put q : = α t − Γ ( t ) X ^ t q:=\alpha_t-\Gamma(t)\widehat X_t q := α t − Γ ( t ) X t ; then q ⋅ ( R ( t ) q ) ≥ 0 q\cdot(R(t)q)\ge0 q ⋅ ( R ( t ) q ) ≥ 0 pointwise with E [ q ⋅ ( R ( t ) q ) ] = 0 \mathbb{E}[q\cdot(R(t)q)]=0 E [ q ⋅ ( R ( t ) q )] = 0 . For every ε > 0 \varepsilon>0 ε > 0 , Markov's inequality (Markov's and Chebyshev's Inequalities ) gives P ( q ⋅ ( R ( t ) q ) ≥ ε ) = 0 P\bigl(q\cdot(R(t)q)\ge\varepsilon\bigr)=0 P ( q ⋅ ( R ( t ) q ) ≥ ε ) = 0 , and taking the union over ε = 1 / n \varepsilon=1/n ε = 1/ n (n ∈ N n\in\mathbb{N} n ∈ N ) with countable subadditivity of the measure yields q ⋅ ( R ( t ) q ) = 0 q\cdot(R(t)q)=0 q ⋅ ( R ( t ) q ) = 0 almost surely. On the event where in addition q q q is defined and q ≠ 0 q\ne0 q = 0 , positive definiteness of R ( t ) R(t) R ( t ) would give q ⋅ ( R ( t ) q ) > 0 q\cdot(R(t)q)>0 q ⋅ ( R ( t ) q ) > 0 ; hence q = 0 q=0 q = 0 almost surely, i.e. α t = Γ ( t ) X ^ t \alpha_t=\Gamma(t)\widehat X_t α t = Γ ( t ) X t componentwise almost surely.
Claim 3. By claim 2 of Existence and Self-Consistency of the Closed-Loop Feedback Control , α ∗ \alpha^{*} α ∗ is admissible and α t ∗ = Γ ( t ) X ^ t ( α ∗ ) \alpha^{*}_t=\Gamma(t)\widehat X_t(\alpha^{*}) α t ∗ = Γ ( t ) X t ( α ∗ ) almost surely for every t t t ; by claim 3 there it is determined by its own observations in the stated sense. By the equality case of claim 2, J [ α ∗ ] = V ∗ J[\alpha^{*}]=V^{*} J [ α ∗ ] = V ∗ . □ \square □