Articulo de referencia

Programación de conos de segundo orden

Un programa de cono de segundo orden ( SOCP ) es un problema de optimización convexa de la forma minimizar F T incógnita {\displaystyle \ f^{T}x\ } sujeto a ‖ A i incógnita ...

Un programa de cono de segundo orden ( SOCP ) es un problema de optimización convexa de la forma

minimizar FTincógnita {\displaystyle \ f^{T}x\ }
sujeto a
Aiincógnita+bi2doiTincógnita+di,i=1,,metro{\displaystyle \lVert A_{i}x+b_{i}\rVert _{2}\leq c_{i}^{T}x+d_{i},\quad i=1,\dots ,m}
Fincógnita=gramo {\displaystyle Fx=g\ }

donde los parámetros del problema sonFRnorte, AiRnortei×norte, biRnortei, doiRnorte, diR, FRpag×norte{\displaystyle f\in \mathbb {R} ^{n},\ A_{i}\in \mathbb {R} ^{{n_{i}}\times n},\ b_{i}\in \mathbb {R} ^{n_{i}},\ c_{i}\in \mathbb {R} ^{n},\ d_{i}\in \mathbb {R} ,\ F\in \mathbb {R} ^{p\times n}}, ygramoRpag{\displaystyle g\in \mathbb {R} ^{p}}.incógnitaRnorte{\displaystyle x\in \mathbb {R} ^{n}}es la variable de optimización. incógnita2{\displaystyle \lVert x\rVert _{2}}es la norma euclidiana yT{\displaystyle ^{T}}indica transposición . [ 1 ]

El nombre "programación de cono de segundo orden" proviene de la naturaleza de las restricciones individuales, cada una de las cuales tiene la forma:

Aincógnita+b2doTincógnita+d{\displaystyle \lVert Ax+b\rVert _{2}\leq c^{T}x+d}

Cada uno de estos define un subespacio que está acotado por una desigualdad basada en una función polinómica de segundo orden definida en la variable de optimización.incógnita{\displaystyle x}Se puede demostrar que esto define un cono convexo , de ahí el nombre de " cono de segundo orden ". [ 2 ] Por definición de conos convexos, su intersección también puede ser un cono convexo, aunque no necesariamente uno que pueda definirse mediante una sola desigualdad de segundo orden. Véase más adelante para un análisis más detallado.

Los SOCP se pueden resolver mediante métodos de punto interior [ 3 ] y, en general, se pueden resolver de manera más eficiente que los problemas de programación semidefinida (SDP) [ 4 ] . Algunas aplicaciones de ingeniería de los SOCP incluyen el diseño de filtros , el diseño de pesos de matrices de antenas, el diseño de estructuras de celosía y la optimización de la fuerza de agarre en robótica [ 5 ] . Las aplicaciones en finanzas cuantitativas incluyen la optimización de carteras ; algunas restricciones de impacto de mercado , debido a que no son lineales, no se pueden resolver mediante programación cuadrática , pero se pueden formular como problemas SOCP [ 6 ] [ 7 ] [ 8 ] .

Conos de segundo orden

El cono estándar o unitario de segundo orden de dimensiónnorte+1{\displaystyle n+1}se define como

donorte+1={[incógnitat]|incógnitaRnorte,tR,incógnita2t}{\displaystyle {\mathcal {C}}_{n+1}=\left\{{\begin{bmatrix}x\\t\end{bmatrix}}{\Bigg |}x\in \mathbb {R} ^{n},t\in \mathbb {R} ,\|x\|_{2}\leq t\right\}}.

El cono de segundo orden también se conoce con los nombres de cono cuadrático , cono de helado [ 9 ] o cono de Lorentz . Por ejemplo, el cono estándar de segundo orden enR3{\displaystyle \mathbb {R} ^{3}}es

{(incógnita,y,z)|incógnita2+y2z}{\displaystyle \left\{(x,y,z){\Big |}{\sqrt {x^{2}+y^{2}}}\leq z\right\}}.

El conjunto de puntos que satisfacen una restricción de cono de segundo orden es la imagen inversa del cono unitario de segundo orden bajo una transformación afín:

Aiincógnita+bi2doiTincógnita+di[AidoiT]incógnita+[bidi]donortei+1{\displaystyle \lVert A_{i}x+b_{i}\rVert _{2}\leq c_{i}^{T}x+d_{i}\Leftrightarrow {\begin{bmatrix}A_{i}\\c_{i}^{T}\end{bmatrix}}x+{\begin{bmatrix}b_{i}\\d_{i}\end{bmatrix}}\in {\mathcal {C}}_{n_{i}+1}}

