Modèle d'état

Représentation d'état

En automatique, une représentation d'état permet de modéliser un système dynamique sous forme matricielle en utilisant des variables d'état. On se place alors dans un espace d'état. Cette représentation, qui peut être linéaire ou non-linéaire, doit rendre compte de l'état du système à n'importe quel instant futur si l'on possède les valeurs initiales. Cette représentation peut être continue ou discrète.

Sommaire

1 Variables d'état
2 Systèmes linéaires
3 Systèmes non linéaires
- 3.1 Multiplicité de la représentation d'état
- 3.2 Exemple du pendule
4 Voir aussi

Variables d'état

Un système peut être entièrement décrit à l'aide d'un ensemble de variables minimal. Les variables d'état sont des grandeurs physiques continues du système (elles doivent être dérivables) et doivent être indépendantes les unes des autres. Elles sont généralement rassemblée dans un vecteur X. La connaissance de toutes les variables d'état à un instant t doit permettre de connaître toutes les valeurs du système à un instant t+dt. Attention, cette représentation n'est pas unique, un même système peut être décrit avec des variables d'état différentes mais leur nombre est toujours le même. Ce nombre, désigné par la lettre n, représente l'ordre du système.

Systèmes linéaires

La représentation d'état la plus générale pour les systèmes linéaires est la suivante en continue:

$\dot X(t) = A(t) \cdot X(t) + B(t) \cdot U(t)$ $Y(t) = C(t) \cdot X(t) + D(t) \cdot U(t)$

$X \in \mathbb{R}^n$ : vecteur qui représente les n variables d'état
$Y \in \mathbb{R}^q$ : vecteur qui représente les q mesures
$U \in \mathbb{R}^p$ : vecteur qui représente les p commandes
$A \in \mathbb{R}^{n\times n}$ : Matrice de dynamique
$B \in \mathbb{R}^{n\times p}$ : Matrice de commande
$C \in \mathbb{R}^{q\times n}$ : Matrice d'observation
$D \in \mathbb{R}^{q\times p}$ : Matrice d'action directe

Ceci est le cas le plus général. Les matrices A, B, C, D sont souvent invariantes selon le temps, elles deviennent alors des matrices constantes et on parle de représentation d'état continue indépendante du temps.

Cette représentation d'état se représente sous forme de schéma-bloc :

Les valeurs propres de la matrices d'état A représentent les pôles du système si ces valeurs propres ne représentent pas des modes cachés du système c’est-à-dire si les valeurs propres sont observables ET gouvernables/commandables. Si ces pôles sont à partie réelles négatives, alors le système est asymptotiquement stable.

L'état du système est un résumé exhaustif du passé du système. En effet, connaissant l'état $x (t i)$ , et la commande $u (t)$ sur l'intervalle fini $[t i; t f]$ , on peut exprimer l'état $x (t f)$ , grâce aux calculs sur l'exponentielle d'une matrice :

$x(t_f) = x(t_i)e^{A(t_f-t_i)} + \int_{t_i}^{t_f}e^{A(t_f-\tau)}B(\tau)u(\tau)d\tau$

Multiplicité de la représentation d'état

Un système donné admet une infinité de formes d'état strictement équivalentes. En effet, considérons le changement de variable $z = P x$ , avec $P \in \mathbb R^{n \times n}$ et inversible. Le système se réécrit donc comme suit :

$\begin{cases} \dot z = PAx + PBu & = PAP^{-1}z + PBu \\ y = Cx + Du & = CP^{-1}z + Du\end{cases}$

Commandabilité et Observabilité

On ne détaillera ici que le cas des systèmes linéaires invariant (SLI). La commandabilité et l'observabilité dans les cas instationnaires et/ou non linéaires sont beaucoup plus difficiles à traiter, et il existe des livres spécialisés dans ce domaine. La commandabilité et l'observabilité sont des propriétés structurelles du système qui n'apparaissent pas dans la représentation par fonction de transfert.

Commandabilité

Un système est dit commandable si quel que soit l'état à l'instant initial $x (t i)$ , il existe une commande $u (t)$ , appliquée sur un intervalle de temps fini $[t i; t f]$ , telle que $x (t f) = 0$ .

Un critère algébrique existe pour la commandabilité, il est dû à Rudolf Kalman et Richard Bucy (1961). Un système est dit commandable si et seulement si :

$rang\begin{bmatrix}B& AB& ...& A^{n-1}B\end{bmatrix} = n$

Cette matrice est communément appelée matrice de commandabilité et ses colonnes se calculent de façon itérative : $A k + 1 B = A * A k B$

