stringtranslate.com

Método de puntos interiores

Ejemplo de búsqueda de una solución. Las líneas azules indican las restricciones y los puntos rojos indican soluciones iteradas.

Los métodos de punto interior (también denominados métodos de barrera o IPM ) son algoritmos para resolver problemas de optimización convexa lineal y no lineal . Los IPM combinan dos ventajas de los algoritmos conocidos anteriormente:

A diferencia del método simplex, que recorre el límite de la región factible, y del método elipsoide, que limita la región factible desde afuera , un IPM alcanza la mejor solución al recorrer el interior de la región factible , de ahí el nombre.

Historia

El matemático soviético II Dikin descubrió un método de punto interior en 1967. [1] El método se reinventó en los EE. UU. a mediados de la década de 1980. En 1984, Narendra Karmarkar desarrolló un método para programación lineal llamado algoritmo de Karmarkar , [2] que se ejecuta en tiempo demostrablemente polinomial ( operaciones en números de L bits, donde n es el número de variables y constantes), y también es muy eficiente en la práctica. El artículo de Karmarkar creó una oleada de interés en los métodos de punto interior. Dos años más tarde, James Renegar inventó el primer método de punto interior de seguimiento de trayectorias , con tiempo de ejecución . El método se extendió más tarde de problemas de optimización lineal a convexos, basándose en una función de barrera autoconcordante utilizada para codificar el conjunto convexo . [3]

Cualquier problema de optimización convexa se puede transformar en minimizar (o maximizar) una función lineal sobre un conjunto convexo mediante la conversión a la forma de epígrafe . [4] : 143  La idea de codificar el conjunto factible utilizando una barrera y diseñar métodos de barrera fue estudiada por Anthony V. Fiacco, Garth P. McCormick y otros a principios de la década de 1960. Estas ideas se desarrollaron principalmente para la programación no lineal general , pero luego se abandonaron debido a la presencia de métodos más competitivos para esta clase de problemas (por ejemplo, programación cuadrática secuencial ).

Yurii Nesterov y Arkadi Nemirovski idearon una clase especial de barreras de este tipo que se pueden utilizar para codificar cualquier conjunto convexo. Garantizan que el número de iteraciones del algoritmo está limitado por un polinomio en la dimensión y la precisión de la solución. [5] [3]

La clase de métodos de seguimiento de trayectorias de punto interior primal-dual se considera la más exitosa. El algoritmo predictor-corrector de Mehrotra proporciona la base para la mayoría de las implementaciones de esta clase de métodos. [6]

Definiciones

Se nos da un programa convexo de la forma: donde f es una función convexa y G es un conjunto convexo . Sin pérdida de generalidad, podemos suponer que el objetivo f es una función lineal . Por lo general, el conjunto convexo G se representa mediante un conjunto de desigualdades convexas e igualdades lineales; las igualdades lineales se pueden eliminar utilizando álgebra lineal, por lo que para simplificar suponemos que solo hay desigualdades convexas, y el programa se puede describir de la siguiente manera, donde g i son funciones convexas: Suponemos que las funciones de restricción pertenecen a alguna familia (por ejemplo, funciones cuadráticas), de modo que el programa se puede representar mediante un vector finito de coeficientes (por ejemplo, los coeficientes de las funciones cuadráticas). La dimensión de este vector de coeficientes se denomina tamaño del programa. Un solucionador numérico para una familia dada de programas es un algoritmo que, dado el vector de coeficientes, genera una secuencia de soluciones aproximadas x t para t = 1, 2, ..., utilizando un número finito de operaciones aritméticas. Un solucionador numérico se denomina convergente si, para cualquier programa de la familia y cualquier ε positivo > 0, existe algún T (que puede depender del programa y de ε ) tal que, para cualquier t > T , la solución aproximada x t es ε-aproximada, es decir: donde es la solución óptima. Un solucionador se denomina polinómico si el número total de operaciones aritméticas en los primeros T pasos es como máximo

poly(tamaño del problema) * log( V / ε ),

