Part01
Part01
Part01
UNIVERSIDAD
NACIONAL
DE COLOMBIA
Sede Palmira
Departamento de Ciencias Básicas
UNIVERSIDAD NACIONAL DE COLOMBIA
Sede Palmira
ISBN : 958-8095-09-3
Derechos reservados .
LISTA DE FIGURAS 7
INTRODUCCiÓN. 9
2. PROGRAMACiÓN LINEAL. 16
3. PROGRAMACiÓN ENTERA. 73
APÉNDICE A. 217
BIBLIOGRAFíA 230
LISTA DE FIGURAS
,,' igura 12. Cambio de flujo de unidades del segundo tablero. lOS
Figura 13. Sistema de vías para el transporte en autobuses de una ciudad. 120
Figura 14. Posibilidades de construir una red eléctrica entre siete municipios. 127
Figura 14-1. Diseño para construir una red eléctrica entre siete municipios. 128
Figura 16. Solución factible para una Red de flujo máximo propuesta. 132
Figura 20. Red sin trayectorias de aumento para el ejemplo de flujo máximo. 135
Figura 22. Hoja de trabajo en EXCEL para el ejemplo de la red de distribución. 141
Figura 23. Cuadro de diálogo de SOLVER en EXCEL. 142
Figura 24. Hoja EXCEL con la solución SOLVER para la red de distribución. 143
Figura 27. Diagrama de red para el proyecto del traslado de las oficinas. 153
Figura 29. Distribución beta para las tres estimaciones dc tiempo dc PERl'. 162
El texto y el material didáctico están divididos en diez secciones que cubren los temas
del curso de Investigación de Operaciones para la carrera de Administración de
Empresas de la Universidad Nacional de Colombia, Sede Palmira.
Cada tema se ilustra con suficientes ejemplos y al final de cada sección se proponen
ejercicios para que sean resueltos por el lector.
9
1. GENERALIDADES SOBRE INVESTIGACiÓN DE OPERACIONES.
10
INVESTIGACION DE OPERAC IONES PARA I NGEN IER IAS y ADM IN ISTRACION DE EMPRESAS
11
LU IS ALBERTO RINCON ABR IL
Las etapas esenciales en el uso de cualquiera de las anteriores técnicas son las
siguientes:
12
INVEST l t'.\C ION DE O I'ER,\C IONES I',\ R,\ I N(;[N IERI ,\ S y , \ I) ~ II N I S ' I R,\C ION DE E~ II'RE S ,\S
1.1.1 Ejemplos.
13
L U IS A L BE RTO RI NCON AB RIL
d) Por aumento en las ventas, una procesadora requiere mayor espacio de bodegas.
Una solución alternativa considera la compra de un depósito a 10 Km de la fábrica,
para almacenar los productos terminados. Este plan requiere el traslado del
producto terminado desde la planta hasta el depósito. La gerencia debe
determinar el número óptimo de camiones que serán alquilados o comprados .
Equivale esto, a estimar:
14
I NYESTIGAC ION DE OPERAC IONES PARA I NGEN I ER I /\S y ADM IN ISTR ACION DE EMPRESAS
15
2. PROGRAMACiÓN LINEAL.
16
IN\ I S-II G\C!ON Dr: Ol'rlU( ' IO NES I',\RA I NGEN I ER I,\ S y ,\D~I I N I S -II {'\C I() N DC UIPRE~AS
J7
LU IS ALI3ERTO R INCON AI3RIL
La función F(Xj ) en las variables X 1, X2, X 3,.... .... ,X n; es lineal , si puede ser
expresada como: F(Xj ) = a1X1 + a2X2 + a3X3 + ........ + anX n
y la condición de no negatividad :
X>O
J- para todo j = 1, 2, 3,.. .... ,n
En donde:
Xj : Variable de decisión asociada con cada actividad j-ésima U=1 ,2,3,... .n) .
18
INVESTI GAC IO N DE O PE RAC ION ES PA R A I NGEN IERI AS y A DMI N ISTR AC ION DE EMPR ESAS
c¡: Coeficiente de efectividad por unidad para la actividad j-ésima U=1,2,... ,n).
b¡: Cantidad limitante del recurso i-ésimo (i=1 ,2,3, .. .m).
a¡¡ : Cantidad de recurso i-ésimo por unidad de actividad j-ésima. Normalmente
reciben el nombre de coeficientes tecnológicos.
Z: Cuantifica la función objeto seleccionada.
Aunque existen varias técnicas para lograr el modelo matemático de cada problema
en particular, se recomienda la siguiente para llegar al propósito.
Es cualquier conjunto de valores positivos para las variables: X 1, X2 , X3 , ... . .... ,X n ; que
cumple cada una y todas las restricciones del modelo matemático de programación
Lineal. En caso contrario, es decir, no cumplen la condición de no negatividad o
algunas de las restricciones, se define como no factible.
19
L U IS ,\ U 3E RTO RI NCON ,\BR IL
Es el conjunto de valores para las variables: X" X2 , X3 , ....... . ,Xn ; que satisfacen el
criterio de factibilidad y optimizan la función objetiva del modelo matenlático de
programación Lineal.
20
INVEST IG¡\C ION DE OPERAC IO NES P;\RA INGEN IER I ,\S y ;\D~1INISTRACION DE H IPR ES /\S
Variables de decisión.
Restricciones.
21
L U I S 1\l.Ilr.RTO RI NCON ,\I1R II .
Producción de salchichas
Alternativa Tipo I Tipo 11
Propuesta Ton/sem Ton/sem
1 O O
2 40 60
3 70 80
4 90 90
5 80 90
6 90 80
22
INVEST Il,AC ION DE OPE RAC IONES PA RA ING EN IERI AS y ADM IN ISTR AC ION DE EMPRESAS
Variables de decisión.
Restricciones.
23
I.lIIS ,\ 1 B !. RTO RINCON ,\I3RIL
X1 + X2 2 50
De acero se mezclarán ::; 70 Ton
X1 ::; 70
De chatarra se mezclarán ::; 40 Ton
X2 ::; 40
El tiempo para la producción ::; 120 horas
2X 1 + 3X 2 ::; 120
La relación chatarra/acero ::; 7/8
7X 1 - 8X 2 2 O
24
INVEST IGAC ION DE OPERAC IONES PARA I NGEN I ER I AS y ,\DM IN IST RAC ION D E EMPRESAS
Variables de decisión.
X 1 : Ton/semana de tomate A.
X2 : Ton/semana de tomate B.
X3 : Ton/semana de tomate C.
Restricciones.
2S
LUIS ALBERTO RINCON ABR IL
26
IN\TS'II(; ,\C!ON DE OPER ,\CIONES PARA INCEN IERI AS y , \I)~I I N I STR /\CION DE E~ IP RES . \S
Variables de decisión.
X, : Número de mesas.
X2 : Número de sillas,
X 3 : Número de escritorios,
X4 : Número de libreros.
Restricciones.
27
LU IS ALBERTO RINCON ABR I L
28
INVEST IGAC ION DE OPERAC I ONES P;\I<A INGEN IERI AS y AD~ II N I STRAC I ON DE EM PRESAS
Ejemplo de solución factible gráfica. Aplicando los conceptos de la sección 2.4 .1,
la región solución para el siguiente conjunto de restricciones se tiene en la figura 1.
2X,+3X 2 ::; 12 R1
3X, + 2X 2 ::; 12 R2
X, + X2 ;c: 2 R3
XI ;c: O para j = 1,2
El gráfico se obtiene trazando las rectas correspondientes a las restricciones R1, R2,
R3 Y ubicando el semiplano solución para cada una de ellas.
29
L U I S I\ L J3ERTO I{ INCON ABR Il.
Uno de los teoremas que se presentará para el método Simplex , considera que la
solución óptima al Programa Lineal coincide con uno de los puntos extremos del
conjunto de soluciones factibles ; por tanto en los programas de dos variables de
decisión , será necesariamente uno de los vértices de la región solución. Así pues,
para encontrar la solución óptima basta reemplazar cada uno de los vértices en la
función objetiva y observar cuál genera el mayor valor en el caso de maximizar o cuál
el menor valor en el caso de minimizar.
30
INVEST IG /\C ION DE OPER /\ClONr.S I',\R ,\ INGEN IER I AS y r\D~ II N I STR ,\C I ON DE H I PRES¡\S
El cual tiene solución para Xl = 2.4 Y X2 = 2.4. Con esto los vértices de la región
están determinados y son: A = (0 ,4) , B = (0 ,2) , C = (2,0) , 0 = (4,0) Y E = (2.4 ,2 .4).
Se tendrá entonces que: l A = 12 , l B = 6 , l c = 8 , l o = 16 Y lE = 16.8. Así que
Max(Z) =16.8, que se obtiene para X1 =2.4 Y X2 =2.4.
31
LU IS 1\ Ll3E RTO RINCO ABRIL
0 11 {l1 2 {I I ~ .. .. . . . . .. 0 1/1
{/ 2 1 {I n {I 23 .......... O 2 /1
A {/ 3 1 0 32 A .1.)
" ....... . .. (l ,
."1
. . . .... . . . . . .. .
XI b,
X2 bo
X = X -', b = b, e - (c
- I Co C,., cJ
X II bJl
= ex
Opt(Z)
sujeto a: AX :s b
x>o
Consiste en escribir el modelo con la Función objetiva para maximizar y todas las
restricciones de la forma S, esto es:
Max(Z) = ex
Sujeto a: AX < b
32
INVESTIGACION DE OPERAC IONES PARA INGEN I ER IAS y ADM IN ISTRAC ION DE EMPRESAS
x > O
Max(Z) = 4X 1 + X2
R1 : X1 + X2 <4
R2 : X1 + 2X 2 ~ 6
Xi 2: O para j = 1,2
Z - 4X 1 - X2 =O
R1: X 1 + X2 + X 3 = 4
R2: X1 + 2X 2 + X4 =6
33
1.1 li S ,\'-BERTO RINCON ABR IL
2.5.4.1 Solución básica factible. Una solución básica factible al programa lineal de
la forma básica con m restricc iones y n+m variables (n de decisión y m de holgura)
es aquella con no más de m componentes positivas que satisface cada una de las
restricciones .
34
IN\'EST IU ,\C ION DE O PER,\C IONES P¡\ R¡\ I N(iEN IERI ,\S y ¡\mdI N ISTR ¡\C ION DE EMPR ESAS
Considerando el programa lineal en su forma básica, éste puede ser escrito como:
Z- ex = o
(A I)X = b
X>O
3S
LU IS ALBE RTO RI NCON ABR IL
36
IN\· I S II (ó .\C ION DE O prR .\C ION F S P,\ R,\ I N(ó l .N II' RI.\S y . \ I) ~ II NIST R . \C I ()N DE E~ II' R I S.\S
Se supone que se inicia con una solución básica factible dada por BX s = b , que
corresponde a un punto extremo de la región de factibilidad , el cual si no es óptimo,
deb e moverse a un punto extremo vecino con objeto de mejorar la función objetiva.
Esto se puede hacer cambiando un vector de la base B por un vector de N. Hay
varias forma s de hacerlo, pero Dantzig 2 fijó la teoría que justifica el cambio que
garantiza el mejor incremento en la función objetiva. Va se había obseNado que
cua lquier columna al de S (no en B) puede escribirse aj = LVkjak con k = 1, .... ,m.
Supón gase que el vector que sale de B es a r y que la componente Vr¡ de VI es
di ferente de cero . Entonces se tiene aj =LVkjak + Vrjar con k = 1 ,.... ,m y diferente de r.
(/ (/ )/,
Esto es, (/ , = )' '- "L....)'¡ , '¡f Y" * O. Por otro lado la solución básica factible puede
1/ ,.,
escribirse LakXbk = b , equivale con arX Sr + LakXbk = b para k=1 ,.... ,m y diferente de
X II Y, {[X II
r. Reemplazando a, en ella se obtiene I (X /Il -
y
"
I}
)u , + '
y
JI
'= iJ, que es una
nueva solución básica. Examinando ahora qué condiciones debe cumplir esta
solución pa ra que sea factible. Ella puede reescribirse como Ix ¡¡¡ + x,(/, = h . En
X Y X
donde se ha definido X, = X ¡¡¡ - 11", X, = 11, las cuales se requieren
Y" Y" '
positivas. Dado que X Br ~ O, se requiere que Vr¡ > O para que X r ~ O. De otra parte si
todas las Y kl S O para k=1 ,2,.... ,m y k # r, se garantiza que X k ~ O. Sin embargo, qué
, OANTZI G G. B .. Compu tationa l algo rit hm of the revised si mp lex meth od. RANO report RM-1 266.
Cali fornia 1953.
37
L U IS ALBERTO RINCON ,\13R II.
x y
ocurri ría si alguna YkJ < O. Analizando se tiene que X ¡¡; - ---.!.' A, ~ O. Y" > O, que
Y"
Xli,
puedo escribir como X ¡¡; ~ O. Y" > O, para garantizar que esto se cumpla se
YA, Y"
requiere que X ¡¡¡ ~ X li, O. Y" > O. De esta condición sólo se puede estar seguro
YA, Y"
Esta regla garantiza que el cambio de un vector de B por otro de N genera una
nueva solución básica factible . Como resultado de este cambio se obtiene una
nueva base B que difiere de la base anterior B en un solo vector. Como cada base se
asocia con un extremo de la región de factibilidad , con el cambio de base el proceso
se ha movido a otro punto extremo XB tal que:
X B = S·1b, Z = CBX B
El siguiente análisis muestra la regla que permite hacer la mejor selección del vector
aJ en N que se debe introducir en B. El valor de Z asociado con el punto extremo XB
de la base es Z = CBXB y el asociado con la nueva base X'B será Z' = C'BX'B. Ya se ha
mencionado que la única diferencia entre las base B y B' es que se ha sacado el
vector ar y se ha reemplazado por el vector aj de N. Entonces la única diferencia entre
CB y C'B se encuentra en la r-ésima componente, es decir:
38
INVESTI G/\C ION DE O l'lcR/\ C IONES I'¡\R ¡\ IN(, EN I E RI ¡\ S Y ¡\D ~ II N I S TR AC I ON D E H IPRES/\ S
X Y e X
·,
es t e t ermlno, se pue d e escn'b'Ir Z '= "L e¡¡¡ X /I! -
"LC m ~Y~-
/JO' "
y- - . con y,,
+ - ,fj, =F O :::::;.
1/ /'1
X
- Z - - n,-
Z '-
y
I ('II!',y + e ,/1,-
X
Y
:::::;. Z '= Z _ X 'J, ~, + e, X n, :::::;. T = Z - ( :: - e ) _
, .1
X
Y
B,
39
L U IS ¡\Lll LRTO RI NCON I\BRIL
El teorema 4 establece las condiciones que debe cumplir una solución para que sea
la solución óptima de un programa lineal.
Zj - Cj ~ 0, para todo j.
Definen el cambio de vector a realizar para encontrar la nueva base cuando alguno
de los requisitos del criterio de optimalidad no se cumple.
2.5.7.1 Criterio Primal. Permite mejorar el valor de la función objetiva. Debe ser
aplicado cuando se tiene una solución básica factible , esto es, Xk ~O, para todo k y
existen algunos ZI - cl < O. Consiste en lo siguiente:
1. Entra en la base aquella variable Xl (vector al) que tiene el ZI - cI más negativo,
(columna de trabajo).
2. Sale de la base aquel a, que genera el X /1, = Mello/"( X ¡¡¡ CO II Y¡¡ > O) , (f'lI a de
Y,¡ ¡ Y¡,
trabajo).
2.5.7.2 Criterio Dual. Permite convertir a factible una solución básica no factible. Por
lo tanto debe ser aplicado cada que la solución básica obtenida sea no factible .
Consiste en lo siguiente:
1. Sale de la base aquel a, para el que se tenga el XB, más negativo. (fila de trabajo) .
40
INVESTI(; /\UON DE OPERAC IONES PARA I NGEN I ER I AS y , \I)~ II N I STR I\CION DE E~lPRESAS
Construir el Tablero
Simplex
Es la solución
óptima? Definir el Elemento Construir el nuevo
Pivot Tablero Simplex
Si
41
LU I S 1\I . JlERTO R INCON MlR IL
z X, X2 X3 Xn X n+, X n+2 X n +3 X n +m
z(c¡ 1 -C l -C2 -C3 -C n O O O O O
a n +, O all a1 2 a1 3 al n 1 O O O b1
a n +2 O a21 a22 a23 a2n O 1 O O b2
a n +3 O a31 a32 a33 a3n O O 1 O b3
En este tablero la base B esta formada por {a n+l , a n+2 , a n+3 ,....... , a n+m }, mientras
que N la forman {al, a2 , a3 ,.... ... , a n }. Equivale a decir que {X n+1 , X n+2 , X n+3 ,... .. .. ,
Xn+m } son las variables básicas y que {Xl , X2 , X 3 ,.... ... , X n } son variables no
básicas . Con esto, cada tablero muestra una solución básica al asignar el valor O a
cada variable no básica , para que las variables básicas tom en como valor el
correspondiente término independi ente . As í pues, en este primer tablero se tiene que :
42
INY L STIG,\C ION DE OPERAC IONES PARA INGEN IER I AS y ,\D~ II N I STRAC I ON DE EMPRESAS
Cada uno de los criterios primal y dual definieron una fila de trabajo (vector de salida)
y una columna de trabajo (vector de entrada) . Esta posición de fila y columna en el
tablero determina el elemento pivote. De tal manera que en el siguiente tablero, este
debe pasar a ser 1 y el resto de la columna conformada por ceros.
Debe ser obtenido sobre el elemento pivote previamente definido usando el método
de eliminación de Gauss que se presenta en el apéndice A .
Max(Z) = 6X 1 + 4X 2
2X 1 + 3X 2 S 12
2X 1 + X 2 S 8
X1 ,2 ~ O
Z - 6X 1 - =O
4X 2
2X 1 + 3X 2 + X3 = 12
43
LUIS I\LBERTO RINCON I\BRII _
2X 1 + X2 + X4 = 8
XI ~ O para j = 1,2,3,4
Z X1 X2 X3 X4 b Tablero
Z¡-c¡ 1 -6 -4 O O O
A3 O 2 3 1 O 12 Inicial
A4 O 2 1 O 1 8
Z--c· 1 O -1 O 3 24
A3 O O 2 1 -1 4 Segundo
A1 O 1 Y2 O Y2 4
Z·-c· 1 O O Y2 5/2 26
A2 O O 1 Y2 Y2 2 Final
A1 O 1 O -1,4 3,4 3
Tablero inicial: Solución básica factible pero no óptima, pues Z1 - C1 = -6, Z2 - C2 = -4,
con variables básicas X3 = 12, X4 = 8 Y variables no básicas: X1 = X2 = o. Se debe
entonces generar un nuevo tablero (nueva base). De acuerdo con el criterio primal , la
variable X1 (vector a1) entra a la base y el vector a4 sale de la base; el elemento
pivote aparece señalado.
Segundo tablero: La nueva base es {a3 ad. Presenta una solución básica factible no
óptima, pues Z2 - C2 = -1 con las variables básicas X3 = 4, X 1 = 4 Y variables no
básicas X2 = X4 = o. Se debe entonces generar un nuevo tablero (nueva base). De
acuerdo con el criterio primal , la variable X2 (vector a2) entra a la base y el vector a3
sale de la base; el elemento pivote aparece señalado.
Tablero final : La nueva base es {a2 ad. Presenta la solución óptima, pues cumple el
criterio de optimalidad, esto es ZI - cI ~ O, para todo j y X Bk ~ O para todo k. Las
variables básicas X2 = 2, X 1 = 3 Y las variables no básicas X3 = X4 = O. Además este
tablero presenta que Max(Z) = 26.
44
l NVrST1(ói\C1()N DE OPER/\C10NES P,\Ri\ l NGEN 1ER1 ¡\S y ¡\D~ II N 1 STR¡\C 10 N DE E~IPRESi\S
Max(Z) = 2X 1 + 5X 2 + 5X 3
X 1 + X2 + X 3 S 60
2X 1 - X 2 S 18
X2 - X3 S 6
Z - 2X 1 - 5X 2 - 5X 3 =O
X 1 + X 2 + X 3 + X 4 = 60
2X 1 - X2 + X5 = 18
X2 - X3 + X6 = 6
Xi?: O para j = 1,2 ,3,4 ,5,6
Z X1 X2 X3 X4 Xs X6 b Tablero
Z -c· 1 -2 -5 -5 O O O O
a4 O 1 1 1 1 O O 60 Inicial
as O 2 -1 O O 1 O 18
a6 O O 1 -1 O O 1 6
Z·-c · 1 -2 O -10 O O 5 30
a4 O 1 O 2 1 O -1 54 Segundo
as O 2 O -1 O 1 1 24
a2 O O 1 -1 O O 1 6
Z¡-c¡ 1 3 O O 5 O 5 300
a3 O % O 1 % O -% 27 Final
as O 5/2 O O % 1 % 51
a2 O Y2 1 O % O % 33
45
LU IS ALBERTO RINCON ABR IL
Tablero inicial: Base {a4 as a6}. Solución básica factible pero no óptima, en donde
las variables básicas X4 = 60 , Xs = 18, X6 = 6 Y las variables no básicas X 1 = X2 = X3 =
O. Se debe generar un nuevo tablero (nueva base). De acuerdo con el criterio primal ,
la variables X2 (vector a2) Y X3 (vector a3) cumplen la condición para entrar a la base;
este empate se rompe arbitrariamente entrando a una de ellas X2, entonces el vector
a6 sale de la base; el elemento pivote aparece señalado.
Tablero final: La base final es {a3 as a2}. Es la solución óptima, pues cumple el
criterio de optimalidad , esto es z¡ - c¡ 2: O, para todo j y X Bk 2: O para todo k. Las
variables básicas X3 = 27, Xs = 51 , X2 = 33 Y las variables no básicas: X 1 = X4 = X6 =
o. Además Max(Z) = 300.
Min(C) = 2X 1 + 2X2
X 1 + X2 S 12
X 1 + 2X 2 2: 10
3X 1 + 2X 2 2: 12
X1.2 2: O
46
INVESTIGACION DE OPERAC IONES PARA INGEN IER IA S y ¡\D~ II N I STRAC I ON DE EMPRESAS
Z X1 X2 X3 X4 Xs b Tablero
z¡-c¡ 1 2 2 O O O O
a3 O 1 1 1 O O 12 Inicial
a4 O -1 -2 O 1 O -10
as O -3 -2 O O 1 -12
z¡-c¡ 1 O 2/3 O O 2/3 -8
a3 O O 1/3 1 O 1/3 8 Segundo
a4 O O -4/3 O 1 -1/3 -6
a1 O 1 2/3 O O -1 /3 4
z·-c· 1 O O O % % -11
a3 O O O 1 1,4 1,4 13/2 Final
a2 O O 1 O -% 1,4 9/2
a1 O 1 O O % -% 1
Tablero inicial: Base inicial {a3 a4 a5}. Solución básica no factible , pues se tiene que
las variables básicas: X 3 = 12, X4 = -10 , X5 = -12 Y variables no básicas: X, = X2 = O.
Debe entonces generarse un nuevo tablero (nueva base). De acuerdo con el criterio
dual, sale de la base a5 Y entra la variable X 1 (vector a1) ; el elemento pivote aparece
señalado.
Segundo tablero: La nueva base es {a3 a4 a,}. Solución básica no factible, pero
menos "infactible" que la anterior, esto es lo que logra el criterio dual. Con variables
básicas X3 = 8, X4 = -6, X, = 4 Y variables no básicas: X2 = X5 = o. Debe generarse un
nuevo tablero (nueva base). De acuerdo con el criterio dual, sale de la base a4 y entra
la variable X2 (vector a2); el elemento pivote aparece señalado.
Tablero final: La base final es {a3 a2 a,}. Solución óptima , pues cumple el criterio de
optimalidad , esto es z¡ - c¡ ~ O, para todo j y X Bk ~ O para todo k. Las variables
47
LlI lS ,\LBrRT O RI NCON ,\13 1<1 1.
Entre los vectores y matrices de los tableros Simplex hay varias re laciones, (algunas
de las cuales se usaron como elementos matemáticos en la demostración del
método) , que se presentarán en esta sección y que resultan ser de gran importancia
en análisis de post-optimización , tema que se trata en el apartado 2.7.
Vector columna y¡: Conformado por los elementos akJ para la variable XJ en el
tablero inicial.
Vector columna Y¡: Conformado por los elementos akJ para la variable XJ en el
tablero final.
Vector fila CN: Conformado con los zJ- cJdel tablero final para las variables de la
base inicial.
• Matriz básica B: Estructurada con los vectores YJ del tablero inicial para la base
final.
• Matriz inversa de la básica B-': Estructurada con los vectores Y¡ del tablero final
para la base inicial.
48
IN\TS11(;/lCION DE OP I i{ ,\C IONES P/lR ,\ I N(;EN I ER I /lS y ,\D~ II N I STR/lCION DE H IPR ESAS
Teorema: Un Programa Lineal carece de solución cuando existiendo algún XB, < O,
todo los correspondientes a" son positivos.
49
LUIS ALBERT O RINCON ABR IL
Min(C) = 3X l + 2X2
Xl + X2 S 3
2X l + 3X 2 ~ 12
Xl ,2 ~ O
El lector podrá fácilmente verificar que este programa carece de solución , trazando
para ello las restricciones en el plano cartesiano y observando que no existe ningún
punto que verifique las restricciones y la condición de no negatividad.
50
INVESTIG ACION DE OPERAC IONES PARA INGEN I ERI AS y ADM IN ISTRAC ION DE EMPRESAS
Max(Z) = 4X 1 + 2X 2
X 1 + X2 ~ 8
2X 1 + X2 ~ 12
X1.2 ~ O
Como este Programa Lineal tiene la forma canónica se construye la forma básica:
Z - 4X 1 - 2X 2 = O
X 1 + X2 + X3 = 8
2X 1 + X2 + X4 = 12
Xj ~ O para j = 1,2,3,4
Z X1 X2 X3 X4 b Tablero
z¡-c¡ 1 -4 -2 O O O
A3 O 1 1 1 O 8 Inicial
A4 O 2 1 O 1 12
z·-c · 1 O O O 2 24
A3 O O 112 1 -% 2 Final 1
A1 O 1 % O Y2 4
z¡-c¡ 1 O O O 2 24
A2 O O 1 2 -1 4 Final 2
A1 O 1 O -1 1 4
51
L U I S A L BERTO RINCON I\BR IL
Para un Máx(Z) = 24 , los tableros final 1 y final 2 presentan soluciones óptimas que
se pueden escribir como:
4
il
.\~ - O 2
d cl .fIll a l I y
IX'1 IX' 1
.1 = d el .fil/u l 2 .
.1 1 2 .\ .1
x" O .\"
Todas las soluciones óptimas pueden ser expresadas como una combinación
convexa de estas dos soluciones encontradas.
Max(Z) = 2X 1 + X2
X1 + X2 ~ 4
X 1 - X2 S 2
X1,2 ~ O
52
INVEST IG¡\C ION DE OPERAC IONES PARA I NGEN I ER I AS y AD~ II N I STR f\C I ON DE EMPRESAS
53
L U I S A L BERTO RI CON AB RIL
Esta nueva solución puede resultar factible o no factible ; si es el primer caso será
la nueva solución óptima . Para el segundo caso , ésta será reemplazada en el
tablero final y se procederá a reestablecer la factibilidad usando el criterio dual.
S4
INVr:ST I(ó ,\C ION DE OPERM.' IONES P;\R ,\ INlóE I ER I,\S y AD~ II N I S TR /\C I ON DE H I PRES ,\S
Solucion.
Variables De Decisión.
55
LU I S ,\L13ERTO RINCON ABRIl.
Algoritmo Simplex
z X, X2 X3 X4 Xs X6 Xl b Tablero
z·-c¡ 1 -4 -6 -8 O O O O O
a4 O -1 -1 -1 1 O O O -120 Inicial
as O -1 O O O 1 O O -30
a6 O O -1 1 O O 1 O O
al O 0.6 0.5 OA O O O 1 80
z·-c· 1 4 2 O -8 O O O 960
a3 O 1 1 1 -1 O O O 120 Segundo
as O -1 O O O 1 O O -30
a6 O -1 -2 O 1 O 1 O -120 Criterio
al O 0.5 0.1 O OA O O 1 32 Dual
z·-c· 1 3 O O -7 O 1 O 840
a3 O Y2 O 1 -Y2 O % O 60 Tercero
as O -1 O O O 1 O O -30
a2 O 0.05 1 O -% O -Y2 O 60 Criterio
a7 O 0.15 O O OA5 O 0.05 1 26 Dual
z·-c· 1 O O O -7 3 1 O 750
a3 O O O 1 -% Y2 Y2 O 45 Cuarto
a, O 1 O O O -1 O O 30
a2 O O 1 O -% % -Y2 O 45 Criterio
a7 O O O O 0.45 0.15 0.05 1 21 .5 Primal
z¡-c¡ 1 O O O O 16/3 16/9 140/9 9760/9
a3 O O O 1 O 2/3 5/9 10/9 620/9 Final
a, O 1 O O O -1 O O 30
a2 O O 1 O O 2/3 -4/9 10/9 620/9
a4 O O O O 1 1/3 1/9 20/9 430/9
S6
INVEST IGAC ION DE OPERAC IONES PARA I NGEN I ER I AS y ADM IN ISTRAC ION DE EMPRESAS
Solución Óptima.
Xl = 30 Ton de concentrado 1.
X2 = 620/9 Ton de concentrado 2.
X3 620/9 Ton de concentrado 3.
X4 = 430/9 Holgura de la restricción 1
Análisis de Post-Optimización.
C u = (8 4 6 O)
a) Cuál es la nueva solución si el total de la producción debe ser al menos 150 Ton ,
no dispone sino de 120 Ton de materia prima A y las demás limitaciones no cambian.
57
LU I S ,\ L BERTO RINCON A BR IL
- 150
-30
En este caso el nuevo vector /J' = y calculando X 'BO = S-1 b ' se obtiene
O
120
.1'¡()
30
X 'lJo= 100 y Z' = CBX'BO = 5120/3. Qu e resulta ser factible , por lo tanto es solución
X1 = 30 Ton de concentrado 1.
X2 = 340/3 Ton de concentrado 2.
X3 = 340/3 Ton de concentrado 3.
X4 = 260/3 Holgura de la restricción 1
El nuevo vector /J '= [=¡~l => X 'BO = S-1 b ' =[~~;: 1 y z' = =
CBX'BO 880/9. Que resulta
40 -~
.'
58
INVr.STI (; ,\C ION DE O I'EI{ ¡\ C IO N ES I',\R J\ I NGEN IER I J\ S y J\ D~ II N I S TR , \ C I ()N DE E~ I PR ES r\S
Z X1 X2 X3 X4 Xs X6 X7 b Tablero
z-c 1 O O O O 16/3 16/9 140/9 880/9
a3 O O O 1 O 2/3 5/9 10/9 -10/3 Nuevo
a1 O 1 O O O -1 O O 30 Tablero
a2 O O 1 O O 2/3 -4/9 10/9 140/3 Final
a4 O O O O 1 1/3 1/9 20/9 -80/3
Para este caso c' 1 =35 - (1 + 0.6*20 + 0.4*30) =10 Y (Z1 - C1)' =- 6
Z X1 X2 X3 X4 Xs X6 X7 b Tablero
z·-c· 1 O O O O 16/3 16/9 140/9 9760/9
a3 O O O 1 O 2/3 5/9 10/9 620/9 Recuperar
O 1 O O O -1 O O 30 La
a2 O O 1 O O 2/3 -4/9 10/9 620/9 Base
a4 O O O O 1 1/3 1/9 20/9 430/9
z¡-c¡ 1 O O O O -2/3 16/9 140/9 11 380/9
a3 O O O 1 O 2/3 5/9 10/9 620/9 Criterio
a1 O 1 O O O -1 O O 30 Primal
a2 O O 1 O O 2/3 -4/9 10/9 620/9
a4 O O O O 1 1/3 1/9 20/9 430/9
z¡-c¡ 1 O O 1 O O 7/3 50/3 4000/3
as O O O 3/2 O 1 5/6 5/3 3 10/3 Final
a1 O 1 O 3/2 O O 5/6 5/3 400/3
a2 O O 1 -1 O O -1 O O
a4 O O O -Y2 1 O -1 /6 5/3 40/3
S9
LU IS ALBERTO RINCON ABRIL
X3 = X6 = X7 = O. Variables no básicas.
En 1950 aparecen los primeros intentos para dar solución a programas lineales
mediante el uso de un computador, pero las primeras soluciones acertadas a un
programa lineal en un computador de alta velocidad , solo se obtuvieron en Enero de
1952 con el uso de la máquina SEAC , del National Bureau of Standard . Desde
entonces, el Algoritmo Simplex, ha sido codificado en varios lenguajes de
programación . El profesional que requiere del uso de la programación lineal , hoy en
día, en verdad no necesita más que del diseño del modelo matemático, pues las
soluciones están al alcance hasta de los más pequeños computadores y la eficiencia
de las respuestas de los programas de computador sólo dependen de un correcto
diseño del modelo matemático al problema . La aplicación de mayor uso actualmente
resulta ser la herramienta SOLVER de EXCEL , la cual se presentará en esta
sección. Hoy en día, muchos libros sobre Investigación de Operaciones vienen
acompañados de Software de aplicación para la solución de varios problemas . Por
ejemplo el texto de Hillier y Lieberman 3 dispone el OR Courseware con
procedimientos didácticos y aplicaciones para la solución de varios problemas.
' HILLlER Frederick y LlEBERMAN Gerald, Introducción a la Investigació n de Ope raciones. Editoria l
Mc Graw Hi11. 1997. México.
60
INVI STI(; ,\C IO N DI' ()I'I RM ' IONI S P,\R ,\ IN(;I.N II I{L\S y \I )~ IINISTR .\C I()N DE E~IPR r· S . \S
1 Este programa puede usarse gratuitamente y el lector podrá hace r una copia de su instalador por
Internet en la direccion
61
L U IS A LllE ln O R INCON .\I3R 11.
Los modelos de Programación Lineal encajan dentro del grupo tres que resue lve la
herramienta Solver, por lo tanto siempre podrá ser usada para este propósi to. Se
ilustrará su uso con un ejemplo para el cual ya se había construído el modelo
matemático en la sección 2.3.4.
Variables de decisión.
62
IN\' I.S1 I(;'\UON D I. O I'I.R ,\ClON I S I'/\ R /\ I N(; EN I ER I,\S y '\ D ~ II N I S T R'\C I ON D E H IP RES r\S
Modelo Matemático
Xl + X2 ~ 50
2X l + 3X 2 ::; 120
7X l - 8X 2 ~ O
Xl ::; 70
X2 ::; 40
Xl , X2 ~ O
Entre las filas uno y 10 se considera una matriz para los parámetros del modelo
matemático del Programa Lineal. Por ejemplo las celdas B5: E5 pueden leerse
como 1 Xl + 1X2 ~ 50 , que resulta ser la primera restricción , mientras que las
celdas B 1O:E 10 pueden leerse como 600X l + 300X 2 , que corresponde con la
función objetiva.
63
I IIl S \1 Ilf Rl() I<I N("() N M1RIL
la función Objetiva . Las celdas B 16:C 16 se han previsto para que Solver escriba
los valores de Xl, X2 , una vez los calcule.
.-
A B I e I oJ E
---- F
MODELO MATEMATICO
-
5 Aleació n 1 1 > 50
6 Capacid ad de Trabajo 2 3 < 120
7 Relac ió n 7 -8 > O
!----- !
8 ,Dispon ibilidad Acero 1 < 70
9 Di spon ibilidad Chatarra 1 < 40
._- J
10 F.O . Min(Costo) 600 300 !
11
12 --
SOLUCION POR SOLVER
13
14 Acero Chatarra Limitación
--
~ X1 X2
I
16 Valores obtenidos .1- I
17 Valores mínimos O O
I
I
---
I
18
19 Aleació n =B5'B$ 16 =C5'C$16 > 50 =SUMA(B19C19)
I 20 -lCapaci dad de Trabajo =B6'B$16 =C6'C$16 < 120 =SUMA(B20 :C20)
121 .Relació n =B7'B$ 16 =C7'C$ 16 > O =S UMA(B21C21 )
122 Dispon ibilidad Ace ro =B8'B$1 6 =C8'C$16 < 70 =SUMA(B22:C22)
23 Dispon ibilidad Chatarra =B9'B$16 =C9'C$ 16 < 40 =S UMA(B23C23)
24 F.O . Min(Costo) =B10' B$16 =C 10'C$ 16 =S UMA(B24C24 )
25
Preparada esta hoja de cálculo se procede a usar la opción SOLVER del menú de
HERRAMIENTAS de EXCEL y se responde el cuadro de diálogo de la siguiente
forma:
64
INVES-I I(, ,\C ION DE OPER /\C IONES P,\Ri\ I NCiEN IER I i\S y i\D~ II N I STR , \CION DE [~IPRES , \S
D E F
Parámetro, de Solver 11
Cerrar
c- t].é xirno
Cam~ando las celdas
Qpoones" ,
Suietas a las si9uientes restricciones:
I IUld .-, ./
I• I
3. Cambiando las Celdas. Se selecciona las celdas donde Solver escribirá los
valores encontrados para las variables .
5. Resolver. Una vez se responde el cuadro de diálogo, se pul sa este botón para
que EXCEL presente los resultados sobre la hoja de trabajo , tal como lo
muestra la siguiente tabla:
6S
l .l ' I S !\ I. HI-.RTO RI NCON ,\131' 11.
A B e J D E F
MODELO MATEMATICO
2
3 Acero Chatarra I Limitación .~
4 X1 X2 ~ ..
5 Aleado n 1 1 ? 50 1
6 Capac idad de Trabaj o 2 3 < 120 -
7 Relacio n 7 -8 > O
8 Dispon ibilidad Acero 1 < 70
9 Dispon ibilidad Chatarra 1 < 40 -
10 F.O. M in (Cos to) 600 300
- l.
.~-
11
I 12 I SO LUCIO N POR SOLVER I
-
13
14 Acero - Chatarra Limitación
15 X1 X2
16 Valore s obtenidos 30 20
17 Valore s mínimos O O
18
19 Al eacio n 30 -
20 > 50 50
20 Capac idad de Trabajo 60 60 < 120 120
21 Relacio n 210 -160 > O 50
- ~
22 Dispon ib il idad Acero 30 O < 70
23 Dispon ibilidad Chatarra O 20 < 40 20
24 F.O. M in (Cos to) 18000 6000
-
24000
1. La fila 16 muestra los valores obtenidos para las variables de decisión; esto es,
30 Ton de Acero y 20 Ton de Chatarra .
66
INVEST I(; ,\CION IlF OPFR ,\C ION lcS !',\RA I N (~E N IERI , \S y , \D~IINISTR;\CI()N DI' r~IPRF S.'\S
1. CONCEPTOS TEÓRICOS.
1.1. Presente una aplicación generalizada de Programación Lin eal para algún caso,
por ejemplo: Diseño industrial de raciones , Mezcla industrial de productos, Política
de préstamos bancarios , Uso y urbanización de la tierra , Asignación de personal,
Asignación de recursos , etc.
1.2. Cómo se construye el programa dual a partir de un programa primal y que
relaciones existen entre éstos?
1.3. Explique la forma en que puede aplicar el análisis de post-optimización con la
herramienta Solver de Excel.
2.2. Una compañía de alquiler de camiones dispone dos tipos de vehículos, el camión A
que tiene 40 pies 3 de espacio refrigerado y 80 pies 3 de espacio no refrigerado; el
camión B tiene 60 pies 3 de cada tipo de espacio. Una procesadora de alimentos
debe transportar 1800 pies 3 de producto refrigerado y 2400 pies 3 de producto no
refrigerado El camión A lo alqui lan a $US 9 el Km y el camión B a $US 12 el Km .
2.2.2. Escriba el modelo matemático, considerando el alqui ler de Camión C con 100
pies 3 refrigerados a $US 8 el Km?
67
L U I S A LB ERTO RI NCON A BRIL
2.3. Una industria de Aceite y Torta de Soya tiene capacidad para procesar hasta 60
Ton/día de Soya , las que adquiere a $55000 la Ton. Debe suministrar por lo
menos 16 Ton/día de Aceite. La industria usa dos procesos con los siguientes
rendimientos por Ton de Soya:
Para completar los pedidos , esta industria puede comprar Aceite a granel con
otros productores a $205000 la Ton , que le cuesta envasarla $8000 la Ton. La
industria vende Ton. de Aceite a $210000. y Ton. de Torta de Soya a $60000.
Presente un modelo matemático que permita planificar la producción en
condiciones óptimas.
2.4. Una imprenta dispone de 1800 tiras de cartulina de 13 pulg de largo. Debe atender
un pedido que le exige cortes de tal manera que disponga al menos 1000 tiras de
7 pulg , 2000 tiras de 5 pulg. La tabla muestra las cantidades de tiras de 5 y 7 pulg
por tira de cartulina de 13 pulg que puede obtener mediante los cortes posibles :
2.4.1. Cuántas tiras de 13 pulg. debe cortar en las formas A y B para cumplir el
pedido minimizando el desperdicio?
2.4 .3. Será solución factible cortar 600 tiras en la forma A y 1100 tiras en la forma
B?
2.5. El Departamento de Servicios de un almacén proporciona servicios de reparación
para los electrodomésticos vendidos. Durante una semana cinco televisores, 12
grabadoras y 18 hornos microondas fueron devueltos para reparación . Dos
técnicos (Juan y Pedro) son contratados temporalmente para ayudar en dicho
68
INVEST IGi\C ION DE O PER AC IONES P¡\RA INGEN I ER I ;\S y c\D~ 1I N I S TR ;\C I ON DE H IPRESAS
departamento. En una jornada de ocho horas Juan puede reparar dos televisores
o 3 grabadoras o 4 hornos , mientras que Pedro puede reparar un televisor o 2
grabadoras o 3 hornos en el mismo tiempo. El salario diario está definido en
$20000 para Juan y $ 15000 para Pedro . Presente un modelo matemática para
calcular el número de horas contratadas para que los costos sean mínimos.
Se supone que los pagos que no se cubren son irrecuperables y, por lo tanto, no
producen ingresos por concepto de intereses. La competencia con otras financieras
del área hace que el banco asigne de los fondos totales al menos el 50% para
actividades agrícolas e industriales y los préstamos para vivienda deben ser al
menos el 50% de los asignados a actividades comerciales y compra de vehículos.
Además , este banco tiene una política establecida que especifica que la relación
global de pagos irrecuperables no puede ser superior a 0.004. Presente un modelo
matemático que le permita al banco calcular las asignaciones a cada uno de los
tipos de préstamos.
3. SOLUCiÓN GRÁFICA.
69
L U IS f\ L BE RTO RINCON ,\BR IL
4X, - 2X 2 S O 6X , + 8X 2 S 120
6X , + 5X 2 S 60 X, + X 2 2: 4
3X, + 4X 2 2: 12 X, S 12
X2 S 8 X 2 S 12
X " X 2 2: O X" X 2 2: O
c) Min(Z) = cX , + cX 2 d) Min(Z) = 6X , + 5X 2
2X , + 3X 2 2: 6 3X , + 2X 2 S 300
X, + X2 s5 2X , + 3X 2 S 300
X, - X2 S 1 X, 2: 50
X" X 2 2: 0 X 2 2: 50
4 . MÉTODO SIMPLEX.
70
INVrST I(ó¡\C ION DE O PER ,\C IONES P,\R ,\ I NCóEN IE RI ¡\S y ¡\D~ II N I ST R ¡\C l üN DE E~ I P R ES ,\ S
71
LU IS ALBERT O RI NCON AB RIL
4.5.1. Qué volumen debe producir de cada vacuna para optimizar utilidades?
4.5.2. Suponga que por una depresión económica el número de empleados debe
reducirse a 5 y el costo de producción a 5000 $/hora. entonces cuál es la nueva
solución?
4.5 .3. Idem anterior para 5 personas y 10000 $/hora.
4.5.4 . Suponga que la utilidad unitaria de la vacuna B se reduce a 1000 $/Iitro.
entonces cuál es la nueva solución?
72