Repaso para la I2: conjuntos y funciones convexas, criterio de la matriz Hessiana, relevancia de la convexidad (mínimo local = global) y existencia de solución óptima (Bolzano-Weierstrass). Con los ejercicios resueltos de la Ayudantía 5.
Un conjunto C es convexo si, al tomar dos puntos cualesquiera de C, todo el segmento que los une queda dentro de C. No tiene "hoyos" ni "entrantes".
La expresión λx1+(1−λ)x2 se llama combinación convexa de x1 y x2.
Son convexos: los semiespacios {x:aTx≤b}, las bolas {x:∥x∥≤r}, los poliedros (intersección de semiespacios) y todo Rn.
Propiedad muy útil: la intersección de conjuntos convexos es convexa. Por eso un poliedro, al ser intersección de semiespacios (convexos), es convexo.
Demostrar que la bola B={x∈Rn:∥x∥≤r} es convexa. Tomamos x1,x2∈B (es decir ∥x1∥≤r, ∥x2∥≤r) y λ∈[0,1], y acotamos la norma de la combinación convexa usando desigualdad triangular y homogeneidad:
Como la combinación convexa también cumple ∥x∥≤r, pertenece a B. Por lo tanto B es convexo. ■
Una función f es convexa si, para todo par de puntos, la cuerda que une (x,f(x)) con (y,f(y)) queda por encima del gráfico de la función:
Si la desigualdad es estricta (<) para x=y y λ∈(0,1), se dice estrictamente convexa.
Demostrar convexidad por la desigualdad de Jensen suele ser engorroso. Para funciones diferenciables hay criterios más rápidos.
Recetario para demostrar convexidad. Para un conjunto: usa la definición (segmento dentro) o que es intersección de convexos. Para una función: primero intenta los atajos (suma/máximo/composición con afín); si no, calcula la Hessiana y verifica d⊤∇2fd≥0.
Si f:Rm→R es convexa y g:Rn→Rm es lineal, entonces h(x)=f(g(x)) es convexa. Como g es lineal, g(λx+(1−λ)y)=λg(x)+(1−λ)g(y); aplicando luego la convexidad de f:
Un problema convexo consiste en minimizar una función convexa sobre un dominio convexo (o maximizar una cóncava). Su importancia está en un teorema potentísimo:
Esto significa que no hay que preocuparse por quedar "atrapado" en un óptimo local: cualquier mínimo local que encuentres ya es la mejor solución posible.
Clave para la I2: como el dominio de un LP (un poliedro) es convexo y el objetivo lineal es convexo, el óptimo que encuentra Simplex es global. La convexidad es lo que justifica que el método funcione.
Si además la función objetivo es estrictamente convexa sobre el dominio convexo, entonces la solución óptima (si existe) es única.
Antes de resolver, conviene saber si el problema tiene óptimo. Hay dos criterios principales.
Si se cumplen estas cuatro condiciones, el problema admite al menos una solución óptima:
Un LP tiene solución si el dominio es no vacío y la función objetivo está acotada en la dirección de optimización. Una variante es la coercitividad: en un problema de maximización, si una variable sin cota superior aparece con coeficiente negativo (penalización), entonces f→−∞ cuando esa variable crece, lo que evita el no acotamiento.
Para min 7x1+4x2−3x3 con x1+3x2+9x3≤12, 2x2−5x3≤5, 3x1−x2≤8, x1,x2≥0, x3∈R: la objetivo es lineal (continua); el dominio es cerrado; es no vacío porque (0,0,0) es factible. Para el acotamiento, de las restricciones se acota x3:
Con x3∈[−1,34] acotado, las restricciones (1) y (3) acotan a x1,x2 (que ya tienen cota inferior 0). El dominio es compacto ⇒ por Bolzano-Weierstrass, existe solución óptima. ■
Sea f(x,y)=ln(ex+ey). Su gradiente y su Hessiana son:
Para verificar que es SDP, tomamos d=(d1,d2) y calculamos la forma cuadrática:
El factor escalar es positivo y (d1−d2)2≥0, así que la Hessiana es SDP en todo R2 y f es convexa. ■