donde V es una constante dependiente de los datos, por ejemplo, la diferencia entre el valor más grande y el más pequeño del conjunto factible. En otras palabras, V / ε es la "precisión relativa" de la solución (la precisión con respecto al coeficiente más grande). log( V / ε ) representa la cantidad de "dígitos de precisión". Por lo tanto, un solucionador es "polinómico" si cada dígito adicional de precisión requiere una cantidad de operaciones que es polinómica en el tamaño del problema.

Tipos

Los tipos de métodos de puntos interiores incluyen:

Métodos de seguimiento de trayectorias

Idea

Dado un programa de optimización convexa (P) con restricciones, podemos convertirlo en un programa sin restricciones añadiendo una función de barrera . En concreto, sea b una función convexa suave, definida en el interior de la región factible G , tal que para cualquier secuencia { x j en interior(G)} cuyo límite esté en la frontera de G : . También suponemos que b no es degenerada, es decir: es definida positiva para todo x en interior(G). Ahora, consideremos la familia de programas:

( P t ) minimizar t * f(x) + b(x)

Técnicamente, el programa está restringido, ya que b está definido solo en el interior de G . Pero prácticamente, es posible resolverlo como un programa sin restricciones, ya que cualquier solucionador que intente minimizar la función no se acercará al límite, donde b tiende al infinito. Por lo tanto, ( P t ) tiene una solución única: denotámosla por x *( t ). La función x * es una función continua de t , que se llama camino central . Todos los puntos límite de x *, cuando t tiende al infinito, son soluciones óptimas del programa original (P).

Un método de seguimiento de trayectorias es un método de seguimiento de la función x * a lo largo de una cierta secuencia creciente t1 , t2 , ..., es decir: calcular una aproximación suficientemente buena x i al punto x *( t i ), de modo que la diferencia x i - x *( t i ) se acerque a 0 cuando i se aproxima al infinito; entonces la secuencia x i se aproxima a la solución óptima de (P). Esto requiere especificar tres cosas:

El principal desafío para demostrar que el método es polinómico es que, a medida que el parámetro de penalización crece, la solución se acerca al límite y la función se vuelve más empinada. El tiempo de ejecución de los solucionadores como el método de Newton se hace más largo y es difícil demostrar que el tiempo de ejecución total es polinómico.

Renegar [7] y Gonzaga [8] demostraron que una instancia específica de un método de seguimiento de trayectoria es politemporal:

Demostraron que, en este caso, la diferencia x i - x *( t i ) permanece como máximo 0,01, y f( x i ) - f* es como máximo 2* m / t i . Por lo tanto, la precisión de la solución es proporcional a 1/ t i , por lo que para agregar un solo dígito de precisión, es suficiente multiplicar t i por 2 (o cualquier otro factor constante), lo que requiere O(sqrt( m )) pasos de Newton. Como cada paso de Newton requiere O( mn 2 ) operaciones, la complejidad total es O( m 3/2 n 2 ) operaciones por dígito de precisión.

Yuri Nesterov extendió la idea de los programas lineales a los no lineales. Observó que la propiedad principal de la barrera logarítmica, utilizada en las demostraciones anteriores, es que es autoconcordante con un parámetro de barrera finito. Por lo tanto, muchas otras clases de programas convexos se pueden resolver en politiempo utilizando un método de seguimiento de trayectorias, si podemos encontrar una función de barrera autoconcordante adecuada para su región factible. [3] : Sec.1 

Detalles

Se nos presenta un problema de optimización convexa (P) en "forma estándar":

minimizar c T x st x en G ,

donde G es convexo y cerrado. También podemos suponer que G está acotado (podemos acotarlo fácilmente añadiendo una restricción | x |≤ R para un R suficientemente grande ). [3] : Sec.4 

Para utilizar el método de puntos interiores, necesitamos una barrera autoconcordante para G . Sea b una barrera autoconcordante M para G , donde M ≥1 es el parámetro de autoconcordancia. Suponemos que podemos calcular de manera eficiente el valor de b , su gradiente y su hessiano , para cada punto x en el interior de G .

