25.06.2013 Views

Observateur à grand gain pour des systèmes non linéaires avec ...

Observateur à grand gain pour des systèmes non linéaires avec ...

Observateur à grand gain pour des systèmes non linéaires avec ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Prépublication n° 40 | Fascicule n° 2<br />

<strong>Observateur</strong> <strong>à</strong> <strong>grand</strong> <strong>gain</strong> <strong>pour</strong> <strong>des</strong><br />

<strong>systèmes</strong> <strong>non</strong> <strong>linéaires</strong> <strong>avec</strong> couplage <strong>non</strong><br />

complètement triangulaire<br />

Fenglong Liu, Mondher Farza, Mohammed M’saad<br />

GREYC, UMR 6072 CNRS, Université de Caen Basse-Normandie, ENSICAEN<br />

6 Bd Maréchal Juin, 14050 Caen Cedex, France<br />

fliu@greyc.ensicaen.fr, mfarza@greyc.ensicaen.fr, msaad@greyc.ensicaen.fr<br />

Résumé :<br />

Un observateur de type <strong>grand</strong> <strong>gain</strong> est synthétisé <strong>pour</strong> une classe de <strong>systèmes</strong> <strong>non</strong>-<strong>linéaires</strong><br />

multi-sorties uniformément observables. Cette classe consiste en <strong>des</strong> sous-<strong>systèmes</strong> en cascade<br />

où chaque sous-système est associé <strong>à</strong> un sous-ensemble <strong>des</strong> sorties. Deux contributions princi-<br />

pales sont <strong>à</strong> souligner : la première est liée <strong>à</strong> la structure considérée qui n’est pas complètement<br />

triangulaire. Autrement dit, la dynamique de certaines variables de chaque sous-système peut<br />

dépendre de <strong>non</strong>-linéarités impliquant tous les états du système. La deuxième contribution ré-<br />

side dans le <strong>gain</strong> de l’observateur dont l’expression est donnée et dont le calibrage s’effectue<br />

<strong>à</strong> travers le choix d’un seul paramètre. Un exemple académique <strong>avec</strong> les résultats de simulation<br />

sont présentés <strong>à</strong> titre d’illustration.<br />

Mots-clés : système <strong>non</strong>-linéaire, observateur <strong>à</strong> <strong>grand</strong> <strong>gain</strong>, couplage <strong>non</strong> triangulaire.<br />

1 Introduction<br />

En dépit d’une activité de recherche intense et continue sur l’observation <strong>des</strong> <strong>systèmes</strong><br />

<strong>non</strong>-<strong>linéaires</strong> (cf. par exemple [1, 2, 3, 4, 5, 6, 7]), ce problème reste encore ouvert <strong>pour</strong> les<br />

<strong>systèmes</strong> multi-sorties.<br />

L’objectif de ce papier est de proposer un observateur exponentiel <strong>pour</strong> une classe de<br />

<strong>systèmes</strong> <strong>non</strong>-<strong>linéaires</strong> multi-sorties qui peuvent se mettre sous la forme suivante :<br />

8<br />

<<br />

:<br />

˙x = Ax + ϕ(u, x)<br />

y = Cx<br />

Fenglong Liu, Mondher Farza, Mohammed M’saad<br />

« <strong>Observateur</strong> <strong>à</strong> <strong>grand</strong> <strong>gain</strong> <strong>pour</strong> <strong>des</strong> <strong>systèmes</strong> <strong>non</strong> <strong>linéaires</strong> <strong>avec</strong> couplage <strong>non</strong> complètement triangulaire »<br />

Schedae, 2007, prépublication n° 40, (fascicule n° 2, p. 275-278).<br />

| 275<br />

(1)


276 |<br />

“<br />

où l’état x = x1 x2 . . . xq ”T<br />

∈ IR n , <strong>avec</strong> xk “<br />

= xk 1 xk 2 . . . xk ”T<br />

∈ IR<br />

λk nk ,<br />

xk i ∈ IRp qX<br />

k , i = 1, . . . , λk, k = 1, . . . , q, nk = n ; la sortie du système<br />

k=1<br />

y =<br />

“<br />

y1 y2 . . . yq<br />

”T<br />

