Projet : META2

Précédent : Fondements scientifiques Remonter : Fondements scientifiques Suivant : Commande Stochastique

Sous-sections

• L'algèbre max-plus
• Systèmes dynamiques max-plus linéaires et graphes d'événements
• Systèmes dynamiques (min, + , x )implicites, réseaux de Petri, programmation dynamique stochastique
• Systèmes dynamiques monotones homogènes
• Ressources partagées, empilements de tâches et automates max-plus
• Vers une approche géométrique des systèmes dynamiques linéaires sur le semianneau (max, +)
• Dualité entre l'optimisation et le calcul des probabilités
• Page WEB

   
Systèmes à événements discrets et algèbre max-plus

Mots clés : Algèbre max-plus, Systèmes à événements discrets, Décision markovienne. .

Participants : Marianne Akian, Guy Cohen, Stéphane Gaubert, Jean-Pierre Quadrat.

Résumé :

L'idée est de substituer au corps des nombres réels le semicorps idempotent $\mathbb{R}$_max des nombres réels auxquels on ajoute - $\infty$ et que l'on munit des opérations max et +. A partir de ce semicorps, on peut développer l'analogue de l'algèbre linéaire classique, une théorie des systèmes dynamiques linéaires, un analogue du calcul des probabilités, l'analogue des espaces de Sobolev etc ...

Les nouveaux systèmes dynamiques linéaires ainsi obtenus ne sont rien d'autre que les systèmes régis par des équations de la programmation dynamique qui jouent un rôle important en ingéniérie et en physique.

L'algèbre max-plus

Le semianneau max-plus, parfois appelé ``algèbre max-plus'', et noté $\mathbb{R}$_max, est l'ensemble $\mathbb{R}$ $\cup$ { - $\infty$}, muni des lois max et +. Traditionnellement, on note $\oplus$ au lieu de max ( 2 $\oplus$ 3 = 3), et $\otimes$ au lieu de + (1 $\otimes$ 1 = 2). L'élément zéro, c'est-à-dire le neutre pour la loi $\oplus$, est noté $\varepsilon$(ici $\varepsilon$ = - $\infty$). L'unité, c'est-à-dire le neutre pour la loi $\otimes$, est notée e(ici e = 0). Les axiomes de structure des semianneaux sont ici satisfaits : $\oplus$ est associative, commutative, a un zéro, $\otimes$ est associative, a une unité, distribue par rapport à $\oplus$, et zéro est absorbant. Le semianneau max-plus est très particulier : il est commutatif (a $\otimes$ b = b $\otimes$ a), idempotent (a $\oplus$ a = a), et les éléments ont un inverse, hormis zéro (on nomme semicorps les semianneaux qui satisfont cette dernière propriété).

Ces nouvelles notations $\oplus$ et $\otimes$sont mieux adaptées que max et + aux réflexes de calcul qui sont liés à l'algèbre ordinaire. Nous écrirons par exemple :

$$\otimes \otimes \otimes \sqrt{3}$$

$$\begin{bmatrix}2 & 0 \\ 4 & \varepsilon \end{bmatrix} \begin{bmatrix}10 \\ 103 \end{bmatrix} \begin{bmatrix}2 \otimes 10\oplus 0 \otimes 103\\ 4 \otimes 10 \oplus \varepsilon \otimes 103 \end{bmatrix} \begin{bmatrix}103 \\ 14 \end{bmatrix}$$

$$\oplus \oplus \oplus \oplus \oplus \oplus$$

Le tableau 1 donne une liste d'autres semianneaux apparentés au semicorps max-plus.

 

Tableau 1: La famille des semianneaux (max, +) et tropicaux