Para cada t > 0, definimos el objetivo penalizado f t (x) := c T x + b( x ) . Definimos la ruta de minimizadores mediante: x*(t) := arg min f t (x) . Aproximamos esta ruta a lo largo de una secuencia creciente t i . La secuencia se inicializa mediante un determinado procedimiento de inicialización de dos fases no trivial. Luego, se actualiza de acuerdo con la siguiente regla: .

Para cada t i , encontramos un mínimo aproximado de f ti , denotado por x i . El mínimo aproximado se elige para satisfacer la siguiente "condición de proximidad" (donde L es la tolerancia de la trayectoria ):

.

Para hallar x i +1 , comenzamos con x i y aplicamos el método de Newton amortiguado . Aplicamos varios pasos de este método, hasta que se satisface la "relación de proximidad" anterior. El primer punto que satisface esta relación se denota por x i +1 . [3] : Sec.4 

Convergencia y complejidad

La tasa de convergencia del método viene dada por la siguiente fórmula, para cada i : [3] : Prop.4.4.1 

Si tomamos , el número de pasos de Newton necesarios para ir de x i a x i +1 es como máximo un número fijo, que depende sólo de r y L . En particular, el número total de pasos de Newton necesarios para encontrar una solución ε -aproximada (es decir, encontrar x en G tal que c T x - c* ≤ ε ) es como máximo: [3] : Teoría 4.4.1 

donde el factor constante O(1) depende únicamente de r y L. El número de pasos de Newton necesarios para el procedimiento de inicialización de dos pasos es como máximo: [3] : Teoría 4.5.1 

[ aclaración necesaria ]

donde el factor constante O(1) depende únicamente de r y L , y , y es algún punto en el interior de G . En general, la complejidad general de Newton para encontrar una solución ε -aproximada es como máximo

, donde V es una constante dependiente del problema: .

Cada paso de Newton requiere O( n 3 ) operaciones aritméticas.

Inicialización: métodos de fase I

Para inicializar los métodos de seguimiento de trayectorias, necesitamos un punto en el interior relativo de la región factible G . En otras palabras: si G se define por las desigualdades g i ( x ) ≤ 0, entonces necesitamos algún x para el cual g i ( x ) < 0 para todo i en 1,..., m . Si no tenemos dicho punto, necesitamos encontrar uno usando el llamado método de fase I . [4] : 11.4  Un método simple de fase I es resolver el siguiente programa convexo: Denotemos la solución óptima por x*, s *.

Para este programa es fácil obtener un punto interior: podemos tomar arbitrariamente x = 0 y tomar s como cualquier número mayor que max( f 1 (0),..., f m (0)). Por lo tanto, se puede resolver utilizando métodos de punto interior. Sin embargo, el tiempo de ejecución es proporcional a log(1/ s *). A medida que s* se acerca a 0, se hace cada vez más difícil encontrar una solución exacta al problema de la fase I y, por lo tanto, más difícil decidir si el problema original es factible.

Consideraciones prácticas

Las garantías teóricas suponen que el parámetro de penalización se incrementa a una tasa de , por lo que el número de pasos de Newton requeridos en el peor de los casos es . En teoría, si μ es mayor (por ejemplo, 2 o más), entonces el número de pasos de Newton requeridos en el peor de los casos está en . Sin embargo, en la práctica, un μ mayor conduce a una convergencia mucho más rápida. Estos métodos se denominan métodos de paso largo . [3] : Sec.4.6  En la práctica, si μ está entre 3 y 100, entonces el programa converge dentro de 20-40 pasos de Newton, independientemente del número de restricciones (aunque el tiempo de ejecución de cada paso de Newton, por supuesto, crece con el número de restricciones). El valor exacto de μ dentro de este rango tiene poco efecto en el rendimiento. [4] : cap.11 

Métodos de reducción de potencial

Para los métodos de reducción de potencial, el problema se presenta en forma cónica : [3] : Sec.5 

minimizar c T x st x en {b+L} ∩ K ,