∈ IR p <strong>avec</strong> yk ∈ IR p qX<br />

k , k = 1, . . . , q et pk = p ; A =<br />

2<br />

3<br />

k=1<br />

h<br />

diag A1 . . . Aq<br />

6<br />

i 6<br />

, Ak = 6<br />

4<br />

0<br />

.<br />

0<br />

Ipk . . .<br />

. ..<br />

0<br />

0<br />

Ipk 7<br />

h<br />

7<br />

7,<br />

C = diag<br />

7<br />

5<br />

C1 . . . Cq<br />

i<br />

, Ck =<br />

h<br />

i<br />

0 . . . 0 0<br />

Ipk 0 . . . 0 et<br />

“<br />

la fonction <strong>non</strong>-linéaire ϕ(u, x) = ϕ1 (u, x) T ϕ2 (u, x) T . . . ϕq (u, x) T<br />

”T<br />

∈ IR n ; ϕk (u, x) =<br />

“<br />

ϕk 1 (u, x)T ϕk 2 (u, x)T . . . ϕk λ (u, x)<br />

k T<br />

”T<br />

∈ IR nk où chaque fonction ϕk i (u, x) ∈ IR pk ,<br />

k = 1, . . . , q, possède la structure suivante :<br />

• <strong>pour</strong> 1 ≤ i ≤ λk − 1 :<br />

• <strong>pour</strong> i = λk :<br />

ϕ k i (u, x) = ϕk i (u; x1 , x 2 , . . . , x k−1 ; x k 1 , . . . , xk i<br />

; xk+1<br />

1<br />

, x k+2<br />

1 , . . . , x q<br />

1 ) (2)<br />

ϕ k λ k (u, x) = ϕ k λ k (u; x 1 , x 2 , . . . , x q ) (3)<br />

Dans ce travail, on montrera, sous <strong>des</strong> hypothèses appropriées, que l’observateur proposé<br />

est global et que sa convergence est exponentielle. La contribution principale de ce tra-<br />

vail réside dans le fait que la classe de <strong>systèmes</strong> considérée renferme toutes les classes<br />

de <strong>systèmes</strong> <strong>pour</strong> lesquelles un observateur <strong>à</strong> <strong>grand</strong> <strong>gain</strong> a été proposé et <strong>pour</strong> lesquelles<br />

l’expression du <strong>gain</strong> a été donnée [1, 4, 8, 6].<br />

2 Synthèse de l’observateur<br />

Comme dans tous les travaux traitant de la synthèse d’observateurs <strong>à</strong> <strong>grand</strong> <strong>gain</strong> [1, 4, 8,<br />

6, 2], nous avons besoin de l’hypothèse suivante :<br />

Hypothèse 1 ϕ(u, x) est une fonction globalement Lipschitzienne en x, uniformément en u.<br />

Avant de donner les équations de l’observateur, nous allons introduire quelques notations<br />

et résultats préliminaires :<br />

• Soit ∆k(θ) une matrice diagonale définie par :<br />

»<br />

∆k(θ) = diag<br />

Ip k ,<br />

1<br />

θδ Ipk , . . .<br />

k<br />

1<br />

θ δ k(λ k−1) Ip k<br />

où θ > 0 est un réel et les δk sont <strong>des</strong> réels positifs définis comme suit :<br />

8<br />

><<br />

>:<br />

δk =<br />

δq = 1<br />

qY<br />

i=k+1<br />

(λi − 1) <strong>pour</strong> 1 ≤ k ≤ q − 1;<br />

• Soit S θ δ k la solution unique de l’équation algébrique de Lyapunov suivante :<br />

Schedae, 2007, prépublication n° 40, (fascicule n° 2, p. 275-278).<br />

θ δ k S θ δ k + A T k S θ δ k + S θ δ k Ak = C T k Ck<br />

–<br />

(4)<br />

(5)<br />

(6)


où les Ak et Ck sont définies dans le système (1). Il a été démontré (cf. par exemple [4]) que<br />

la solution de (6) est SDP (Symétrique Définie Positive) <strong>pour</strong> θ > 0 et que l’on a<br />

Sθδk (i, j) = (−1)(i+j) C j−1<br />