y por lo tanto es convexa.

El cono de segundo orden puede incrustarse en el cono de las matrices semidefinidas positivas ya que

||incógnita||t[tIincógnitaincógnitaTt]0,{\displaystyle ||x||\leq t\Leftrightarrow {\begin{bmatrix}tI&x\\x^{T}&t\end{bmatrix}}\succcurlyeq 0,}

Es decir, una restricción de cono de segundo orden es equivalente a una desigualdad matricial lineal . La nomenclatura aquí puede ser confusa; aquíMETRO0{\displaystyle M\succcurlyeq 0}medioMETRO{\displaystyle M}es una matriz semidefinida: es decir

incógnitaTMETROincógnita0 a pesar de incógnitaRnorte{\displaystyle x^{T}Mx\geq 0{\text{ para todo }}x\in \mathbb {R} ^{n}}

lo cual no es una desigualdad lineal en el sentido convencional.

De manera similar, también tenemos,

Aiincógnita+bi2doiTincógnita+di[(doiTincógnita+di)IAiincógnita+bi(Aiincógnita+bi)TdoiTincógnita+di]0{\displaystyle \lVert A_{i}x+b_{i}\rVert _{2}\leq c_{i}^{T}x+d_{i}\Leftrightarrow {\begin{bmatrix}(c_{i}^{T}x+d_{i})I&A_{i}x+b_{i}\\(A_{i}x+b_{i})^{T}&c_{i}^{T}x+d_{i}\end{bmatrix}}\succcurlyeq 0}.

Relación con otros problemas de optimización

Jerarquía de problemas de optimización convexa. (LP: programa lineal, QP: programa cuadrático, SOCP: programa de cono de segundo orden, SDP: programa semidefinido, CP: programa de cono).

CuandoAi=0{\displaystyle A_{i}=0}parai=1,,metro{\displaystyle i=1,\dots ,m}, el SOCP se reduce a un programa lineal . Cuandodoi=0{\displaystyle c_{i}=0}parai=1,,metro{\displaystyle i=1,\dots ,m}El SOCP es equivalente a un programa lineal convexo con restricciones cuadráticas.

Los programas cuadráticos con restricciones cuadráticas convexas también pueden formularse como SOCP reformulando la función objetivo como una restricción. [ 5 ] La programación semidefinida engloba a los SOCP, ya que las restricciones de los SOCP pueden escribirse como desigualdades matriciales lineales (LMI) y pueden reformularse como una instancia de programa semidefinido. [ 5 ] Sin embargo, lo contrario no es válido: existen conos semidefinidos positivos que no admiten ninguna representación de cono de segundo orden. [ 4 ]

Cualquier conjunto semialgebraico convexo cerrado en el plano puede escribirse como una región factible de un SOCP. [ 10 ] Sin embargo, se sabe que existen conjuntos semialgebraicos convexos de dimensión superior que no son representables por SDP; es decir, existen conjuntos semialgebraicos convexos que no pueden escribirse como la región factible de un SDP (ni, con mayor razón , como la región factible de un SOCP). [ 11 ]

Ejemplos

Restricción cuadrática

Consideremos una restricción cuadrática convexa de la forma

incógnitaTAincógnita+bTincógnita+do0.{\displaystyle x^{T}Ax+b^{T}x+c\leq 0.}

Esto es equivalente a la restricción SOCP.

A1/2incógnita+12A1/2b(14bTA1bdo)12{\displaystyle \lVert A^{1/2}x+{\frac {1}{2}}A^{-1/2}b\rVert \leq \left({\frac {1}{4}}b^{T}A^{-1}bc\right)^{\frac {1}{2}}}

Programación lineal estocástica

Consideremos un programa lineal estocástico en forma de desigualdad.

minimizar doTincógnita {\displaystyle \ c^{T}x\ }
sujeto a
PAG(aiTincógnitabi)pag,i=1,,metro{\displaystyle \mathbb {P} (a_{i}^{T}x\leq b_{i})\geq p,\quad i=1,\dots ,m}

donde los parámetrosai {\displaystyle a_{i}\ }son vectores aleatorios gaussianos independientes con mediaa¯i{\displaystyle {\bar {a}}_{i}}y covarianzaΣi {\displaystyle \Sigma _{i}\ }ypag0,5{\displaystyle p\geq 0.5}Este problema puede expresarse como el SOCP.

