Observateur d'état : Estimation et Synthèse
Observateur d'état : Estimation et Synthèse
1. Observateurs
La totalité des variables d’état d’un système ne sont souvent pas accessibles à la mesure, car
dans certains cas soit il est impossible de les mesurer physiquement, soit le système deviendrait
trop cher et complexe si on y ajoute des capteurs pour la mesure de chaque état. Dans ce cas,
l’implémentation directe de la commande u Kx ne peut pas se faire directement. Il est donc
nécessaire dans ce cas de reconstruire l’état x à partir des informations disponibles sur l’entrée
u et la sortie y . On utilise pour cela un Observateur, qui est constitué d’un système dynamique
qui permet d’obtenir une estimation x̂ de x Cet observateur est souvent connu sous le nom de
reconstructeur, d’estimateur, de filtre.
Le schéma de la commande par retour d’état doit être donc modifié comme celui présenté
sur la figure 4.1.
Où u(t), y(t) sont respectivement l’entrée et la sortie du système. Remarquons que ce sont
également les deux entrées de l’observateur. La sortie de l’observateur est la variable x̂ .
Si z et x ont même dimension, l’observateur est dit complet (tout l’état est estimé). Dans ce cas
T I ; donc z x et ẑ xˆ .
Si dim( z ) dim(x ) , alors l’observateur est dit d’ordre réduit.
Un observateur doit satisfaire aux minimums les deux conditions suivantes :
1. Un observateur doit être stable.
2. Un observateur doit assurer la convergence de ẑ vers z (estimation sans biais) :
limt ( zˆ(t ) z(t )) limt e(t ) 0, u, x(t0 ) (4.2)
Où e(t ) est dénommée erreur de reconstruction.
45 | P a g e
Chapitre IV : Observateur d’état
L’observateur identité est un observateur complet et sans biais ( zˆ(t ) z, si, t ) avec z x
xˆ Fxˆ Ly Ju
x Ax Bu (4.3)
y Cx
En considérant la dérivée de l’erreur d’estimation :
e x xˆ ( Ax Bu ) ( Fxˆ Ly Ju )
e ( Ax Bu ) ( Fxˆ LCy Ju ) ( A LC ) x ( B J )u Fxˆ
Comme x̂ x e
e ( A LC F ) x ( B J )u Fe (4.4)
Pour une estimation sans biais il faut que e 0 et que e 0, x, u pour cela il faut donc
satisfaire les équations :
A LC F 0 F A LC
B J 0 J B
F , stable ( A LC ), stable
xˆ ( A LC ) xˆ Ly Bu (4.5)
xˆ Axˆ Bu L( y yˆ ) (4.7)
Cette relation, qui caractérise l’Observateur peut être représenté par le schéma fonctionnel :
46 | P a g e
Chapitre IV : Observateur d’état
Exemple 4.1 :
Soit le système linéaire x Ax Bu; y Cx défini par les matrices :
2 3 0
A B C 1 0
1 4 1
On suppose ici que seule la sortie y x1 est mesurable, il faut donc un observateur pour
déterminer l’état non mesurable x2 . Même si on a besoin de déterminer seulement un des deux
états, on va construire un observateur qui permet de déterminer tous les états, un observateur
complet. Il est possible cependant de construire un observateur réduit, mais souvent il est
convenable de construire un observateur d’ordre complet et utiliser les états ou l’information
est redondante pour corriger un éventuel effet des bruits de mesures ou autres perturbations. Un
exemple, est le Filtre de Kalman (observateur optimal à temps variant) qui permet de résoudre
le problème d’observations avec bruit de mesure.
47 | P a g e
Chapitre IV : Observateur d’état
Il faut initialement déterminer si le système est observable, car si ce n’est pas le cas il n’est pas
possible de construire un observateur.
C 1 0
rang (O) rang 2n
CA 2 3
Le système est donc observable, il est possible de construire un observateur pour déterminer les
états du système.
Pour l’observateur on va choisir une dynamique telle que le comportement de l’observateur soit
équivalent à un système de deuxième ordre avec équation caractéristique :
L1 6 16
4 L1 3L2 11 100
Donc la matrice L est :
L 22
L 1
L2 59
L’observateur du système est donc :
2 3 0 22
xˆ Axˆ Bu L( y Cxˆ )
ˆ
x 1 u 59 y 1 0 xˆ
1 4
48 | P a g e
Chapitre IV : Observateur d’état
Où les coefficients son choisis pour obtenir des spécifications désirées pour l’observateur.
La matrice de Gain de l’Observateur
L L1 Ln
T
Ln (4.9)
Peut être obtenue avec :
L p( A)O1 0 0 1
T
(4.10)
Où O est la matrice d’observabilité du système et :
p( A) An n1 An1 2 A2 1 A 0 I (4.11)
Exemple 4.2 :
Pour le même système antérieur.
2 3 0
A B C 1 0
1 4 1
Et le même choix de l’équation caractéristique désiré :
p( ) 2 16 100
On obtient :
2
2 3 2 3 1 0 133 66
p( A) A 16 A 100 I
2
16 100
1 4 1 4 0 1 22 177
49 | P a g e
Chapitre IV : Observateur d’état
commandables restent stables et que les non observables ne soient pas nécessaires pour la
commande. Le problème de la commande se résout ensuite en trois grandes étapes :
u v Kxˆ (4.12)
xˆ ( A LC ) xˆ Ly Bu (4.13)
50 | P a g e
Chapitre IV : Observateur d’état
xˆ ( A LC BK ) xˆ Ly
(4.14)
u Kxˆ
e Ax Bu ( A LC) xˆ Ly Bu (4.15)
Ce qui revient à :
Ceci montre que la dynamique de l’erreur ne dépend pas de l’entrée. D’un autre coté si on
revient au système à commander x Ax Bu, y Cx et, on substitue la loi de commande dans
l’équation dynamique on obtient :
x Ax BKxˆ (4.17)
x ( A BK ) x BKe
(4.19)
e ( A LC )e
Ou sous forme Matricielle :
x A BK BK x
e 0 A LC e
(4.20)
Les valeurs propres de la nouvelle matrice dynamique du système sont donc les racines de :
I ( A BK ) ( A LC ) 0 (4.21)
Ce qui veut dire que la dynamique du système en boucle fermée est fixée d’une part par les
valeurs propres assignées à ( A BK ) , et par les valeurs propres assignées à ( A LC ) d’autre
part. Théoriquement ces valeurs propres se fixent de manière indépendante. C’est le principe
de la séparation.
51 | P a g e
Chapitre IV : Observateur d’état
3 2 2
x(t ) x(t ) u (t )
4 5 0
y (t ) 2 1 x(t )
x(0) x0
n 2 p m 1
Cahier de charges
1- Amélioration de la dynamique.
2- Valeur finale de la sortie y 1 .
Analyse du système
Analyse de la stabilité
La matrice d’état A possède deux valeurs propres strictement négatives 1 1 2 7 . Le
système est par conséquent stable.
Analyse transitoire
La matrice d’état A possède deux valeurs propres réelles auxquelles on peut associer deux
1 1
constantes de temps 1 1 2 0.14 . En négligeant la plus petite constante de
1 2
temps, le temps de réponse peut être estimé à 31 3 .
2 6
Qc B AB
0 8
det(QC ) 16 0
rang (QC ) 2 (n)
Système commandable
C 2 1
Q0
CA 10 9
det(QO ) 8 0
rang (QO ) 2 (n)
Système observable
52 | P a g e
Chapitre IV : Observateur d’état
Equation de l’observateur
La loi de commande
u(t ) v(t ) 0.5z1 (t ) 0.25z2 (t )
53 | P a g e
Chapitre IV : Observateur d’état
Simulation 4.1 :
Les résultats de la simulation sont illustrés ci-dessous. Le système est supposé dans l’état
initial x0 1 2 alors que celui de l’observateur est z0 0 0 .
T T
54 | P a g e
Chapitre IV : Observateur d’état
Fig. 4.5 Résultats de la simulation de la commande par retour d’état avec observateur
x(t ) Ax(t ) Bu (t )
Sboucle ouverte (4.22)
y (t ) Cx(t )
55 | P a g e
Chapitre IV : Observateur d’état
x1 (t ) A11 A12 x1 (t ) B1
x (t ) A
A22 x2 (t ) B2
u (t )
2 21 (4.23)
y (t ) x1 (t )
d’où :
x2 (t ) A22 x2 (t ) u (t )
(4.25)
w(t ) A12 x2 (t )
Cette représentation fait apparaitre un système d’ordre réduit dont la matrice d’état est A22 , et
la matrice de sortie est A12 .
Pour éviter la dérivée x1 (t ) dans le second membre de la première équation, on introduit une
nouvelle variable notée s(t ) telle que :
Ces dernières équations représentent bien la structure d’un observateur excité par les signaux
disponibles, à savoir, d’entrée u (t ) et la sortie y (t ) . La sortie z (t ) est de dimension réduite par
rapport à l’observateur classique.
56 | P a g e
Chapitre IV : Observateur d’état
On montre par ailleurs que si la paire (C , A) est observable alors il en est de même pour la paire
( A12 , A22 ) , ce qui permet d’assurer l’existence de la matrice L . Celle-ci sera calculée en fixant
une dynamique appropriée à l’observateur d’ordre réduit.
Cette erreur transitoire tend asymptotiquement vers zéro quelle que soit l’erreur initiale e0 et
ce, selon une dynamique fixée par les valeurs propres attribuées à la matrice ( A22 LA12 ) .
Exemple 4.4 :
On considère le système suivant qui est naturellement partitionné :
x1 (t ) 1 1 x1 (t ) 1
x (t ) 2 2 x (t ) 1 u (t )
2 2
y (t ) 1 0 x(t )
x1 (t ) y(t ) , est accessible et donc seule x2 (t ) est à reconstruire. Selon la notation adoptée, on
a:
Il s’en suit :
57 | P a g e
Chapitre IV : Observateur d’état
Fig. 4.6 Schéma de principe pour la réalisation d’un observateur d’ordre réduit
u v Kxˆ (4.31)
xˆ ( A LC ) xˆ Ly Bu
(4.32)
y Cx
58 | P a g e
Chapitre IV : Observateur d’état
e Ax Bu ( A LC) xˆ Ly Bu (4.33)
Ce qui revient à :
x Ax Bv BKxˆ (4.35)
x ( A BK ) x BKe Bv
e ( A LC )e (4.37)
y Cx
x A BK BK x B
e 0
A LC e 0
v
(4.38)
y Cx
Si yc est constante, une fois atteint le régime permanent pour que la valeur de y tende a yc on a :
0 A BK BK x B
0 0
A LC e 0
v
(4.39)
y yc Cx
Puisque A LC 0 car la matrice L a été choisi pour que l’observateur soit stable, donc les
pôles strictement inférieurs à zéro, donc forcément e 0 (erreur d’observation nul en régime
stationnaire).
L’équation devient :
( A BK ) x Bv 0 (4.40)
Donc :
x ( A BK )1 Bv 0 (4.41)
59 | P a g e
Chapitre IV : Observateur d’état
donc finalement :
1
v C (( A BK )1 B) yc Nyc (4.43)
Identique à l’équation pour la commande par retour d’état directe. Donc la sortie y tends ver la
1
consigne yc si la matrice N est calculé avec : N C (( A BK )1 B) .
60 | P a g e
Chapitre IV : Observateur d’état
Exercices : Observateur d’état et commande par retour de sortie des systèmes linéaires
Exercice 01
Soit un oscillateur décrit par l'équation (t ) 02 (t ) u(t ) on demande :
1- De mettre le système sous sa forme d'état (matrices ( A, B, C, D) ), sachant que la variable
observée est .
2- D'effectuer la synthèse d'un contrôleur d'état de telle sorte que le système en boucle
fermée ait une fréquence naturelle deux fois supérieure et un facteur d'amortissement
égal à l'unité.
3- D'effectuer la synthèse d'un observateur d'état cinq fois plus rapide que le contrôleur
d'état.
Exercice 02
On considère un système régi par l’équation d’état :
1 5 4
x(t ) Ax(t ) Bu (t ) avec A et B
2 8 1
Exercice 03
On considère un système régi par l’équation d’état :
1 2 1
x(t ) Ax(t ) Bu (t ) avec A et B
2 2 2
Montrer que ce système n’est pas commandable. Déterminer l’expression du mode non
commandable de ce système et étudier la possibilité de commander ce système par un retour
d’état.
Exercice 04
Un navire de 100m de longueur en mouvement latérale avec une vitesse constante de 10m/s est
décrit par les équations d'état suivantes
0.0895 0.286 0 0.0145
r 0.0439 0.272 0 r 0.0122
0 1 0
0
1- Déterminer la fonction de transfert entre et , ainsi que les pôles du système en
boucle ouverte, le système est-il stable ?
61 | P a g e
Chapitre IV : Observateur d’état
Exercice 05
Soit le système linéaire représenté par sa fonction de transfert :
9
G( s)
s 92
Exercice 06
Soit le système décrit par la fonction de transfert en boucle ouverte suivante :
1
G( s)
( s 1)( s 0.7 s 2)
2
62 | P a g e
Chapitre IV : Observateur d’état
63 | P a g e