$\mathbb{R}$max	($\mathbb{R}$ $\cup$ { - $\infty$}, max, +)	semianneau (max, +)	semicorps idempotent
$\overline{\mathbb{R}}_{\max}^{}$	($\mathbb{R}$ $\cup$ { $\pm$ $\infty$}, max, +)	semianneau (max, +) complété	- $\infty$ + (+ $\infty$) = - $\infty$, car $\varepsilon$ $\otimes$ a = $\varepsilon$
$\mathbb{R}$max, x	($\mathbb{R}$+, max, x )	semianneau (max, x )	isomorphe à $\mathbb{R}$max ( x $\mapsto$ log x)
$\mathbb{R}$min	($\mathbb{R}$ $\cup$ { + $\infty$}, min, +)	semianneau (min, +)	isomorphe à $\mathbb{R}$max ( x $\mapsto$ - x)
$\mathbb{N}$min	($\mathbb{N}$ $\cup$ { + $\infty$}, min, +)	semianneau tropical	fameux en théorie des langages
$\mathbb{R}$max, min	($\mathbb{R}$ $\cup$ { $\pm$ $\infty$}, max, min)	algèbre goulot, ou floue	n'est pas un semicorps.
$\mathbb{B}$	({false,true},or,and)	semianneau de Boole	isomorphe à ({$\varepsilon$, e}, $\oplus$ , $\otimes$ ), pour n'importe lequel des semianneaux ci-dessus
$\mathbb{R}$h	($\mathbb{R}$ $\cup$ { - $\infty$}, $\oplus_{h}^{}$ , +) a $\oplus_{h}^{}$ b = hlog(ea/h + eb/h)	semianneaux ``de Maslov''	déformation de ($\mathbb{R}$+, + , x ) $\lim_{h\to 0^+}^{}$$\mathbb{R}$h = $\mathbb{R}$0 = $\mathbb{R}$max

 

Systèmes dynamiques max-plus linéaires et graphes d'événements

