Théorie mathématique de la viabilité (détails)

La théorie mathématique de la viabilité, introduite par Jean-Pierre Aubin, fournit un cadre formel pour étudier la compatibilité entre un système dynamique et des contraintes dans l'espace d'état.

Les problèmes étudiés se veulent les plus larges possibles : ce qui est demandé à l'ensemble de contraintes, en général, c'est d'être un fermé. Cela permet de prendre en compte des contraintes qui soient, par exemple, des unions de convexes (qui ne sont pas nécessairement convexes) ; ou des intersections de sous-variétés régulières (qui ne sont pas nécessairement régulières).

On peut écrire les problèmes de viabilité comme des problèmes d'optimisation où la fonction objectif est discontinue (elle vaut, par exemple, $0$ dans l'ensemble de contraintes et $-\mathrm{\infty }$ en dehors. Les stratégies de contrôle viable peuvent faire appel à des opérateurs comme la vitesse de norme minimale. De ce fait les résultats et théorèmes peuvent être différents de ce qui existe en théorie du contrôle ou en contrôle optimal, en particulier pour les systèmes en temps continu.

On considère :

• un système dynamique $S$ défini sur un ouvert $\mathrm{\Omega }$ d'un espace vectoriel de dimension finie (${ℝ}^n$ dans la suite),
• un ensemble de contraintes $K$ dans l'espace d'état : $K\subset \mathrm{\Omega }$.

Soit $F$ une correspondance de ${ℝ}^n$ dans ${ℝ}^n$. On définit $S$ par :

• en temps continu : $x^{\prime }(t)\in F(x(t))$ for almost all $t\ge 0$ (1)
• en temps discret : $x_{k+1}\in F(x(t))$ for all $k\in \mathbb{N}$ (2)

On note $S(x)$ l'ensemble des évolutions gouvernées par $S$ issues de l'état $x$.

Un domaine de viabilité pour $S$ défini sur ${ℝ}^n$ est un sous-ensemble $D\subset {ℝ}^n$ tel que, à partir de tout $x\in D$, il existe au moins une évolution gouvernée par $S(x)$ qui reste dans $D.$

Notation : dans la suite on utilisera $x(t),t\in \mathbb{N}$ ou $x_k,k\in \mathbb{N}$ suivant le cas pour des raisons de simplification des formules en temps discret.

Noyau de viabilité

Par définition, le noyau de viabilité est le plus grand des domaines de viabilité.

Le noyau de viabilité est un ensemble fermé moyennant certaines conditions sur le système dynamique, pour lesquelles nous rappelons les définitions suivantes :

Soit $G$ une correspondance de $A$ vers $B$ ($A$ et $B$ étant des espaces vectoriels normés de dimension finie), $G:A⇝B$.

• Le domaine de $G$ est $Dom(G)=\{x\in A,G(x)\ne \mathrm{\varnothing }\}$
• Le graphe de $G$ est $Graph(G)=\{(x,y)\in A\times B,y\in G(x)\}$.
• On note $\Vert G(x)\Vert =\underset{y\in G(x)}{\sup }\Vert y\Vert$
• $G$ est à croissance linéaire s'il existe une constante $c>0$ telle que $\mathrm{\forall }x\in Dom(G)$, $\Vert G(x)\Vert \le c(\Vert x\Vert +1)$.

Définition : Soit $F$ une correspondance de $A$ vers $B$. $F$ est Marchaud si :

• Le graphe et le domaine de $F$ sont des fermés non vides
• $\mathrm{\forall }x\in Dom(F)$, $F(x)$ est convexe
• $F$ est à croissance linéaire

Le système dynamique $S$ défini par (1) sera dit Marchaud si la correspondance $F$ est Marchaud.

Modèle:Théorème Il existe des variantes de ces théorèmes avec des conditions plus ou moins faibles pour l'ensemble de contraintes $K$ et le système dynamique $S$.

Caractérisation par le cône contingent

Il est possible de caractériser le noyau de viabilité comme le plus grand des domaines de viabilité de la dynamique $S$.

Un domaine de viabilité peut être caractérisé par le cône contingent.

Soit $x\in K$, $T_K(x)$, le cône contingent à $K$ en $x$ est le cône (fermé) des éléments $v$ tels que :

$\underset{h\to 0+}{\mathrm{lim\; inf}}\frac{d(x+hv,K)}{h}=0$

où $d(x,K)$ est la distance à un ensemble induite par la métrique sur $K$.

Concrètement, cela veut dire que si $v\in T_K(x)$, il y a une suite $(h_n>0)$ convergeant vers $0$, et une suite $(v_n)$ convergeant vers $v$, telle que $x+h_n{v}_n\in K$ pour tout $n$.

On peut montrer l'équivalence entre cette définition et la définition par les évolutions, car le cône contingent caractérise les évolutions non sortantes.

Cette caractérisation permet à la théorie de disposer d'algorithmes généraux de calcul d'approximation du noyau de viabilité avec preuve de convergence.

Systèmes contrôlés

