x = A () x + Bu W ( s, ) y = Cx rational with respect to - - PowerPoint PPT Presentation

x a x bu
SMART_READER_LITE
LIVE PREVIEW

x = A () x + Bu W ( s, ) y = Cx rational with respect to - - PowerPoint PPT Presentation

Robust passification via static output feedback - LMI results Dimitri PEAUCELLE & Alexander FRADKOV & Boris ANDRIEVSKY LAAS-CNRS - Toulouse, FRANCE IPME-RAS - St Petersburg, RUSSIA Problem statement Passification and


slide-1
SLIDE 1

Robust passification via static output feedback - LMI results

Dimitri PEAUCELLE† & Alexander FRADKOV‡ & Boris ANDRIEVSKY‡

† LAAS-CNRS - Toulouse, FRANCE ‡ IPME-RAS - St Petersburg, RUSSIA

slide-2
SLIDE 2

Problem statement

Passification and Passivity-based techniques :

➞ linear and nonlinear control ➞ simplicity and physical meaning ➞ robustness ➞ applications to adaptive control, control of partially linear composite systems,

flight control, process control... Passification of LTI systems :

➞ SISO and MIMO ➞ SOF for Strict Positive Real ⇔ hyper-minimum-phaseness ➞ Proof of robustness w.r.t. parametric uncertainty (norm-bounded) ➞ Passification of non square systems: G-passification

& 1 IFAC’05, 4-8 July 2005, Prague

slide-3
SLIDE 3

Problem statement

Let an LTI uncertain system:

W(s, ∆) ∼

    

˙ x = A(∆)x + Bu y = Cx

rational with respect to ∆

A(∆) = A + B∆∆(I − D∆∆)−1C∆

uncertain constant real or complex norm-bounded:

∆ ∆C = {∆ ∈ Cm∆×l∆ : ∆∗∆ ≤ I} , ∆ ∆R = {∆ ∈ Rm∆×l∆ : ∆T∆ ≤ I} .

And let G ∈ Cm×p be given, where B ∈ Cn×m and C ∈ Cp×n

& 2 IFAC’05, 4-8 July 2005, Prague

slide-4
SLIDE 4

Problem statement

Robust G-Hyper-Minimum-Phaseness The system is robustly G-HMP if ∀∆ ∈ ∆

∆ φ(s, ∆) = det(sI − A(∆)) det GW(s, ∆) = det

   sI − A(∆)

−B GC O

   is Hurwitz and the high-frequency gain of GW(s, ∆) is a square symmetric positive definite matrix: GCB = B∗C∗G∗ > O.

➞ Generalizes HMP to non-square systems. ➞ Robustness: infinite number of conditions to test.

& 3 IFAC’05, 4-8 July 2005, Prague

slide-5
SLIDE 5

Problem statement

Parameter-dependent and unique static output-feedback PD-SOF :u = K(∆)y + v SOF :u = Ky + v Robust G-Passive control The closed-loop system is robustly strictly G-passive if ∀∆ ∈ ∆

there exists a quadratic PD storage function V (x, ∆) = x∗H(∆)x > 0 and a scalar ρ(∆) > 0 such that

V (x(t), ∆) ≤ V (x(0), ∆) +

t

  • v(θ)∗Gy(θ) − ρ(∆)|x(θ)|2

dθ ➞ Generalizes strict passivity for non-square systems ➞ G-passification : find K that makes the closed-loop G-passive ➞ G-passification of W(s) = Passification of GW(s)

& 4 IFAC’05, 4-8 July 2005, Prague

slide-6
SLIDE 6

Problem statement

Theorem 1 : [Fradkov 1976-2003] Equivalence of

① W(s, ∆) is robustly G-HMP ② W(s, ∆) is robustly G-passifiable by PD-SOF K(∆) ③ ∃K unique that robustly G-passificates W(s, ∆)

Proof (Sketch)

G-HMP ⇒ High gain control for any ∆: K(∆) = −k(∆)G : k(∆) > 0 , sufficiently large

Well-posedness of uncertain modeling: K = − max∆∈∆

∆ k(∆)G

Outline

① LMI results for robust G-HMP analysis ③ LMI results for robust G-passifying SOF design ➞ Numerical example : cruise missile model.

& 5 IFAC’05, 4-8 July 2005, Prague

slide-7
SLIDE 7

Robust G-HMP analsyis

Theorem 1 Let the following matrices

