Algorithme du lagrangien augmenté
Les algorithmes du lagrangien augmenté sont une certaine classe d'algorithmes pour résoudre des problèmes d'optimisation sous contraintes. Elles présentent des similitudes avec les méthodes de pénalité dans le sens où elles remplacent un problème d'optimisation sous contraintes par une série de problèmes sans contrainte et ajoutent un terme de pénalité à l'objectif ; la différence est qu'une méthode du lagrangien augmenté ajoute encore un autre terme, conçu pour agir comme un multiplicateur de Lagrange. Le Lagrangien augmenté est apparenté, mais non identique à la méthode des multiplicateurs de Lagrange.
Vu différemment, l'objectif sans contrainte est le lagrangien du problème contraint, avec un terme de pénalité supplémentaire (l'augmentation).
La méthode était à l'origine connue sous le nom de méthode des multiplicateurs et a été beaucoup étudiée dans les années 1970 et 1980 comme une bonne alternative aux méthodes de pénalité. Il a été discuté pour la première fois par Magnus Hestenes[1] et par Michael J. D. Powell en 1969[2]. La méthode a été étudiée par R. Tyrrell Rockafellar en relation avec la dualité de Fenchel (en), en particulier en relation avec les méthodes des points proximaux, la régularisation Moreau-Yosida et les opérateurs monotones maximaux : ces méthodes ont été utilisées en optimisation structurelle. La méthode a également été étudiée par Dimitri Bertsekas, notamment dans son livre de 1982[3], ainsi que des extensions impliquant des fonctions de régularisation non quadratiques, telles que la régularisation entropique, qui donne naissance à la « méthode exponentielle des multiplicateurs », une méthode qui gère les contraintes d'inégalité avec une fonction lagrangienne augmentée deux fois dérivable.
Depuis les années 1970, la programmation quadratique séquentielle (PQS) et les méthodes de points intérieurs (MPI) ont suscité une attention croissante, en partie parce qu'elles utilisent plus facilement des sous-programmes de matrice creuse à partir de bibliothèques de logiciels numériques, et en partie parce que les MPI ont prouvé des résultats de complexité via la théorie de fonctions auto-concordantes (en). La méthode du Lagrangien augmenté a été rajeunie par les systèmes d'optimisation LANCELOT, ALGENCAN [4] et AMPL, qui ont permis d'utiliser des techniques de matrices creuses sur des problèmes apparemment denses mais "partiellement séparables". La méthode est toujours utile pour certains problèmes. Vers 2007, il y a eu une résurgence des méthodes du lagrangien augmenté dans des domaines tels que le débruitage à variation totale et la détection compressée. En particulier, une variante de la méthode lagrangienne augmentée standard qui utilise des mises à jour partielles (similaire à la méthode de Gauss-Seidel pour résoudre des équations linéaires) connue sous le nom d'algorithme des directions alternées des multiplicateurs ou MDAM a attiré l'attention.
Méthode générale
modifierOn cherche à résoudre le problème sous contraintes suivant :
sous conditions
où désigne les indices des contraintes d'égalité. Ce problème peut être résolu comme une série de problèmes de minimisation sans contrainte. Pour référence, on liste d'abord la ke étape de l'approche de la méthode de pénalité :
La méthode de pénalité résout ce problème, puis à l'itération suivante, elle résout à nouveau le problème en utilisant une plus grande valeur de (et en utilisant l'ancienne solution comme estimation initiale ou "démarrage à chaud").
La méthode du lagrangien augmenté utilise l'objectif non contraint suivant :
et après chaque itération, en plus de la mise à jour des , la variable est également mise à jour selon la règle
où est la solution du problème sans contrainte à la ke étape, c'est-à-dire
La variable est une estimation du multiplicateur de Lagrange, et la précision de cette estimation s'améliore à chaque étape. L' avantage majeur de la méthode est que contrairement à la méthode de pénalisation, il n'est pas nécessaire de prendre afin de résoudre le problème sous contraintes initial. Au lieu de cela, en raison de la présence du terme de multiplicateur de Lagrange, peut rester beaucoup plus petit, évitant ainsi le mauvais conditionnement. Néanmoins, il est courant dans les implémentations pratiques de projeter les estimations des multiplicateurs dans un grand ensemble borné (sauvegardes), évitant les instabilités numériques et conduisant à une forte convergence théorique.
La méthode peut être étendue pour gérer les contraintes d'inégalité .
Méthode des directions alternées
modifierLa méthode des directions alternées est une variante du schéma du lagrangien augmenté qui utilise des mises à jour partielles pour les variables duales. Cette méthode est souvent appliquée pour résoudre des problèmes tels que
Ceci est équivalent au problème contraint
Bien que ce changement puisse sembler anodin, le problème peut maintenant être envisagé en utilisant des méthodes d'optimisation sous contraintes (en particulier, la méthode du lagrangien augmenté), et la fonction objectif est séparable en x et y. La double mise à jour nécessite de résoudre une fonction de proximité en x et y en même temps ; la technique ADA permet de résoudre approximativement ce problème en résolvant d'abord pour x avec y fixe, puis en résolvant pour y avec x fixe. Plutôt que d'itérer jusqu'à avoir convergé (comme la méthode de Jacobi), l'algorithme procède directement à la mise à jour de la variable duale, puis à la répétition du processus. Ce n'est pas équivalent à la minimisation exacte, mais on peut quand même montrer que cette méthode converge vers la bonne réponse sous certaines hypothèses. Du fait de cette approximation, l'algorithme se distingue de la méthode du lagrangien augmenté pure.
L'ADA peut être considéré comme une application de l'algorithme de Douglas-Rachford, et l'algorithme de Douglas-Rachford est à son tour une application de l'algorithme du point proximal[5]. Il existe plusieurs progiciels modernes qui résolvent la poursuite de base et se variantes et utilisent l'ADA ; ces packages incluent YALL1 [6] (2009), SpaRSA [7] (2009) et SALSA [8] (2009). Il existe également des packages qui utilisent l'ADA pour résoudre des problèmes plus généraux, dont certains peuvent exploiter plusieurs cœurs de calcul SNAPVX [9] (2015), parADMM [10] (2016).
Optimisation stochastique
modifierL'optimisation stochastique considère le problème de la minimisation d'une fonction de perte avec accès à des échantillons bruités (du gradient) de la fonction. Le but est d'avoir une estimation du paramètre optimal (minimiseur) pour chaque nouvel échantillon. L'ADA est à l'origine une méthode groupée. Cependant, avec quelques modifications, elle peut également être utilisée pour l'optimisation stochastique. Étant donné que dans le cadre stochastique, on n'a accès qu'à des échantillons bruités de gradient, on utilise une approximation inexacte du lagrangien comme
où est une taille de pas temporel variable[11].
La méthode des directions alternées est une méthode populaire pour l'optimisation en ligne et distribuée à grande échelle[12], et est utilisée dans de nombreuses applications [13],[14]. L'ADA est souvent appliqué pour résoudre des problèmes régularisés, où l'optimisation et la régularisation de la fonction peuvent être effectuées localement, puis coordonnées globalement via des contraintes. Les problèmes d'optimisation régularisés sont particulièrement pertinents dans le régime de haute dimension puisque la régularisation est un mécanisme naturel pour surmonter le caractère mal-posé et pour encourager la parcimonie dans la solution optimale, par exemple, pour une matrice creuse et de faible rang. En raison de l'efficacité de l'ADA dans la résolution de problèmes régularisés, il a un bon potentiel d'optimisation stochastique en grandes dimensions.
Approches alternatives
modifier- Programmation quadratique séquentielle
- Programmation linéaire séquentielle
- Programmation séquentielle linéaire-quadratique
Articles connexes
modifierRéférences
modifier- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Augmented Lagrangian method » (voir la liste des auteurs).
- Hestenes, « Multiplier and gradient methods », Journal of Optimization Theory and Applications, vol. 4, no 5, , p. 303–320 (DOI 10.1007/BF00927673, S2CID 121584579)
- M. J. D. Powell, Optimization, New York, Academic Press, , 283–298 p. (ISBN 0-12-260650-7), « A method for nonlinear constraints in minimization problems »
- Dimitri P. Bertsekas, Constrained optimization and Lagrange multiplier methods, Athena Scientific, (1re éd. 1982)
- Andreani, Birgin, Martínez et Schuverdt, « On Augmented Lagrangian Methods with General Lower-Level Constraints », SIAM Journal on Optimization, vol. 18, no 4, , p. 1286–1309 (DOI 10.1137/060654797, lire en ligne)
- Eckstein et Bertsekas, « On the Douglas—Rachford splitting method and the proximal point algorithm for maximal monotone operators », Mathematical Programming, vol. 55, nos 1–3, , p. 293–318 (DOI 10.1007/BF01581204, S2CID 15551627, CiteSeerx 10.1.1.141.6246)
- « YALL1: Your ALgorithms for L1 », yall1.blogs.rice.edu
- « SpaRSA », www.lx.it.pt
- « (C)SALSA: A Solver for Convex Optimization Problems in Image Recovery », cascais.lx.it.pt
- « SnapVX », snap.stanford.edu
- « parADMM/engine »,
- Ouyang, H., He, N., Tran, L. et Gray, A. G, « Stochastic alternating direction method of multipliers », Proceedings of the 30th International Conference on Machine Learning, , p. 80–88
- Boyd, S., Parikh, N., Chu, E. et Peleato, B., « Distributed optimization and statistical learning via the alternating direction method of multipliers », Foundations and Trends{\textregistered} in Machine Learning, vol. 3, no 1, , p. 1–122 (DOI 10.1561/2200000016, CiteSeerx 10.1.1.360.1664)
- Esser, E., Zhang, X. et Chan, T., « A general framework for a class of first order primal-dual algorithms for convex optimization in imaging science », SIAM Journal on Imaging Sciences, vol. 3, no 4, , p. 1015–1046 (DOI 10.1137/09076934X)
- Mota, J. FC, Xavier, J. MF, Aguiar, P. MQ et Puschel, M., « Distributed ADMM for model predictive control and congestion control », Decision and Control (CDC), 2012 IEEE 51st Annual Conference O, , p. 5110–5115 (ISBN 978-1-4673-2066-5, DOI 10.1109/CDC.2012.6426141, S2CID 12128421)
Bibliographie
modifier- Dimitri P. Bertsekas, Nonlinear Programming, Belmont, Mass, Athena Scientific, (ISBN 978-1-886529-00-7)
- E. G. Birgin et J. M. Martínez, Practical Augmented Lagrangian Methods for Constrained Optimization, Philadelphia, Society for Industrial and Applied Mathematics, (ISBN 978-1-611973-35-8, DOI 10.1137/1.9781611973365
)
- Jorge Nocedal et Stephen J. Wright, Numerical Optimization, Berlin, New York, Springer-Verlag, (ISBN 978-0-387-30303-1)