Ejercicios Resueltos Optimizacion PDF
Ejercicios Resueltos Optimizacion PDF
Ejercicios Resueltos Optimizacion PDF
Mtodos de Optimizacin
Introduccin V
I Programacin Lineal 1
1. Introduccin a la Programacin Lineal 2
1.1. Solucin Grca . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2. Modelacin de problemas simples . . . . . . . . . . . . . . . . . . . . . 7
1.3. Modelacin de problemas ms complejos . . . . . . . . . . . . . . . . . 10
1.4. Ejercicios Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3. El Mtodo Simplex 37
3.1. Simplex Analtico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.2. Forma Tabular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.3. Ejercicios Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
4. Dualidad 58
4.1. Teorema de Holguras Complementarias . . . . . . . . . . . . . . . . . . 61
4.2. Mtodo Dual-Simplex . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
5. Post Optimizacin 72
6. Optimizacin en Redes 76
6.1. Transporte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
6.2. Problemas de Rutas Mnimas . . . . . . . . . . . . . . . . . . . . . . . 81
6.3. Digrafos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
6.4. Algoritmo del camino ms corto de Dijsktra . . . . . . . . . . . . . . . 82
ii
NDICE GENERAL iii
II Programacin No Lineal 90
7. Condiciones de Opimalidad 91
7.1. Sin Restricciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.2. Con Restricciones de Igualdad . . . . . . . . . . . . . . . . . . . . . . . 99
7.3. Con Restricciones de Desigualdad . . . . . . . . . . . . . . . . . . . . . 102
Prefacio
iv
Introduccin
v
vi INTRODUCCIN
Parte I
Programacin Lineal
1
Captulo 1
Introduccin a la Programacin
Lineal
0 1
mn z = ct x
@ s:a : Ai x bi ; i = 1; : : : ; m A
x 0
ct x = k, con k 2 R
2
1.1. SOLUCIN GRFICA 3
Fig. 1
El conjunto factible es el poliedro limitado por las restricciones, que es no acotado (ver
Fig. 2):
Fig. 2
Repitiendo el mismo procedimiento se llega a que el mximo, cuyo valor es 24, se obtiene
en el punto (7; 4), como lo indica la curva de nivel que se traz (linea punteada, ver Fig.
2) en direccin rz y a que no existe mnimo, puesto que en la direccin rz, siempre
podemos encontrar un punto que haga decrecer la funcin objetivo y pase por una curva
de nivel, este es el caso en que la funcin objetivo no est acotada inferiormente.
Ejercicio 3 0 1
max z1 = 3x1 + 2x2
B m n z2 = 2x1 2x2 C
B C
B s:a: C
B C
B 2x1 x2 2 C
B C
@ x1 + 2x2 8 A
x 1 ; x2 0
Solucin.
1.1. SOLUCIN GRFICA 5
x1 x2 x1 x2
R1 : 0 2 ; R2 : 0 4
3 8 8 0
4 18
(0; 2) ; (0; 0) ; (8; 0) ; ;
5 5
El ltimo punto sale de: R1 \ R2 : Adems tenemos que rz1 = (3; 2) y rz2 =
(2; 2) : De aqui sale que el mximo de z1 se alcanza en (8; 0) y el mnimo de z2
en (4=5; 18=5) :
Una fbrica est a desarrollando una nueva lnea de galletas de avena y chocolate.
Las galletas deben tener al menos 5 gramos de protenas, pero no ms de 5 gramos de
carbohidratos, y 3 gramos de grasas saturadas. Construye un modelo lineal que permita
determinar la cantidad de cada ingrediente que se debe utilizar para que se satisfagan
6 CAPTULO 1. INTRODUCCIN A LA PROGRAMACIN LINEAL
Solucin. Para este caso lo primero que observamos que las variable x2 y x3 son
irrestrictas, por ende hacemos el cambio de variable x2 = x+
2 x2 y x3 = x+3 x3
1.2. MODELACIN DE PROBLEMAS SIMPLES 7
donde x+ +
2 ; x2 ; x3 ; x 3 0 y reemplazamos en nuestro (P L) obteniendo
0 1
max z = 2x1 x+2 x2 + 4 x+ 3 x3
B s:a : C
B C
]
(P L) : B
B 5x1 + 2 x+ 2 x2 3 x+ 3 x3 x4 = 7 C
C
@ +
2x1 2 x2 x2 + x3 x3 + x5 = 8 + A
x 1 ; x+ +
2 ; x 2 ; x3 ; x3 ; x4 ; x5 0
Solucin.
a. 0 1
max z = 3 10 5 xA + 2 10 4 xB
B s:a C
B C
B 50 xA + xB 150 C
B C
@ xA xB A
x A ; xB 0
Aqu las funcin objetivo maximiza las vitaminas. Explicacin: Debido a que
1000 gr contienen 30 mg de vitaminas, xA gramos (del producto A) contienen
30mg
1000gr
xA = 3 10 5 xA vitaminas, anlogo para las vitaminas de tipo B.
1.2. MODELACIN DE PROBLEMAS SIMPLES 9
b. 0 1
450 150
max z = 1000 xA + 1000 xB
B s:a C
B C
B 50 xA + xB 150 C
B C
@ xA xB A
x A ; xB 0
Aqu las funcin objetivo minimiza las caloras. Explicacin: Debido a que 1000 gr
20mg
contienen 450 mg de caloras, xA gramos (del producto A) contienen 100gr xA =
4
2 10 xA vitaminas, anlogo para las caloras del tipo B.
0;8x1 + x2 98;
10 CAPTULO 1. INTRODUCCIN A LA PROGRAMACIN LINEAL
x2 26;
La materia prima D est asociada a la materia prima C, lo que signica que por cada
Kg que se compre de C se debe adquirir 1.5 Kg de D.
1.3. MODELACIN DE PROBLEMAS MS COMPLEJOS 11
Escriba un modelo de programacin lineal que permita determinar cuanta materia prima
de cada clase debe comprarse de modo que se obtenga una tonelada de mezcla (1000
Kg.) con un contenido del facor F de por lo menos 18 %.
Solucin. Denamos las siguientes variables:
Ejercicio 12 Una Granja dispone de 10.000 ha. para el cultivo de maz, papas y frijoles,
adems de que una parte de esa supercie puede utilizarse para la cra de cerdos. El
manager de la granja debe decidir cuantas h. debe dedicar a cada cultivo y cuntas
a la cra de cerdos para maximizar la ganancia. El costo de la mano de obra es de $5
la hora/hombre para cualquiera de las tareas de cultivo y $7 la hora/hombre para la
crianza de cerdos. El nmero mximo de horas/hombre que se puede contratar es de
48.000 y las necesidades de horas/hombre por h. para cada cultivo y para los cerdos
se muestran en la siguiente tabla:
horas / hombre
Maz 3 por h
Papa 2 por h
Frijol 2 por h
Cerdo 4 por c/ 10 cerdos
El costo de cada saco de semillas de maz, papa y frijol y el rendimiento de cada saco
de semillas por hectrea viene dado por:
Adems, se sabe que para sembrar una h. de maz, papa o frijoles se requiere de 2,5
y 3 sacos de semillas respectivamente. La disponibilidad de semillas en el mercado es
12 CAPTULO 1. INTRODUCCIN A LA PROGRAMACIN LINEAL
ilimitada. Si se decide criar cerdos, estos sern alimentados con maz, el cual puede
comprarse al precio del mercado y/o puede utilizarse parte del que se cultive en la
granja. Las necesidades de maz para alimentar a los cerdos hasta que alcancen un peso
de 0,3 Ton., adecuado para la venta, es de 3 ton. por cerdo y se sabe que 10 h. es lo
que se necesita para criar 20 cerdos. El Precio de venta de cada producto(por Ton.) en
el mercado est dado en la tabla siguiente:
Precio de Mercado / ton
Maz $25
Papa $15
Frijol $10
Cerdo $200
Denicin de variables:
Con lo que respecta a las restricciones, tenemos 3 limitaciones: con las hora/hombre,
el terreno disponible, y el nmero de
1. Terreno disponible:
4
3 (xM 1 + xM 2 ) + 2xP + 2xF + 2 (xC1 + xC2 ) 48;000 (h=h)
10
20
3 xC1 5xM 1
10
4
CHH = 5 (3xM 1 + 3xM 2 + 2xP + 2xF ) + 7 2 (xC1 + xC2 )
10
14 CAPTULO 1. INTRODUCCIN A LA PROGRAMACIN LINEAL
Ejercicio 13 El dueo de una empresa tiene que ir al banco a sacar dinero para pagar
el sueldo de sus empleados. El nmero de empleados que tiene junto con el sueldo
respectivo est descrito en la siguiente tabla:
El pago de los sueldos puede efectuarse utilizando billetes de 1.000, 5.000 y 10.000
pesos, pero el banco pone las siguientes restricciones para el uso de estas denomina-
ciones:
1. Los billetes de 1.000 deben ser, a lo ms, el doble de los billetes de 5.000.
2. Los billetes de 5.000 deben ser al menos, la cuarta parte de los de 10.000.
Modelo: 0 P P 1
mn z= xij
B i2f1;5;10g j2fA;B;Cg C
B C
B s:a C
B C
B 1000x1A + 5000x5A + 10000x10A = 555 25 C
B C
B 1000x1B + 5000x5B + 10000x10 = 351 10 C
B C
B 1000x1C + 5000x5C + 10000x10C = 225 25 C
B C
B x1A + x1B + x1C 2 (x5A + x5B + x5C ) C
B C
B x5A + X5b + X5c 14 (x10A +
Px10B +P x10C ) C
B 40 C
B x10A + x10B + x10C 100
xij C
@ i2f1;5;10g j2fA;B;Cg A
xij 0; i 2 f1; 5; 10g ; j 2 fA; B; Cg
Restricciones de mezcla:
Scotch no ms del 20 % de C :
Kilt no ms del 80 % de C :
Tartan no ms del 50 % de C :
No ms de 2000 lts. de A :
No ms de 2500 lts. de B :
No ms de 1200 lts. de C :
Restricciones de positividad:
xij 0; i = A; B; C; j = S; K; T:
Funcin objetivo:
1.3. MODELACIN DE PROBLEMAS MS COMPLEJOS 17
Total de ingresos:
Total de gastos:
Total de ganancia:
z=I G:
Ejercicio 15 Una empresa fabricante de tostadoras elctricas debe tomar una decisin
sobre la produccin de un nuevo modelo. La empresa tiene la posibilidad de emplear
3 tcnicas alternativas de produccin: manula, semi-automtica y robotizada. Los re-
querimientos de cada tcnica para fabricar una tostadora se resumen en el siguiente
cuadro:
TCNICA DE ENSAMBLADO
Manual Sem-Aut. Robot.
Mano de Obra Especializada 1 min 4 min 8 min
Mano de Obra No-Especializada 40 min 30 min 20 min
Tiempo de Taller de Ensamblado 3 min 2 min 4 min
Ejercicio 16 Bajo las condiciones de trabajo normales una fbrica produce sobre 100
unidades de cierto producto en cada uno de los 4 periodos de tiempo consecutivos a un
costo que vara de periodo en periodo de tiempo como se muestra a continuacin.
Unidades adicionales pueden ser producidas en horas extras de trabajo. La cantidad
mxima y el costo son mostradas en la siguiente tabla, junto con el pronstico de
demanda para el producto en cada uno de los 4 periodos de tiempo
Es posible manejar 70 unidades del producto en stock desde un periodo hasta el siguiente
a un costo de $1;5 por unidad por periodo.
Formule(Sin resolver) un Problema de programacin lineal que permita determinar
el calendario de produccin y almacenamiento que satisfagan las demandas durante los
cuatro perodos de tiempo a un costo mnimo, sabiendo que al comienzo del perodo 1
tenemos 15 unidades en stock.
m n z = (6x1 +4x2 +8x3 +9x4 )+(8y1 +6y2 +10y3 +11y4 )+(1;5z1 +1;5z2 +1;5z3 +1;5z4 )
1.4. EJERCICIOS PROPUESTOS 19
4. A una persona que quiere adelgazar se le ofrecen dos productos, A y B para que
tome una mezcla de ambos con las siguientes recomendaciones: No debe tomar
ms de 150 gr. ni menos de 50 gr., la cantidad de A debe ser igual o superior a la
de B. Si 1000 gr. de A contiene 30 mg de vitaminas y 450 caloras; y 100 gr de B
contienen 20 gr. y 150 caloras:
Precio Venta
Marca Especicacin
por litro
No menos del 60 % de A
A $ 680
No ms del 20 % de C
No menos del 15 % de A
B $ 570
No ms del 80 % de C
C No ms del 50 % de C $ 450
mercado indican que hay que producir por lo menos 26 unidades de E2 y la can-
tidad a producir de E1 debe ser por lo menos 3 veces la cantidad que se produzca
de E2 . Formule un modelo de programacin lineal que, respetando las condiciones
de produccin y de mercado, permita maximizar las ganancias de la Empresa en
la prxima semana.
7. Una empresa fabricante de tostadoras elctricas debe tomar una decisin sobre
la produccin de un nuevo modelo. La empresa tiene la posibilidad de emplear 3
tcnicas alternativas de produccin: manula, semi-automtica y robotizada. Los
requerimientos de cada tcnica para fabricar una tostadora se resumen en el sigu-
iente cuadro:
TCNICA DE ENSAMBLADO
Manual Sem-Aut. Robot.
Mano de Obra Especializada 1 min 4 min 8 min
Mano de Obra No-Especializada 40 min 30 min 20 min
Tiempo de Taller de Ensamblado 3 min 2 min 4 min
La materia prima D est asociada a la materia prima C, lo que signica que por
cada Kg que se compre de C se debe adquirir 1.5 Kg de D.
10. Un comerciante compra azcar a granel y vende al detalle. Para venderla tiene
dos alternativas: envases de 1 kg y envases de 5 kg. El precio de venta es $300
y $250 por kg respectivamente, y en el mercado del azcar al detalle se pueden
vender 20.000 kg en envases de 1 kg y 17.000 en envases de 5 kg. Debido a un
contrato anterior se deben entregar 5.000 kg en envases de 5 kg a un determinado
cliente. El comerciante se puede abastecer de azcar desde dos proveedores. El
primero le puede vender hasta 15.000 kg a un precio de $90 por kg, y el segundo
le ofrece la cantidad de azcar que el comerciante desee, pero a un precio de $110
por kg y debido a requerimientos de sus distribuidores el comerciante debe vender
menos del tercio del azcar en envases de 1 kg. Adems, suponga que el precio
de los envases y el proceso de envasado son nulos,y que el comerciante no tiene
azcar almacenada y vende todo el azcar que compra.
Formule un problema de programacin lineal que permita al comerciante decidir
cual es el plan de abastecimiento y ventas de modo de obtener el mayor benecio
en su negocio.
12. Una empresa fabrica y tambin importa cierta mercanca (por ejemplo televisores)
y la distribuye entre sus clientes. Las cantidades disponibles de la mercanca en
los almacenes suministradores extranjeros y los pedidos de los clientes para los
primeros 6 meses del ao son conocidos y vienen dados en la tabla siguiente:
La empresa tambin puede fabricar las mercancas, pero por su tecnologa el costo
le sale a $200 cada una y tiene un lmite de produccin de 40 unidades mensuales.
Los costos por unidad de la mercanca importada y el precio de venta varan segn
el mes de ao y estn dados por:
Suponga que las mercancas fabricadas en un mes slo pueden venderse en el mes
siguiente, que las compras se realizan a principios de cada mes y las ventas a los
clientes al nal de cada mes. Entonces es posible que haya que tener mercancas
almacenadas, con lo cual se incurrira en un costo de almacenamiento de $5 por
unidad. A principios de Enero la empresa dispone de 10 unidades de la mercanca
en almacn.
Plantee un modelo de PL que resuelva el problema de cmo planicar la produc-
cin y las compras de cada mes de modo de satisfacer las demandas de los clientes
y obtener una ganancia mxima.
26 CAPTULO 1. INTRODUCCIN A LA PROGRAMACIN LINEAL
13. Una empresa que fabrica computadoras debe satisfacer las demandas de los prx-
imos cuatro trimestres, que vienen dadas por:
Trimestre Demanda
1 700
2 1500
3 1000
4 800
Ejercicio 17 Un empresario cuenta con tres millones de pesos para invertir en una
microempresa de produccin de mermelada que tiene en mente. Para empezar, fabricar
a mermelada de frutilla y frambuesa, ya que son las de mayor demanda en la zona. Un
estudio de mercado le indica que, para satisfacer la demanda de la localidad deber a
producir las siguientes cantidades mensuales de mermelada, expresada en kilos en la
siguiente tabla:
Adems, ha descubierto que la calidad de la fruta var a dependiendo la poca del ao,
lo cual hace que el kilo de fruta rinda una cierta cantidad de kilos de mermelada, cosa
que se aprecia en la siguiente tabla:
P = x 2 R2 : x 1 + x 2 4; x2 4; x1 1
a) x1 = (2; 3)T
b) x2 = (2; 2)T
c) x3 = (1; 3)T
d) x4 = (1; 2)T
a) d1 = (0; 1)T
b) d2 = (1; 1)T
c) d3 = (1; 2)T
d) d4 = (2; 0)T
Por ltimo, use el teorema de Caratheodory, para expresar el punto x = (3; 4) co-
mo una combinacin convexa de los puntos extremales y una combinacin positiva
de las direcciones extremales de P.
29
30 CAPTULO 2. CONJUNTOS CONVEXOS Y POLIEDROS
Solucin.
luego
x1 = (2; 3)T ! x1 = (2; 3; 1; 1; 1)T
x1 no es extremal, ya que al menos n m=5 3 = 2 de sus componentes
deberian ser nulas.
b) x2 = (2; 2) : Anlogo al caso anterior, busquemos las componentes que faltan:
det (B) = 1 6= 0.
2. Recordemos que :
d 2 Pe1 () Ad = 0; d 0y
B 1 aj
d 2 Pe1
e
() d 2 Pe1 y 9B, A = [B; N ] : d = ; 0, det (B) 6= 0;
ej
B 1 aj 0, aj 2 N
a) d1 = (0; 1)T . Busquemos el equivalente en Pe; esto es, debemos buscar las
componentes en el poliedro homogneo:
8 9
< d1 + d2 d3 = 0 =
d 2 R5 : d2 d4 = 0 ; d 0 ;
: ;
d1 d5 = 10
0+1 d3 = 0 =) d3 = 1
1 d4 = 0 =) d4 = 1
0 d5 = 0 =) d5 = 0
luego obtenemos el punto (0; 1; 1; 1; 0)T 0; debemos vericar que al menos
tiene n m 1 = 1 ceros, los cual es el caso, adems haciendo e4 = (1; 0) y
tomando B como
0 1
1 1 1
B = @ 0 1 0 A =) det (B) = 1
1 0 0
e
luego el punto d1 2 P1 :
32 CAPTULO 2. CONJUNTOS CONVEXOS Y POLIEDROS
Solucin.
2 1 2 0
B1 = y B2 =
1 0 1 1
2D11 2D12 + D2 + D3 = 0;
D11 D12 + 2D2 D4 = 0:
4 + 1 + D3 = 0;
2 + 2 D4 = 0;
2 1
B= ; det(B) = 1:
1 0
P = x 2 R2 : Ax b; x 0 :
Solucin.
a) Una posible grca en R2 del poliedro se muestra en la gura. Los puntos ex-
tremales estn marcados y las rectas paralelas y = 0; y = 3 obligan a que la nica
direccin extremal sea d = (1; 0)t : Pueden construirse multitud de ejemplos de
esta forma, poniendo una recta paralela al eje X e intersectando otras dos rectas
en un punto que est entre el eje X y la recta paralela al eje.
F ig:
b) Una recta intersecta los ejes en los puntos (0; 3); (1; 0) y la otra en los puntos
(3; 0); (0; 1). Las respectivas ecuaciones se escriben:
y = 3(x 1);
1
y = (x 3):
3
Por lo tanto, el poliedro se dene por las desigualdades:
8 9
>
> 3x + y 1 >
>
< =
x + 3y 3
P = (x; y) : ;
>
> 0 y 3 > >
: ;
x 0:
y las matrices que lo representan en forma cannica sern:
0 1 0 1
3 1 1
A= @ 1 3 A; b = @ 3 A
0 1 3
P = x 2 R2 : Ax b; x 02 :
a) Dibuje el grco de P y halle, con ayuda del grco, todos los puntos extremales
y todas las direcciones extremales de P:
b) Utilice el teorema que caracteriza los puntos extremales para vericar que el vector
(x1 ; x2 ) = (0; 5)t es un punto extremal de P:
c) Utilice el teorema que caracteriza las direcciones extremales para vericar que el
vector (d1 ; d2 )t = (2; 5)t es una direccin extremal de P:
Solucin.
a) En el grco:
GRAFICO 1
se encuentran fcilmente los puntos extremales x1 = (0; 3)t y x2 = (0; 5)t . Las
direcciones extremales son los vectores directores de las rectas 5x1 + 2x2 = 10 y
3x1 5x2 = 15, en la direccin del cuadrante positivo. Son entonces d1 = (2; 5)t
y d2 = (5; 3), o cualquier mltiplo positivo de ellas.
que tiene:
3 5 1 0 15
A= ; b= :
5 2 0 1 10
Con respecto a P~ , los puntos extremales hallados tienen coordenadas:
5 0
B=
2 1
es no singular, det(B) = 5 6= 0.
36 CAPTULO 2. CONJUNTOS CONVEXOS Y POLIEDROS
B 1b xB
x= =
0 0
llamando z = cTB B 1 b; zN
T
= cTB B 1 N y Y = B 1 N nos queda
0 T
1
mn z=z zN cTN xN
B s:a C
(P L)F E :B
@
C
A
x B = x B Y xN
(xB ; xN ) 0
37
38 CAPTULO 3. EL MTODO SIMPLEX
El punto equivalente ]
a (0; 0) en (P L) es x = (0; 0; 150; 200; 175)T . Los datos del
problema son
0 1 0 1
1 2 1 0 0 150
A= @ 3 4 0 1 0 A ; b = @ 200 A ; c = (10; 11; 0; 0; 0)T
6 1 0 0 1 175
Aqu 0 1
1 2
N1 = @ 3 4 A
6 1
Puesto que z2 c2 = 11 es mayor que cero, la variable bsica que debe ingresar a la base
es x2 : Y la que sale, la determinamos de la siguiente manera: Introducimos la variable
3.1. SIMPLEX ANALTICO 39
500 225
Observacin 3.1 Si partimos en el punto ; 7
21 0
obtenemos
1 el punto equivalente en
1 2 1
]
(P L) : 500 ; 225 ; 1300 ; 0; 0 con base asociada B = @ 3 4 0 A luego IB = f1; 2; 3g ; IN 0 =
21 7 21
6 1 0
f4; 5g luego tenemos
0 1 4
10 1 0 1 4
1
0 21 21
1 2 1 0 0 1 0 0 21 21
@ 0 2 1 A@
3 4 0 1 0 A = @ 0 1 0 72 1 A
,
7 7 7
11 2 11 2
1 21 21
6 1 0 0 1 0 0 1 21 21
luego el punto x = 45 ; 18 ]
; 0; 0 es una SFB de (P L): Expresemos el problema en forma
5
explcita con respecto a esta base:
2 1 2 1
5 5
2 1 1 0 1 0 5 5
1 2 = 1 2
5 5
1 2 0 1 0 1 5 5
luego
4 2 1
x1 = + x3 + x4
5 5 5
18 1 2
x2 = x3 + x4
5 5 5
Reemplazando en la funcin objetivo, obtenemos:
4 2 1 18 1 2
z = 3 + x3 + x4 + 2 x3 + x4
5 5 5 5 5 5
4 7 48
= x3 x4
5 5 5
El problema es no acotado debido a que
1 7
Y4 = 5 0 y z4 c4 = >0
2
5 5
Explicar esto ltimo con ms detalles. Tenemos adems una direccin extremal visible:
1 2
d= Y 4T ; e4 = ; ; 0; 1
5 5
puesto que
T
2 1 1 0 1 2 0
; ; 0; 1 =
1 2 0 1 5 5 0
con
2 1
B= ; y e4 = (0; 1) .
1 2
3.2. FORMA TABULAR 43
x1 x2 x3 x4 x5 x6 b
x4 h2i 3 1 1 0 0 3
x5 1 2 1 0 1 0 4
x6 1 3 0 0 0 1 9
z 2 5 1 0 0 0 0
x1 x2 x3 x4 x5 x6 b
x1 1 3=2 1=2 1=2 0 0 3=2
x5 0 7=2 h1=2i 1=2 1 0 11=2
x6 0 3=2 1=2 1=2 0 1 15=2
z 0 8 2 1 0 0 3
x1 x2 x3 x4 x5 x6 b
x1 1 5 0 1 1 0 7
x3 0 7 1 1 2 0 11
x6 0 2 0 1 1 1 2
z 0 22 0 3 4 0 25
con IB2 = f2; 4; 1g ; IN2 = f3; 5g y x2 = (12=7; 15=7; 0; 73=7; 0)T , como observamos el
la tabla simplex anterior, los zj cj 0, 8j 2 IN2 = f3; 5g, luego tambin x2 es un
ptimo para el problema. Concluimos que estamos en precencia de solucin mltiple, y
de hecho cualquier punto que tomemos del segmento que une x1 con x2 es un ptimo,
es decir, cualquier punto que se pueda escribir como combinacin convexa de dichos
puntos:
x1 + (1 ) x2 ; con 2 [0; 1] :
3.2. FORMA TABULAR 45
x1 x2 x3 x4 x5 b
x4 h2i 3 2 1 0 3
x5 1 1 1 0 1 5
z 3 2 1 0 0 0
x1 x2 x3 x4 x5 b
x1 1 3=2 1 1=2 0 3=2
x5 0 1=2 2 1=2 1 13=2
z 0 13=2 2 3=2 0 9=2
2 3 2 1 0
A=
1 1 1 0 1
y 0 1
3=2
B 1 C
2 3 2 1 0 B C 0
Ad = B 0 C= :
1 1 1 0 1 B C 0
@ 0 A
1=2
46 CAPTULO 3. EL MTODO SIMPLEX
^
a) Escriba el problema (P L1 ) en la forma estndar (P L1 ):
^
c) Escriba el problema (P L1 ) en forma explicita con respecto a la base obtenida en
b).
Solucin.
1 1 2 4
B x0B = b () = ;
1 3 2 4
^
se deduce que X0 es un punto extremal de (P L1 ) y por lo tanto X es un punto
extremal de (P L1 ):
3.2. FORMA TABULAR 47
^
c) Para obtener la forma explcita de (P L1 ) debemos calcular la inversa de B :
t
1 cofactores(B)t 1 3 1
B = =
det(B) 4 1 1
3=4 1=4
=) B 1 = :
1=4 1=4
Entonces:
3=4 1=4 1 1 1 0 1 0 3=4 1=4
B 1A = = ;
1=4 1=4 1 3 0 1 0 1 1=4 1=4
3=4 1=4 4 2
B 1b = = ;
1=4 1=4 4 2
y el sistema en forma explcita se escribe:
x1 + 43 x3 + 14 xh4 = 2;
x2 + 41 x3 41 xh4 = 2:
Eliminando de z = 2x1 6x2 + 3x3 las variables bsicas, tendremos:
3 1 h 1 1
z = 2 2 x3 x 6 2 x3 + xh4 + 3x3 ;
4 4 4 4 4
=) z= 8 + 3x3 2xh4 :
El problema en forma explcita con respecto a la base B ser:
0 1
m n z = 8 + 3x3 2xh4
B s:a : x1 + 3 x3 + 1 xh4 = 2 C
B 4 4 C
@ x2 + 14 x3 14 xh4 = 2: A
x 0
]
d) Resolvamos el problema (P L1) a partir de la solucin factible bsica X0 : La tabla
Simplex correspondiente ser:
x1 x2 x3 xh4 b
x1 1 0 3=4 1=4 2
x2 0 1 1=4 1=4 2
z 0 0 3 2 8
xh4 es la variable que entra a la base y x1 es la variable que sale. Transformando,
nos queda la siguiente tabla:
x1 x2 x3 xh4 b
xh4 4 0 3 1 8
x2 1 1 1 0 4
z 8 0 9 0 24
48 CAPTULO 3. EL MTODO SIMPLEX
]
La solucin X1 = (0; 4; 0; 8)t es ptima para (P L1) y por tanto Xopt = (0; 4; 0)t
es solucin ptima de (P L1), con valor zopt = 24. Como no hay (zj cj ) = 0
para los ndices j no bsicos, se tiene que NO hay solucin mltiple.
b) Halle una solucin factible bsica inicial donde aparezca x3 como variable bsica,
adems de las variables de holgura y articiales adecuadas.
2do mtodo:
Solucin.
b) Aadimos una variable articial xa6 en la primera ecuacin para obtener la base
inicial unitaria: 0 1
1 0 0
B=@ 0 1 0 A
0 0 1
correspondiente a los ndices bsicos IB = f6; 5; 3g, en ese orden. La SFB inicial
ser:
x1 = x2 = xh4 = 0;
x3 = 10; xh5 = 15; xa6 = 3:
que todava no es ptima. La variable que entra es x1 y la que sale es xh5 . Trans-
formando, obtenemos la tabla:
x1 [x2 ] x3 xh4 xh5 b
x2 0 1 1 0 1 5
x1 1 0 3=5 0 2=5 0
xh4 0 0 2=5 1 3=5 2
z 0 0 7=2 0 1=5 5
que es ptima. La solucin, en las variables originales, es el vector (x1 ; x2 ; x3 ) =
(0; 5; 0)t con valor z = 5.
d) Eliminando la variable x3 utilizando la ecuacin x3 = 10 + 5x1 2x2 , pero re-
spetando la restriccin x3 0, se obtiene el problema equivalente:
0 1
m n z = 22x1 9x2 + 40
B s:a : x1 + x2 3 C
B C
(P L2 ) : B
B 5x 1 + 3x 2 15 C
C
@ 5x1 + 2x2 10 A
x1 ; x2 0:
a) Escriba el problema en forma explcita con respecto a la base formada por los
ndices IB = f1; 3; 6g
2 1 3
1. Demuestre que d = ; ; ; 1; 0
5 5 5
es una direccin extremal de P.
2 1
B=
5 1
S = (x; y) 2 R2 j x y 4; 2x y+3 0, x y 3
Para qu valores del parmetros A es el punto (4; 2) una solucin ptima del problema
anterior?
3.3. EJERCICIOS PROPUESTOS 53
Ejercicio 35 Resolver el siguiente programa lineal segn los valores del parmetro real
a. 0 1
max z = ax1 + 3x2
B s:a : x1 + x 2 C
B C
B x + x 4 C
B 1 2 C
@ x1 + 2x2 1 A
x 1 ; x2 0
x1 x2 x3 x4 x5 b
1 1
x3 2
0 2
3
x4 2
0 0 1
x2 0
z 0 0
c) Resolver el programa.
0 1
max z = x1 + 9x2 + x3
B s:a : x1 + 2x2 + 3x3 9 C
1. B
@
C
3x1 + 2x2 + 2x3 15 A
x 1 ; x 2 ; x3 0
0 1
m n z = 2x1 + 5x2 x3
B s:a : 2x1 + 3x2 x3 3 C
B C
2. B
B x 1 + 2x 2 + x 3 4 C
C
@ x1 + 3x2 9 A
x 1 ; x 2 ; x3 0
0 1
m n z = x1 + x2 4x3
B s:a : x1 + x2 + 2x3 9 C
B C
3. B
B x1 + x2 x3 2 C
C
@ x1 + x2 + x3 4 A
x 1 ; x2 ; x3 0
0 1
m n z = 5x1 3x2
B s:a : x1 + 2x2 6 C
B C
4. B
B 2x1 + x2 4 C
C
@ 5x1 + 3x2 15 A
x 1 ; x2 0
0 1
max z = 3x1 + 2x2 + x3
B s:a : 2x1 3x2 + 2x3 3 C
5. B
@
C
x1 + x2 + x3 5 A
x 1 ; x2 ; x 3 0
0 1
max z = 4x1 6x2 18x3
B s:a : 2x1 + 3x3 3 C
B C
6. B
B 3x 2 + 2x3 5 C
C
@ x1 + 3x2 9 A
x 1 ; x 2 ; x3 0
0 1
max z = 2x1 + 3x2 + 5x3 + 6x4
B s:a : x1 + 2x2 + 3x3 0 C
7. B
@
C
x1 x2 + x3 3x4 3 A
x 1 ; x2 ; x3 ; x 4 0
3.3. EJERCICIOS PROPUESTOS 55
0 1
m n z = x1 2x2
B s:a : x1 + x2 2 C
8. B
@
C
x1 + x2 1 A
x 1 ; x2 0
0 1
m n z = 3x1 + 4x2
B s:a : x1 + x2 4 C
B C
9. B
B 2x 1 + 3x2 18 C
C
@ x2 3 A
x 1 ; x2 0
a) Aplique la Fase I del mtodo Simplex para hallar una solucin factible bsica
del problema.
b) Al nal de la Fase I deben quedar variables articiales en la base con valor
0. Determine si existen restricciones redundantes y elimnelas.
c) Aplique la Fase II del mtodo Simplex y resuelva el problema completamente.
Diga si existe solucin nica.
56 CAPTULO 3. EL MTODO SIMPLEX
13. La tabla nal ptima del mtodo Simplex (incluyendo los costos bsicos cB ) que
resolvi cierto problema de PL es el siguiente:
donde xh5 es una variable de holgura y xa6 , xa7 son variables articiales, las cuales
integraban la solucin bsica inicial(en ese orden).
1
a) Cul es la matriz inversa B de la base nal ptima B?
b) Escriba en forma estndar el problema inicial cuya ltima tabla es la anterior.
x1 x2 x3 xh4 xh5 b
xh4 1 0 6
xh5 1 2 0 1 1
z 1 3 0 0 0
donde xh4 , xh5 son variables de holgura positivas. La tabla que se obtiene a partir
de ella despus de una iteracin del mtodo Simplex es la siguiente:
x1 x2 x3 xh4 xh5 b
x1 2=3 2=3 1=3 0
xh5 1=3 1=3 1 3
z 1 4
3.3. EJERCICIOS PROPUESTOS 57
En ambas tablas se observa que faltan nmeros de las casillas marcadas con
asterisco( ).Halle los nmeros faltantes en cada tabla teniendo en cuenta las fr-
mulas del mtodo Simplex.
Captulo 4
Dualidad
Recordemos:
58
59
Solucin.
a) Tenemos IB = f3; 1g, adems, puesto que A 4 y A 5 son las columnas asociadas
a las variables de holgura, debemos tener que
1 1=2 0 2 0
B = =) B = :
1=6 1=3 1 3
Llamemos T a la primeras dos las de la tabla nal:
0 1=4 1 1=2 0 5=2
T =
1 1=2 0 1=6 1=3 5=2
60 CAPTULO 4. DUALIDAD
luego
1
2 0 0 1=4 1 1=2 0 5=2 0 2
2 1 0 5
= 5 ;
1 3 1 1=2 0 1=6 1=3 5=2 3 4
1 0 1 10
1 2 0
(c3 ; c1 ) B = ( 3; 2) =) (c3 ; c1 ) = ( 3; 2) = 8 6
1 3
1
z2 = 8 6 4 = 1 =) z2 c2 = 2 =) c2 = 1 + 2 = 3
1
2
luego el problema original es
0 1
m n z = 6x1 + 3x2 8x3
B s:a C
B C
B 1
x + 2x3 5 C
B 2 2 C
@ 3x1 54 x2 + x3 10 A
x 1 ; x2 ; x 3 0
b) El dual es 0 1
max w = 5y1 + 10y2
B s:a : C
B C
B 3y2 6 C
B C
B 1
y 5
y 11 C
B 2 1 4 2 2 C
@ 2y1 + y2 2 A
y 1 ; y2 0
2. Con la tabla ptima: Puesto que como dato, nos dicen que las variables
x4 y x5 son variables de holguras, basta mirar las los zj cj correspondientes
a estas variables en la tabla ptima, esto es:
y T = (z4 c4 ; z5 c5 ) = ( 3; 2) :
4.1. TEOREMA DE HOLGURAS COMPLEMENTARIAS 61
8 1
(1; 0) 2 H1 ; ; 2 H1 \ H3 ; (0; 5=3) 2 H3
5 5
De los cuales
w (1; 0) = 2;
8 1 26
w ; = ;
5 5 5
w (0; 5=3) = 5:
8 1
Es claro que el punto ptimo es y T = ;
5 5
:
c) Debido a que y 2 H1 \ H3 =)
T
c AT y = (0; ; 0; )
T
y debido a que c AT y x = 0; necesariamente debemos tener que x2 = 0; x4 =
0: Adems tambin y T (Ax b) = 0, pero y > 0 =) Ax b = 0 esto es
xT = (7=5; 0; 1=5; 0) :
c) Usando la informacin del inciso a), encuentre la solucin ptima del problema
dual (P L)0 :
Solucin.
x1 + x2 = 1;
x1 + 2x2 = 4:
2 5 t
La solucin del sistema de ecuaciones lineales es Xmax = ;
3 3
y el valor
ptimo es zmax = 43 10 3
= 143
:
y ponemos y2 = y2 ; y 3 = y3 :
0 1
max w = y1 + y2 + 4y3
B s:a : y1 + y2 + y3 2; C
(P L)0 : B C ()
@ y1 + y2 + 2y3 2; A
y1 0; y2 0; y3 2 R
0 1
m n w = y1 y2 4y3
B s:a : y1 + y2 + y3 2; C
() (P L)0 : B @
C
A
y1 y2 2y3 2;
y1 0; y2 0; y3 2 R
0 1
m n w = y1 + y2 + 4y3
B s:a: y1 y2 y3 2; C
() (P L)0 : B @
C
A
y1 + y2 + 2y3 2;
y1 0; y2 0; y3 2 R
En conclusin:
0 1
max z = 2x1 2x2
B s:a : x1 + x2 1 C
B C
(P L) : BB x 1 + x 2 1 C =)
C
@ x1 + 2x2 = 4 A
x1 0; x2 0
0 1
m n w = y1 + y2 + 4y3
B s:a : y1 y2 y3 2; C
=) (P L)0 : B
@
C
A
y1 + y2 + 2y3 2;
y1 0; y2 0; y3 2 R
Ntese que los datos del problema (P L) son:
0 1 0 1
1 1 1
2
A = @ 1 1 A; b = @ 1 A; c = ;
2
1 2 4
y entonces, si (P L2) se escribe:
0 1
0
max z = ct x; s:a : Ax @ A b; con x ;
0
=
el dual (P L2)0 se escribe con los mismos parmetros A; b; c:
0 1
0
m n w = bt y; s:a : At y c; con y @ 0 A:
2R
4.1. TEOREMA DE HOLGURAS COMPLEMENTARIAS 65
y la inversa ser:
0 1 0 1
2 1 1 2=3 0 1=3
cof (B) = @ 0 0 3 A =) B 1
= @ 1=3 0 1=3 A :
1 1 2 1=3 1 2=3
ctB B 1 A ct =
0 1t
0 1t 0 10 1 2
2 2=3 0 1=3 1 1 1 0 B 2 C
=@ 2 A @ 1=3 0 1=3 A @ 1 1 0 1 A B C
@ 0 A =
0 1=3 1 2=3 1 2 0 0
0
0 1
1 0 2=3 0
= ( 2; 2; 0) @ 0 1 1=3 0 A ( 2; 2; 0; 0)
0 0 1=3 1
2 2
= 2; 2; ; 0 (2; 2; 0; 0) = 0; 0; ; 0 ;
3 3
66 CAPTULO 4. DUALIDAD
Solucin.
b) La primera tabla del mtodo Simplex (primera fase), incluyendo variables arti-
ciales es:
x1 x2 x3 x4 xh5 xh6 xa7 b
xh5 4 2 0 1 1 0 0 8
xa7 3 1 2 0 0 1 1 3
za 3 1 2 0 0 1 0 3
y donde los zj cj = ctB Y j cj de za fueron calculados por la denicin.
Entra la variable x1 y sale la variable articial xa7 y obtenemos una tabla nal
ptima de la primera fase:
Hemos conservado la columna de la variable articial porque ella nos servir para
hallar la segunda componente de la solucin ptima del dual. Es claro que la
solucin factible bsica que tenemos X1 = (1; 0; 0; 0; 4; 0)t es ptima si 8=3
pues en tal caso todos los zj cj seran negativos o nulos. Por lo tanto, X1 es
ptima para todo 8=3: La solucin ptima del dual ser Y1 = (0; 4=3)t : El
valor ptimo = 4 coincide en ambas funciones primal y dual.
mn cT x
s:a : x 0; x 2 Rn
Condiciones de positividad: x1 ; x2 ; x3 0.
x1 x2 x3 x4 x5 x6 x7 x8 b
x1 1 0 0 1 0 0 0 0 180
x5 0 0 0 1 1 1 1 0 5
x2 0 1 0 1 0 1 0 0 45
x3 0 0 1 0 0 0 1 0 75
x8 0 0 0 4 0 4 1 1 105
z 0 0 0 0;03 0 0;12 0;12 0 41;4
a) Cul es el precio mximo que estara dispuesta a pagar la enlatadora por duraznos
adicionales de categora A?, Cunto comprara a ese precio?.
e) Cunto tendra que variar el precio de las mitades de duraznos para que cambiara
la solucin ptima?.
x1 x2 x3 x4 x5 x6 x7 b
x6 1 2 3 1 0 1 0 6
x7 1 2 2 1 1 0 1 3 !
z 2 3 5 6 0 0 0 0
"
x1 x2 x3 x4 x5 x6 x7 b
x6 5=2 1 0 1=2 3=2 1 3=2 3=2
x3 1=2 1 1 1=2 1=2 0 1=2 3=2
z 9=2 8 0 7=2 5 0 5 15=2
Y estamos en el ptimo.
Captulo 5
Post Optimizacin
x1 x2 x3 xh4 xh5 b
x1 1 1 1 1 0 6
xh5 0 3 1 1 1 10
z 0 3 1 2 0 12
a) Encuentre el problema original (P L) sabiendo que xh4 y xh5 son la base inicial
(positivas).
b) Reemplazar en (P L) c2 por e
c2 = 5:
c) Reemplazar en (P L) c1 por e
c1 = 0:
Solucin.
a) Se tiene que
1 1 0 1 0
B = =) B =
1 1 1 1
Luego llamando a la tabla nal (correspondientes a las primeras las) Tf in a
1 1 1 1 0 6 1 1 1 1 0 6
Tf in = =) Tini = BTf in =
0 3 1 1 1 10 1 2 0 0 1 4
72
73
c1 c2 = 3
c1 c3 = 1
c1 = 2
=) c1 = 2; c2 = 1; c3 = 1:
Luego el problema original es
0 1
m n z = 2x1 + x2 x3
B s:a : C
B C
B x1 + x2 + x3 6 C
B C
@ x1 + 2x2 4 A
x 1 ; x2 ; x3 0
z1 e
c1 = (0; 0)T Y1 e
c1 = 0
z2 c2 = c2 = 1
z3 c3 = c3 = 1
z4 c4 = c4 = 0
z5 c5 = 0
x1 x2 x3 xh4 xh5 b
x1 1 1 1 1 0 6
xh5 0 3 1 1 1 10
z 0 1 1 2 0 0
se observa que perdimos optimalidad. Hay que aplicar simplex, entrando la vari-
able x3 y sacando x1 de la base.
74 CAPTULO 5. POST OPTIMIZACIN
e) Cambiar el vector b por el nuevo vector b0 = (3; 2; 4), si la solucin deja de ser
factible, obtenga la nueva solucin.
6.1. Transporte
Resuelva los siguientes problemas usando el mtodo de transporte.
O1 O2 O3 O4
A1 10 10 9 11
A2 7 11 12 13
A3 8 14 4 9
b) Determine una solucin factible bsica inicial por el mtodo del costo mnimo
y aplique el algoritmo stepping-stone para hallar la solucin ptima del prob-
lema y el costo total. Considerando que el aeropuerto A4 debe ser siempre
satisfecho.
76
6.1. TRANSPORTE 77
D1 D2 D3 D4
A1 9 7 12 8
A2 15 12 12 15
A3 8 9 6 12
A4 14 12 11 12
D1 D2 D3 D4 Ai
A1 4 14 18
A2 4 4
A3 2 4 6
A4 7 5 12
Dj 6 14 15 5
a) Sabiendo que los vectores columna de la base optimal B del problema equi-
librado son los siguientes:
0 1
1 1 1 0 0 0 0
B 0 0 0 1 0 0 0 C
B C
B 0 0 0 0 1 1 0 C
B C
B=B B 1 0 0 1 0 0 0 C
C
B 0 0 0 0 1 0 0 C
B C
@ 0 1 0 0 0 0 0 A
0 0 1 0 0 1 1 7 7
78 CAPTULO 6. OPTIMIZACIN EN REDES
halle los valores de las variables bsicas asociadas a la base B y el costo total
de transporte.
b) Compruebe que la SFB correspondiente a B es ptima, calculando los zij cij
no bsicos. Solucin nica? De no serlo calcule otra solucin ptima.
(1; 1); (1; 3); (1; 4); (2; 1); (3; 2); (3; 4)
y por tanto:
9 6 4 7
? ? ?
2 4 6 3
?
8 1 8 6
? ?
Para hallar los valores de las variables bsicas utilizamos las restricciones del
modelo de transporte; ellas indican que la suma por la y por columna debe
coincidir con las existencias y las demandas. Por consiguiente, tenemos una nica
SFB posible:
C1 C2 C3 C4 Ei
9 6 4 7
O1 10 10 15 35
2 4 6 3
O2 20 20
8 1 8 6
O3 40 5 45
Dj 30 40 10 20
b) Calculemos los ui ; vj y los zij cij :
C1 C2 C3 C4 Ei ui
9 6 4 7
O1 10 4 10 15 35 0
2 4 6 3
O2 20 9 9 3 20 -7
8 1 8 6
O3 0 40 5 5 45 -1
Dj 30 40 10 20
vj 9 2 4 7
y como todos son 0, ello indica que la SFB dada es ptima. La solucin no es
6.1. TRANSPORTE 79
3 0 0 4
a) Utilice el mtodo de costo mnimo, para encontrar una SFB con WinQSB, y
compruebe el Teorema de Dualidad Dbil en este punto utilizando MATLAB.
b) Resuelva el problema, utilizando WinQSB y MATLAB, calcule la solucin
Dual.
c) Encuentre la base asociada al ptimo.
d) Haga un anlisis de sensibilidad que diga cuanto puede variar los costos, de
tal forma que el ptimo se mantenga.
Ejemplo:
La siguiente gura muestra un ejemplo de un grafo dirigido sobre V = fa; b; c; d; eg
con A = f(a; a) ; (a; b) ; (a; d) ; (b; c)g : La direccin de una arista se indica al colocar
una echa dirigida sobre ella, como se muestra aqu. Para cualquier arista, como (b; c),
decimos que la arista es incidente con los vrtices b; c; b es adyacente hacia c, mientras
que c es adyacente desde b. Adems el vrtice c es el origen, o fuente de la arista (b; c)
y el vrtice c es el trmino, o vrtice terminal. La arista (a; a) es un ejemplo de un lazo,
y el vrtice e que no tiene aristas es un vrtice aislado.
Cuando no importa la direccin de las aristas, la estructura G = (V; A), donde A
es ahora un conjunto de pares no ordenados sobre V , es un grafo no dirigido. La gura
siguiente muestra un grafo no dirigido sobre V = fa; b; c; dg.
11
7
6
3 4
9
6
11
4 4 9 4 5
1
3
5
5 2
6
11
17
1. Seleccionamos un vrtice vi+1 , tal que L (vi+1 ) sea mnimo (para todo v de
este tipo). Puede haber varios de estos vrtices, en cuyo caso podemos elegir
cualquiera de los posibles candidatos. El vrtice vi+1 es un elemento de Si
que es el ms cercano a v0 .
2. Asignamos Si [ fvi+1 g a Si+1 .
a) Incrementamos el contador i en 1.
Si i = n 1, el grafo etiquetado contiene la informacin deseada. Si i < n 1:
Regresamos al paso 2.
i = 0; n = 6; S0 = fcg ; S0 = fa; b; h; f; gg ;
Iteracin 1:
S0 = fa; b; h; g; f g ; i = 0
L (a) = m n fL (a) ; L (c) + p (c; a)g = m n f1; 0 + 1g = 1
L (b) = m n fL (b) ; L (c) + p (c; b)g = m n f1; 0 + 1g = 1
L (h) = m n fL (h) ; L (c) + p (c; h)g = m n f1; 0 + 11g = 11
L (g) = 1
L (f ) = 6
As etiquetamos el vrtice f con (6; c) y el vrtice h con (11; c). Los dems vrtices de
S0 siguen etiquetados con (1; ). En el paso 3 vemos que f es el vrtice v1 en el S0 ms
cercano a v0 por lo que asignamos a S1 el conjunto S0 [ ff g = fc; f g e incrementamos
el contador i a 1. Como i = 1 < 5(= 6 1); regresamos al paso 2.
Iteracin 2:
S1 = fc; f g ; S1 = fa; b; h; gg , i = 1
L (a) = m n fL (a) ; L (c) + p (c; a) ; L (f ) + p (f; a)g
u2S1
= m n f1; 0 + 1; 6 + 11g = 17
L (b) = m n fL (b) ; L (c) + p (c; b) ; L (f ) + p (f; b)g
u2S1
= m n f1; 0 + 1; 6 + 1g = 1
L (g) = m n fL (g) ; L (c) + p (c; g) ; L (f ) + p (f; g)g
u2S1
= m n f1; 0 + 1; 6 + 9g = 15
L (h) = m n fL (h) ; L (c) + p (c; h) ; L (f ) + p (f; h)g
= m n f11; 0 + 11; 6 + 4g = 10
Sin mayores detalles, tenemos que
S2 = S1 [ fhg = fc; f; hg
Iteracin 3:
S2 = fc; f; hg , S2 = fa; b; gg ; i = 2
L (a) = m n fL (a) ; L (c) + p (c; a) ; L (f ) + p (f; a) ; L (h) + p (h; a)g
u2S2
= m n f17; 0 + 1; 6 + 1; 10 + 11g = 17
u2S2
L (b) = m n fL (b) ; L (c) + p (c; b) ; L (f ) + p (f; b) ; L (h) + p (h; b)g
u2S2
= m n f1; 0 + 1; 6 + 1; 10 + 1g
L (g) = m n fL (g) ; L (c) + p (c; g) ; L (f ) + p (f; g) ; L (h) + p (h; b)g
= m n f15; 0 + 1; 6 + 9; 10 + 4g = 14
u2S2
6.4. ALGORITMO DEL CAMINO MS CORTO DE DIJSKTRA 85
S4 = fc; f; h; g; bg ; S4 = fbg
d (c; f ) = L (f ) = 6
d (c; g) = L (g) = 14
d (c; b) = L (b) = 22
d (c; h) = L (h) = 10
d (c; a) = L (a) = 17
Ahora bien, para determinar una camino dirigido ms corto desde c a b, partimos del
vrtice b, que tiene la etiqueta (22; a). El vrtice que lo precede es (17; f ) y despus
(6; c). As el camino es:
matlab.
11
7
6
3 4
9
6
11
4 4 9 4 5
1
3
5
5 2
6
11
17
Fig. 1
>> W=[6,5,17,4,3,7,6,11,7,11,4,9,5,11,9,4];
2. Relacionar cada arista con los pesos, cada vrtice del cual sale una arista se ingresa
en un vector y cada componente de este vector esta relacionada con los vertices
de incidencia. Por ejemplo la arista a12 y a34 se relaciona a travs de los vectores
[1; 3] y [2; 4] donde cada componente del primer vector son los nodos de los cuales
salen las aristas que inciden en 2 y 4 respectivamente. Para generar la estructura
del grafo se utiliza la funcin sparce a continuacin mostramos la rutina:
>> DG=sparse([1,1,1,2,2,2,3,3,4,4,4,4,5,6,6,6],[3,2,5,3,1,6,4,6,3,1,6,5,6,1,4,5],W);
6.4. ALGORITMO DEL CAMINO MS CORTO DE DIJSKTRA 87
DG =
(2; 1) 3
(4; 1) 11
(6;1) 11
(1; 2) 5
(1; 3) 6
(2; 3) 4
(4; 3) 7
(3; 4) 6
(6; 4) 9
(1; 5) 17
(4; 5) 9
(6; 5) 4
(2; 6) 7
(3; 6) 11
(4; 6) 4
(5; 6) 5
DistanciaMin = 14
Camino = 3 4 6 5
88 CAPTULO 6. OPTIMIZACIN EN REDES
pred = 4 1 0 3 6 4
Node1
5
3
11
Node2
6
4
11
Node6
9 17
4
Node4 4
4 11 5
7 9
6
Node3 Node5
Fig. 2
>> set(edges,LineWidth,1.5);
6.4. ALGORITMO DEL CAMINO MS CORTO DE DIJSKTRA 89
Node 1
5
3
11
Node 2
6
7
11
Node 6
9 17
4
Node 4 4
5
4 11
7 9
6
Node 3 Node 5
Parte II
Programacin No Lineal
90
Captulo 7
Condiciones de Opimalidad
P
n
Ejercicio 52 Encuentre los puntos crticos de f (x) = x4i y analice si los puntos
i=1
obtenidos cumplen condiciones para mximo o mnimo.
2 5 3
f (x) = xT x + xT + 7:
1 1 4
b) Encontrar todos los puntos que satisfacen las CN1O para (m n f (x) ; x 2 R2 ). f
tiene un mnimo? Si tiene, encuntrelos todos; en caso contrario explicar por qu
no existen.
Solucin.
91
92 CAPTULO 7. CONDICIONES DE OPIMALIDAD
1.
a) Tenemos:
@
f (x1 ; x2 ) = 4x31 + 6x21 + 4x1 2x2 ;
@x1
@
f (x1 ; x2 ) = 2x2 2x1
@x2
@
(a) : @x1
f (x1 ; x2 ) = 4x31 + 6x21 + 4x1 2x2 = 0
rf (x; y) = 0 () @
(b) : @x2
f (x1 ; x2 ) = 2x2 2x1 = 0
de (b) ; se tiene x1 = x2 , luego reemplazando en (b):
12x21 + 12x1 + 4 2
H (x1 ; x2 ) = r2 f (x1 ; x2 ) = ;
2 2
luego
4 2
H (0; 0) = H ( 1; 1) =
2 2
=) det H1 (0; 0) = 4; det H2 (0; 0) = 4
1 1 1 2
H ; =
2 2 2 2
1 1 1 1
=) det H1 ; = 1; det H2 ; = 1
2 2 2 2
1 1
Concluimos que los puntos (0; 0)y ( 1; 1) son puntos de mnimos y 2
; 2
es un punto de silla.
Ejercicio 56 Supongamos que c es el costo unitario por orden, co el costo por sobre
stock y cu el costo por quiebres de stock, luego el costo total es:
(x; D) = cx + c0 (x D)+ + cu (D x)+
donde D es la demanda aleatoria del producto cuya funcin de distribucin la denotamos
por F . As se requiere encontrar el pedido ptimo x que minimice los costos:
m n (x; D)
x 0
1 c0 c
x =F .
co + cu
94 CAPTULO 7. CONDICIONES DE OPIMALIDAD
d2 (x )
= (c0 + cu ) f (x ) > 0;
dx2
por lo que se satisfacen condiciones sucientes de segundo orden.
3. Realice una iteracin del mtodo de Newton modicado con BLE a partir del
punto xt0 = (0; 1; 0) considerando = 0.
2. Para qu valores del parmetro los vectores X ( ) son mximos y para cules
son mnimos?
6x2 + 6y 12x 6
rf (x; y) = ; r2 f (x; y) =
3 y 2 + 6x 6 6 y
7.1. SIN RESTRICCIONES 95
6x2 + 6y 0 y= x2
= ()
3 y 2 + 6x 0 x= 2
y2
r s !T
2 4
)X( )= 3
; 3
2
0 6
r2 f (0; 0) =
6 0
Luego
r
3
2
H1 = 12
r s ! s
3
2 3
4 3
8
H2 = 12 6 2 36 = 72 3 36 = 108
Una ltima nota es que si X1 ; X2 ; :::; Xn es una muestra aleatoria, entonces las Xi
son iid, as
f (x1 ; x2 ; :::; xn j ) = f (x1 j ) f (x2 j ) f (xn j ) .
a) Pruebe que max 2Rk f (x1 ; x2 ; :::; xn j ) se alcanza en el mismo punto ^ que en el
problema max 2Rk ln f (x1 ; x2 ; :::; xn j ).
Xi exp ( ) ; 8i = 1; :::; n
g ( ) := ln f (x1 ; x2 ; :::; xn j )
X
n
1 2
g( ) = ln p e
1
2 ( xi )
i=1
2
!
Xn
1 xi
2 p
= ln 2
i=1
2
1X
n 2
xi p
) g( ) = n ln 2
2 i=1
Busquemos los puntos que satisfacen las CN1O, es decir, los puntos tales que
0 1
@g (^)
0
rg ^ = @ @g@ ^ A =
() 0
@
@g ( ) X
n
xi ^ 1X
n
) = =0)^= xi
@ =^ i=1
^ n i=1
98 CAPTULO 7. CONDICIONES DE OPIMALIDAD
!1=2
@g ( ) X
n
(xi ^ )2 n 1X
n
) = =0)^= (xi ^ )2
@ =^ i=1
^3 ^ n i=1
En conclusin, los estimadores de mxima verosimilitud, son
1X
n
^= xi (media muestral)
n i=1
!1=2
1X
n
^= (xi ^ )2 (desviacin estandar muestral)
n i=1
Probemos ahora, que efectivamente son mximos para max 2R2 g ( ), para ello probemos
que (^ ; ^ ) satisface CN2O. En efecto,
1 X
n
@2g ( ) n @g ( )
= ; = (xi )
@ 2 @ @ 2
i=1
@g ( ) 1 X
n
@2g ( ) 3 X
n
n
= 2
(xi ); = 4 (xi )2 +
@ @ i=1
@ 2 i=1
2
0 P
n 1
n 1
2 (xi )
B C
r2 g ( ) = @ P
n P
n
i=1
A
1
2 (xi ) 3
4 (xi )2 + n
2
i=1 i=1
Ejercicio 62 Suponga que usted conoce que la cantidad de individuos que entraron en
un da cualquier a un supermercado se distribuye Poisson a tasa , es decir,
n
e
P (X = n) = ; n = 0; 1; 2; :::
n!
7.2. CON RESTRICCIONES DE IGUALDAD 99
A usted le encargan estimar el parmetro para hacer inferencia del futuro. Para
ello usted toma una muestra iid fx1 ; x2 ; :::; xK g.
a) Utilice el mtodo de mxima verosimilitud para encontrar ^ en funcin de la
muestra.
b) Pruebe que el estimador obtenido en a) es efectivamente el mximo de la funcin
de verosimilitud. Utilice condiciones sucientes de segundo orden.
Ejercicio 63 Otra aplicacin en estadstica es para ajustar una curva lineal o polino-
mial a un conjunto de datos extraidos de una poblacin.
Suponga que se van a observar las respuestas Yi que corresponden a los valores de
entrada xi , i = 1; :::; n; y que se van a usar para estimar y en un modelos de
regresin lineal simple, es decir se quiere ajustar una curva de la forma y = + x
al conjunto de datos (x1 ; y1 ) ; (x2 ; y2 ) ; :::; (xn ; yn ). Para determinar los valores de y
razonamos como sigue: Si ^ y ^ son los estimadores de y , entonces el estimador de
la respuesta correspondiente a la variable de entrada xi ser ^ + ^ xi . Como la verdadera
2
respuesta es Yi , el cuadrado de la diferencia es Yi ^ ^ xi , y de esta manera si ^
y ^ son los estimadores de y respectivamente, entonces la suma de los cuadrados
de las diferencias entre los valores de la respuesta estimada y de la respuesta verdadera,
llammosle SS, est dada por
X
n
2
SS = Yi ^ ^ xi :
i=1
optimizar f (x; y) = x + y
s:a : x2 + y 2 = 1
Ejercicio 67 Encuentre todos los puntos crticos del siguiente problema y compruebe
por medio de las condiciones sucientes de segundo orden, si es un mximo o un mnimo.
@g (y; z)
= z 2yz z2 = 0
@y
@g (y; z)
= y y2 2yz = 0
@z
Aqu salen 4 puntos crticos: (0; 0) ; (1; 0) ; (0; 1) y (1=3; 1=3) : El jacobiano de la funcin
objetivo es el siguiente:
2z 1 2y 2z
H (y; z) = r2 g (y; z) = ;
1 2y 2z 2y
7.2. CON RESTRICCIONES DE IGUALDAD 101
0 1
H (0; 0) = ; det H1 (0; 0) = 0; det H2 (0; 0) = 1
1 0
0 1
H (1; 0) = ; det H1 (1; 0) = 0; det H2 (1; 0) = 1
1 2
2 1
H (0; 1) = ; det H1 (0; 1) = 2; det H2 (0; 1) = 1
1 0
1 1 2 1 1 1 2 1 1 1
H ; = 3 3 ; det H1 ; = ; det H2 ; =
1 2
3 3 3 3 3 3 3 3 3 3
De aqu sale que (0; 0) ; (1; 0) ; (0; 1) son puntos de ensilladura y (1=3; 1=3) es un mximo,
puesto que H 31 ; 13 es denida negativa.
con a y b 2 R. Determinar:
a) Las condiciones que deben cumplir los parmetros a y b para que el punto (1; 1)
sea el punto critico del programa.
b) Las condiciones que deben cumplir los parmetros a y b para que en el punto
(1; 1) halla un mximo local estricto del programa.
a) Encuentre todos los puntos que satisfacen condiciones necesarias de 1er orden
para mnimo local.
b) Usando condiciones Sucientes de 2do orden, pruebe que el o los puntos obtenidos
en a) son mnimos locales.
@
L (x; ) = 2x1 + 2 1 x1 + 2 2 x2 = 0
@x1
@
L (x; ) = 2 (x2 + 2) + 1 2 =0
@x2
@
L (x; ) = x21 + x2 1 0
@ 1
@
L (x; ) = x21 x2 1 0
@ 1
1 x21 + x2 1 = 0
2 x21 x2 1 = 0
2x1 (1 + 1 + 2) = 0; =) x1 = 0
1 2 = 2 (x1 + x2 )
Esta matriz no es denida positiva ni denida negativa en todo R2 . Por lo tanto, esta-
mos obligados a calcular el subespacio tangente (en este caso, el hiperplano tangente
apliado):
^ = fd 2 Rn : rgi (x) d = 0; i = 1; :::; p; i > 0g
M
como 1 > 0; y 2 = 0, tenemos que el hiperplano tangente ampliado es:
^ = d 2 R2 : rg1 (0; 1) d = 0
M (no se considera rg2 (0; 1) d = 0, porque 2 = 0)
= f( ; 0) : 2 Rg
8 0
( ; 0) r2 L (x; ) = ( ; 0) =8 2
> 0; 6= 0:
0 0 2 0
y as el punto satisface condiciones sucientes de segundo orden para mnimo local.