Commande H-infini

Dans le domaine de l'automatique et de la théorie du contrôle, la synthèse H-infini ou H_∞ est une méthode qui sert à la conception de commandes optimales, imposant des contraintes sur la norme H-infini d'un système dynamique (minimisation, ou imposition d'une borne supérieure). La norme H-infini d'un système peut s'interpréter comme l'amplification énergétique maximale qu'il peut exercer en sortie par rapport à l'énergie de ses signaux d'entrée. Généralement réservée aux systèmes linéaires temps invariant, la commande H-infini a néanmoins été généralisée dans des cadres non-linéaires. Tout comme pour la commande H2, le terme H-infini provient de l'espace fonctionnel de Hardy du même nom.

Présentation

La synthèse H-infini est une méthode qui sert à la conception d'un contrôleur capable de délivrer une commande optimale particulière, prenant en compte des contraintes sur le comportement énergétique attendu d'un système muni de ce contrôleur. Comme pour toute commande optimale, le mot « optimal » est ici utilisé dans un sens strictement mathématique, l'optimisation étant toujours relative à un objectif recherché.

• En synthèse H-infini optimale, on cherche à obtenir un contrôleur qui minimisera le transfert d’énergie maximal des entrées du système vers les sorties du système.
• En synthèse H-infini suboptimale, on cherche simplement un contrôleur pour lequel ce transfert d'énergie sera majoré par une valeur choisie à l'avance.

La commande H-infini a pour principal avantage de synthétiser des concepts liés :

• à la commande classique, typiquement par son interprétation fréquentielle
• au contrôle moderne, par son aspect MIMO, son interprétation temporelle, et sa capacité à tenir compte d'incertitudes modèle.

C'est une méthodologie très utilisée pour assurer la robustesse d'un contrôleur.

Définitions de la norme H-infini

Définition fréquentielle

Si on considère la matrice de transfert $\mathbf {T} (s)$ d'un système dynamique linéaire stable qui à une entrée $w(s)$ associe une sortie $z(s)=\mathbf {T} (s)w(s)$, la norme ${\mathcal {H}}_{\infty }$ du système dynamique est donnée par le maximum de la valeur singulière maximale (ou plus rigoureusement, par le supremum essentiel de la valeur singulière maximale) de $\mathbf {T} (j\omega )$ avec $\omega \in \mathbb {R}$ (mais on peut se restreindre au sup pour $\omega \geq 0$). Il s'agit de la définition de la norme de l'espace des matrices $n_{z}\times n_{w}$ dont les coefficients appartiennent à l'espace de Hardy ${\mathcal {H}}_{\infty }$.

$\lVert \mathbf {T} \rVert _{{\mathcal {H}}_{\infty }}\triangleq {\underset {\omega \in \mathbb {R} }{\operatorname {ess} \sup }}\;{\bar {\sigma }}(\mathbf {T} (j\omega ))=\sup _{\omega \geq 0}{\bar {\sigma }}(\mathbf {T} (j\omega ))$

Cette définition a le mérite d'être visuellement claire dans le cas SISO, car la norme ${\mathcal {H}}_{\infty }$ d'un système peut alors directement se déduire de la valeur de gain maximale de $|\mathbf {T} (j\omega )|$ sur un diagramme de Bode. Dans ce cas, la norme ${\mathcal {H}}_{\infty }$ d'un système peut également s'interpréter comme le rapport d'amplitude maximal entre un signal sinusoïdal d'entrée et le signal sinusoïdal de sortie associé (en régime permanent). La valeur $\omega ^{*}$ où ce maximum est atteint fournit la fréquence de résonance maximale du système.

Définition temporelle

Si on considère un système dynamique stable $\mathbf {T}$ (de condition initiale $x(0)=0$) comme un opérateur entre des espaces fonctionnels de Lebesgue $L^{2}$, qui à une entrée $w:[0,+\infty )\rightarrow \mathbb {R} ^{n_{w}}$ dans $L^{2}$ associe une sortie $z=\mathbf {T} w:[0,+\infty )\rightarrow \mathbb {R} ^{n_{z}}$ dans $L^{2}$, alors la norme ${\mathcal {H}}_{\infty }$ du système dynamique est définie comme le suprémum du rapport entre la norme $L^{2}$ des signaux de sorties et la norme $L^{2}$ des signaux d'entrées. Formellement, il s'agit de la norme d'opérateur subordonnée aux normes $L^{2}$ des signaux.

$\lVert \mathbf {T} \rVert _{{\mathcal {H}}_{\infty }}\triangleq \sup _{w\in L^{2*}}{\frac {\lVert \mathbf {T} w\rVert _{L^{2}}}{\lVert w\rVert _{L^{2}}}}=\sup _{w\in L^{2*}}{\sqrt {\frac {\int _{0}^{+\infty }z^{\top }(t)z(t)dt}{\int _{0}^{+\infty }w^{\top }(t)w(t)dt}}}$

Le carré de la norme $L^{2}$ d'un signal pouvant s'interpréter comme l'énergie de ce signal, il est ici clair que les critères ${\mathcal {H}}_{\infty }$ sont des critères d'amplification maximale de l'énergie des signaux de sortie par rapport à l'énergie de leurs signaux d'entrée associés. Cette définition temporelle est équivalente à la définition fréquentielle dans le cas des systèmes linéaires, ce qui peut se démontrer à l'aide de l'égalité de Parseval. Elle a également le mérite de pouvoir se généraliser aux systèmes non-linéaires, ainsi qu'à des rapports de normes d'espaces fonctionnels autres que $L^{2}$ (on applique parfois des normes d'espaces de Sobolev). Bien que le terme ${\mathcal {H}}_{\infty }$ soit toujours employé en contrôle pour parler de ces généralisations, il n'a plus de rapport avec l'espace de Hardy qu'il désignait originellement.

