Théorie du transport

En mathématiques et en économie, la théorie du transport est le nom donné à l'étude du transfert optimal de matière et à l'allocation optimale de ressources. Le problème a été formalisé par le mathématicien français Gaspard Monge en 1781[1]. D'importants développements ont été réalisés dans ce domaine pendant la Seconde Guerre mondiale par le mathématicien et économiste russe Léonid Kantorovitch[2]. Par conséquent, le problème dans sa forme actuelle est parfois baptisé problème (du transport) de Monge-Kantorovitch.

Motivation modifier

Mines et usines modifier

 
Transport de minerai

On se donne un ensemble de   mines d'où est extrait un minerai de fer, et un ensemble de   usines utilisant ce minerai comme matière première. Ces mines et ces usines ont une certaine aire. On suppose donc pour la clarté du propos que ces mines et ces usines constituent deux sous-ensembles disjoints et limités par une courbe fermée   et   du plan euclidien  . On suppose également qu'on dispose d'une fonction, que l'on appelle le coût, à savoir  , telle que   soit le coût de transport d'un transfert de minerai du site   au site  . Par souci de simplicité, on ignore le temps mis lors de ce transfert. On considère également que chaque mine ne peut fournir en minerai qu'une seule usine (pas de partage de minerai durant le transfert) et que la quantité transférée doit être intégralement distribuée à une usine donnée pour que celle-ci soit opérationnelle (les usines ne peuvent fonctionner ni en sur-capacité ni en sous-capacité).

Ayant fait ces hypothèses, un plan de transport est une bijection   — c.-à-d. un arrangement où chaque mine   alimente précisément une usine  . On désire trouver le plan de transport optimal, le plan   dont le coût total

 

est minimal vis-à-vis de tous les plans de transport possibles de   vers  .

Déplacement de livres : l'importance du type de la fonction de coût modifier

L'exemple suivant illustre l'importance de la fonction coût dans la détermination du plan de transport optimal. On suppose qu'on a   livres d'égale épaisseur sur une étagère (la droite réelle), arrangés sans espaces entre eux. On désire les réarranger toujours sans espace entre eux, mais selon un décalage égal à l'épaisseur d'un livre sur la droite. Deux candidats pour le plan de transport optimal se présentent ainsi :

  1. déplacer l'ensemble des   livres (en commençant par le plus à droite) d'une distance égale à l'épaisseur d'un livre vers la droite ; (« beaucoup de petits déplacements »)
  2. déplacer le livre le plus à gauche d'une distance égale à l'épaisseur de   livres vers la droite et laisser les autres livres en place. (« peu de grands déplacements »)

Si la fonction de coût est proportionnelle à la distance euclidienne ( ) alors ces deux candidats sont tous deux optimaux. Si, au contraire, on choisit une fonction de coût strictement convexe proportionnelle au carré de la distance euclidienne ( ), alors l'option « beaucoup de petits déplacements » devient l'unique minimiseur.

De manière intéressante, alors que les mathématiciens préfèrent travailler avec des fonctions de coût convexes, les économistes préfèrent les concaves. La justification intuitive de ce constat est que, une fois que des biens ont été chargés, disons, sur un train de marchandises, le transport de ces biens sur une distance de 200 kilomètres coûtera bien moins que le double du coût de transport sur une distance de 100 kilomètres. Les fonctions de coûts concaves représentent cette économie d'échelle.

Histoire modifier

La formulation de Monge[3] modifier

Le problème du transport optimal est formalisé pour la première fois par le mathématicien français Gaspard Monge en 1781, dans un Mémoire sur la théorie des déblais et des remblais destiné à l'Académie des Sciences[4]. Dans cet ouvrage, il considère la question du déplacement d'une quantité de terre, appelée les déblais, vers un espace qu'ils doivent occuper après leur transport, appelé les remblais. En considérant un coût de transport proportionnel à la masse déplacée ainsi qu'à la distance parcourue, il conjecture que tous les déplacements ne résulteront pas en un même coût total. Monge affirme de plus que parmi tous ces transports possibles, il en existe un, appelé transport optimal, pour lequel le coût total sera minimum. Contrairement à ce qu'il annonce, la preuve de l'existence et de l'unicité d'un tel transport n'est pas donnée dans son mémoire. Malgré quelques avancées réalisées par le mathématicien français Charles Dupin dans un mémoire de 1822[5], le problème ne connaît pas de développements importants pendant presque cent ans. En 1884, l'Académie des Sciences propose un prix à quiconque obtiendra des avancées majeures sur la question du transport optimal. C'est le mathématicien français Paul Appell qui remporte 2000 des 3000 francs proposés en 1886. Il publie en 1887 un mémoire sur ses travaux[6], dans lesquels la question de l'existence n'est toujours pas clairement formulée ni résolue.