Certains systèmes à événements discrets (les graphes d'événements) se modélisent par des systèmes dynamiques sur le semianneau max-plus. Considérons par exemple un système comprenant n tâches répétitives (on peut penser à des assemblages dans un atelier, à des rendez-vous dans un réseau). L'instant de démarrage d'une tâche est conditionné par l'exécution d'autres tâches devant être effectuées antérieurement. Supposons que la k-ième occurrence de la tâche i (tir de la transition i) ne puisse commencer que $\tau_{ij}^{}$ unités de temps après que la dernière des occurences (k - $\nu_{ij}^{}$) des tâches j n'ait eu lieu. Alors, clairement, les dates au plus tôt d'occurrence des k-ième tâches sont données par la récurrence

 

$$\max_{j}^{} \left[\vphantom{ \tau_{ij} + x_j(k-\nu_{ij})}\right. \tau_{ij}^{} \nu_{ij}^{} \left.\vphantom{ \tau_{ij} + x_j(k-\nu_{ij})}\right]$$ (1)

qui n'est autre qu'un système linéaire stationnaire sur le semianneau max-plus. On peut l'écrire en termes matriciels

$$\bigoplus_{\nu \in \mathsf{F}}^{} \nu \nu$$ (2)

où $\mathsf {F}$ est un ensemble fini, et {A_$\nu$}_{$\nu$ $\in$ $\mathsf {F}$}une famille de matrices que l'on écrira aisément.

On a obtenu un modèle linéaire pour l'évaluation de performance d'un système dans lequel les préconditions d'une tâche sont fixées a priori et invariantes en temps. En termes de système de production cela signifie à peu près qu'un ordonnancement périodique a été défini.

L'algèbre max-plus permet d'analyser complètement cette classe de systèmes. Un des résultats intéressants obtenus consiste à remarquer que la relation : instants d'entrée -instants de sortie des pièces, est max-plus linéaire et shift invariante. C'est une inf-convolution de la suite des instants d'entrée et d'une suite caractéristique du système (sa réponse impulsionnelle au sens max-plus). On a donc une bonne notion de fonction de transfert. Ce sont les séries formelles rationnelles (si le graphe d'événement n'a qu'un nombre fini de transitions) à coefficients dans l'algèbre max-plus. Ces résultats ont été largement expliqués dans [[6],[3],[9]].

Cette classe est trop restreinte (dans la majorité des applications, trouver un bon ordonnancement fait partie du problème). Elle peut être élargie au prix de la perte de la linéarité max-plus.

Systèmes dynamiques (min, + , x )implicites, réseaux de Petri, programmation dynamique stochastique

La modélisation de réseaux de Petri temporisés généraux ne peut pas se faire dans le semianneau $\mathbb{R}$_max ou $\mathbb{R}$_min, la multiplication est indispensable.

On note $\mathcal {P}$ l'ensemble des places, $\mathcal {Q}$ l'ensemble des transitions, M les multiplicités des arcs, m le marquage initial, $\tau$ les temps de séjour. À chaque place p $\in$ $\mathcal {P}$ on associe un compteur Z_p : $\mathbb{R}$ $\rightarrow$ $\mathbb{Z}$, (Z_p(t) est le numéro du dernier jeton entré avant l'instant t, en incluant les jetons du marquage initial) ; de même, pour chaque transition q $\in$ $\mathcal {Q}$, Z_q(t) note le numéro du dernier tir de q arrivé avant l'instant t. Ces fonctions sont croissantes, par définition. On a les bilans suivants (on note xⁱⁿ et x^out les ensembles des prédecesseurs et successeurs d'un noeud x, respectivement):

 

$$\sum_{q\in p^{\text{\rm in}}}^{}$$ (3)

 

$$\tau_{p}^{} \geq \sum_{q'\in p^{\text{\rm out}}}^{}$$ (4)

Par exemple, pour le réseau de Petri de la Figure 1, Z_q3(t) + 2Z_q2(t) $\leq$ Z_p2(t - $\tau_{p_2}^{}$)traduit le fait qu'un jeton en p₂ peut contribuer au tir d'une transition q₃, et qu'il en faut deux pour tirer q₂.

  

Figure 1: Un réseau de Petri temporisé

Dans un fonctionnement au plus tôt, au moins une inégalité de type (4) est saturée (i.e. au moins une place est indisponible) en amont de chaque transition. On a donc:

 

$$\forall \in \mathcal {Q} \min_{p\in q^{\text{\rm in}}}^{} \left\{\vphantom{m_p+\sum_{q''\in p^{\text{\rm in}}}M_{pq''}Z_{q''}(t) - \sum_{q'\in p^{\text{\rm out}}} Z_{q'}(t) M_{q'p}}\right. \sum_{q''\in p^{\text{\rm in}}}^{} \sum_{q'\in p^{\text{\rm out}}}^{} \left.\vphantom{m_p+\sum_{q''\in p^{\text{\rm in}}}M_{pq''}Z_{q''}(t) - \sum_{q'\in p^{\text{\rm out}}} Z_{q'}(t) M_{q'p}}\right\}$$ (5)

Ce problème posé dans [Lib97] est un système dynamique (min, + , x ) implicite sur lequel on ne sait à peu près rien dire. Résoudre ce genre de problème, même dans le cas statique, revient à étudier les variétés algébriques max-plus. La programmation linéaire devrait être d'une aide précieuse pour la résolution de ces systèmes.

On peut néanmoins simplifier ce problème en forçant des routages stationnaires : on notera alors $\rho_{qp}^{}$ la proportion de fluide (cas où on ne se préoccupe pas de l'intégrité des jetons) routé vers q par la place p, avec $\sum_{q\in p^{\text{\rm out}}}^{}$$\rho_{qp}^{}$ = 1. En posant,

$$\mu_{pq}^{}$$$$\;\stackrel{\text{\rm d\'ef}}{=}\;$$M_pq,    $$\mu_{qp}^{}$$$$\;\stackrel{\text{\rm d\'ef}}{=}\;$$M_qp^-1,  $$\mu{^\prime}_{qp}$$$$\;\stackrel{\text{\rm d\'ef}}{=}\;$$$$\mu_{qp}^{}$$$$\rho_{qp}^{}$$  

l'évolution du système devient cette fois :

 

$$\min_{p\in q^{\text{\rm in}}}^{} \Bigl[ \mu{^\prime}_{qp} \bigl( \sum_{q''\in p^{\text{\rm in}}}^{} \mu_{pq''}^{} \tau_{p}^{} \bigr) \Bigr]$$ (6)

On reconnait une équation de la programmation dynamique, pour un processus de décision semi-Markovien, avec coût additif actualisé général. On peut alors calculer l'asymptotique de Z par des algorithmes du type itération sur les politiques [CGQ95].

Systèmes dynamiques monotones homogènes

La modélisation de systèmes à événements discrets plus généraux peut se faire en remplaçant (1), (6) par des dynamiques f : $\mathbb{R}$ⁿ $\rightarrow$ $\mathbb{R}$ⁿqui satisfont les trois axiomes :

Homogénéité	f ($$\lambda$$ + x)	=	$$\lambda$$ + f (x)	(H)
Monotonie	x $$\leq$$ y	$$\Rightarrow$$	f (x) $$\leq$$ f (y)	(M)
Non-expansivité	| f (x) - f (y)|$\infty$	$$\leq$$	| x - y|$\infty$ ,	(N)

pour tout $\lambda$ $\in$ $\mathbb{R}$, x, y $\in$ $\mathbb{R}$ⁿ, avec $\lambda$ + x$\;\stackrel{\text{\rm def}}{=}\;$($\lambda$ + x₁,...,$\lambda$ + x_n), et où $\leq$ désigne l'ordre partiel usuel sur $\mathbb{R}$ⁿ.

Sous l'hypothèse H, on sait que les propriétés M et N sont équivalentes. Gunawardena et Keane ont proposé d'appeler ``topicales'' les applications qui vérifient ces trois propriétés dont l'importance est bien connue [Nus90,CT80].

On peut dire que H, N sont le minimum vital pour faire une théorie de Perron-Frobenius, la propriété H donne du sens au problème spectral f (x) = $\lambda$ + x, la propriété N garantit que la limite $\chi$(f )= $\lim_{k\to\infty}^{}$f^k(x)/k, lorsqu'elle existe, est indépendante du point de départ x $\in$ $\mathbb{R}$ⁿ.

On essaie donc dans ce cadre de généraliser la théorie de Perron-Frobenius, et de comprendre géométriquement les asymptotiques des itérées de f (qui donnent par exemple, les performances asymptotiques de systèmes à événements discrets : ateliers, circuits digitaux; la quantité $\chi$(f ) donne alors l'inverse du taux de production, ou le temps de cycle).

Ressources partagées, empilements de tâches et automates max-plus

Un système à ressources partagées peut être vu comme un empilement de tâches (représentées par les contraintes temporelles imposées sur les ressources par ces tâches) du type jeu de Tetris.

Dans l'espace temps-ressource $\mathbb{R}$ x $\mathcal {R}$ avec $\mathcal {R}$ = {1,..., n}. Une tâche a est une classe d'équivalence de formes géométriques (deux formes sont équivalentes si elles sont superposables par une translation temporelle). D'un point de vue physique, une tâche mobilise des ressources de façon synchronisée. Elle peut être retardée ou avancée. La tâche a utilise les ressources R(a). Une tâche en position $\lambda$ immobilise la ressource r à la date d_r(a) + $\lambda$et la libère à la date f_r(a) + $\lambda$. d (a), f (a) : R(a) $\rightarrow$ $\mathbb{R}$ sont donnés.

Un travail w est un ensemble ordonné de tâches à accomplir w = a₁ ^... a_k. On peut alors vouloir calculer les instants de libération des ressources au plus tôt x(w) pour un travail w donné, sachant que les ressources étaient disponibles au départ aux instants g ( g : $\mathcal {R}$ $\rightarrow$ $\mathbb{R}$). Cela correspond au calcul du profil haut d'un empilement de pièces si l'on donne une interprétation spatiale à l'axe temporel. Le vecteur g est alors le profil du sol. L'exécution au plus tôt des différentes tâches d'un travail w revient à laisser tomber les pièces a₁,..., a_k sur le sol g, comme dans le jeu de tétris. Ce profil haut se calcule récursivement par la programmation dynamique.

 

  

Figure 2: Un jeu de Tetris à trois types de pièces et quatre ressources

R(c) = {2, 4}, d (c) = [ . , 0, . , 0], f (c) = [ . , 2, . , 2] R(b) = {1, 2}, d (b) = [0, 0, . , . ], f (b) = [2, 2, . , . ] R(a) = {1, 2, 3}, d (a) = [0, 0, 0, . ], f (a) = [1, 1, 3, . ]

En termes algébriques, à chaque pièce a prise dans un ensemble de pièces $\mathcal {T}$, on associe la matrice M(a) $\in$ ($\mathbb{R}$_max)^{$\mathcal {R}$ x $\mathcal {R}$} M(a)_{r, s} = f (a)_s - d (a)_r si r, s $\in$ R(a), et M(a)_{r, r} = e pour les coefficients diagonaux hors R(a)(les autres coefficients valent $\varepsilon$). On a :

 

$$\mathcal {R}$$ (10)

où e_{$\mathcal {R}$} désigne un vecteur de zéro et donc y(w)représente la date de fin du travail (la hauteur du tas).

Autrement dit, la série génératrice $\bigoplus_{w\in \mathcal{T}^*}^{}$y(w)w $\in$ $\mathbb{R}$_max$\langle$$\langle$$\Sigma$$\rangle$$\rangle$ est reconnue par un automate à coefficients dans le semianneau max-plus. D'après le théorème de Kleene-Schützenberger, c'est une série rationnelle.

Les classiques problèmes d'ordonancement reviennent alors à minimiser la vitesse de croissance du tas par rapport à l'ordre des lettres dans le mot w. Ils se ramènent à l'étude des semigroupes de matrices à plusieurs générateurs M(a₁), ^... , M(a_k). Ce travail est présenté dans [[8]].

Vers une approche géométrique des systèmes dynamiques linéaires sur le semianneau (max, +)

La solution du système

x(k) = Ax(k - 1) $$\oplus$$ Bu(k)  , y(k) = Cx(k)  ,

où x(- 1) = - $\infty$, est donnée par :

y(k) = $$\bigoplus_{l=0}^{k}$$CA^{k - l}Bu(l )  .

Dans un atelier, y(k) s'interprète comme la date d'obtention de la k-ième pièce (on suppose, afin d'alléger la notation, que le système a une seule entrée u(k) $\in$ $\mathbb{R}$ et une seule sortie y(k) $\in$ $\mathbb{R}$).

Les problème de suivi de trajectoire consiste à rechercher u tel que y(k) $\leq$ $\eta$(k) , pour k = 0, 1,... ($\eta$ représente un carnet de commandes à suivre). Le problème d'atteignabilité en horizon k consiste à rechercher u tel que x(k) $\leq$ $\xi$($\xi$ est un vecteur des dates de disponibilité des ``ressources'' internes, machines, palettes, personnel, etc.).

Ces deux problèmes sont de la forme Fu $\leq$ v, ou F est un opérateur linéaire. E.g., dans le deuxième cas, Fu = $\mathcal {C}$_ku_k, où $\mathcal {C}$_k = [B, AB,..., A^kB]et u_k = [u(k),..., u(0)]^T. Il n'est pas en général possible de résoudre exactement Fu = v pour toute valeur de v, car génériquement, un opérateur linéaire (max, +) n'est ni surjectif ni injectif. Cependant, il est facile de voir que Fu $\leq$ v a toujours une solution maximale, F$\backslash$v, donnée en termes matriciels par

(F$$\backslash$$v)_i = $$\inf_{j}^{}$$(- F_ji + v)  .

(il s'agit d'un cas particulier de résiduation [BJ72], [[3]]). On peut ainsi calculer la solution au plus tard satisfaisant le carnet de commande, où l'état atteignable maximal en dessous de $\xi$. En outre, on sait tester l'atteignabilité exacte de $\xi$ en temps k : il suffit de vérifier que l'égalité $\mathcal {C}$_k($\mathcal {C}$_k$\backslash$$\xi$) = $\xi$ est vraie.

La version la plus achevée de ces problèmes d'atteignabilité dans le cas classique a été donnée dans [Won79]. Elle est basée sur une discussion des espaces invariants d'un opérateur linéaire. L'extension de cette théorie au cas max-plus passe par une meilleure compréhension des semi-modules (analogue sur un semianneau des espaces vectoriels). Même lorsque les semimodules sont finiment engendrés, très peu de choses sont connues. Plusieurs notions d'indépendance existent, conduisant, chacunes, à des notions, distinctes, de rang. De plus un sous semimodule n'a pas toujours de supplémentaire. Ces difficultés sont répertoriées dans la théorie des modules. Cependant, une nouvelle difficulté surgit ici : l'analogue d'un noyau doit être posé dans un espace double ker C = {(x, y) | Cx = Cy}. En effet, dans cette structure tous les termes d'une équation ne peuvent pas être ramenés dans un même membre, puisque, en général, un élément n'a pas d'opposé. Par contre, la résiduation donne des résultats intéressants que nous sommes en train d'explorer [CGQ97].

Dualité entre l'optimisation et le calcul des probabilités

Dans de nombreux domaines (grandes déviations en probabilité, entropie en mécanique statistique et théorie de l'information, dualité filtrage-commande en commande, maximum de vraisemblance en statistique, passage de la mécanique quantique à la mécanique classique, méthode de Hopf pour la résolution explicite d'équations d'Hamilton Jacobi etc.) une dualité entre le calcul des probabilités et l'optimisation est utilisé implicitement.

Dans un certain nombre de travaux relativement récents, cette dualité a été explicitée par divers groupes indépendamment [Mas87,Pap95,dTRS90,Qua90,AQV94]. On dispose maintenant d'une théorie trés satisfaisante [[21]]. Cette dualité, a pour intérêt, en dehors de son caractère esthétique et unifiant, de faire ressortir les dissymétries existantes dans le développement des deux domaines. Le calcul des probabilités contient beaucoup de résultats analytiques et conceptuels. L'optimisation est beaucoup plus développée au niveau algorithmique.

Un problème d'optimisation est souvent présenté comme le calcul de

$$\inf_{x\in A}^{}$$c(x)  ,

pour des c et des A suffisamment simples pour qu'il existe un x optimal calculable algorithmiquement. Des notions importantes sont apparues, comme la dualité, surtout pour des raisons algorithmiques.

Dans le calcul des probabilités on s'intéresse au calcul de la probabilité d'événements, c.a.d.

$$\int_{A}^{}$$p(x)dx

pour des A et p pouvant être trés compliqués. On appelle P(A)le résultat et on fait la théorie des applications A $\mapsto$ P(A). On fait apparaître toutes sortes de notions intéressantes : moyenne, variance, fonction caractéristique, etc. On fait aussi des calculs, mais surtout, dans des cas où il est possible d'obtenir des résultats explicites.

Il suffit d'appeler C(A)$\;\stackrel{\text{\rm d\'ef}}{=}\;$$\inf_{x\in A}^{}$c(x)   pour être dans une situation complètement analogue au calcul des probabilités au remplacement de l'algèbre ordinaire par $\mathbb{R}$_min près. Pratiquement tous les concepts des probabilités ont une notion duale, utile en optimisation. Citons en trois. Le dual de la loi Gaussienne $\mathcal {N}$_m,$\sigma$ est la forme quadratique

$$\mathcal {Q}$$_m,$\sigma$(x)$$\;\stackrel{\text{\rm d\'ef}}{=}\;$$$${\textstyle\frac{1}{2}}$$($${\frac{x-m}{\sigma}}$$)²  .

Le dual de la transformée de Laplace est la transformée de Fenchel. Le dual de la convolution des lois est l'inf-convolution des coûts.

A l'opposé, aucun algorithme sérieux n'est donné, pour calculer la probabilité d'un événement un peu compliqué de $\mathbb{R}$ⁿ, exceptée la méthode de Monte-Carlo.

Il existe un moyen systématique de transférer les résultats du calcul des probabilités. C'est l'utilisation de la transformée de Cramer $\mathcal {C}$ définie par

$$\mathcal {C}$$ = $$\mathcal {F}$$ologo$$\mathcal {L}$$  ,

où $\mathcal {L}$ désigne la transformée de Laplace, $\mathcal {F}$ la transformée de Fenchel. On a par exemple :

$$\mathcal {C}$$($$\mathcal {N}$$_m,$\sigma$) = $$\mathcal {Q}$$_m,$\sigma$  .

La transformée de Cramer a de nombreuses propriétés. Par exemple, elle transforme les convolutions de mesures en inf-convolutions de coûts convexes. Ses propriétes de continuité pour des topologies adéquates ont été étudiées dans [Aki96]. Elles permettent dans des cas assez restreints (loi log concave) le transfert immédiat des théorèmes asymptotiques du calcul de probabilités à des théorèmes analogues en optimisation.

Page WEB

On pourra consulter les pages web http://amadeus.inria.fr/TROPICAL/ pour obtenir des informations supplémentaires.

Footnotes

... rationnelles

Elles sont caractérisés par leurs accroissements ultimement périodiques.

Précédent : Fondements scientifiques Remonter : Fondements scientifiques Suivant : Commande Stochastique