Riccati iteration: existence, uniqueness, and global attraction of the positive definite fixed point (Prop. 4.4.1, parts 1–3)
ProvedBertsekasDP.riccati_psd_fixed_point_exists_unique_convergentcontrol-theorylinear-quadratic-regulatormatrix-analysisriccati-equation
Let , , with , and . Assume that is controllable and is observable in the sense of Definition 4.1.1 of Bertsekas, and let
be the discrete-time Riccati operator. Then there exists a positive definite matrix such that:
- solves the algebraic Riccati equation ;
- is the unique positive semidefinite fixed point: every with equals ;
- the Riccati iteration converges from every positive semidefinite start: as for every .
Proof sketch (Bertsekas, Vol. I, \S4.4): is order-preserving on the symmetric positive semidefinite cone; controllability gives a uniform upper bound for the iterates , which are nondecreasing since , hence they converge to some , and continuity of gives . Observability shows (an with forces for all , contradicting observability). Minimality of among positive semidefinite fixed points follows from , and the attraction property from comparing the iterates from with those from and from any fixed point.
Preamble
import Mathlib import Definitions.Def_BertsekasRiccatiMap open Matrix
Formal statement
namespace BertsekasDP
theorem riccati_psd_fixed_point_exists_unique_convergent {n m q : ℕ}
(A : Matrix (Fin n) (Fin n) ℝ) (B : Matrix (Fin n) (Fin m) ℝ)
(Q : Matrix (Fin n) (Fin n) ℝ) (R : Matrix (Fin m) (Fin m) ℝ)
(C : Matrix (Fin q) (Fin n) ℝ)
(hQ : Q = Cᵀ * C) (hQpsd : Q.PosSemidef) (hR : R.PosDef)
(hctrb : BertsekasControllablePair A B)
(hobs : BertsekasObservablePair A C) :
∃ P : Matrix (Fin n) (Fin n) ℝ, P.PosDef ∧
BertsekasRiccatiMap A B Q R P = P ∧
(∀ P' : Matrix (Fin n) (Fin n) ℝ, P'.PosSemidef →
BertsekasRiccatiMap A B Q R P' = P' → P' = P) ∧
(∀ P₀ : Matrix (Fin n) (Fin n) ℝ, P₀.PosSemidef →
Filter.Tendsto (fun k => (BertsekasRiccatiMap A B Q R)^[k] P₀)
Filter.atTop (nhds P)) := by sorry
end BertsekasDPSource
D. P. Bertsekas, Dynamic Programming and Optimal Control, Vol. I, Athena Scientific, 3rd ed., 2005, Section 4.4, Proposition 4.4.1 (parts 1–3); cf. J. C. Willems, Least squares stationary optimal control and the algebraic Riccati equation, IEEE Trans. Automat. Control 16 (1971), no. 6, 621–634.