La relaxation de Kantorovitch[7] modifier

Le mathématicien et économiste russe Leonid Kantorovitch est un des inventeurs du concept d'optimisation linéaire dans les années 1930. Cette théorie le conduit à s'intéresser au problème de transport optimal. Durant la Seconde Guerre Mondiale, il propose une formulation plus générale du problème de Monge, en n'imposant pas que le transport des masses soit nécessairement décrit par une fonction, mais plutôt par une mesure appelée plan de transport. Dans ce cadre, Kantorovich parvient à obtenir l'existence et l'unicité d'un plan de transport optimal, avec des hypothèses très faibles sur le coût de transport et la répartition des masses à transporter et des espaces les accueillant[8]. Le mathématicien russe V. N. Sudavok utilise les outils développés par son prédécesseur pour revenir au problème original de Monge, auquel il propose une solution en 1976[9]. Le mathématicien italien Luigi Ambrosio repère et corrige une erreur fondamentale dans l'argument de Sudakov en 2003[10], apportant une réponse jugée satisfaisante par la communauté mathématique à la question de l'existence et de l'unicité dans le problème de Monge.

Caractérisation et régularité[11] modifier

Dans les années 1990, le mathématicien français Yann Brenier réalise une avancée déterminante dans la caractérisation du problème pour le coût quadratique. Sous de faibles hypothèses de régularité, il montre que le transport se fait nécessairement via une fonction de transport qui est le gradient d'une fonction convexe[12]. Le mathématicien italien Luigi Cafferelli poursuit les travaux sur le coût quadratique et obtient en 2002 des résultats importants sur la régularité de la fonction de transport[13]. Poussant encore plus loin l'analyse du problème de transport, Ma, Trudinger et Wang établissent des propriétés de régularité de la fonction de transport en prenant des coûts   en 2005[14].

Formulation abstraite du problème modifier

Les formulations de Monge et de Kantorovitch modifier

La forme donnée à l'énoncé du problème du transport pourra varier quelque peu dans la littérature technique moderne en fonction des développements en géométrie riemannienne et en théorie de la mesure. L'exemple mines-usines, par sa simplicité, pourra être vu comme un bon point de départ pour aborder l'étude abstraite. Dans ce cadre, on s'autorise à considérer que les mines et les usines ne sont pas toutes nécessairement ouvertes pendant les transactions, que les mines peuvent alimenter plus d'une usine, et que chaque usine peut être ravitaillée en fer par plus d'une mine.