(le rang d'une matrice est le nombre maximal de vecteurs lignes (ou colonnes) linéairement indépendants)

Observabilité

Un système est dit observable si l'observation de ses entrées et sorties pendant un intervalle de temps fini $[t i; t f]$ permet de retrouver l'état initial $x (t i)$ . En fait, puisqu'il est possible pour les SLI d'avoir une solution analytique, l'observabilité nous permet donc de connaître l'état $x (t)$ à tout instant compris dans l'intervalle $[t i; t f]$ . De manière duale, Kalman et Bucy ont proposé un critère algébrique pour la commandabilité. Une représentation d'état continue indépendante du temps est observable si et seulement si :

$rang\begin{bmatrix}C\\ CA\\ ...\\ CA^{n-1}\end{bmatrix} = n$

Cette matrice est communément appelée matrice d'observabilité et ses lignes se calculent de façon itérative : $C A k + 1 = C A k * A$

Remarques

Un système à la fois commandable et observable est dit minimal.
Les propriétés de commandabilité et d'observabilité sont dites duales.
Les propriétés de commandabilité et d'observabilité sont des propriétés structurelles fortes sur le système. Il existe également des extensions plus faibles.
Voir : Commandabilité et Observabilité.

Cas des systèmes sous forme de fonction de transfert

Comme dit ci-dessus, les propriétés de commandabilité et d'observabilité ne sont pas apparentes sur les représentations par fonction de transfert. Voici un exemple illustratif.

Considérons les systèmes $S 1$ et $S 2$ correspondants aux équations différentielles suivantes :

$S_1 : \dot y + y = u$ $S_2 : \ddot y - y = \dot u - u$

Le système $S 1$ possède la forme d'état triviale suivante :

$S_1 : \begin{cases} \dot x = - x + u \\ y = x \end{cases}$

On remarque trivialement que le système est sous forme minimale.

Le système $S 2$ est légèrement plus difficile à exprimer sous forme d'état :

$S_2 : \begin{cases} \begin{bmatrix} \dot x \\ \ddot x - \dot u \end{bmatrix} = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix} \begin{bmatrix} x \\ \dot x - u \end{bmatrix} + \begin{bmatrix} 1 \\ -1 \end{bmatrix} u \\ y = \begin{bmatrix} 1 & 0 \end{bmatrix} \begin{bmatrix} x \\ \dot x - u \end{bmatrix} \end{cases}$

On peut alors calculer la matrice de commandabilité $\mathcal{C}$ :

$\begin{bmatrix} B & AB \end{bmatrix} = \begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix}$

On remarque que le rang de cette matrice est égal à 1 alors que le système est d'ordre $n = 2$ , le système n'est donc pas commandable.

Malgré tout cela, on peut noter que les deux systèmes ont la "même" représentation sous forme de fonction de transfert (avec $s$ la variable de Laplace), alors qu'ils sont fondamentalement différents :

$S_1 : \xrightarrow{\mathcal{L}} : \frac{Y(s)}{U(s)} = \frac{1}{s+1}$ $S_2 : \xrightarrow{\mathcal{L}} : \frac{Y(s)}{U(s)} = \frac{s-1}{s^2 - 1} = \frac{s-1}{(s+1)(s-1)} = \frac{1}{s+1}$

En fait, ce sont les simplifications des pôles et zéros non commandables, qui fait que la représentation sous forme de fonction de transfert n'est pas adaptée pour refléter toutes les propriétés structurelles du système.

Passage Représentation d'état-Fonction de transfert

Une représentation d'état linéaire continue indépendante du temps peut être transformée en fonction de transfert. Chacune de ces deux représentations dynamiques contient les mêmes informations sur le système. Soit la représentation d'état suivante :

$\displaystyle \dot{x}(t) = A x(t) + B u(t)$

Après un passage en transformée de Laplace (à conditions initiales nulles) :

$\begin{array}{rcl} sX(s) &=& A X(s) + B U(s) \\ (sI - A)X(s) &=& BU(s) \\ X(s) &=& (sI - A)^{-1}BU(s) \end{array}$

on substitue alors $X (s)$ dans l'équation de sortie :

$\begin{array}{rcl} Y(s) &=& CX(s) + DU(s) \\ Y(s) &=& C((sI - A)^{-1}BU(s)) + DU(s) \\ Y(s) &=& G(s)U(s) \\ \end{array}$

$G (s)$ représente la fonction de transfert du système :

$\displaystyle G(s) = C(sI - A)^{-1}B + D$

La fonction de transfert $G (s)$ doit avoir une dimension $q \times p$ , et ainsi un total de qp éléments. On obtient donc $k$ fonctions de transfert différentes, une pour chaque couple entrée/sortie. On comprend aisément que la représentation d'état est plus simple à manipuler lorsqu'il y a plusieurs sorties.

Représentation discrète

Pour passer d'un espace continu à un espace discret de période d'échantillonnage $Δ$ il suffit de faire une approximation de la dérivée. Il existe plusieurs méthodes dont l'approximation de Padé. En général, on utilise toujours $\dot X = \frac{X_{k+1}-X_k}{\Delta}$ . On représente le système de la manière suivante (l'indice k représente l'échantillonnage) :

$X_{k+1} = A_k \cdot X_k + B_k \cdot U_k$ $Y_{k} = C_k \cdot X_k + D_k \cdot U_k$

Pour qu'un système reste observable et commandable par discrétisation, il faut que les valeurs propres $λ i$ de la matrice d'état A du système continue vérifie la condition :

$Re(\lambda _i)=Re(\lambda _j) \Rightarrow Im(\lambda _i - \lambda _j) \ne \frac{2k \pi}{\Delta}$

Les valeurs propres de la matrices d'état A représentent les pôles du système. Si ces pôles possèdent un module inférieur à 1, alors le système est asymptotiquement stable.

Systèmes non linéaires

La forme la plus générale d'une représentation d'état est composée de deux fonctions:

$\begin{cases} \dot{x}(t) &= f(t, x(t), u(t)) \\ y(t) &= h(t, x(t), u(t)) \end{cases}$

x est toujours appelé le vecteur d'état, y vecteur de sortie et u vecteur d'entré. x évolue généralement dans un sous espace J de $\mathbb{R}^n$ .

La première équation représente l'équation d'évolution et la seconde l'équation d'observation. La représentation d'état linéaire décrite préalablement est un cas particulier de cette forme si les fonctions f et h sont linéaires (on représente alors les fonctions sous forme matricielle).