Articulo de referencia

Método de subgradiente

Los métodos de subgradiente son métodos de optimización convexa que utilizan subderivadas . Desarrollados originalmente por Naum Z. Shor y otros en las décadas de 1960 y 1970, l...

Los métodos de subgradiente son métodos de optimización convexa que utilizan subderivadas . Desarrollados originalmente por Naum Z. Shor y otros en las décadas de 1960 y 1970, los métodos de subgradiente convergen incluso cuando se aplican a una función objetivo no diferenciable. Cuando la función objetivo es diferenciable, los métodos de subgradiente para problemas sin restricciones utilizan la misma dirección de búsqueda que el método de descenso de gradiente .

Los métodos de subgradiente son más lentos que el método de Newton cuando se aplican para minimizar funciones convexas dos veces continuamente diferenciables. Sin embargo, el método de Newton no converge en problemas que presentan discontinuidades no diferenciables.

En los últimos años, se han propuesto algunos métodos de punto interior para problemas de minimización convexa, pero los métodos de proyección de subgradiente y los métodos de descenso de haces relacionados siguen siendo competitivos. Para problemas de minimización convexa con un número muy grande de dimensiones, los métodos de proyección de subgradiente son adecuados, ya que requieren poco almacenamiento.

Los métodos de proyección de subgradientes se aplican con frecuencia a problemas a gran escala mediante técnicas de descomposición. Estos métodos de descomposición suelen permitir un método distribuido sencillo para resolver un problema.

Reglas clásicas de subgradiente

DejarF:RnorteR{\displaystyle f:\mathbb {R} ^{n}\to \mathbb {R} }sea ​​una función convexa con dominioRnorte.{\displaystyle \mathbb {R} ^{n}.} Un método clásico de subgradiente itera incógnita(k+1)=incógnita(k)αkgramo(k) {\displaystyle x^{(k+1)}=x^{(k)}-\alpha _{k}g^{(k)}\ } dóndegramo(k){\displaystyle g^{(k)}}denota cualquier subgradiente deF {\displaystyle f\ }enincógnita(k), {\displaystyle x^{(k)},\ }yincógnita(k){\displaystyle x^{(k)}}es elkth{\displaystyle k^{th}}iterar deincógnita.{\displaystyle x.} SiF {\displaystyle f\ }Si es diferenciable, entonces su único subgradiente es el vector gradiente.F{\displaystyle \nabla f}mismo. Puede suceder quegramo(k){\displaystyle -g^{(k)}}no es una dirección decente paraF {\displaystyle f\ } enincógnita(k).{\displaystyle x^{(k)}.} Por lo tanto, mantenemos una lista.Fbmist {\displaystyle f_{\rm {mejor}}\ }que registra el valor más bajo de la función objetivo encontrado hasta el momento, es decir Fbmist(k)=min{Fbmist(k1),F(incógnita(k))}.{\displaystyle f_{\rm {best}}^{(k)}=\min\{f_{\rm {best}}^{(k-1)},f(x^{(k)})\}.}

Reglas del tamaño del paso

Los métodos de subgradiente utilizan muchos tipos diferentes de reglas de tamaño de paso. Este artículo describe cinco reglas clásicas de tamaño de paso para las que se conocen demostraciones de convergencia:

  • tamaño de paso constante,αk=α.{\displaystyle \alpha _{k}=\alpha .}
  • longitud de paso constante,αk=γ/gramo(k)2,{\displaystyle \alpha _ {k}=\gamma /\lVert g^{(k)}\rVert _ {2},}lo cual daincógnita(k+1)incógnita(k)2=γ.{\displaystyle \lVert x^{(k+1)}-x^{(k)}\rVert _{2}=\gamma .}
  • Tamaño de paso cuadrado sumable pero no sumable, es decir, cualquier tamaño de paso que satisfagaαk0,k=1αk2<,k=1αk=.{\displaystyle \alpha _{k}\geq 0,\qquad \sum _{k=1}^{\infty }\alpha _{k}^{2}<\infty ,\qquad \sum _{k=1}^{\infty }\alpha _{k}=\infty .}
  • Decreciente no sumable, es decir, cualquier tamaño de paso que satisfagaαk0,límitekαk=0,k=1αk=.{\displaystyle \alpha _{k}\geq 0,\qquad \lim _{k\to \infty }\alpha _{k}=0,\qquad \sum _{k=1}^{\infty }\alpha _{k}=\infty .}
  • Longitudes de paso decrecientes no sumables, es decirαk=γk/gramo(k)2,{\displaystyle \alpha _{k}=\gamma _{k}/\lVert g^{(k)}\rVert _{2},}dóndeγk0,límitekγk=0,k=1γk=.{\displaystyle \gamma _{k}\geq 0,\qquad \lim _{k\to \infty }\gamma _{k}=0,\qquad \sum _{k=1}^{\infty }\gamma _{k}=\infty .}