donde b es un vector en R n , L es un subespacio lineal en R n (por lo que b + L es un plano afín ), y K es un cono convexo puntiagudo cerrado con un interior no vacío. Todo programa convexo se puede convertir a la forma cónica. Para utilizar el método de reducción de potencial (específicamente, la extensión del algoritmo de Karmarkar a la programación convexa), necesitamos las siguientes suposiciones: [3] : Sec.6 

Los supuestos A, B y D son necesarios en la mayoría de los métodos de puntos interiores. El supuesto C es específico del enfoque de Karmarkar; se puede aliviar utilizando un "valor objetivo deslizante". Es posible reducir aún más el programa al formato de Karmarkar :

minimizar s T x st x en M ∩ K y e T x = 1

donde M es un subespacio lineal de en R n , y el valor objetivo óptimo es 0. El método se basa en la siguiente función potencial escalar :

v ( x ) = F ( x ) + M ln ( s T x )

donde F es la barrera M -autoconcordante para el cono factible. Es posible demostrar que, cuando x es estrictamente factible y v ( x ) es muy pequeña (- muy negativa), x es aproximadamente óptima. La idea del método de reducción de potencial es modificar x de modo que el potencial en cada iteración caiga al menos en una constante fija X (específicamente, X =1/3-ln(4/3)). Esto implica que, después de i iteraciones, la diferencia entre el valor objetivo y el valor objetivo óptimo es como máximo V * exp(- i X / M ), donde V es una constante dependiente de los datos. Por lo tanto, el número de pasos de Newton necesarios para una solución ε -aproximada es como máximo .

Obsérvese que en los métodos de seguimiento de trayectorias la expresión es en lugar de M , lo que es mejor en teoría. Pero en la práctica, el método de Karmarkar permite dar pasos mucho mayores hacia el objetivo, por lo que puede converger mucho más rápido que las garantías teóricas.

Métodos primal-duales

La idea del método primal-dual es fácil de demostrar para la optimización no lineal restringida . [9] [10] Para simplificar, considere el siguiente problema de optimización no lineal con restricciones de desigualdad:

Este problema de optimización restringido por desigualdades se resuelve convirtiéndolo en una función objetivo sin restricciones cuyo mínimo esperamos encontrar de manera eficiente. En concreto, la función de barrera logarítmica asociada con (1) es

Aquí hay un pequeño escalar positivo, a veces llamado "parámetro de barrera". Como converge a cero, el mínimo de debería converger a una solución de (1).

El gradiente de una función diferenciable se denota por . El gradiente de la función barrera es

Además de la variable original ("primal"), introducimos una variable dual inspirada en el multiplicador de Lagrange.

La ecuación (4) a veces se denomina condición de "complementariedad perturbada", por su parecido con la "holgura complementaria" en las condiciones KKT .

Intentamos encontrar aquellos para los que el gradiente de la función barrera es cero.

Sustituyendo (4) en (3), obtenemos una ecuación para el gradiente: donde la matriz es el jacobiano de las restricciones .

La intuición detrás de (5) es que el gradiente de debe estar en el subespacio abarcado por los gradientes de las restricciones. La "complementariedad perturbada" con (4) pequeña puede entenderse como la condición de que la solución debe estar cerca del límite o que la proyección del gradiente sobre la normal del componente de restricción debe ser casi cero.

Sea la dirección de búsqueda para actualizar iterativamente . Aplicando el método de Newton a (4) y (5), obtenemos una ecuación para :

donde es la matriz hessiana de , es una matriz diagonal de , y es la matriz diagonal de .

Debido a (1), (4) la condición

Deben aplicarse en cada paso. Esto se puede hacer eligiendo lo apropiado :

Trayectoria de los iterados de x utilizando el método del punto interior.

Tipos de programas convexos que se pueden resolver mediante métodos de puntos interiores

A continuación se presentan algunos casos especiales de programas convexos que pueden resolverse de manera eficiente mediante métodos de puntos interiores. [3] : Sección 10 

Programas lineales

Consideremos un programa lineal de la forma: Podemos aplicar métodos de seguimiento de trayectorias con la barrera La función es autoconcordante con el parámetro M = m (el número de restricciones). Por lo tanto, el número de pasos de Newton necesarios para el método de seguimiento de trayectorias es O( mn 2 ), y la complejidad total del tiempo de ejecución es O( m 3/2 n 2 ). [ aclaración necesaria ]

