Ejercicios Resueltos Optimización v3
Ejercicios Resueltos Optimización v3
Ejercicios Resueltos Optimización v3
Este documento es material de estudio del curso ICS1113 y fue financiado por fondos de la Escuela de Ingenierı́a UC.
1.2 Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Capı́tulo 2: Convexidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1 Definiciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2 Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.2 Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
4.1 Definiciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
4.2 Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
I
CAPÍTULO 1
EXISTENCIA DE SOLUCIÓN ÓPTIMA
3. D es no vacı́o,
2. D es cerrado,
3. D es no vacı́o,
Definición 1.1. Una función f (x) es coercitiva sobre un dominio D ⊂ Rn si lı́m||x||→∞ f (x) → ∞.
1. D es no vacı́o
1
1.2 Ejercicios
(P ) mı́n f (x1 , x2 )
s.a. (x1 , x2 ) ∈ D.
Indique para cuáles de los siguientes casos se puede garantizar existencia de solución óptima.
(a) f (x1 , x2 ) = 4x1 + 5x2 ; D = {(x1 , x2 ) ∈ R2 : 2x1 + 5x2 ≤ 20, x1 + 3x2 ≥ 3, 0 ≤ x1 ≤ 10}
(c) f (x1 , x2 , x3 ) = 3x1 + 2x2 − 3x3 ; D = {2x1 + x2 + 2x3 ≤ 5, x1 − x3 ≤ 10, 3x2 − 2x3 ≥ 0, x1 , x3 ≥ 0}
(d) f (x1 , x2 , x3 , x4 ) = 7x1 − x2 − 6x3 + 2x4 ; D = {3x1 + 2x2 − 6x3 ≤ 4, 5x1 − 2x2 ≤ 7, x2 − 3x4 ≥ 10,
x1 , x4 ≥ 0}
¿Admite solución óptima? Demuestre.
Solución
De R1:
2x1 + 5x2 ≤ 20
⇒ 5x2 ≤ 20 − x1
⇒ 5x2 ≤ 20
⇒ x2 ≤ 4
2
De R2:
x1 + 3x2 ≥ 3
⇒ 3x2 ≥ 3 − x1
⇒ 3x2 ≥ −7
7
⇒ x2 ≥ −
3
Finalmente, como se cumplen las 4 condiciones del teorema de BW, Pa admite solución óptima.
Función objetivo continua: f (x1 , x2 , x3 ) = 3x1 + 2x2 − 3x3 es una función lineal ✓.
Dominio cerrado: todas las restricciones del dominio son de desigualdad irrestricta ✓.
Dominio no vacı́o: (0, 0, 0) es factible ✓.
Dominio acotado: debemos acotar las variables por arriba y por abajo.
Por NV:
x1 ≥ 0 z≥0
Por R3:
3x2 − 2x3 ≥ 0
3
2x3 toma su mı́nimo valor cuando z = 0:
⇒ x2 ≥ 0
Por R1:
2x1 + x2 + 2x3 ≤ 5
Esto nos permite acotar superiormente las 3 variables. Primero, acotamos x1 :
2x1 ≤ 5 − x2 − 2x3
5 − x2 − 2x3 toma su máximo valor cuando x2 = x3 = 0:
⇒ x1 ≤ 2,5
Análogamente, acotamos superiormente x2 y x3 :
x2 ≤ 5 − 2x1 − 2x3 2x3 ≤ 5 − 2x1 − x2
⇒ x2 ≤ 5 ⇒ x3 ≤ 2,5
Problema lineal: la función objetivo es lineal y las restricciones son todas lineales ⇒ (P ) es lineal ✓.
Dominio cerrado: todas las restricciones del dominio son de desigualdad irrestricta (se incluyen los
bordes del dominio) ⇒ D es cerrado ✓.
Dominio no vacı́o: se puede encontrar un punto perteneciente al dominio (que cumpla todas las restric-
ciones), por ejemplo, (x1 , x2 , x3 , x4 ) = (0, 10, 3, 0) ∈ D ⇒ D ̸= ∅ ✓.
Función objetivo acotada superiormente: al ser un problema de maximización, interesa acotar la f.o.
superiormente, es decir, debemos encontrar un k ∈ R : f (x) ≤ k ∀x ∈ D. Para esto, debemos
formar la función objetivo usando las restricciones. Reescribimos las restricciones para que todas
las desigualdades sean “≤” y las sumamos:
3x1 + 2x2 − 6x3 ≤ 4 (R1)
5x1 − 2x2 ≤ 7 (R2)
− x2 + 3x4 ≤ −10 (R3)
− x1 ≤ 0 (NV)
− x4 ≤ 0 (NV)
Finalmente, como se cumplen todas las condiciones, queda demostrado por teorema de existencia de
solución para problemas lineales que (P ) admite solución óptima.
4
Problema 1.2. I1 2021-1
(a) En un problema de optimización con un dominio abierto y acotado: Explique por qué podrı́a no existir solu-
ción. De un ejemplo.
(c) Sea (P ) mı́n{f (x) : x ∈ D}, con D ⊂ Rn convexo, cerrado y no vacı́o, y f : D → R continua y
estrictamente convexa. ¿Es verdad que (P ) tiene solución óptima? Demuestre o refute con un contraejemplo.
Solución
(a) Podrı́an existir sucesiones de soluciones factibles cuyos valores objetivos convergen al valor óptimo del
problema, pero no lo alcanzan nunca y, por ende, no existe solución. Por ejemplo, con mı́n{x : x ∈
(0, 1]}, cuyo valor óptimo es 0, pero no es alcanzado por ningún x factible.
(b) Queremos demostrar que f alcanza su mı́nimo sobre los reales, es decir, que mı́nx∈Rn f (x). Se definen
los problemas:
Por teorema de BW, sabemos que (Q) admite solución óptima: f continua, D = L compacto (cerrado y
acotado) y no vacı́o.
Además, sabemos que (Q) abarca un menor espacio de soluciones que (P ): L ⊆ Rn , pero como L solo
toma las soluciones x tales que f (x) ≤ α, entonces para todo x ∈ Rn \ L se cumple que f (x) > α. Es
decir, las soluciones en Rn que no se encuentran en L son peores o iguales en objetivo que las de (Q). Por
lo tanto, el problema (P ) admite solución óptima (f alcanza su mı́nimo sobre Rn ) y se cumple que:
v(P ) = v(Q)
(P ) mı́n ex
s.a : x ≤ 0
En este ejemplo, D y f cumplen las condiciones descritas, pero como f (x) decrece a medida que x
disminuye y D es no acotado, entonces (P ) no tiene solución óptima.
P ) max f (x)
s.a gj (x) ≤ bj ∀j = 1, ..., n
x ∈ Rn
5
Con f (x) una función continua. Usted sabe que el siguiente problema (Pk ) admite solución óptima con cualquier
k = 1, ..., n y cualquier valor escalar de c:
Pk ) max c · xk
s.a gj (x) ≤ bj ∀j = 1, ..., n
x ∈ Rn
Solución
Debemos verificar si P ) cumple las condiciones del teorema de Bolzano-Weierstrass:
Para las siguientes condiciones, P tiene el mismo dominio que los problemas Pk : D = {x ∈ Rn : gj (x) ≥
bj ; ∀j = 1, ..., n}.
Dominio cerrado: todas las restricciones del dominio son de desigualdad irrestricta (se incluyen los bordes del
dominio) ⇒ D es cerrado ✓.
Dominio no vacı́o: los problemas Pk admiten solución óptima, es decir, existe al menos una solución pertene-
ciente al dominio D ✓.
Dominio acotado: podemos considerar los siguientes problemas que son casos particulares de los problemas
Pk :
Pk1 ) máx −1 · xk ∼ Pk1 ) mı́n xk
s.a x∈D s.a x ∈ D
Pk2 ) máx 1 · xk
s.a x∈D
Como ambos problemas admiten solución, entonces cada una de las variables está acotada superior e
inferiormente para dicho dominio. Es decir, D está acotado inferiormente, ya que Pk encuentra solución
en el sentido de minimización, y también está acotado superiormente, ya que Pk encuentra solución en el
sentido de maximización. Por lo tanto, D es acotado ✓.
Finalmente, como se cumplen las 4 condiciones del teorema de BW, el problema P ) admite solución óptima.
6
(b) Considere el problema
min ax21 + 2bx1 x2 + cx22
s.a: x ∈ D ⊂ R2
donde a, b y c son escalares reales no nulos. En cada uno de los siguientes casos, indique si el problema tiene
solución óptima o no:
i) D = {(x1 , x2 ) ∈ R2 : −1 ≤ x1 ≤ 1, −1 ≤ x2 ≤ 1}
ii) D = R2 y se tiene que ax21 + 2bx1 x2 + cx22 ≥ 0 para todo x = (x1 , x2 ) ∈ R2
iii) D = R2 y existe y = (y1 , y2 ) ∈ R2 tal que ay12 + 2by1 y2 + cy22 < 0.
Solución
(a) Notemos que el dominio es el conjunto Rn y este no es acotado, la función objetivo es de segundo grado,
por lo que el teorema de existencia de solución para problemas lineales no aplica. Revisemos entonces las
condiciones del teorema de función coercitiva:
Función objetivo continua: es la suma de funciones exponenciales, por lo que no se indefine en ningún
punto ⇒ la f.o. es continua ✓.
Dominio cerrado: el dominio son los reales y este es un conjunto cerrado (se incluyen los bordes del
dominio) ✓.
Dominio no vacı́o: todo punto perteneciente a Rn es parte del dominio, por ejemplo el origen ⇒ D ̸= ∅
✓.
Función objetivo coercitiva sobre el dominio: debemos demostrar que:
7
Función objetivo continua: es un polinomio, por lo que no se indefine en ningún punto del domi-
nio ⇒ es continua ✓.
Dominio cerrado: el dominio son los reales y este es un conjunto cerrado (se incluyen los bordes
del dominio) ✓.
Dominio no vacı́o: todo punto perteneciente a R2 es parte del dominio, por ejemplo el origen ⇒
D ̸= ∅ ✓.
Función objetivo acotada inferiormente: por enunciado, la función objetivo está acotada inferior-
mente por 0: ax21 + 2bx1 x2 + cx22 ≥ 0 ✓
Finalmente, como se cumplen todas las condiciones, queda demostrado por teorema de acotamiento
inferior que (P ) admite solución óptima.
iii) El problema es de minimización, por lo que busca una solución (x1 , x2 ) que entregue el menor valor
posible para f (x1 , x2 ) = ax2 + 2bx1 x2 + 2x22 . Sin embargo, el enunciado indica que existe una
solución y = (y1 , y2 ) tal que f (y) < 0. Luego, si ponderamos esta solución y por un escalar distinto
de cero λ ∈ R \ {0}, nos queda:
Es decir, las soluciones de la forma λ · y son factibles y su valor objetivo tiende a −∞ si |λ| tiende
a infinito. Luego el problema es no acotado y no tiene solución.
8
CAPÍTULO 2
CONVEXIDAD
2.1 Definiciones
Nota:
D es convexo.
2. Si además (P ) admite solución óptima y f es estrictamente convexa, entonces la solución óptima es única.
9
Definición 2.7. Gradiente de una función
∂f (x)
∂f∂x(x)
1
∂x
∇f (x) = . 2
..
∂f (x)
∂xn
Definición 2.8. Matriz Hessiana de una función
∂ 2 f (x) ∂ 2 f (x) ∂ 2 f (x)
2 ∂x1 ∂x2 ··· ∂x1 ∂xn
∂ 2∂x 1
f (x) ∂ 2 f (x) ∂ 2 f (x)
∂x22
···
H(f (x)) = ∇2 f (x) =
∂x2 ∂x1 ∂x2 ∂xn
.. .. ... ..
. . .
∂ 2 f (x) ∂ 2 f (x) ∂ 2 f (x)
∂xn ∂x1 ∂xn ∂x2 ··· ∂x2n
Definición 2.9. Matriz (semi) definida positiva
La matriz A simétrica n × n es . . .
Dada la matriz:
a11 a12 ··· a1n
a21 a22 ··· a2n
A= .
.. .. ..
.. . . .
am1 am2 · · · amn
Los menores principales superiores son:
∆1 = |a11 |
a11 a12
∆2 =
a21 a22
∆n = |A|
10
A es definida positiva λi >0 ∀i
A es definida negativa λi <0 ∀i
A es semidefinida positiva λi ≥0 ∀i
A es semidefinida negativa λi ≤0 ∀i
2.2 Ejercicios
(a) Para cada una de las siguientes afirmaciones, determine si es verdadera o no. En caso de ser verdadera, justi-
fique por qué lo es. En caso de no serlo, exhiba un contraejemplo.
(b) [I1 2020-2] Sea (P ) mı́n{f (x) : x ∈ D} un problema de optimización convexo. Demuestre que si x1 , x2 ∈ D
son soluciones óptimas para (P ), entonces toda combinación convexa de ellas también es solución óptima para
(P ).
(c) [I1 2018-2] Sea f (x) = máx{f1 (x), ..., fk (x)}, donde f1 (x), ..., fk (x) : Rn → R son k funciones convexas
en un conjunto convexo S. Demuestre formalmente que f es una función convexa en S.
(d) [Examen 2020-2] Sean f : Rm → R una función convexa y g : Rn → Rm una transformación lineal.
Demuestre que h : Rn → R dada por h(x) = f (g(x)) es convexa.
(e) [I2 2019-1] Argumente matemáticamente que para un problema convexo de minimización todo óptimo local
es global.
Solución
(a) (i) Falso. Sean f (x) = x y g(x) = −x funciones convexas. Sin embargo, su multiplicación no es
convexa: h(x) = f (x) · g(x) = −x2 .
(ii) Falso. Sean f (x) = 0 y g(x) = x2 funciones convexas. Sin embargo, su diferencia no es convexa:
h(x) = f (x) − g(x) = −x2 .
(iii) Verdadero, pues la función es convexa. Es “un” minimizador y no “el único” minimizador porque la
función es convexa, no estrictamente convexa.
(b) Sea v ∗ = f (x1 ) = f (x2 ) el valor óptimo de (P ). Se tiene que v ∗ ≤ f (x0 ) para todo x0 ∈ D. Dado
λ ∈ [0, 1], se define z como una combinación convexa de las soluciones óptimas x1 y x2 :
z = λx1 + (1 − λ)x2
11
Por la convexidad del dominio, sabemos que la combinación convexa de dos puntos en D también perte-
nece a D. Por lo tanto, z ∈ D, lo que significa que z es solución factible para (P ) y que v ∗ ≤ f (z).
Además, por la convexidad de f tenemos:
f (z) = f (λx1 + (1 − λ)x2 ) ≤ λf (x1 ) + (1 − λ)f (x2 ) = λv ∗ + (1 − λ)v ∗ = v ∗
Como f (z) ≤ v ∗ y f (z) ≥ v ∗ , concluimos entonces que f (z) = v ∗ y por ende, z es solución óptima de
(P ).
(c) Sean x1 , x2 ∈ S y λ ∈ [0, 1]. Para simplificar el desarrollo, se define z = λx1 + (1 − λ)x2 .
Por convexidad de fi , con i ∈ {1, ..., k} un ı́ndice cualquiera, se cumple que:
fi (z) ≤ λfi (x1 ) + (1 − λ)fi (x2 )
Luego, como f (x) = máx{f1 (x), ..., fk (x)}, se tiene que fi (x) ≤ f (x), ∀x ∈ S, ∀i ∈ {1, ..., k}.
Entonces:
fi (z) ≤ λfi (x1 ) + (1 − λ)fi (x2 ) ≤ λf (x1 ) + (1 − λ)f (x2 ) ∀i ∈ {1, ..., k}
Como toda función fi cumple con la condición:
f1 (z) ≤ λf (x1 ) + (1 − λ)f (x2 )
f2 (z) ≤ λf (x1 ) + (1 − λ)f (x2 )
..
.
fk (z) ≤ λf (x1 ) + (1 − λ)f (x2 )
Finalmente, se tiene:
f (λx1 + (1 − λ)x2 ) ≤ λf (x1 ) + (1 − λ)f (x2 )
Lo que demuestra que f = máx{f1 , ..., fk } es convexa.
(d) Sean x, y ∈ Rn y λ ∈ [0, 1]. Como g es lineal, entonces es cóncava y convexa al mismo tiempo, es decir:
g(λx + (1 − λ)y) = λg(x) + (1 − λ)g(y). Por lo tanto:
h(λx + (1 − λ)y) = f (g(λx + (1 − λ)y)) = f (λ · g(x) + (1 − λ) · g(y))
12
Problema 2.2. Convexidad de una función
Considere la siguiente función f (x1 , x2 ) = ax21 + 2x22 − bx1 x2 + 4 ¿Qué condiciones deben cumplir los parámetros
a y b para que la función sea estrictamente convexa?
Solución
Para que la función sea estrictamente convexa el Hessiano debe ser definido positivo. Tenemos, entonces que:
2ax1 − bx2
∇f (x1 , x2 ) =
4x2 − bx1
2 2a −b
Hf (x1 , x2 ) = ∇ (f (x1 , x2 )) =
−b 4
Para que el Hessiano sea definido positivo basta que los dos determinantes principales sean > 0. Esto lleva a las
condiciones:
2a > 0
8a − b2 > 0
luego, las condiciones que deben cumplir los parámetros para que la función sea convexa es que a > 0 y b2 < 8a
n
αi e−xi es convexa para cualquier valor α ∈ Rn .
P
(a) [Examen 2022-1] La función g(x) =
i=1
x21
(b) [Examen 2017-2] La función f (x1 , x2 ) = es convexa para el dominio definido por x2 > 0.
x2
Solución
(a) Falso. Si αi < 0 para cierto i ∈ {1, ..., n}, entonces g no es necesariamente convexa. Por ejemplo si
n = 2, α1 = 1 y α2 = −2 entonces g(x) = e−x1 − 2e−x2 . La cual no es convexa.
13
Revisemos si la matriz es (semi) definida positiva:
2 −2x1
x2
x22
d1
T
d Hd = d1 d2 · ·
2x21 d2
−2x1
x22 x32
−x1
1
x2
2 d1
= · d1 d2 · ·
x2 −x1 x21 d2
x2 x22
−d1 · x1 d2 · x21
2 d2 · x1 d
= · d1 − + 2
· 1
x2 x2 x2 x2 d2
d1 · d2 · x1 d22 · x21
2 2
= · d1 − 2 · +
x2 x2 x22
d2 x1 2
2
= · d1 − ≥0 ∀x2 > 0, ∀d ∈ R2
x2 x2
(a) Defina la siguiente función unidimensional f (x) = − ln(100 − |t|). Muestre que f (t) es convexa para todo
t ∈ (−100, 100).
3
3
P
(b) Sea D = x ∈ R : xi = 1 ; 0 ≤ xi ≤ 1, i = 1, 2, 3 . Usando el resultado de la pregunta anterior,
i=1
muestre que la función f cumple:
f (c1 · x1 + c2 · x2 + c3 · x3 ) ≤ x1 · f (c1 ) + x2 · f (c2 ) + x3 · f (c3 )
para todo x ∈ D. Suponga que c es tal que f (·) estará siempre bien definida, es decir, |c1 | + |c2 | + |c3 | < 100.
Solución
(a) Notemos que |t| es una función convexa para todo t ∈ (−100, 100), por lo que 100 − |t| es cóncava. Sa-
bemos que ln(·) es una función creciente, por lo tanto mantiene la propiedad de concavidad o convexidad
de la función interior, por lo que ln(100 − |t|) es cóncava y − ln(100 − |t|) es convexa.
(b) Sabemos que f (x) es convexa y que x1 , x2 , x3 ∈ [0, 1]. Entonces, dada la convexidad de f , podemos
14
separar el lado izquierdo de la desigualdad usando x1 como parámetro:
(1 − x1 )
f (c1 · x1 + c2 · x2 + c3 · x3 ) = f c1 · x1 + (c2 · x2 + c3 · x3 ) ·
(1 − x1 )
c2 · x2 c3 · x3
= f c1 · x1 + + · (1 − x1 )
1 − x1 1 − x1
x2 x3
≤ x1 · f (c1 ) + (1 − x1 ) · f c2 · + c3 ·
1 − x1 1 − x1
x2 x3
Nos queda la desigualdad: f (c1 ·x1 +c2 ·x2 +c3 ·x3 ) ≤ x1 ·f (c1 )+(1−x1 )·f c2 · + c3 ·
1 − x1 1 − x1
x2 x3
Como x1 + x2 + x3 = 1, entonces + = 1, lo que nos permite aplicar nuevamente las
1 − x1 1 − x1
propiedades de convexidad sobre f y ordenar los términos:
x2 x3
f (c1 · x1 + c2 · x2 + c3 · x3 ) ≤ x1 · f (c1 ) + (1 − x1 ) · f c2 · + c3 ·
1 − x1 1 − x1
x2 x3
≤ x1 · f (c1 ) + (1 − x1 ) · · f (c2 ) + · f (c3 )
1 − x1 1 − x1
= x1 · f (c1 ) + x2 · f (c2 ) + x3 · f (c3 )
√ √
a) ¿Es {(x1 , x2 ) ∈ R2+ : x1 + x2 ≤ 1} un conjunto convexo? Justifique.
b) Demuestre que la función f (x) : R → R definida por f (x) = e−|x| no es convexa y tampoco es cóncava.
c) Sean g : Rn → R una función convexa y f (x) = eg(x) . Demuestre utilizando la definición de convexidad que
f (x) es convexa. Hint: ex es convexa.
Solución
a) Para que el conjunto sea convexo debe cumplirse que la combinación convexa de cualquier par de puntos
9 9
del conjunto también debe pertenecer al conjunto. Tomando: y1 = ( 10 , 0), y2 = (0, 10 ) y y3 = 12 y1 + 12 y2 .
√ √
Como y1 y y2 cumplen con x1 + x2 ≤ 1 pero y3 no, el conjunto no es convexo.
15
f (λx1 + (1 − λ)x2 ) = f (−1) = 0,3679
λf (x1 ) + (1 − λ)f (x2 ) = 0,5 · 0,1353 + 0,5 · 1 = 0,56765
f (λx1 + (1 − λ)x3 ) = f (0) = 1
λf (x1 ) + (1 − λ)f (x3 ) = 0,5 · 0,1353 + 0,5 · 0,1353 = 0,1353
Como f (λx1 + (1 − λ)x2 ) < λf (x1 ) + (1 − λ)f (x2 ) ←→ 0, 3679 < 0, 56765, la función no es cóncava.
Como f (λx1 + (1 − λ)x3 ) > λf (x1 ) + (1 − λ)f (x3 ) ←→ 1 > 0, 1353, la función no es convexa.
Es decir:
eg(λx1 +(1−λ)x2 ) ≤ λeg(x1 ) + (1 − λ)eg(x2 )
Solución
16
Reemplazamos la función y desarrollamos:
1 1 1
≤ λ· + (1 − λ) ·
λx + (1 − λ)y x y
1 λy + (1 − λ)x
⇔ ≤
λx + (1 − λ)y xy
(b)
2 2 1 2
C = {(x1 , x2 ) ∈ R+ : x1 · x2 ≥ 1} = {(x1 , x2 ) ∈ R+ : x2 ≥ } = {(x1 , x2 ) ∈ R+ : x2 ≥ f (x1 )}
x1
Es decir, C es el epigrafo de f . Como f es convexa, concluimos que C es convexo.
(a) f1 (x)
Pn
(b) i=1 fi (x)
Pn/2
(c) (Asuma que n es par) i=1 f2i
Solución
(a) La función fi (x) evaluada en i = 1 entrega como resultado f1 (x) = x − a1 , la cual es una función lineal
y, por lo tanto es convexa
17
(b) No es convexa, si evaluamos para n = 3 obtenemos
3
X
fi (x) = (x − a1 ) + (x − a2 )2 + (x − a3 )3
i=1
f1 (x) · g(x) = (x − a1 ) · x2 = x3 − a1 x
* Se pueden redefinir las funciones impares de esta forma porque jamás serán máximo cuando su valor sea
menor a 0 (las funciones pares siempre serán mayores o iguales a 0). Ahora cada función fi es convexa.
Como fue explicado en el ejercicio 2.1 (c), el máximo de funciones convexas es una función convexa, por
lo cual queda demostrado lo pedido.
18
CAPÍTULO 3
MODELOS EQUIVALENTES
Se dice que (P1 ) y (P2 ) son modelos equivalentes si ambos tienen el mismo conjunto de soluciones óptimas o
si existe una transformación simple que permite construir el conjunto de soluciones óptimas de uno a partir de
las soluciones óptimas del otro. En palabras simples, dos modelos son equivalentes si, y solo si, al resolver uno,
resolvemos el otro.
Notación: (P1 ) ∼ (P2 ) equivale a señalar “(P1 ) y (P2 ) son equivalentes”.
A continuación, se presentan algunas equivalencias útiles.
Si el problema es de maximización:
n
P n
P
(P1 ) mı́n fi (x) (P2 ) mı́n ui
x∈D i=1 ∼ x∈D, u∈Rn i=1
s.a: fi (x) ≤ ui ∀i ∈ {1, ..., n}
Equivalencia 3.1.4. Equivalencia 4: composición de funciones
Sea g una función estrictamente creciente sobre la imagen de f .
(P1 ) mı́n f (x) ∼ (P2 ) mı́n g f (x)
x∈D x∈D
19
Equivalencia 3.1.5. Equivalencia 5: funciones recı́procas
Si la función objetivo es positiva para todo x en el dominio, f (x) > 0 ∀x ∈ D, entonces:
1
(P1 ) mı́n f (x) ∼ (P2 ) máx
x∈D x∈D f (x)
Sin embargo, si |f (x)| ≥ b, las restricciones resultantes son disjuntas y no son poliedralmente representables.
20
3.2 Ejercicios
Solución
La afirmación es verdadera. Como −g es estrictamente creciente, podemos usar la equivalencia de composición
de funciones para obtener (P0 ) equivalente a (P ). Luego, ocupamos la equivalencia mı́n / máx y obtenemos
(P ′ ), equivalente a (P0 ) y, por ende, a (P ):
Solución
Usando la equivalencia mı́n / máx, tenemos que (P1 ) es equivalente a máx{−f (x) : x ∈ D}.
Luego, como −f (x) > 0 para todo x∈ D, podemos ocupar la equivalencia de funciones recı́procas, por lo que
1
este problema es equivalente a: mı́n : x∈D .
−f (x)
1
Como > 0 para todo x ∈ D y la función cuadrática es creciente sobre R+ , podemos usar la equivalencia
−f (x) ( )
2
1
de composición de funciones y obtenemos: mı́n : x∈D
−f (x)
Por ende, la afirmación es verdadera.
21
Solución
Primero, juntamos la funciones exponenciales de la función objetivo. Podemos ver que, si descomponemos la
pitatoria, los exponentes se suman:
n
− máx {3xj −xj−1 } Y
f (x) = e j∈{2,...,n}
· e3xi +4
i=1
− máx {3xj −xj−1 }
=e j∈{2,...,n}
· e3x1 +4 · e3x2 +4 · · · e3xn +4
n
P
− máx {3xj −xj−1 } (3xi +4)
=e j∈{2,...,n}
· ei=1
n
P
− máx {3xj −xj−1 }+ 3xi +4n
j∈{2,...,n}
=e i=1
Como la función exponencial es estrictamente creciente, podemos usar la equivalencia de composición de fun-
ciones y obtenemos el problema (P ′ ), esto por que la función de logaritmo natural es una función creciente sobre
la imagen de e:
n
(P ) ∼ (P ′ ) máx −
P
máx {3xj − xj−1 } + 3xi + 4n
x j∈{2,...,n} i=1
Omitiendo la constante 4n y multiplicando la f.o. por −1 para cambiar el problema a uno de minimización,
obtenemos (P ′′ ):
n
(P ) ∼ (P ′′ ) mı́n
P
máx {3xj − xj−1 } − 3xi
x j∈{2,...,n} i=1
Luego, linearizamos la función objetivo usando la equivalencia mı́nimo de máximos y linearizamos la función
22
“mı́n” de la primera restricción. Obtenemos (P ′′′ ):
n
(P ) ∼ (P ′′′ ) mı́n u −
P
3xi
x,u i=1
23
donde c es tal que la función objetivo está siempre bien definida, es decir: |c1 | + |c2 | + |c3 | < 100. Muestre que el
problema (P ) es equivalente al siguiente problema (Q) de optimización lineal:
(Q) máx w
s.a: w + c1 · x1 + c2 · x2 + c3 · x3 ≤ 100
w − c1 · x1 + c2 · x2 + c3 · x3 ≤ 100
x1 + x2 + x3 = 1
0 ≤ xi ≤ 1 i = 1, 2, 3
Solución
Sabemos que la función exponencial es estrictamente creciente, por lo que podemos ocupar la equivalencia de
composición de funciones para eliminar el logaritmo de la siguiente forma, obteniendo (P ′ ):
(P ) ∼ mı́n e− ln(100−|c1 ·x1 +c2 ·x2 +c3 ·x3 |) ∼ (P(1) ) mı́n (100 − |c1 · x1 + c2 · x2 + c3 · x3 |)−1
s.a: x1 + x2 + x3 = 1 s.a: x1 + x2 + x3 = 1
0 ≤ xi ≤ 1 i = 1, 2, 3 0 ≤ xi ≤ 1 i = 1, 2, 3
Si nos fijamos en (Q), su función objetivo es una única variable, por lo que ocupamos la equivalencia linealiza-
ción de objetivo y obtenemos (P(2) ):
(P ) ∼ (P(2) ) mı́n u
s.a: u ≥ (100 − |c1 · x1 + c2 · x2 + c3 · x3 |)−1
x1 + x2 + x3 = 1
0 ≤ xi ≤ 1 i = 1, 2, 3
u∈R
Como sabemos que |c1 |+|c2 |+|c3 | < 100, entonces u ≥ 100−1 > 0, por lo que podemos invertir la desigualdad:
Reordenamos los términos para que sea más simple linearizar el valor absoluto:
c1 · x1 + c2 · x2 + c3 · x3 ≤ 100 − u−1
|c1 · x1 + c2 · x2 + c3 · x3 | ≤ 100 − u−1 ∼
−(c1 · x1 + c2 · x2 + c3 · x3 ) ≤ 100 − u−1
−1
u + c1 · x1 + c2 · x2 + c3 · x3 ≤ 100
⇒
u−1 − c1 · x1 − c2 · x2 − c3 · x3 ≤ 100
(P ) ∼ (P(3) ) mı́n u
s.a: u−1 + c1 · x1 + c2 · x2 + c3 · x3 ≤ 100
u−1 − c1 · x1 − c2 · x2 − c3 · x3 ≤ 100
x1 + x2 + x3 = 1
0 ≤ xi ≤ 1 i = 1, 2, 3
u ≥ 100−1
24
Luego, hacemos un cambio de variable: u−1 = w. Como u = w−1 ≥ 100−1 > 0, entonces 0 < w ≤ 100.
Reemplazamos y obtenemos (P(4) ):
Usamos la equivalencia de funciones recı́procas para pasar de mı́n a máx e invertir la función objetivo. Obtene-
mos (P(5) ):
(P ) ∼ (P(5) ) máx w
s.a: w + c1 · x1 + c2 · x2 + c3 · x3 ≤ 100
w − c1 · x1 − c2 · x2 − c3 · x3 ≤ 100
x1 + x2 + x3 = 1
0 ≤ xi ≤ 1 i = 1, 2, 3
0 < w ≤ 100
Por último, dado que w = u−1 ≤ 100 − |c1 · x1 + c2 · x2 + c3 · x3 | ≤ 100, y esto está descrito por las
dos primeras restricciones, provenientes de la descomposición del valor absoluto, se puede prescindir de la
restricción w ≤ 100. Además, como el objetivo de este problema es maximizar w, esta siempre va a buscar ser
lo más grande poible, por lo que se puede prescindir de la restricción w ≥ 0. Obtenemos (P(6) ):
(P ) ∼ (P(6) ) máx w
s.a: w + c1 · x1 + c2 · x2 + c3 · x3 ≤ 100
w − c1 · x1 − c2 · x2 − c3 · x3 ≤ 100
x1 + x2 + x3 = 1
0 ≤ xi ≤ 1 i = 1, 2, 3
25
Solución
Empecemos con la primera restricción, donde podemos sumar la variable de holgura h1 ≥ 0, y para la segunda
restricción, podemos restar la variable de exceso e1 ≥ 0. El problema queda de la siguiente forma:
mı́n z = −x1 + x2 + x3 + 5
s.a 2x1 + 3x2 + x3 + h1 = 5
− x1 + x2 − e 2 = 3
x1 ≤ 0
x2 ≥ 1
h1 , e1 ≥ 0
En cuanto a la naturaleza de las variables, notemos que x1 ≤ 0, o bien, x′1 = −x1 ≥ 0. Por otro lado, vemos
que x2 ≥ 1, o también, x′2 = x2 − 1 ≥ 0. En cuanto a x3 , esta es libre de signo, por lo que podemos asumir que
x3 = x′3 − x′′3 con x′3 ≥ 0 y x′′3 ≥ 0. Reemplazamos en el modelo, y obtenemos,
Falta la función objetivo, donde z = −x1 + x2 + x′3 − x′′3 + 6, o más bien, z − 6 = −x1 + x2 + x′3 − x′′3 .
Podemos hacer el cambio de variable z ′ = z − 6. Finalmente, la formulación estándar es la siguiente:
26
usa parámetros tradicionales, indique explı́citamente su valor.
Solución
Para linearizar la función objetivo, primero debemos descomponer la sumatoria y la pitatoria para simplificar
términos:
K
ekx1 − e(k+1)x1
P
k=1 (ex1 − e2x1 ) + (e2x1 − e3x1 ) + (e3x1 − e4x1 ) + · · · + (eKx1 − e(K+1)x1 )
=
K+1 ex1 · 2ex1 · 3ex1 · · · (K + 1)ex1
kex1
Q
k=1
ex1 − e(K+1)x1
=
(K + 1)! · e(K+1)x1
ex1 eKx1 · ex1
= −
(K + 1)! · eKx1 · ex1 (K + 1)! · eKx1 · ex1
1 1
= Kx
−
(K + 1)! · e 1 (K + 1)!
−1
Eliminamos la constante y reemplazamos en (P ); nos queda el problema equivalente (P(1) ):
(K + 1)!
1
(P ) ∼ (P(1) ) máx
(K + 1)! · eKx1
Usando la equivalencia de función recı́proca, nos queda un problema que minimiza f (x1 ) = (K + 1)! · eKx1 .
Eliminamos la constante (K + 1)! y nos queda el problema (P(2) ):
Debemos sacar la función exponencial para terminar de linearizar la función objetivo. Para esto, podemos usar la
equivalencia de composición de funciones, ya que la función logarı́tmica es estrictamente creciente y ln(eKx1 ) =
Kx1 :
mı́n eKx1 ∼ mı́n ln(eKx1 ) ∼ mı́n Kx1 ∼ mı́n x1
27
Nos queda el problema equivalente (P(3) )
(P ) ∼ (P(3) ) mı́n x
En este caso, como el valor absoluto es ≥ 4, no es tan simple de linearizar como en otros casos, ya que se separa
en dos restricciones disjuntas:
x1 − x2 ≥ 4 (R2a )
máx{x1 + x2 , 5x1 + x2 − 32} ≤ 10 ∼ ó
x1 − x2 ≤ −4 (R2b )
28
Dado que introdujimos un nuevo parámetro, debemos calcular su valor. En este caso, dada la primera restricción,
el mayor valor que puede tomar cualquiera de las variables es 10, solo cuando la otra es 0. Por lo tanto: −10 ≤
x1 − x2 ≤ 10. Revisemos nuevamente los casos:
x1 − x2 ≤ −4 + 14 · 1 = 10
(P ) ∼ (P ′ ) mı́n x1
s.a: x1 + x2 ≤ 10 (R1)
5x1 + x2 − 32 ≤ 10 (R2)
x1 − x2 ≤ −4 + 14x3 (R3)
x1 − x2 ≥ 4 − 14(1 − x3 ) (R4)
x1 , x2 ≥ 0
x3 ∈ {0, 1}
(1) mı́n(x1 ,x2 )∈R2 {máx{|2x1 + 3x2 |, |x1 − x2 |}} s.a |x1 | + 2 máx{x1 , x2 } ≤ 1
(2) mı́n(x1 ,x2 )∈R2 {máx{|2x1 + 3x2 |, |x1 − x2 |}} s.a |x1 | − |x2 | ≥ 1
(3) máx(x1 ,x2 )∈R2 {2 mı́n{|x1 + x2 |, 2x2 } − |x1 − x2 |} s.a (x21 + x22 ) ≤ 2, |x1 | ≤ 1, |x2 | ≤ 1
Solución
mı́n u
(x1 , x2 )∈R2 , u∈R
s.a: |x1 | + 2 máx{x1 , x2 } ≤ 1
máx{|2x1 + 3x2 |, |x1 − x2 |} ≤ u
Luego, separamos la primera restricción en todos los casos posibles y eliminamos la función máximo de
29
la segunda restricción:
mı́n u
(x1 , x2 )∈R2 , u∈R
s.a: x1 + 2x1 ≤ 1
x1 + 2x2 ≤ 1
−x1 + 2x1 ≤ 1
−x1 + 2x2 ≤ 1
|2x1 + 3x2 | ≤ u
|x1 − x2 | ≤ u
Por último, sacamos los valores absolutos de las últimas dos restricciones:
(Q) mı́n u
(x1 , x2 )∈R2 , u∈R
s.a: x1 + 2x1 ≤ 1
x1 + 2x2 ≤ 1
−x1 + 2x1 ≤ 1
−x1 + 2x2 ≤ 1
2x1 + 3x2 ≤ u
−(2x1 + 3x2 ) ≤ u
x1 − x2 ≤ u
−(x1 − x2 ) ≤ u
El problema (1) puede ser reformulado en un problema lineal equivalente dado por (Q).
(2) Este problema tiene la misma función objetivo que (1), pero diferente dominio. Si eliminamos los valores
absolutos, nos quedan 4 restricciones disjuntas:
x1 − x2 ≥ 1
x1 + x2 ≥ 1
−x1 − x2 ≥ 1
−x1 + x2 ≥ 1
Para transformar el dominio en uno equivalente sin valor absoluto, es necesario agregar variables binarias,
obteniendo un problema lineal entero mixto. Esto se debe a que el dominio no es convexo: (1, 0) y (−1, 0)
son soluciones factibles, pero (0, 0) no lo es. Como el problema no es convexo, no puede ser reformulado
como un problema lineal equivalente.
Usamos la equivalencia de suma de funciones para linealizar la función objetivo, agregando dos variables
30
al problema:
∼ máx u1 + u2
(x1 ,x2 )∈R2 , u1 ,u2 ∈R
s.a: 2 mı́n{|x1 + x2 |, 2x2 } ≥ u1
−|x1 − x2 | ≥ u2
x21 + x22 ≤ 2
|x1 | ≤ 1
|x2 | ≤ 1
∼ máx u1 + u2
(x1 ,x2 )∈R2 , u1 ,u2 ∈R
s.a: 2(x1 + x2 ) ≥ u1 (R1)
−2(x1 + x2 ) ≥ u1 (R2)
2x2 ≥ u1 (R3)
−(x1 − x2 ) ≥ u2 (R4)
x1 − x2 ≥ u2 (R5)
x21 + x22 ≤ 2 (R6)
−1 ≤ x1 ≤ 1 (R7)
−1 ≤ x2 ≤ 1 (R8)
Lo único que falta linealizar es R6. Si analizamos el problema, podemos notar que R6 es redundante
debido a R7 y R8: restringen que (x)2 puede ser máximo 1, por lo que x21 + x22 ≤ 2 gracias a estas dos
restricciones. Por lo tanto, podemos eliminar R6. Finalmente, (3) puede ser reformulado en un problema
lineal equivalente dado por:
máx u1 + u2
(x1 ,x2 )∈R2 , u1 ,u2 ∈R
s.a: 2(x1 + x2 ) ≥ u1
−2(x1 + x2 ) ≥ u1
2x2 ≥ u1
−(x1 − x2 ) ≥ u2
x1 − x2 ≥ u2
−1 ≤ x1 ≤ 1
−1 ≤ x2 ≤ 1
31
CAPÍTULO 4
TEORÍA POLIEDRAL Y PROGRAMACIÓN LINEAL
4.1 Definiciones
P = {x ∈ Rn : Ax ≤ b}
C ⊂ Rn C = {x ∈ Rn : Ax ≤ 0}
Contiene al origen
Si x, y ∈ C ⇒ x + y ∈ C
32
Nota: Si un cono C solo contiene al origen se denomina cono “trivial”. De lo contrario, se denomina “regular” (es
decir, si el cono contiene un x ̸= 0).
C(X) = {y ∈ Rn : Ay ≤ 0}.
Nota: La definición Envoltura Convexa aplica para cualquier lista de vectores L, incluyendo listas infinitas e incon-
tables de vectores como un poliedro o un conjunto cualquiera.
p
X
x = x′ + y, donde x′ = λ i vi .
i=1
(i) (P ) es no acotado
(ii) (P ) tiene solución óptima y al menos uno de los vértices de su dominio alcanza el valor óptimo.
33
4.2 Ejercicios
(i) X := {x ∈ Rn : | ni=1 xi | ≥ 1}
P
(ii) X := {x ∈ Rn : −2 ≤ ni=1 xi ≤ 2}
P
(iv) X := x ∈ R2 : máx(x1 , x2 ) ≤ 3, x1 ≥ 0, x2 ≥ 0
(iv) X := {x ∈ R2 : x1 − x2 ≥ 1 y x1 + x2 ≥ 1}.
(v) X := {x ∈ R3 : max [x3 − min(x1 + 3, 2x2 ), |2x1 + 3x” |] ≤ 5}
Solución
(iii) Este conjunto no es un poliedro, puesto que posee desigualdades estrictas, lo que implica que no es
un conjunto cerrado.
(iv) La función máx(x1 , x2 ) ≤ 3 se puede descomponer en dos desigualdades lineales:
X := x ∈ R2 : x1 ≤ 3, x2 ≤ 3, x1 ≥ 0, x2 ≥ 0
Al poder representarse mediante una cantidad finita de restricciones lineales, el conjunto es un po-
liedro.
(v) El único punto factible es x = 0. Por lo tanto, X := {x ∈ Rn : xi ≤ 0, −xi ≤ 0, ∀i ∈ {1, . . . , n}}
representa al mismo conjunto.
(vi)
X := {x ∈ R3 : max [x3 − min(x1 + 3, 2x2 ), |2x1 + 3x2 |] ≤ 5}
34
es equivalente a:
x ∈ R3 : x3 − 5 ≤ min(x1 + 3, 2x2 )
: |2x1 + 3x2 | ≤ 5
x ∈ R3 : x3 − 5 ≤ x1 + 3
; x3 − 5 ≤ 2x2
; 2x1 + 3x2 ≤ 5
; −2x1 − 3x2 ≤ 5
Por lo tanto, como se puede representar mediante una cantidad finita de restricciones lineales, el
conjunto sı́ es un poliedro.
x ∈ R3 : x3 − 5 ≤ x1 + 3
; x3 − 5 ≤ 2x2
; 2x1 + 3x2 ≤ 5
; −2x1 − 3x2 ≤ 5
x3 ≤ x1 + 8
35
Determine si las siguientes afirmaciones son verdaderas o falsas. Justifique su respuesta demostrando la afirmación
si es verdadera o dando un contraejemplo si es falsa.
(b) Si los conjuntos C1 y C2 en Rn son conos poliedrales, entonces su intersección C1 ∩ C2 es un cono poliedral.
(c) El poliedro
X = {x ∈ Rn+1 : x0 ≥ 0; x0 + xk ≤ bk , ∀k ∈ {1, ..., n}}
definido por un parámetro b ∈ Rn tiene un único vértice dado por (0, b1 , ..., bn ).
Solución
(a) Falso. En palabras, lo que nos dice la afirmación es que todo poliedro no acotado contiene una recta que
ponderada por cualquier λ ∈ R seguirá perteneciendo a X. Se puede encontrar un contraejemplo simple.
Tomemos como poliedro el primer cuadrante en R2 :
X = {x ∈ R2 : x ≥ 0}
Se puede ver que para toda recta en X, si λ < 0 y decrece lo suficiente, esta se saldrá del poliedro.
Veámoslo gráficamente. Sea x1 = (0, 0) ∈ X, x2 = (1, 1) y λ ∈ R tales que x′ (λ) = (0, 0) + λ(1, 1).
Graficamos:
x2
x′
x1
En este caso, x′ (λ) se escapa de X para todo λ < 0. Si revisamos las otras rectas en el cuadrante X,
podemos notar que para todo x′ = x1 + λx2 , con x1 ∈ X, x2 ̸= 0, existe algún λ ∈ R tal que x′ ∈
/ X.
(b) Verdadero. Como C1 y C2 son conos poliedrales en Rn , entonces existen matrices A y B tales que:
C1 = {x ∈ Rn : Ax ≤ 0}
C2 = {x ∈ Rn : Bx ≤ 0}
A
Entonces, definiendo la matriz D = , se tiene que:
B
C1 ∩ C2 = {x ∈ Rn : Ax ≤ 0; Bx ≤ 0} = {x ∈ Rn : Dx ≤ 0}
36
(c) Verdadero. Un vértice es un punto en el dominio en el que se encuentra una cantidad de restricciones l.i.
activas igual a la dimensión del dominio. En este caso, como P está en Rn+1 , un vértice de P debe estar
definido por n+1 restricciones l.i. activas. El poliedro P está definido por exactamente n+1 restricciones,
por lo que el vértice, si existe, es tal que activa todas las restricciones del poliedro, es decir:
x0 + x1 = b1
x0 + x2 = b2
..
.
x0 + xn = bn
x0 = 0
¿Cuál de las siguientes direcciones d contenidas en el cono recesivo del poliedro certifica que el problema no está
acotado? Marque la(s) alternativa(s) correcta(s).
(a) d = (0, 0, 1)
(b) d = (−1, −1, −1)
(c) d = (0, −1, 0)
(d) d = (0, 0, 0)
(e) d = (−2, 0, 1)
Solución
d debe estar contenido en el cono recesivo del dominio de (P ), donde:
Y para certificar que el problema sea no acotado, debemos recordar el teorema de problema acotado.
Es decir, si d ∈ C(P ) y, además, la función objetivo evaluada en d da un valor negativo, entonce el problema
es no acotado. Esto significa que podemos ponderar d por un λ ∈ R tal que λ · d es una solución factible del
cono. Luego cT · (λd) → ∞ cuando λ → ∞, por lo que d es una dirección de escape.
Vamos a analizar cada alternativa de forma independiente. Debemos evaluar si d cumple con las 3 restricciones
de C(P ) y si f (d) > 0.
37
(a) d = (0, 0, 1):
(d) d = (0, 0, 0)
g1 (d) = 3x1 − 2x2 + 4x3 = 0 ≤ 0
g2 (
d) = x1 − 9x2 + 2x3 = 0 ≥ 0
g3 (d) = −2x1 − 5x2 + 3x3 = 0 ≥ 0
cT d = 3x1 − 2x2 − 3x3 = 0 ≮ 0
(e) d = (−2, 0, 1)
d es una dirección en el cono recesivo del poliedro que garantiza que (P ) sea no acotado.
38
(P ) mı́n f (x)
x∈X
Solución
Supongamos que la premisa es falsa. Entonces existe una solución óptima x∗ con valor óptimo v ∗ que puede ser
escrita como el punto medio de otras dos soluciones del dominio x1 , x2 ∈ X, es decir: x∗ = 0,5x1 + 0,5x2 .
Por concavidad estricta de f , se cumple:
Luego, 0,5 · f (x1 ) + 0,5 · f (x2 ) ≥ mı́n{f (x1 ), f (x2 )}, por lo que:
En otras palabras, al menos una solución factible en el dominio, x1 o x2 , posee mejor valor objetivo que el valor
óptimo → contradicción.
Lo anterior muestra que toda solución óptima de (P ) no puede ser escrita como punto medio de otras dos
soluciones factibles, por ende, toda solución óptima de (P ) es necesariamente un vértice de X.
Solución
(i) Primero, el poliedro no es vacı́o pues x0 = (1, 1, 1)T es tal que Ax0 ≤ b. Por otra parte, existe un vector
y = (0, 1, 1)T ̸= 0 que cumple Ay ≤ 0. Luego, si nos paramos en x0 y nos movemos en la dirección de
39
escape y una cantidad arbitraria t ≥ 0, nos encontraremos siempre en el espacio de soluciones factibles:
(ii) Se debe encontrar un y en el cono con producto positivo con el gradiente del objetivo. Esto se cumple
para
y = (0, 1, 1)T
pues
cT y = 2 > 0,
es decir, la función objetivo crece a medida que nos movemos en la dirección y y podemos movernos
infinitamente en esta dirección quedando dentro del espacio factible del dominio. Por ende, el problema
es no acotado.
Solución
(a) Hay 10 vértices en el poliedro D. En este poliedro, para que un punto sea considerado un vértice, debe
tener 4 restricciones l.i. en igualdad.
(0, 0, 0, 0) Activa las 4 restricciones de no negatividad.
(3, 0, 0, 0), (0, 3, 0, 0), (0, 0, 3, 0), (0, 0, 0, 3). Activan 3 no negatividades y una restricción l.i. junto
a ellas.
40
( 32 , 32 , 23 , 0), ( 32 , 23 , 0, 32 ), ( 32 , 0, 23 , 23 ), (0, 32 , 32 , 32 ). Activan una no negatividad y tres restricciones l.i.
adicionales.
( 32 , 32 , 23 , 32 ). Activan todas las restricciones adicionales a las no negatividades. Son 4 restricciones
l.i.
(b) (i) x0 = (0, 0, 0, 0) es un vértice del poliedro D, por lo cual tiene que existir dicho c. Por ejemplo
cT = (−1, −1, −1, −1) (la suma de los gradientes de restricciones activas en ese punto).
(ii) x0 = ( 23 , 0, 32 , 32 ). Al no ser un vértice del poliedro D, no existe un c tal que sea solución óptima. Si
es que fuera solución óptima, esta no se serı́a única.
(iii) x0 = ( 32 , 0, 23 , 23 ). Es un vértice del poliedro D. La suma de los gradientes de sus restricciones activas
logran un vector objetivo cuyo problema asociado solo se maximiza en ( 32 , 0, 32 , 32 ). Por ejemplo, la
suma de los gradientes c = (0, −1, 0, 0) + (1, 0, 1, 0) + (1, 0, 0, 1) + (0, 0, 1, 1) = (2, −1, 2, 2)
(c) El dominio es no vacı́o y compacto. Además el objetivo es continuo, por lo que existe solución. Es nece-
sario evaluar los vértices del poliedro D especificados en la parte (a), y especificar cuál obtiene el máximo
valor.
(0, 32 , 23 , 3 3 3 3 3
2 ), ( 2 , 2 , 2 , 2 ) → 27
2
27
De los vértices, los dos que alcanzan el valor óptimo 2 son x1 = (0, 32 , 32 , 32 ), x2 = ( 32 , 32 , 23 , 32 ).
Luego, el conjunto de soluciones óptimas, es toda aquella solución formada por:
Es decir: x∗ = ( 32 · λ, 32 , 23 , 32 )
2023
( )
2023
X 3
X := x∈R : xi ≤ 2; 0 ≤ xi ≤ , ∀i{1, · · · , 2023}
2
i=1
(a) Sea ei ∈ R2023 un vector con valor 1 en la componente i y 0 en el resto. ¿Cuáles de los siguientes puntos son
vértices de X? Explique.
(0, 0, 0, · · · , 0)
41
3
2 e1 es decir, ( 32 , 0, 0, · · · , 0)
3
2 e1 + 23 e2 , es decir, ( 32 , 32 , 0, 0, · · · , 0)
3
2 e1 + 21 e2
e1 + e2
(b) Obtenga una solución óptima para máx{ 2023
P
i=1 i · xi s.t. x ∈ X}.
Solución
(b) De la pregunta anterior es fácil deducir que los vértices del poliedro son: 0, 23 ei para todo i ∈ {1, · · · , 2023},
3 1 3 1
2 ei + 2 ej para todo i, j ∈ {1, · · · , 2023} : i ̸= j. Luego, el vértice óptimo debe ser 2 e2023 + 2 e2022 . El
3·2023+2022
valor óptimo es 2
(a) Suponga que al tratar de resolver el problema (P ) mı́n{cT x : Ax ≥ b}, detectó que existe un vector y ̸= 0
tal que Ay ≥ 0 y cT y < 0. Explique por qué (P ) no puede tener solución óptima.
(b) Sea X = {x ∈ Rn : Ax ≥ 0}.
(i) Demuestre que si x̄ ∈ X es tal que x̄ ̸= 0, entonces x̄ no es vértice de P
(ii) ¿Es verdad que x = 0 es necesariamente vértice de P ? Demuestre o refute con un contraejemplo.
Solución
(a) Si (P ) es infactible, entonces claramente no tiene solución óptima. Si (P ) es factible, las condiciones
sobre r implican que (P ) es no acotado, puesto que r ̸= 0, lo que implica que no hay cota inferior, por lo
que tampoco tiene solución óptima.
(b) (i) Si x̄ ̸= 0 pertenece a X, entonces el origen y 2x̄ pertenecen a X. Luego, x̄ es el promedio de dos
vectores distintos pertenecientes a X por lo que no puede ser vértice.
(ii) No necesariamente. Por ejemplo, tomando X = {(x1 , x2 ∈ R2 : x1 ≥ 0)} vemos que el origen no
es vértice de X, esto porque en el origen solo se tendrı́a una restricción l.i. activa y no dos como es
necesario en dos dimensiones.
42
Sea X = {x ∈ Rn : Ax ≤ b} = ̸ ∅. Decimos que X contiene una lı́nea (o recta) si existe y ∈ Rn no nulo, tal que x
+ δy pertenece a X para todo x ∈ X y para todo escalar δ (positivo o negativo).
Solución
Ya que para todo punto en el poliedro podemos sumarle o restarle y tanto como queramos.
x2
6
5
4
3
2
1
−6 −5 −4 −3 −2 −1 1 2 3 4 5 6 x1
−1
−2
−3
−4
−5
−6
(b) Para determinar que no tiene vértices basta demostrar que para todo punto x ∈ X existen dos puntos x1 y
x2 tales que x es su punto medio. Podemos notar que basta usar x1 = x − y y x2 = x + y (sabemos que
tanto x1 y x2 pertenecen a X ya que y es lı́nea).
(c) Supongamos que y define una lı́nea en X. Esto implica que x + δy es factible para cualquier x ∈ X y
todo δ ∈ R. Esto implica que A(x + δy) ≤ b donde Ax ≤ b. Luego, se cumple que δ · Ay ≤ b − Ax ≤ 0
para todo δ ∈ R. En particular, si δ = 1, se cumple que Ay ≤ 0. Por otro lado, si δ = −1, se cumple que
−Ay ≤ 0. Luego, Ay = 0.
43
Problema 4.10. [I1 2020-2] Método gráfico y geometrı́a poliedral
Sea (PE ) un problema lineal en forma estándar, donde h1 , e2 y h3 son las holguras/excesos de las restricciones R1 ,
R2 y R3 , respectivamente:
(a) Grafique el dominio factible del problema original (P ) (es decir, del problema antes de haber sido transfor-
mado a su equivalente estándar), identificando cada una de las soluciones básicas factibles del problema con
sus respectivas bases (variables que conforman la base). Considere que, en el problema original, todas las
variables eran de naturaleza no negativa.
(b) Determine la(s) solución(es) óptima(s) para las siguientes funciones objetivos:
(c) Determine el rango de valores para δ de modo que la dirección d = (δ, 1) corresponda a un rayo de escape
para el problema original.
(d) Demuestre que agregando la restricción 3x1 + 2x2 ≤ 24 al dominio del problema, este siempre poseerá
solución óptima para toda función objetivo continua.
Solución
(a) Primero, construimos el problema original (sin holguras) para que nos quede un dominio en 2 dimensio-
nes:
(P ) ??? f (x)
s.a x1 − x2 ≤ 5 (R1 )
3x1 + x2 ≥ 6 (R2 )
−2x1 + x2 ≤ 1 (R3 )
x1 , x2 ≥ 0
El siguiente gráfico contiene en azul las restricciones y en celeste el dominio factible, además de marcar
los tres vértices factibles:
44
x2
R3
6
5 R1
3 (1, 3)
1
(2, 0) (5, 0)
1 2 3 4 5 6 7 8 9 10 x1
R2
Para identificar las bases de un vértice, se puede proceder de forma visual o numérica. Visualmente,
se pueden reconocer las restricciones y holguras que están activas (es decir, las que forman parte de
la base) y aquellas que no lo están. De manera numérica, se pueden sustituir los valores del vértice en
las restricciones: aquellas que no se anulen deben estar en la base. Esto aplica tanto para las variables
originales como para las variables de holgura.
Las soluciones básicas factibles se encuentran en los vértices, estos son:
Vértice (2, 0) al evaluarlo en las restricciones se obtienen los siguientes valores de las variables
x1 = 2, x2 = 0, h1 = 3, h2 = 0, h3 = 5. Por lo tanto, la base correspondiente es aquella que
contiene los valores no nulos: B = {x1 , h1 , h3 }.
Vértice (5, 0) al evaluarlo en las restricciones se obtienen los siguientes valores de las variables
x1 = 5, x2 = 0, h1 = 0, h2 = 9, h3 = 11. Por lo tanto, la base correspondiente es aquella que
contiene los valores no nulos: B = {x1 , h2 , h3 }.
Vértice (1, 3) al evaluarlo en las restricciones se obtienen los siguientes valores de las variables
x1 = 1, x2 = 3, h1 = 7, h2 = 0, h3 = 0. Por lo tanto, la base correspondiente es aquella que
contiene los valores no nulos: B = {x1 , x2 , h1 }.
(b) El dominio se mantiene. Debemos agregar la función objetivo al gráfico. Para esto, lo más fácil es encon-
trar la curva de nivel f (x) = 0 (ası́ encontramos la pendiente) y trasladarla en la dirección de ascenso o
de descenso, dependiendo del objetivo del problema. La dirección de descenso está dada por el gradiente
negativo de f (x) y la dirección de ascenso es el gradiente de f (x), es decir: si se minimiza debemos mo-
ver la función objetivo en el sentido de −∇f (x), mientras que si se maximiza debemos mover la función
objetivo en el sentido de ∇f (x).
45
x2 f (x) = −x
R3
6
5 R1
3 (1, 3)
2
−∇f
1
(2, 0) (5, 0)
1 2 3 4 5 6 7 8 9 10 x1
R2
x2
R3
6
5 R1
3 (1, 3)
2 f (x) = x1 − 2x2
1
∇f
(2, 0) (5, 0)
1 2 3 4 5 6 7 8 9 10 x1
R2
46
x2
R3
6
−∇f f (x) = −4x1 + 2x2
5 R1
3 (1, 3)
1
(2, 0) (5, 0)
1 2 3 4 5 6 7 8 9 10 x1
R2
El problema tiene múltiples soluciones. Corresponde al vértice (1, 3) y todo el rayo extremo deter-
minado por la restricción R3 .
(c) Llamemos DP al dominio del problema (P ). Cualquier vector perteneciente al cono recesivo de DP es
un rayo de escape de DP . En particular, los rayos extremos son r1 = (1, 2) y r2 = (1, 1), determinados
por las restricciones R1 y R3 , respectivamente. Es decir, son rayos de escape del dominio los vectores
pertenecientes al cono:
Entonces, d = (δ, 1) debe cumplir las restricciones de C(S) para que corresponda a un rayo de escape:
δ−1≤0⇒δ ≤1
−2δ + 1 ≤ 0 ⇒ δ ≥ 0,5
(d) Agregamos la nueva restricción al gráfico. Esta contiene al origen. El nuevo dominio escrito de la forma
y = mx + n es:
x2 = x1 − 5 (R1 )
x2 = −3x1 + 6 (R2 )
x2 = 2x1 + 1 (R3 )
x2 = −3/2 · x1 + 12 (R4 )
x1 , x2 ≥ 0
Graficamos:
47
x2
R3 R4
6
5 R1
3 (1, 3)
1
(2, 0) (5, 0)
1 2 3 4 5 6 7 8 9 10 x1
R2
Por lo tanto, si se tiene una función objetivo continua, por teorema de Bolzano-Weierstrass, el problema
admite solución óptima.
48