Para las cinco reglas, los tamaños de paso se determinan "fuera de línea", antes de iterar el método; los tamaños de paso no dependen de las iteraciones previas. Esta propiedad "fuera de línea" de los métodos de subgradiente difiere de las reglas de tamaño de paso "en línea" utilizadas para los métodos de descenso para funciones diferenciables: muchos métodos para minimizar funciones diferenciables satisfacen las condiciones suficientes de Wolfe para la convergencia, donde los tamaños de paso generalmente dependen del punto actual y de la dirección de búsqueda actual. Una discusión extensa de las reglas de tamaño de paso para métodos de subgradiente, incluidas las versiones incrementales, se encuentra en los libros de Bertsekas [ 1 ] y de Bertsekas, Nedic y Ozdaglar. [ 2 ]

Resultados de convergencia

Para subgradientes de longitud de paso constante y escalados con norma euclidiana igual a uno, el método de subgradiente converge a una aproximación arbitrariamente cercana al valor mínimo, es decir

límitekFbmist(k)F<ϵ{\displaystyle \lim _{k\to \infty }f_{\rm {best}}^{(k)}-f^{*}<\epsilon }por un resultado de Shor . [ 3 ]

Estos métodos clásicos de subgradiente tienen un rendimiento deficiente y ya no se recomiendan para uso general. [ 4 ] [ 5 ] Sin embargo, todavía se utilizan ampliamente en aplicaciones especializadas porque son sencillos y se pueden adaptar fácilmente para aprovechar la estructura especial del problema en cuestión.

Métodos de proyección de subgradiente y de haces

Durante la década de 1970, Claude Lemaréchal y Phil Wolfe propusieron los "métodos de haces" de descenso para problemas de minimización convexa. [ 6 ] El significado del término "métodos de haces" ha cambiado significativamente desde entonces. Kiwiel proporcionó versiones modernas y un análisis completo de convergencia. [ 7 ] Los métodos de haces contemporáneos a menudo utilizan reglas de " control de nivel " para elegir tamaños de paso, desarrollando técnicas a partir del método de "proyección de subgradiente" de Boris T. Polyak (1969). Sin embargo, hay problemas en los que los métodos de haces ofrecen poca ventaja sobre los métodos de proyección de subgradiente. [ 4 ] [ 5 ]

Optimización con restricciones

Subgradiente proyectado

Una extensión del método del subgradiente es el método del subgradiente proyectado , que resuelve el problema de optimización con restricciones.

minimizarF(incógnita) {\displaystyle f(x)\ }sujeto aincógnitado{\displaystyle x\in {\mathcal {C}}}

dóndedo{\displaystyle {\mathcal {C}}}es un conjunto convexo . El método del subgradiente proyectado utiliza la iteración incógnita(k+1)=PAG(incógnita(k)αkgramo(k)){\displaystyle x^{(k+1)}=P\left(x^{(k)}-\alpha _{k}g^{(k)}\right)} dóndePAG{\displaystyle P}es proyección endo{\displaystyle {\mathcal {C}}}ygramo(k){\displaystyle g^{(k)}}es cualquier subgradiente deF {\displaystyle f\ }enincógnita(k).{\displaystyle x^{(k)}.}

Restricciones generales

El método del subgradiente puede extenderse para resolver el problema con restricciones de desigualdad.

minimizarF0(incógnita) {\displaystyle f_{0}(x)\ }sujeto aFi(incógnita)0,i=1,,metro{\displaystyle f_{i}(x)\leq 0,\quad i=1,\ldots ,m}