Formulation du problème

Considérant le système représenté par la forme standard suivante:

Le système P possède deux entrées:

• w représente les entrées extérieures, notamment le vecteur de consigne de commande, les perturbations et les bruits,
• u représente le vecteur de commande ;

et deux sorties :

• z les sorties à optimiser pour avoir un bon comportement de la commande
• v les mesures disponibles utilisées par le contrôleur pour calculer la commande.

Toutes ces données sont d’une façon générale des vecteurs et P et K sont des matrices de transfert.

Le système précédent peut s’écrire en fréquentiel sous la forme:

${\begin{bmatrix}z\\v\end{bmatrix}}=\mathbf {P} (s)\,{\begin{bmatrix}w\\u\end{bmatrix}}={\begin{bmatrix}P_{11}(s)&P_{12}(s)\\P_{21}(s)&P_{22}(s)\end{bmatrix}}\,{\begin{bmatrix}w\\u\end{bmatrix}}$

avec le contrôleur :

$u=\mathbf {K} (s)\,v$

Il est alors possible d’exprimer la relation entre z et w par:

$z=F_{\ell }(\mathbf {P} ,\mathbf {K} )\,w$

Où $F_{\ell }$ est appelé la ‘’Transformation Linéaire Fractionnaire’’ (LFT) définie par:

$F_{\ell }(\mathbf {P} ,\mathbf {K} )=P_{11}+P_{12}\,\mathbf {K} \,(I-P_{22}\,\mathbf {K} )^{-1}\,P_{21}$

L’objectif de la synthèse ${\mathcal {H}}_{\infty }$ est de trouver le contrôleur $\mathbf {K}$ tel que la norme ${\mathcal {H}}_{\infty }$ de $F_{\ell }(\mathbf {P} ,\mathbf {K} )$ soit minimum. La norme ${\mathcal {H}}_{\infty }$ de la matrice de fonctions de transfert $F_{\ell }(\mathbf {P} ,\mathbf {K} )$ est donnée par :

$||F_{\ell }(\mathbf {P} ,\mathbf {K} )||_{{\mathcal {H}}_{\infty }}=\sup _{\omega \geq 0}{\bar {\sigma }}(F_{\ell }(\mathbf {P} ,\mathbf {K} )(j\omega ))$

où ${\bar {\sigma }}$ est la valeur singulière maximale de la matrice $F_{\ell }(\mathbf {P} ,\mathbf {K} )(j\omega )$.

Il existe plusieurs méthodes permettant le calcul et l’optimisation de la norme ${\mathcal {H}}_{\infty }$ :

• méthode basée sur la résolution des équations de Riccati ;
• méthode(s) basée sur les Inégalités matricielles linéaires ;
• méthode utilisant la paramétrisation de Youla.

Applications

La synthèse de commande ${\mathcal {H}}_{\infty }$ est depuis le début des années 2000 utilisée dans le cadre du contrôle de structures flexibles des systèmes spatiaux comme les satellites ou les lanceurs. Ces systèmes possèdent des modes de flexions proches de la bande passante du contrôleur et dont les caractéristiques sont incertaines (fréquence, gain, amortissement). La méthode ${\mathcal {H}}_{\infty }$ permet de garantir théoriquement une certaine robustesse de la stabilité du système (via la marge de module) tout en assurant les performances de réjection de perturbation.

Sur Ariane 5, le pilotage basé ${\mathcal {H}}_{\infty }$ a ainsi permis de garantir la stabilité du lanceur tout en optimisant la commande de braquage de la tuyère, faisant économiser plusieurs centaines de kilos sur le liquide hydraulique de commande.

Limitations

Comme toutes les techniques de contrôle, la synthèse H-infini possède quelques désavantages. Il peut être par exemple difficile d'imposer un temps de réponse particulier à la sortie du système, ou encore un taux de dépassement maximal.

Articles connexes

• Commande H2
• Commande LQ
• Commande LQG
• Commande optimale

Notes et références

Références

• Daniel Alazard et al, Robustesse et Commande Optimale, Toulouse, Cepadues, 1999, 348 p. (ISBN 2-85428-516-6)
• Gilles Duc et Stéphane Font, Commande Hinfini et mu-analyse -- des outils pour la robustesse, Paris, Hermès-Science, 1999, 121 p. (ISBN 2-7462-0041-4)