minimizar doTincógnita {\displaystyle \ c^{T}x\ }
sujeto a
a¯iTincógnita+Φ1(pag)Σi1/2incógnita2bi,i=1,,metro{\displaystyle {\bar {a}}_{i}^{T}x+\Phi ^{-1}(p)\lVert \Sigma _{i}^{1/2}x\rVert _{2}\leq b_{i},\quad i=1,\dots ,m}

dóndeΦ1() {\displaystyle \Phi ^{-1}(\cdot )\ }es la función de distribución acumulativa normal inversa . [ 1 ]

Programación estocástica de conos de segundo orden

Nos referimos a los programas de cono de segundo orden como programas de cono deterministas de segundo orden, ya que los datos que los definen son deterministas. Los programas de cono estocásticos de segundo orden son una clase de problemas de optimización que se definen para manejar la incertidumbre en los datos que definen los programas de cono deterministas de segundo orden. [ 12 ]

Otros ejemplos

Otros ejemplos de modelado están disponibles en el libro de recetas de modelado de MOSEK. [ 13 ]

Solucionadores y lenguajes de programación (scripting)

Véase también

Referencias

  1. 1 2 Boyd, Stephen; Vandenberghe, Lieven (2004). Optimización convexa (PDF) . Cambridge University Press. ISBN 978-0-521-83378-3. Consultado el 15 de julio de 2019 .
  2. Jibrin, Shafiu; Swift, James W. (2024). "Sobre funciones de cono de segundo orden" . Journal of Optimization . 2024 (1) 7090058. doi : 10.1155/2024/7090058 . ISSN 2314-6486 . 
  3. Potra, Lorian A.; Wright, Stephen J. (1 de diciembre de 2000). "Métodos de punto interior". Journal of Computational and Applied Mathematics . 124 ( 1– 2): 281– 302. Bibcode : 2000JCoAM.124..281P . doi : 10.1016/S0377-0427(00)00433-7 .
  4. 1 2 Fawzi, Hamza (2019). "Sobre la representación del cono semidefinido positivo utilizando el cono de segundo orden". Mathematical Programming . 175 ( 1– 2): 109– 118. arXiv : 1610.04901 . doi : 10.1007/s10107-018-1233-0 . ISSN 0025-5610 . S2CID 119324071 .  
  5. 1 2 3 Lobo, Miguel Sousa; Vandenberghe, Lieven; Boyd, Stephen; Lebret, Hervé (1998). "Aplicaciones de la programación de cono de segundo orden" . Álgebra lineal y sus aplicaciones . 284 ( 1–3 ): 193–228 . doi : 10.1016/S0024-3795(98)10032-0 .
  6. "Resolviendo SOCP" (PDF) .
  7. "Optimización de cartera" (PDF) .
  8. Li, Haksun (16 de enero de 2022). Métodos numéricos con Java: para ciencia de datos, análisis e ingeniería . APress. pp. Capítulo 10. ISBN  978-1-4842-6796-7.
  9. Chen, Ying; Bolzern, Paolo; Colaneri, Patrizio (16 de mayo de 2019). "Estabilidad, rendimiento ℒ1 y diseño de retroalimentación de estado para sistemas lineales en conos de helado" . International Journal of Control . 94 (3): 784. doi : 10.1080/00207179.2019.1616825 . hdl : 11311/1090203 . Recuperado el 19 de febrero de 2026 .
  10. Scheiderer, Claus (2020-04-08). "Representación de cono de segundo orden para subconjuntos convexos del plano". arXiv : 2004.04196 [ math.OC ].
  11. Scheiderer, Claus (2018). "Sombras espectrales" . SIAM Journal on Applied Algebra and Geometry . 2 (1): 26– 44. doi : 10.1137/17M1118981 . ISSN 2470-6566 . 
  12. Alzalg, Baha M. (2012-10-01). "Programación estocástica de cono de segundo orden: modelos de aplicaciones" . Modelado matemático aplicado . 36 (10): 5122– 5134. doi : 10.1016/j.apm.2011.12.053 . ISSN 0307-904X . 
  13. "Manual de modelado MOSEK - Optimización cuadrática cónica" .
  14. "Solucionador de programación de conos de segundo orden - MATLAB coneprog" . MathWorks . 1 de marzo de 2021. Consultado el 15 de julio de 2021 .
  15. "Algoritmo de programación de cono de segundo orden - MATLAB y Simulink" . MathWorks . 1 de marzo de 2021. Consultado el 15 de julio de 2021 .
  16. "Manual de modelado MOSEK: los conos de potencia" .