All products are matrix products, (⋅)⊤ is the transpose, and we use freely the transpose rules of claim 3 of Componentwise Estimates, Transpose Identities, and Indefinite Riemann Integrals, the involutivity (S⊤)⊤=S (immediate from the definition of the transpose), and associativity (Associativity of the Matrix Product).
Step 1: reversed data. For s∈[0,T] define
Aˇ(s):=A(T−s)⊤−V(T−s)R(T−s)−1B(T−s)⊤,Dˇ(s):=B(T−s)R(T−s)−1B(T−s)⊤,Cˇ(s):=Q(T−s)−V(T−s)R(T−s)−1V(T−s)⊤.
The map s↦T−s is continuous from [0,T] to [0,T] (Sums and Products of Continuous Real-Valued Functions), t↦R(t)−1 has continuous entries by Invertibility of Symmetric Positive Definite Matrices and claim 1 of Continuity of the Inverse of a Continuous Matrix Function, and precomposition with a continuous map preserves continuity (Composition of Continuous Euclidean Maps); hence all entries of Aˇ,Dˇ,Cˇ are continuous, by Sums and Products of Continuous Real-Valued Functions.
Each Dˇ(s) is symmetric positive semidefinite: symmetry follows from the transpose rules and the symmetry of R−1 (Invertibility of Symmetric Positive Definite Matrices), and for x∈Rl, with y:=B(T−s)⊤x, the transpose rule gives x⋅(Dˇ(s)x)=y⋅(R(T−s)−1y)≥0, since R(T−s)−1 is positive definite (Invertibility of Symmetric Positive Definite Matrices), positive semidefiniteness including the value 0 at y=0. Each Cˇ(s) is symmetric positive semidefinite by hypothesis (evaluated at T−s; symmetry again by the transpose rules and symmetry of Q and R−1), and F is symmetric positive semidefinite.
Step 2: an algebraic identity. For every symmetric real l×l matrix S and every r∈[0,T]:
Aˇ(T−r)S+SAˇ(T−r)⊤−SDˇ(T−r)S+Cˇ(T−r)=A(r)⊤S+SA(r)−(SB(r)+V(r))R(r)−1(SB(r)+V(r))⊤+Q(r).(†)
Indeed, (SB+V)⊤=B⊤S⊤+V⊤=B⊤S+V⊤ by symmetry of S, so
(SB+V)R−1(SB+V)⊤=SBR−1B⊤S+SBR−1V⊤+VR−1B⊤S+VR−1V⊤,
and the right-hand side of (†) equals A⊤S+SA−SBR−1B⊤S−SBR−1V⊤−VR−1B⊤S−VR−1V⊤+Q (all at r), which regroups exactly as Aˇ(T−r)S+SAˇ(T−r)⊤−SDˇ(T−r)S+Cˇ(T−r), since Aˇ(T−r)=A(r)⊤−V(r)R(r)−1B(r)⊤, Aˇ(T−r)⊤=A(r)−B(r)R(r)−1V(r)⊤ (symmetry of R−1), Dˇ(T−r)=B(r)R(r)−1B(r)⊤, and Cˇ(T−r)=Q(r)−V(r)R(r)−1V(r)⊤.
Step 3: existence. By Global Existence and Uniqueness for the Kalman Covariance Riccati Equation applied on [0,T] with data Aˇ, Cˇ, Dˇ, and initial matrix F, there is exactly one assignment P with continuous entries such that
P(s)=F+∫0s(Aˇ(r)P(r)+P(r)Aˇ(r)⊤−P(r)Dˇ(r)P(r)+Cˇ(r))dr(0≤s≤T),
and every P(s) is symmetric positive semidefinite. Define Z(t):=P(T−t); its entries are continuous (composition with the continuous map t↦T−t) and every Z(t) is symmetric positive semidefinite. Let g(r) denote the integrand of the forward equation, with continuous entries, so that Z(t)=P(T−t)=F+∫0T−tg(ρ)dρ. By claim 2 of Product Rule and Reflection for Indefinite Riemann Integrals on [0,T] applied entrywise to the continuous function r↦g(T−r),
∫tTg(T−r)dr=∫0T−tg(T−(T−ρ))dρ=∫0T−tg(ρ)dρ,
so Z(t)=F+∫tTg(T−r)dr. By (†) with the symmetric matrix S=P(T−r)=Z(r),
g(T−r)=A(r)⊤Z(r)+Z(r)A(r)−(Z(r)B(r)+V(r))R(r)−1(Z(r)B(r)+V(r))⊤+Q(r),
so Z satisfies the backward Riccati equation.
Step 4: every continuous solution is symmetric. Let Z′ be any assignment with continuous entries satisfying the backward equation. Its integrand g′(r):=A(r)⊤Z′(r)+Z′(r)A(r)−(Z′(r)B(r)+V(r))R(r)−1(Z′(r)B(r)+V(r))⊤+Q(r) has continuous entries. Transposing the backward equation entrywise (transposition commutes with entrywise integration, being a relabeling of entries) and using F⊤=F, Q⊤=Q, the transpose rules, and the fact that MR−1M⊤ is symmetric for any l×k matrix M (symmetry of R−1), we find that Z′⊤ satisfies
Z′⊤(t)=F+∫tT(A(r)⊤Z′⊤(r)+Z′⊤(r)A(r)−(Z′(r)B(r)+V(r))R(r)−1(Z′(r)B(r)+V(r))⊤+Q(r))dr,
with the same quadratic term. Subtracting, the difference Δ:=Z′−Z′⊤ satisfies the homogeneous linear equation Δ(t)=∫tT(A(r)⊤Δ(r)+Δ(r)A(r))dr. Reflect: by claim 2 of Product Rule and Reflection for Indefinite Riemann Integrals, Δˇ(s):=Δ(T−s) satisfies Δˇ(s)=∫0s(A(T−ρ)⊤Δˇ(ρ)+Δˇ(ρ)A(T−ρ))dρ. Set δ(s):=∑i,j∣Δˇij(s)∣, a continuous function. The entries of A are bounded on [0,T] by Extreme Value Theorem on a Compact Interval; by the product entry bound (claim 2 of Componentwise Estimates, Transpose Identities, and Indefinite Riemann Integrals) and the triangle inequality for Riemann integrals of continuous functions (claims 4-5 of the same lemma give the needed splitting and bounding of entrywise integrals), there is a real C≥0 with δ(s)≤C∫0sδ(ρ)dρ for all s. By Gronwall's Lemma (Integral Form) with zero constant term, δ≡0, i.e. Z′(t) is symmetric for every t.
Step 5: uniqueness. Let Z′ be as in Step 4, so every Z′(t) is symmetric, and define P′(s):=Z′(T−s), with continuous entries. Evaluating the backward equation at T−s and applying claim 2 of Product Rule and Reflection for Indefinite Riemann Integrals to the continuous g′ itself, ∫T−sTg′(r)dr=∫0sg′(T−ρ)dρ, whence
P′(s)=Z′(T−s)=F+∫0sg′(T−ρ)dρ.
By (†) with the symmetric matrix S=Z′(T−ρ)=P′(ρ), read from right to left, g′(T−ρ)=Aˇ(ρ)P′(ρ)+P′(ρ)Aˇ(ρ)⊤−P′(ρ)Dˇ(ρ)P′(ρ)+Cˇ(ρ), so P′ satisfies the forward equation of Step 3. By the uniqueness assertion of Global Existence and Uniqueness for the Kalman Covariance Riccati Equation, P′=P, hence Z′=Z. □