Algebraic Properties of Riccati equations. Ruth Curtain University of Groningen, The Netherlands

Transcription

1 Algebraic Properties of Riccati equations Ruth Curtain University of Groningen, The Netherlands

2 Special Recognition Peter Falb Jan Willems Tony Pritchard Hans Zwart

3 Riccati equation Ṗ(t) + A(t) P(t) + P(t)A(t) = P(t)B(t)B(t) P(t) + C(t) C(t), P(T) = Q. W.T. Reid: Riccati Differential Equations, Academic Press, 1972.

4 Algebraic Riccati equation Let A, B, C be complex matrices. If (A, B) is stabilizable and (A, C) is detectable, then A P + PA PBB P + C C = 0 has a unique stabilizing solution P. Stabilizing solution: P = P 0 and A BB P is stable. (A, B) is stabilizable if F s.t. A + BF is stable. (A, C) is detectable if L s.t. A + LC is stable.

5 Algebraic question Let A, B, C A, a Banach algebra with involution. When will A P + PA PBB P + C C = 0 have a unique stabilizing solution P A n n? Stabilizing solution: P = P and A BB P is exponentially stable, i.e., the spectrum of A is contained in the open left half-plane.

6 Hilbert space Riccati eqns Suppose that H is a Hilbert space and A, B, C L(H). If (A, B) is exponentially stabilizable and (A, C) is exponentially detectable, then a unique exponentially stabilizing solution P L(H), P = P 0 of A P + PA PBB P + C C = 0. A Positive result: L(H) is a Banach algebra with the involution the adjoint operation.

7 Gelfand Representation Theorem for commutative Banach algebras Every commutative semi-simple Banach algebra is isomorphic to an algebra of continuous functions on its maximal ideal space M(A) (a compact Hausdorff space, equipped with the weak * topology). Denote the Gelfand transform of a A by â. Then â(ϕ) = ϕ(a), ϕ M(A), a A.

8 Example (Even-Weighted Wiener algebras on unit circle T) Let α = (α n ) n Z be any sequence of positive real numbers satisfying α n+m α n α m, α n = α n (even). An even-weighted Wiener algebra is { W α (T) := f : f (φ) = f n φ n (φ T) and } α n f n < +, n Z n Z with norm f Wα(T) := n Z α n f n. This is a commutative semi-simple Banach algebra and the maximal ideal space is isomorphic to the annulus around T: A(ρ) = {φ C : 1/ρ φ ρ}, ρ = inf n>0 n αn = lim n n αn. Gelfand transform: ˆf (ϕ) = f (ϕ) for ϕ A(ρ).

9 Chris Byrnes, CDC 1980 Let A be a commutative, unital, complex, semisimple Banach algebra, with an involution. Let A A n n, B A n m, C A p n be such that for all ϕ M(A) (the maximal ideal space) (Â(ϕ), B(ϕ)) is controllable and (Â(ϕ), Ĉ(ϕ)) is observable. Then there exists a solution P A n n such that PA + A P PBB P + C C = 0, and A BB P is stable.

10 Example (Counterexample to Chris Byrnes claim) Let â(ϕ) = 2 + ϕ, ˆb(ϕ) = 1, ĉ(ϕ) = 1 W α (T) for arbitrary even weights α n = α n. Now (â(ϕ), ˆb(ϕ)) is controllable and and (â(ϕ), ĉ(ϕ)) is observable for all ϕ C. According to Byrnes we should have a unique solution to the Riccati equation P W α (T) for for arbitrary even weights α k. Now ( 4 + ϕ + 1/ϕ + ) (4 + ϕ + 1/ϕ) P(ϕ) = ˆP(ϕ) = ( 5+1 Note: Singularity at ϕ = 1 ± + i 1 ± 2 2 But elements of W α (T) are analytic in the interior of A(ρ), where ρ = lim n n αn. Thus P / W α (T) for for arbitrary symmetric weights α k. ).

11 Correct result with Amol Sasane, SIAM 2011 Let A be a commutative, unital, complex, semisimple Banach algebra. For all ϕ M(A) let A A n n, B A n m, C A p n satisfy: (A )(ϕ) = (Â(ϕ)) (BB )(ϕ) = B(ϕ)( B(ϕ)) (C C)(ϕ) = (Ĉ(ϕ)) Ĉ(ϕ) (Â(ϕ), B(ϕ)) is stabilizable (Â(ϕ), Ĉ(ϕ)) is detectable. Then there exists a P A n n such that PBB P PA A P C C = 0, A BB P is exponentially stable, and P = P.

12 Essential, but restrictive condition on the involution: Gelfand transform of involution must match complex conjugation: (A )(ϕ) = (Â(ϕ)). (1) Symmetric Banach algebras always satisfy (1). Two possible involutions for W α (T) are f (φ) := f (φ), f (φ) := f (1/φ). In our counter example neither satisfies (1): (2 + z) = 2+z, (A + z) = 2 + z, (2 + z) = 1+ 1 z.