Soit   et   deux espaces métriques séparables tels que toute mesure de probabilité sur   (ou sur  ) soit une mesure de Radon (i.e. ce sont des espaces de Radon). Soit   une fonction mesurable au sens des boréliens. Étant donné des mesures   sur   et   sur  , la formulation de Monge du problème de transport optimal est issue de la recherche du plan de transport   qui réalise l'infimum

 

  représente la mesure image de   par  . Un plan   qui atteint cet infimum (i. e. l'infimum doit alors être appelé minimum) est appelé un « plan de transport optimal ».

Le problème de transport optimal selon la formulation de Monge peut être mal posé (car il n'y a parfois pas de   satisfaisant   : ceci se produit, par exemple, quand   est une mesure de Dirac et que   n'en est pas une).

On peut alors améliorer cette formulation en adoptant la formulation de Kantorovich du problème de transport optimal, qui consiste à trouver la mesure   sur   qui atteint l'infimum

 

  représente l'ensemble des mesures définies sur   de marginales données par   sur   et par   sur  . On peut montrer[15] qu'un minimiseur de ce problème existe toujours quand la fonction coût   est semi-continue inférieurement et que   est un ensemble de mesures à espérances et variances bornées (ce qui est assuré par la nature des espaces de Radon   et  ). (Comparer cette formulation avec la définition de la métrique de Wasserstein   sur l'espace des mesures de probabilité.)

Formulation duale modifier

Le minimum du problème de Kantorovitch est égal à

 

où la borne supérieure est à prendre parmi toutes les paires de fonctions continues bornées   et   telles que

 

Solution du problème modifier

Transport optimal sur la droite réelle modifier

Pour  , soit   l'ensemble des mesures de probabilité sur   de   moment fini. Soient   et soit  , où   est une fonction convexe.

  1. Si   n'a pas d'atome, i.e., si la fonction de répartition   de   est une fonction continue, alors   est un plan de transport optimal. Ce dernier est unique si   est strictement convexe.
  2. On a
 

Dans le cas de distributions de probabilités uniformes discrètes, et pour   une version très élémentaire du problème de transport optimal est l'inégalité de réarrangement.

Espaces de Hilbert séparables modifier

Soit   un espace de Hilbert séparable. Soit   l'ensemble des mesures de probabilité sur   de   moment fini ; soit   les   qui sont réguliers au sens gaussien : pour toute mesure gaussienne   strictement positive sur   avec  , alors   aussi.

Soit  ,  ,   pour  ,  . Alors le problème de Kantorovitch a une solution unique  , et cette solution est induite par le plan de transport optimal : i.e., il existe une carte   au sens de Borel telle que

 

De surcroît, si   est à support borné, alors

  pour  -presque tout  

pour un certain potentiel maximal de Kantorovitch   localement lipschitzien et  -concave. (Ici   représente la dérivée de Gateaux de  .)

Régularisation entropique modifier

En pratique, on considère souvent une variante du problème du transport optimal, faisant intervenir un terme de régularisation entropique, celle-ci étant plus simple à résoudre numériquement[16]. Le nouveau problème s'écrit alors, pour un certain   donné :

 

  est la divergence de Kullback-Leibler, aussi appelée entropie relative, du fait de son lien avec l'entropie de Shannon. Ce problème peut alors être résolu par l'algorithme de Sinkhorn.

Références modifier

  1. G. Monge. Mémoire sur la théorie des déblais et de remblais. Histoire de l’Académie Royale des Sciences de Paris, avec les Mémoires de Mathématique et de Physique pour la même année, pages 666–704, 1781.
  2. L. Kantorovich. On the translocation of masses. C.R. (Doklady) Acad. Sci. URSS (N.S.), 37:199–201, 1942.
  3. Etienne Ghys, « Gaspard Monge, le mémoire sur les déblais et les remblais », sur Images des mathématiques, (consulté le )
  4. Gaspard Monge, Mémoire sur la théorie des déblais et de remblais. Histoire de l’Académie Royale des Sciences de Paris, avec les Mémoires de Mathématique et de Physique pour la même année, pages 666–704, 1781.
  5. Charles Dupin, Applications de géométrie et de méchanique : à la marine aux ponts et chaussées, etc., pour faire suite aux développements de géométrie,
  6. Paul Appell, Mémoire sur les déblais et les remblais de systèmes continus ou discontinus, Mémoires présentés par divers savants à l’Académie royale des sciences de l’Institut de France... Sciences mathématiques et physiques. 1827-1914 (2e s. I-XXXV),
  7. Yann Brenier, Thierry Viéville, « La brouette de Monge ou le transport optimal », sur Images de mathématiques, (consulté le )
  8. Leonid Kantorovitch, On the transfer of masses. Dokl. Acad. Nauk. USSR, 37, 7–8, 1942.
  9. V. N. Sudakov, Geometric problems in the theory of infinite-dimensional probability distributions. Cover to cover translation of Trudy Mat. Inst. Steklov 141 (1976). Proc. Steklov Inst. Math. 2, i–v, 1–178,
  10. L. Ambrosio, Lecture Notes on Optimal Transport Problems, Mathematical Aspects of Evolving Interfaces, Springer Verlag, Lecture Notes in Mathematics (1812), 1–52,
  11. F. Santambrogio, Optimal Transport for Applied Mathematicians,
  12. Yann Brenier, Décomposition polaire et réarrangement monotone des champs de vecteurs. (French) C. R. Acad. Sci. Paris. I Math. 305 (19), 805–808, 1987
  13. Luigi Caffarelli, Constructing optimal maps for Monge’s transport problem as a limit of strictly convex costs J. Amer. Math. Soc. 15, 1–26,
  14. Ma, Trundinger, Wang, Regularity of potential functions of the optimal transportation problem. Arch. Ration. Mech. Anal., 177 (2), 151–183,
  15. L. Ambrosio, N. Gigli & G. Savaré. Gradient Flows in Metric Spaces and in the Space of Probability Measures. Lectures in Mathematics ETH Zürich, Birkhäuser Verlag, Basel. (2005)
  16. Marco Cuturi, « Sinkhorn Distances: Lightspeed Computation of Optimal Transport », Advances in Neural Information Processing Systems, Curran Associates, Inc., vol. 26,‎ (lire en ligne, consulté le )

Liens internes modifier

Liens externes modifier