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 ...
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