13 Example ( Symmetric even-weighted Wiener algebras W α (T)) The Gelfand-Raikov-Shilov condition: ρ = 1 Involution reduces to = A(ρ) = T. f (φ) := f (φ) (φ T), f W α (T). With this involution W α (T) is a symmetric Banach algebra and (1) is always satisfied. Exponential weights: ρ > 1. Subexponential weights: ρ = 1, where α n = e α n β, α > 0, 0 β < 1.

14 Corollary, Curtain & Sasane, SIAM 2011 Let A be a commutative, unital, complex, symmetric Banach algebra. For all ϕ M(A) let A A n n, B A n m, C A p n satisfy: (Â(ϕ), B(ϕ)) is stabilizable (Â(ϕ), Ĉ(ϕ)) is detectable. Then there exists a P A n n such that PBB P PA A P C C = 0, A BB P is exponentially stable, and P = P. Example : A, B, C are matrices with entries from W α (T) where α n are subexponential weights.

15 Platoons of vehicles

16 Spatially Invariant Systems (Bamieh 2002) Σ(A, B, C) with A, B, C matrices whose entries are bounded convolution operators on the Hilbert space l 2 (Z; C n ) = {x = (x r ) r Z : x r 2 < }. r=z Example ż r (t) = z r (t) + u r (t) + 2u r 1 (t), y r (t) = z r (t) z r 1 (t), r.

17 L(l 2 (Z; C)) is a Banach algebra with the adjoint operation as involution. Convolution operators (Tx) l = r Z T r lx r form a subalgebra: Under the Fourier transform F : l 2 (Z) L 2 (T) and spatially invariant systems Σ(A, B, C) are isometrically isomorphic to systems on L 2 (T; C n ): Σ(FAF 1, FBF 1, FCF 1 ) := Σ(Â, ˆB, Ĉ). Since Â, ˆB, Ĉ are multiplication operators on L (T; C n ) they are much easier to handle mathematically.

18 Example (Fourier transform of example) ż r (t) = z r (t) + u r (t) + 2u r 1 (t), y r (t) = z r (t) z r 1 (t), r. ẑ(φ) = (10 φ 1/φ)û(φ), ŷ(φ) = ẑ(φ), φ T. The Riccati equation can be solved pointwise: δ = Â(φ) ˆP(φ) + ˆP(φ)Â(φ) + Ĉ(φ) Ĉ(φ) = ˆP(φ)ˆB(φ)ˆB(φ) ˆP(φ), φ T. 1 ˆP(φ) = 10 φ 1/φ = k Z δ k φ k,

19 Control of Platoon-type spatially invariant systems Spatially invariant operators are isomorphic to L(L 2 (T; C n )) = L (T; C n n ). Existence of Riccati solutions If (Â(φ), ˆB(φ), Ĉ(φ)) is stabilizable and detectable for all φ T the Riccati equation has a unique stabilizing solution P L (T; C n n ). Even for simple examples ˆP(φ) = r Z p rφ r. But for an implementable control law you need to truncate and the truncation should approximate ˆP, i.e., r T α rp r < for some (α r ), i.e., ˆP must be in a Wiener algebra W α (T) L (T).

20 Spatially distributed systems Introduced in Bamieh et al IEEE AC Trans 2002 and further studied in Motee & Jadbabaie in IEEE AC Trans Unlike spatially invariant systems, spatially distributed systems allow for the interaction of an array of infinitely many distinct linear systems. The requirement of implementable control laws corresponds to the algebraic properties of the Riccati equation for noncommutative Banach algebras. Claim in Motee & Jadbabaie in IEEE AC Trans was false. Bunce 1985 had proven a positive result for C -algebras, but this does not cover the spatially distributed case.

21 Algebraic properties for noncommutative Banach algebras (Curtain 2011) Suppose A is a unital symmetric Banach algebra and A is a Banach *-subalgebra of L(Z), where Z is a Hilbert space; A has the inverse-closed property: D A, D 1 L(Z) = D 1 A; A A M A L(Z). If A, B, C A and (A, B, C) is exponentially stabilizable and detectable wrt Z, then P A, where P L(Z) is the unique nonnegative solution to the control Riccati equation: A P + PA PBB P + C C = 0.

22 Applications Using results by Gröchenig & Leinert (2006) it can be shown that a subclass of the spatially distributed systems studied in Motee & Jadbabaie in 2008 satisfies the conditions of the previous theorem which leads the way to designing implementable control laws. Again, only subexponential weights are allowed and not exponential ones.

23 Conclusions For commutative Banach algebras we have an elegant result for the LQR Riccati equation; in particular, for symmetric algebras. Application to platoon-type spatially invariant systems: design of implementable control laws. There are generalizations to other Riccati equations, including H type equations. Algebraic properties of the LQR Riccati equation for inverse-closed noncommutative algebras have direct applications to spatially distributed systems. Engineers provide motivation for interesting mathematical problems.

24 Open: spatially invariant systems with unbounded A Linear p.d.e. s on infinite domains are examples of spatially invariant systems with unbounded A. Under the Fourier transform these corresponds to systems with multiplicative operators on L (ir). Semigroup properties have been well studied by the Nagel school. An algebraic property of the Riccati equation is shown in Curtain System theoretic properties need closer investigation.