dóndeFi{\displaystyle f_{i}}son convexas. El algoritmo toma la misma forma que el caso sin restricciones. incógnita(k+1)=incógnita(k)αkgramo(k) {\displaystyle x^{(k+1)}=x^{(k)}-\alpha _{k}g^{(k)}\ } dóndeαk>0{\displaystyle \alpha _{k}>0}es un tamaño de paso, ygramo(k){\displaystyle g^{(k)}}es un subgradiente de la función objetivo o de una de las funciones de restricción enincógnita. {\displaystyle x.\ } Llevar gramo(k)={F0(incógnita) si Fi(incógnita)0i=1metroFj(incógnita) para algunos j de tal manera que Fj(incógnita)>0{\displaystyle g^{(k)}={\begin{cases}\partial f_{0}(x)&{\text{ if }}f_{i}(x)\leq 0\;\forall i=1\dots m\\\partial f_{j}(x)&{\text{ for some }}j{\text{ such that }}f_{j}(x)>0\end{cases}}} dóndeF{\displaystyle \partial f}denota el subgradiente deF. {\displaystyle f.\ } Si el punto actual es factible, el algoritmo utiliza un subgradiente objetivo; si el punto actual no es factible, el algoritmo elige un subgradiente de cualquier restricción violada.

Véase también

Referencias

  1. Bertsekas, Dimitri P. (2015). Algoritmos de optimización convexa (Segunda  edición). Belmont, MA: Athena Scientific. ISBN 978-1-886529-28-1.
  2. Bertsekas, Dimitri P.; Nedic, Angelia; Ozdaglar, Asuman (2003). Análisis convexo y optimización (Segunda ed.). Belmont, MA.: Athena Scientific. ISBN  1-886529-45-0.
  3. La convergencia aproximada del método de subgradiente de paso constante (escalado) se indica en el Ejercicio 6.3.14(a) de Bertsekas (página 636): Bertsekas, Dimitri P. (1999). Nonlinear Programming (Segunda edición). Cambridge, MA.: Athena Scientific. ISBN  1-886529-00-0.En la página 636, Bertsekas atribuye este resultado a Shor: Shor, Naum Z. (1985). Minimization Methods for Non-differentiable Functions . Springer-Verlag . ISBN 0-387-12763-1.
  4. ^ Lemaréchal , Claude (2001). "Relajación lagrangiana". En Michael Jünger y Denis Naddef (ed.). Optimización combinatoria computacional: artículos de la escuela de primavera celebrada en Schloß Dagstuhl, del 15 al 19 de mayo de 2000 . Apuntes de conferencias sobre informática. vol. 2241. Berlín: Springer-Verlag. págs. 112-156 . doi : 10.1007/3-540-45586-8_4 . ISBN   3-540-42877-1. MR 1900016 . S2CID 9048698 .  
  5. 1 2 Kiwiel, Krzysztof C.; Larsson, Torbjörn; Lindberg, P. O. (agosto de 2007). "Relajación lagrangiana mediante métodos de subgradiente de paso de bola" (PDF) . Matemáticas de la Investigación Operativa . 32 (3): 669– 686. doi : 10.1287/moor.1070.0261 . MR 2348241 .   
  6. Bertsekas, Dimitri P. (1999). Programación no lineal (Segunda edición). Cambridge, MA: Athena Scientific. ISBN  1-886529-00-0.
  7. Kiwiel, Krzysztof (1985). Métodos de descenso para la optimización no diferenciable . Berlín: Springer Verlag . pág. 362. ISBN  978-3540156420. SR 0797754 . 

Lecturas adicionales

  • Bertsekas, Dimitri P. (1999). Programación no lineal . Belmont, MA: Athena Scientific. ISBN 1-886529-00-0.
  • Bertsekas, Dimitri P.; Nedic, Angelia; Ozdaglar, Asuman (2003). Análisis convexo y optimización (Segunda  edición). Belmont, MA: Athena Scientific. ISBN 1-886529-45-0.
  • Bertsekas, Dimitri P. (2015). Algoritmos de optimización convexa . Belmont, MA: Athena Scientific. ISBN 978-1-886529-28-1.
  • Shor, Naum Z. (1985). Métodos de minimización para funciones no diferenciables . Springer-Verlag . ISBN 0-387-12763-1.
  • Ruszczyński, Andrzej (2006). Optimización no lineal . Princeton, NJ: Princeton University Press . págs.  xii+454. ISBN 978-0691119151MR 2199043 .​ 
  • EE364A y EE364B , la secuencia de cursos de optimización convexa de Stanford.