Il existe une formulation plus explicite de ces théorèmes pour les systèmes contrôlés.

$U$ est la correspondance qui à chaque état $x$ associe l'ensemble des contrôles admissibles en $x$. $f$ est la fonction qui à un état et un contrôle associe un nouvel état.

Le système dynamique contrôlé $S_{(f,U)}$ en temps continu est Marchaud si :

• f est continue
• Le graphe de U est fermé
• Les ensembles $f(x,u),u\in U(x)$ sont convexes pour tout $x\in Dom(F)$
• f et U sont à croissance linéaire

Cartes de régulation

Dans le cas des systèmes contrôlés, quand le noyau de viabilité n'est pas vide, on peut définir la correspondance des contrôles viables $\tilde{U}:{ℝ}^n⇝{ℝ}^p$ qui à chaque état du noyau associe le sous-ensemble des contrôles admissibles qui permet de garder les évolutions dans le noyau : ${\text{viab}}_{S_{(f,U)}}(K)$ est un ensemble invariant pour le système contrôlé $S_{(f,\tilde{U})}$.

En temps continu, si l'état $x$ est dans l'intérieur du noyau de viabilité, tous les contrôles de $U(x)$ sont viables : $\tilde{U}(x)=U(x)$. C'est sur la frontière du noyau que l'ensemble des contrôles viables peut être strictement plus petit que l'ensemble des contrôles admissibles.

À partir des cartes de régulation viable, on peut définir des stratégies de régulation en boucle fermée. On peut citer en particulier :

• La stratégie du contrôle lourd qui consiste à sélectionner un contrôle qui minimise $\Vert \frac{du}{dt}\Vert$
• La stratégie du contrôle lent qui sélectionne le contrôle de norme minimale en tout état $x$ : ${\text{argmin}}_{\tilde{U}(x)}\Vert u\Vert$

L'inconvénient de ces stratégies est que les trajectoires qu'elles gouvernent passent beaucoup de temps sur la frontière du noyau de viabilité, ce qui les expose à quitter le noyau de viabilité en cas de perturbation. Des heuristiques de trajectoires "prudentes" ont été proposées pour remédier à ce problème.

Bassin de capture

La théorie de la viabilité s'intéresse à la possibilité intrinsèque de capturer une cible en temps fini tout en respectant des contraintes de viabilité. Le bassin de capture d'une cible rassemble les états de l'ensemble de contraintes à partir desquels il existe une évolution viable atteignant la cible en temps fini.

Viabilité garantie

La propriété de viabilité garantie est recherchée dans le cas des systèmes dynamiques contrôlés qui prennent en compte des incertitudes de nature ensembliste (tychastique).

$U$ est la correspondance qui à chaque état $x$ associe l'ensemble des contrôles admissibles en $x$. $V$ est la correspondance qui à chaque état $x$ associe l'ensemble des perturbations admissibles en $x$. $f$ est la fonction qui à un état et un contrôle associe un nouvel état.

Aubin définit la propriété de viabilité garantie d'un ensemble comme la possibilité de trouver une carte de contrôle définie sur cet ensemble, pour laquelle l'ensemble est invariant, c'est-à-dire que pour le système défini avec cette carte de contrôle, toutes les évolutions restent dans l'ensemble, quelles que soient les perturbations (admissibles).

Contrairement au cas de la viabilité contrôlé, ou de la viabilité tychastique sans contrôle, le plus grand des domaines de viabilité garantie dans un ensemble de contraintes $K$ fermé n'a pas, en général, cette propriété pour un système dynamique en temps continu. Il est nécessaire d'exiger des conditions supplémentaires

On rappelle la définition d'une correspondance lipschitzienne. Soit $G:A⇝B$ (espaces vectoriels de dimension finie), $G$ est lipschitzienne de constante $\lambda >0$ si pour tout $x_1,x_2$ dans $A$, $G(x_1)\subset G(x_2)+\lambda \Vert x_1-x_2\Vert B$ (où $B$ est la boule unité).

Dans le cas des systèmes continus, on peut montrer l'existence du noyau de viabilité garantie dans le cas où la carte de contrôle est lipschitzienne de constante $\lambda >0$. Le noyau de viabilité n'est donc pas absolu, il est défini en fonction de la constante de Lipschitz.

Dans le cas discret, la fermeture du noyau de viabilité garantie est réalisée dans les mêmes conditions que pour le noyau de viabilité. Dans le cas continu, on peut montrer qu'elle est réalisée si $f$ lipschitzienne, ainsi que $U$ et $V$.

Bibliographie

Les ouvrages et documents suivants ont servis à rédiger ce texte :

Viability Theory, Aubin J.P, Birkhauser, 1991.

Viability Theory: New Directions, Aubin, J.-P., Bayen, A., & Saint-Pierre, P. (2011). Springer.

Modélisation, exploration de modèles, théorie de la viabilité pour l’aide à la gestion des espèces végétales invasives : application au cas des renouées asiatiques le long d’un cours d’eau. Lavallée, F. Mathématiques [math]. Thèse de Sorbonne Université, 2020.