Programas cuadráticos con restricciones cuadráticas

Dado un programa cuadrático restringido cuadráticamente de la forma: donde todas las matrices A j son matrices semidefinidas positivas . Podemos aplicar métodos de seguimiento de trayectorias con la barrera La función es una barrera autoconcordante con parámetro M = m . La complejidad de Newton es O( (m+n)n 2 ), y la complejidad total en tiempo de ejecución es O( m 1/2 (m+n) n 2 ).

yopagaproximación de norma

Consideremos un problema de la forma donde cada uno es un vector, cada uno es un escalar y es una norma L p con Después de convertir a la forma estándar, podemos aplicar métodos de seguimiento de trayectorias con una barrera autoconcordante con parámetro M = 4 m . La complejidad de Newton es O( (m+n)n 2 ), y la complejidad total en tiempo de ejecución es O( m 1/2 (m+n) n 2 ).

Programas geométricos

Considere el problema

Existe una barrera autoconcordante con parámetro 2 k + m . El método de seguimiento de trayectoria tiene una complejidad de Newton O( mk 2 + k 3 + n 3 ) y una complejidad total O(( k+m ) 1/2 [ mk 2 + k 3 + n 3 ]).

Programas semidefinidos

Los métodos de puntos interiores se pueden utilizar para resolver programas semidefinidos. [3] : Sec.11 

Véase también

Referencias

  1. ^ Dikin, II (1967). "Solución iterativa de problemas de programación lineal y cuadrática". Dokl. Akad. Nauk SSSR . 174 (1): 747–748.
  2. ^ Karmarkar, N. (1984). "Un nuevo algoritmo de tiempo polinomial para programación lineal" (PDF) . Actas del decimosexto simposio anual de la ACM sobre teoría de la computación – STOC '84 . p. 302. doi : 10.1145/800057.808695 . ISBN . 0-89791-133-4. Archivado desde el original (PDF) el 28 de diciembre de 2013.
  3. ^ abcdefghijklm Arkadi Nemirovsky (2004). Métodos de tiempo polinomial de punto interior en programación convexa.
  4. ^ abc Boyd, Stephen; Vandenberghe, Lieven (2004). Optimización convexa . Cambridge: Cambridge University Press . ISBN 978-0-521-83378-3.Sr. 2061575  .
  5. ^ Wright, Margaret H. (2004). "La revolución del punto interior en la optimización: Historia, desarrollos recientes y consecuencias duraderas". Boletín de la American Mathematical Society . 42 : 39–57. doi : 10.1090/S0273-0979-04-01040-7 . MR  2115066.
  6. ^ Potra, Florian A.; Stephen J. Wright (2000). "Métodos de puntos interiores". Revista de Matemática Computacional y Aplicada . 124 (1–2): 281–302. doi : 10.1016/S0377-0427(00)00433-7 .
  7. ^ ab Renegar, James (1 de enero de 1988). "Un algoritmo de tiempo polinomial, basado en el método de Newton, para programación lineal". Programación matemática . 40 (1): 59–93. doi :10.1007/BF01580724. ISSN  1436-4646.
  8. ^ ab Gonzaga, Clovis C. (1989), Megiddo, Nimrod (ed.), "Un algoritmo para resolver problemas de programación lineal en operaciones O(n3L)", Progress in Mathematical Programming: Interior-Point and Related Methods , Nueva York, NY: Springer, págs. 1–28, doi :10.1007/978-1-4613-9617-8_1, ISBN 978-1-4613-9617-8, consultado el 22 de noviembre de 2023
  9. ^ Mehrotra, Sanjay (1992). "Sobre la implementación de un método de punto interior dual primario". Revista SIAM sobre optimización . 2 (4): 575–601. doi :10.1137/0802028.
  10. ^ Wright, Stephen (1997). Métodos de puntos interiores primal-dual . Filadelfia, PA: SIAM. ISBN 978-0-89871-382-4.