N = (GC)⊥ , M = (NN ∗ + BB∗)−1 , ˜ A = N ∗MAN . W(s, ∆) , ∆ ∈ ∆ ∆C is robustly G-HMP

if and only if GCB > O and ∃P > O ∈ C    P ˜

A + ˜ A∗P PN ∗MB∆ B∗

∆MNP

−I

   +    N ∗C∗

D∗

      N ∗C∗

D∗

  

< O

where N = (GC)⊥ and M = (NN ∗ + BB∗)−1. In case ∆ ∈ ∆

∆R, P ∈ R ; LMI conditions are only sufficient.

Proof Robust G-HMP is reformulated as the robust Hurwitz stability of a reduced order system.

& 6 IFAC’05, 4-8 July 2005, Prague

slide-8
SLIDE 8

Robust G-passifying design

Theorem 2 W(s, ∆) is uniformly robustly strictly G-passifiable via SOF if and only if ∃H > O ∈ C, ∃K ∈ C:

HB = C∗G∗

   HA + A∗H + C∗(G∗K + K∗G)C

HB∆ B∗

∆H

−I

   +    C∗

D∗

      C∗

D∗

  

< O

Proof Classical LMI results for ’quadratic’ stability

➞ Uniform storage function V (x, ∆) = V (x) = x∗Hx.

& 7 IFAC’05, 4-8 July 2005, Prague

slide-9
SLIDE 9

Robust G-passifying design

Theorem 3 W(s, ∆) is uniformly robustly strictly G-passifiable via SOF if and only if ∃H > O ∈ C, ∃K ∈ C:

HB = C∗G∗

   HA + A∗H + C∗(G∗K + K∗G)C

HB∆ B∗

∆H

−I

   +    C∗

D∗

      C∗

D∗

  

< O

Remarks

✪ Thm 2 ⇒ Thm 1 with P = N ∗HN (conjecture : converse also holds) ✪ PB : design K and G simultaneously ?

LMI problem if ∃S such that PB = BS (conservative)

✪ Always possible to take K = −kG if feasible. ✪ Possible to add LMI constraints on K, e.g. find K with minimum norm.

& 8 IFAC’05, 4-8 July 2005, Prague

slide-10
SLIDE 10

Numerical example : cruise missile

Model definiion

➞ 4th order model of lateral dynamics for cruise missile + actuator dynamics ➞ Dynamics depend on altitude h ∈ [h ¯ h] ⊂ R+ (converted into ∆ ∈ ∆ ∆R) ➞ Measured outputs:

yaw angle ϕ(t), yaw angular rate r(t) and the rudder deflection angle δr(t)

➞ Control input: rudder servo command signal ➞ G is chosen a priori such that GCB > O.

& 9 IFAC’05, 4-8 July 2005, Prague

slide-11
SLIDE 11

Numerical example : cruise missile

Robust G-HMP analysis

✪ For h = 0 and ¯ h = 10km : feasible ✪ For h = 9.9925km and ¯ h = 10.2105km : feasible ✪ For h = 0 and ¯ h = 10.2105km : infeasible ➞ h = 10.1 + 0.5i makes system non G-HMP

.

➞ Conservatism for real-valued uncertainty.

0.98 1 1.02 1.04 !0.02 0.02 0.04 0.06 0.08 0.1 0.12

Real(!) Imaginary(!)

➞ Exists a SOF for h ∈ [0 10.2105], cannot be found with Thm 2.

& 10 IFAC’05, 4-8 July 2005, Prague

slide-12
SLIDE 12

Numerical example : cruise missile

Robust G-passifying SOF design Assume h ∈ [0 10]

✪ Thm 2 → K1 = −

  • 79.28

50.34 11.92

  • ✪ Thm 2, min K:

→ K2 = −

  • 60.75

34.47 10.67

  • ✪ Thm 2, minK=−kG k, :

→ K3 = −

  • 118.53

44.45 14.82

  • Yaw angle and rudder deflection

for control K2 and for h = 0.1 , 5 , 9

& 11 IFAC’05, 4-8 July 2005, Prague

slide-13
SLIDE 13

Conclusion ➞ Non-conservative (complex case) LMI conditions of robust strict G-passification ➞ Conservative LMI design method ✪ Design simultaneously K and G ✪ Design of robust G-passifying adaptive control u(t) = K(t)y(t)

& 12 IFAC’05, 4-8 July 2005, Prague