i+j−2<br />

θδk(i+j−1) De plus, on a :<br />

Ip k <strong>pour</strong> 1 ≤ i, j ≤ nk, où C p n =<br />

n!<br />

(n − p)!p!<br />

Sθδk = 1<br />

θδ ∆k(θ)S1k∆k(θ) (8)<br />

k<br />

où S1k = S θ δ k |θ=1. En particulier, S −1<br />

θ δ k CT k = (θδ k C 1 n k Ip k , θ 2δ k C 2 n k Ip k , . . . , θ n k δ k C n k<br />

nk Ip k ) T<br />

Soit maintenant le système dynamique suivant :<br />

˙ˆx = Aˆx + ϕ(u, ˆx) − S −1<br />

Θ CT (C ˆx − y) (9)<br />

où<br />

“<br />

• ˆx = ˆx 1 ˆx 2 . . . ˆx q<br />

”T<br />

∈ IR n , ˆx k “<br />

= ˆx k 1 ˆx k 2 . . . ˆx k ”T<br />

∈ IR<br />

λk nk ,ˆx k<br />

i ∈<br />

IR p qX<br />

k , i = 1, . . . , λk, k = 1, . . . , q, nk = n.<br />

•<br />

k=1<br />

ˆx k 1 = xk 1 <strong>pour</strong> k = 1, . . . , q (injection de la sortie).<br />

• u et y sont respectivement les entrées et les sorties du système (1).<br />

h<br />

i<br />

• SΘ = diag Sθδ1 . . . S δq θ<br />

Sθδk sont données par (8) <strong>pour</strong> k = 1, . . . , q.<br />

Nous é<strong>non</strong>çons le théorème suivant :<br />

est une matrice diagonale en blocs et les matrices<br />

Théorème 1 Supposons que système (1) satisfait l’hypothèse (1), alors :<br />

∃θ0 > 0; ∀θ ≥ θ0; ∃λθ > 0; ∃µθ > 0,<br />

ˆx(t) − x(t) ≤ λθe − µ θ t<br />

2 ˆx(0) − x(0)<br />

<strong>pour</strong> toute entrée bornée. De plus, lim<br />

θ→∞ µθ = +∞. Autrement dit, le système (9) est un<br />

observateur exponentiel <strong>pour</strong> le système (1) <strong>pour</strong> <strong>des</strong> entrées bornées.<br />

3 Exemple<br />

Soit le système dynamique suivant,<br />

8<br />

><<br />

>:<br />

˙z1 = z3 − 0.01u1tanh(z1z2z5) − z1<br />

˙z2 = z4 − 0.01u2tanh(z1z2z5) − z2<br />

˙z3 = −z3 3 −<br />

u1z6<br />

1 + (z1z2z3z4z5z6z7) 2<br />

u2z7<br />

˙z4 = −z4 +<br />

2 + tanh(z1z2z3z4z5z6z7)<br />

˙z5 = z6 + 0.01z1z2z3z4z5u1<br />

˙z6 = z7 − z6 + (1 − z2 1 )z2z5<br />

u1 tanh(u2)<br />

+<br />

+<br />

1 + (z3z4) 2 1 + z2 ˙z7 =<br />

6<br />

−z7 + (1 − z2) 2z1z5 + tanh(u1) u2<br />

+<br />

1 + (z6z7) 2 1 + (z3z4) 2<br />

y =<br />

“<br />

”T<br />

z1 z2 z5<br />

Il est clair que le système (10) est sous la forme (1). En considérant les valeurs initiales<br />

0<br />

zi(0) = 1, i = 1, . . . , 7 et les entrées @ u1<br />

1 0<br />

A = @<br />

u2<br />

sin(0.3t)<br />

1<br />

A, les trajectoires du système<br />

cos(0.3t)<br />

(10) sont bornées et l’hypothèse (1) est donc satisfaite (les plans de phase ne sont pas<br />

reproduits par manque de place). En conséquence, un observateur de type (9) permettant<br />

| 277<br />

(7)<br />

(10)<br />

Schedae, 2007, prépublication n° 40, (fascicule n° 2, p. 275-278).


278 |<br />

d’estimer les états <strong>non</strong> mesurés du système (10) peut être synthétisé. Pour la raison de<br />

manque <strong>des</strong> espaces, nous ne montrons pas les équations d’observateur.<br />

La simulation a été effectuée <strong>avec</strong> les conditions initiales zi(0) = 1, et ˆzi(0) = −1, i =<br />

1, . . . , 7. La valeur de paramètre θ utilisée au cours de la simulation est égale <strong>à</strong> √ 5. Les ré-<br />

sultats obtenus (cf. figures 3 montrent bien la convergence rapide <strong>des</strong> trajectoires estimées<br />

vers leurs vraies trajectoires issues de la simulation du modèle.<br />

2.5<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

−0.5<br />

−1<br />

−1.5<br />

0 5 10 15<br />

4<br />

3<br />

2<br />

1<br />

0<br />

−1<br />

−2<br />

−3<br />

Estimation<br />

Simulation<br />

−4<br />

0 5 10 15<br />

4 Conclusion<br />

Estimation<br />

z 3<br />

z 6<br />

Simulation<br />

3<br />

2.5<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

−0.5<br />

−1<br />

Estimation<br />

Simulation<br />

−1.5<br />

0 5 10 15<br />

4<br />

2<br />

0<br />

−2<br />

−4<br />

−6<br />

Estimation<br />

−8<br />

0 5 10 15<br />

z 7<br />

Simulation<br />

Fig. 1: Comparaison <strong>des</strong> états estimés et simulés <strong>pour</strong> zi,i = 3, 4, 6, 7.<br />

Dans cette contribution, nous avons proposé un observateur de type <strong>grand</strong> <strong>gain</strong> <strong>pour</strong> une<br />

classe de <strong>systèmes</strong> <strong>non</strong>-<strong>linéaires</strong> multi-sorties uniformément observables. Cette classe in-<br />

clut toutes les classes de <strong>systèmes</strong> considérées jusqu’ici et <strong>pour</strong> lesquelles le <strong>gain</strong> de l’ob-<br />

servateur est explicitement donné. Toutefois, beaucoup d’autres <strong>systèmes</strong> uniformément<br />

observables n’appartiennent pas <strong>à</strong> cette classe et nous travaillons <strong>pour</strong> élargir la synthèse<br />

de l’observateur <strong>à</strong> de tels <strong>systèmes</strong>.<br />

Références<br />

[1] G. Bornard and H. Hammouri. A high <strong>gain</strong> observer for a class of uniformly observable systems. In<br />

Proc. 30th IEEE Conference on Decision and Control, volume 122, Brighton, England, 1991.<br />

[2] M. Farza, M. M’Saad, and L. Rossignol. Observer <strong>des</strong>ign based on triangular form generated by<br />

injective map. Automatica, 40 :135–143, 2004.<br />

[3] J.P. Gauthier and G. Bornard. Observability for any u(t) of a class of <strong>non</strong>linear systems. IEEE Trans.<br />

on Aut. Control, 26 :922–926, 1981.<br />

[4] J.P. Gauthier, H. Hammouri, and S. Othman. A simple observer for <strong>non</strong>linear systems - application<br />

to bioreactors. IEEE Trans. on Aut. Control, 37 :875–880, 1992.<br />

[5] J.P. Gauthier and I. A. K. Kupka. Observability and observers for <strong>non</strong>linear systems. SIAM J. Control.<br />

Optim., 32 :975–994, 1994.<br />

[6] H Hammouri and M. Farza. Nonlinear observers for locally uniformly observable systems. ESAIM J.<br />

on Control, Optimisation and Calculus of Variations, 9 :353–370, 2003.<br />

[7] M. Hou and A. C. Pugh. Observer with linear error dynamics for <strong>non</strong>linear multi-output systems.<br />

Syst. Contr. Lett., 37 :1–9, 1999.<br />

[8] G. Bornard and H. Hammouri. A graph approach to uniform observability of <strong>non</strong>linear multi output<br />

systems. In Proc. of the 41st IEEE Conference on Decision and Control, Las Vegas, Nevada, USA,<br />

December 2002.<br />

Schedae, 2007, prépublication n° 40, (fascicule n° 2, p. 275-278).<br />

z 4

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!