Técnicas de Optimación
Técnicas de Optimación
FACULTAD DE INGENIERIA
Tecnicas De Optimaci6n
iii
EOOOl9
r:; "4) 1 0 0 0 3
~'!!
PRO LOGO
Marzo de 1994
'.j
111-12.1- Minim1zaci6n 88
111-12.2.- Desigualdad con sentido invertido 91
111-12.3.- Valores negativos para b 95
111-12.4.- lgualdades en las restricciones 102
111-12.5.- Variables no restringidas en signa 107
111-12.6.- Empate para entrar como variable basica 108
111-12.7.- Degeneraci6n 108
111-12.8.- Soluciones multiples 112
!11-12.9.- Ausencia de soluciones factibles 113
111-12.10.- Soluci6n optima sin limite 114
111-13.- METODO DE LAS DOS FASES 114
111-14.- PROBLEMAS 125
viii
\11-2.- METODO SIMPLEX REVISADO USANDO FASE I Y
FASEII 238
:x
Vill-5. 7.- Consideraciones 328
Vlll-6.- SISTEMAS DE REDES CON ACTIVIDADES EN EL
NODO 329
CAPITULO I
INTRODUCCION
1.1.· ANTECEDENTES.
Es por todos conocido el rapido desarrollo que han experimentado,
durante el presente siglo, el tamario y Ia complejidad de las organizaciones
humanas. Par ejemplo, el tamario de las empresas modernas, implica que
las decisiones administrativas pueden tener un efecto sabre grandes can-
tidades de capital y gran numero de personas. Los errores pueden ser
tremendamente costosos y una sola decision equivocada puede requerir
alios para rectificarse.
El cambia revolucionario sufrido par las organizaciones humanas, ha
traido como consecuencia Ia division del trabajo y Ia segmentacion de las
responsabilidades administrativas en estas organizaciones. Gran parte de los
componentes de una organizacion han tendido a crecer dentro de un plan
relatrvamente autonomo; debido a que cada uno de elias tiene sus propias
metas, han perdido, por lo tanto.la vision de como sus actividades y objetivos
se comparan con los de Ia organizacion . Lo que es buena para un com-
ponente, frecuentemente es perjudicial para algun otro. Es asi que, en una
empresa, por ejemplo, existen las siguientes tendencias:
Producci6n.- Pocas lfneas de productos y grandes corridas
de produccion (grandes inventarios).
V e n t a s _ - Abundantes y variados inventarios.
F ina nza s.- Minimizar inventarios
Person a I.- Producir inventarios solamente durante los
perfodos flojos (mantener producci6n constante).
(.Cual es en realidad Ia mejor politica para Ia organizaci6n? Es aqui
donde entra en acci6n Ia investigaci6n de operaciones.
La investigaci6n de operaciones siempre toma el punta de vista general
de Ia organizaci6n e intenta llegar a soluciones 6ptimas. o sea Ia mejor
soluci6n para Ia organizaci6n. Esta soluci6n no siempre es Ia optima para
sus componentes par separado. En algunos casas, Ia ejecuci6n de un
proyecto ocasiona efectos laterales, que si no son analizados adecuada-
mente, disminuyen su beneficia (p ej. problemas eco16gicos). par lo cual es
importante considerar Ia magnitud real del sistema considerado.
La investigaci6n de operaciones involucra estudios de grupo, es decir,
agrupa a toda clase de expertos en sus estudios. Esto es debido al aspecto
tan general que toma Ia misma y al hecho de que serfa imposible que un
especialista en investigaci6n de operaciones cuente con todos los co-
nocimientos necesarios para atacar los diversos problemas que se presentan
en el anal isis de sistemas de gran magnitud, par Ia que requiere de los demas
expertos, a quienes el coordina. Si bien Ia investigaci6n de operaciones no
sustituye los conocimientos particulares en una rama de Ia ciencia 6 tecnica,
si en cambia ayuda a integrar esos conocimientos con los de otras ramas.
La metodologfa de Ia investigaci6n de operaciones de ninguna manera
sustituye los conocimientos especfficos en ramas particulares.
Otro de los aspectos que favoreci6 el desarrollo de Ia investigaci6n de
operaciones, es el heche de que el ritmo de Ia empresa moderna es tal que
las decisiones se requieren mas rapidamente que nunca; el posponer Ia
acci6n puede poner en una decidida ventaja a un competidor.
No es sorprendente que el aumento en Ia dificultad de tamar las
decisiones haya requerido de esfuerzos para dar a esta actividad una base
mas objetiva y rutinaria. La investigaci6n de operaciones forma parte de este
esfuerzo.
Desde el punta de vista de Ia investigaci6n de operaciones, una
decision es una recomendaci6n de que se lleve a cabo un curso de acci6n
particular que afecta al sistema. Quien toma Ia decision, intenta que el sistema
sea mas "elective" para alcanzar las metas de Ia organizaci6n.
Generalmente, Ia medida de efectividad tiene unidades caracterfsticas
asociadas con ella (pesos, etc.), pero a veces son adimensionales (pro-
babilidad de destruir un blanco durante un ataque aereo).
Los orfgenes de Ia investigaci6n de operaciones datan de hace muchos
arias, cuando se hizo el intento de usar un acercamiento cientffico en Ia
administraci6n de las organizaciones. Sin embargo, estos orfgenes fueron
aislados y faltos de coordinaci6n; no fue sino hasta Ia segunda guerra
mundial cuando Ia complejidad de las operaciones militares hizo surgir lo
2
que hoy se conoce por investigacion de operaciones. A causa de Ia guerra
hubo Ia necesidad de distribuir. de manera eticiente. materiales escasos a las
distintas operaciones militares. y dentro de estas. a las actividades que !as
campanian. Para el efecto se reuni6 un gran numero de cientificos para que
reaiizaran una investigaci6n de las operaciones militares.
Debido al exito lagrada en el aspecta militar, Ia industria se tue inte-
resanda en este campo. Los investigadores descubrieron que Ia industria
sutria de los mismas problemas basicas que los militares.
Una vez iniciadas las tecnicas de Ia investigaci6n de operacianes,
atrajeron Ia atenci6n de numerosos investigadores, quienes lograron desa-
rrollarlas a grandes pasas. Las tecnicas desarrolladas contribuyeran a em-
pujar ia rapida Carrera que llevaba Ia investigaci6n de operaciones, de aquf
que para 1950 muchas de las tecnicas hubiesen alcanzada un grado de
desarrollo extraordinario.
En Ia misma epoca en que se perfeccionaban las tecnicas de Ia
investigaci6n de operaciones, tambien se lograron grandes avances en el
diserio de computadoras, lo que tacilit6 Ia utilizaci6n de las mismas, debido
a Ia gran cantidad de camputaci6n requerida par Ia investigaci6n de
operacianes.
Aunque no existe una definicion correcta para Ia investigaci6n de
operacianes, podemos decir que se aplica a problemas que canciernen a Ia
conducci6n y coordinaci6n de operaciones o actividades dentro de una
organizaci6n.
La investigaci6n de aperaciones se aplica extensamente en el campo
de los negocios, ei militar, Ia industria, el gobierno, haspitales, co-
municacianes, transparte, etc.
El metoda de analisis empleada es el del metoda cientffico, esta es, Ia
investigaci6n de aperacianes introduce el razanamienta cientffica creativo
dentro de las prapiedades tundamentales de las operaciones. Dewey, a
principios de siglo, indic6 que Ia salucion de un problema consiste en dar
respuesta a las siguientes preguntas:
1)i cual es el problema?
2)i cuales son las alternativas?
3}i cual aiternativa es mejor?
Las herramientas con las que cuenta Ia investigaci6n de aperaciones,
son las siguientes:
a) La programac16n lineal ha sida usada con exito en problemas de
asignaci6n de personal, mezciado de materiales, distribucion y
transportaci6n, asf como en paliticas de inversion
3
b) La programacion dinamica ha tenido exito en Ia planeacion de gastos de
propaganda, distribucion del esfuerzo de ventas y planeacion de !a
produccion.
c) Los fenomenos de espera se han aplicado a problemas de conges-
tionamiento de tratico, reparacion de maquinaria, planeaci6n de trafico
aereo y diseno de presas.
Otras de las herramientas de Ia investigac16n de operaciones son:
d) Teorfas de juegos.
e) Simulacion.
f) Teoria de inventarios.
g)Teorfa de redes.
En Ia actualidad existen dos organizaciones en Estados Unidos que
agrupan gente interesada en Ia investigacion de operaciones:
The Operations Research Society of America, que edita Ia revista
Operations Research.
The Institute of Management Science, que edita Ia revista Management
Science.
1.2.- METODOLOGIA.
Aunque no existe una secuencia fija de pasos a seguir en un estudio
de investigacion de operaciones, podemos analizar el procedimiento dentro
de un problema tfpico.
a) Formulaci6n del problema. En Ia practica, los problemas surgen
de una manera vaga, raz6n par Ia cual es necesario estudiar el sistema
bajo consideraci6n y definir exactamente el problema, indicando clara-
mente cuales son los objetivos, las restricciones y que es lo que se
puede hacer: interrelaciones existentes, alternativas posibles,
limitaci6n del tiempo, etc. Seve que Ia formulaci6n es cruciaL pues es
imposible obtener buenos resultados partiendo de un mal planteamien-
to. De aquf que en cuanto surjan nuevas aspectos durante el desarrollo
del estudio, Ia formulaci6n debe ser revisada para ver si no ha cam-
biado.
Se debe tener cuidado, ya que aunque se mencion6 que el objetivo
principal era Ia soluci6n optima para toda Ia organizaci6n, existen casas
en que el problema afecta a una sola secci6n y el incluir todas las de mas
complica terriblemente el problema.
4
b) Construcci6n del modelo matematico (muchas veces es este ei
paso mas sencillo). Consiste en Ia reformulaci6n del problema. de
manera tal, que sea conveniente para analizarlo. En este paso repre-
sentamcis matematicamente ei sistema bajo estudio.
Un modelo matematico es una representacion idealizada de Ia
realidad,expresada en terminos de sfmbolos matematicos.
Par ejemplo,en programaci6n lineal podemos representar las n va-
rables de decision, cuyo valor se va a determinar, como: x,, x 2 •.. ,x" .
La medida compuesta de efectividad (llamada funcion objetivo) (p. ej.
ganancia), se expresa entonces como una funcion matematica asf:
Z = 5x. ·- 3x 2 + 7x 3 ~ .. .. -r 2x"
Las restricciones a ias variables de decision tambu!m se expresan
matematicamente:
3x. + 8x2 + .... - 3X 0
:S 100
5
Muchas veces Ia so!uci6n obtenida se emplea como datos para
formu!ar nuevamente el problema y asi, mediante una serie de ciclos,
se obtiene una mejor solucion.
d) Probar el modelo y su soluci6n. No importa que tan buena parezca
nuestro modelo, siempre se debe probar. El primer paso consiste en
verificar errores obvios o deta!les pasados par alto.
Prueba retrospectiva. Cuando es aplicable, cansiste en analizar Ia
historia de Ia campania y encontrar que hubiese pasado de haberse
aplicado este modele. Este metoda puede mostrar fallas en el modelo
y asi modificarlo. Una gran desventaja de este metoda, es que utiliza
los mismos datos con ayuda de los cuales se desarrollo el modele. Otra
falla es que las condiciones pueden haber cambiado y que los datos
pasados ya no representen las condiciones futuras.
Otra prueba consiste en decir que no se aplique el modelo hasta dentro
de seis meses y observar que hubiese pasado en este lapso de haberse
aplicado el modele, pero para entonces las condiciones pueden haber
cambiado.
e) Establecer control sobre Ia soluci6n. Este paso se usa solo cuando
el modelo desarrollado se va a usar repetidamente, debido a que las
condiciones (datos) pueden estar variando constantemente.
Este paso puede incluir el descubrir los parametro criticos, (aquellos
cuyo cambia puede afectar significativamente Ia soluci6n) Esto se
efectua mediante el analisis de sensibilidad.
Una vez descubiertos los parametres criticos, se puede establecer un
metoda estadfstico para detectar cambios significativos en los mismos.
Cuando se descubre un cambia en un pan3metro critico, es necesario
ajustar Ia solucion al nuevo valor.
f) lmplementaci6n.- El grupo de investigadores debe participar activa-
mente en esta fase. con objeto de supervisar que Ia soluci6n obtenida
sea convertida con exactitud en un procedimiento operative.
Los pasos a seguir son:
1) Se comunica a Ia administraci6n Ia soluci6n abtenida.
2) La administracion y los investigadares comparten, conjunta-
mente, Ia responsabilidad deponer esta solucion en operacion.
3) Adiestramiento del personal involucrado en Ia solucion.
4) Retroalimentacion.
6
CAPITULO II
METODOS
MATEMATICOS
DE CALCULO
11-1.- MATRICES.
Definicion 11-1.-
Una matriz es un arreglo rectangular de numeros escalares.
Denominaremos por "m" el numero de renglones y par ' n" el nurnero
de columnas de una rnatriz.
Una matriz con "m" rer.giones y ''n" columnas. se conoce como una
matriz m x n.
Una matriz n x n se dice que es de orden n.
Tenemos entonces que:
!4 7
I 8
A= '~
:~ 2 6 14
I 1 5 12 3
es una matriz 2 x 4 y
4 5
B = 8 3
7 9
es una martiz 3 x 2
A los numeros que integran el arreglo rectangular, se les denomina
elementos de Ia matriz ·. Es asi que los numeros: 4, 3. 1, ?, 2, 5, 8, 6. 12, 1,
14. 3, son los elementos de Ia matriz A, mientras que: 4. 8, 7, 5. 3, 9, son los
elementos de Ia matriz B.
Las matnces las denominaremos con letras mayusculas y a los elemen-
tos correspondientes de Ia matriz con Ia misma letra, pero minuscula, seguida
de 2 subindices, el primero de los cuales representa el renglon al que el
elemento corresponde, mientras que el segundo indica Ia columna.
Es asi que los elementos de Ia matriz A seran a j , por ejemplo: 1
a 24
= 14;
A=
rI .
a,, a. 2 ... a.. . .. a.n
I
matriz. N6tese que mientras los elementos aii y ak 1 pueden no ser iguales,
se considera que las matrices i a1j ii y if ak 1 !I son identicas.
Una matriz no posee ningun valor numerico (a diferencia de su deter-
minante).
Sin embargo, es conveniente realizar ciertas manipulaciones con los
arreglos de numeros; es asi que se han desarrollado reglas que permiten
realizar operaciones con matrices, las cuales son aniilogas a las operaciones
aritmeticas.
8
Definicion 11-2
Sean: A = , a 'I:, y
dos matrices. Se dice que A y B son iguales , esto es A = B. sf y solo sf todos
y cada una de los elementos correspondientes son iguales. Es decir, A y B
tienen exactamente los mismos componentes (a, = b, para todos los val ores
de i y de j) ·
De Ia anterior es clara que para que A = B. es necesario que ambas
matrices tengan el mismo numero de reng!ones, asi como el mismo numero
de columnas entre sf.
Definicion. 11-3
La diagonal principal de una matriz :; aii i es el conjunto de elementos
{a 11 .. atd donde t = min {m, n}.
Definicion. 11-4
Una matriz diagonal es una matriz cuadrada de arden n, en Ia cuallos
elementos que no pertenecen a Ia diagonal principal son cero.
Definicion. 11-5
La operaci6n de multiplicar una matriz por un numero k, se efectua
multiplicanda cada una de los elementos de Ia matriz par k.
Asf:
rka,
k~,
ka, 2
ka22
ka,,
ka 2,,
ka,n
ka2n
l
I
kA II ka; 1 il
Ik~, ka 2 "' ka 11 kain
rs
51
4 13 7 l
I
140 20 65 35 l
11
1..
0 9 2
J
I IL 5 0 45 10 l
J
9
Definicion. 11-6
Para 'sumar dos matrices A y B (que deben tener el mismo numero de
renglones y de columnas entre s0, sola mente es necesario sumar los elemen-
tos correspondientes de ambas, lo cual puede quedar representado en Ia
siguiente forma:
A+ B 'a
' ,, + b
Ejemplo. 11-2
3
H
4
1
0
-3l
~J + r
L
~
-1
7
2
~~ "
=
lw
l-~ 8
2
n
Definicion. 11-7
La resta de dos matrices A y B podemos interpretarla como Ia suma de
Ia matriz A con Ia matriz C donde C = kB = ( -1) B .
Asf, tenemos:
A- B = A+ C = A + (-1) B = II a ,1 - b,i il
4 -3 l [ 7 -1 10 l -4
1 8
0 2
J- 3 7 -4
-2 2 3
j = [ -1~
Definicion. 11-8 (Multiplicaci6n)
Para encontrar el elemento cii en el rengl6n i y Ia columna j de Ia matriz
resultante de multiplicar A par B, es necesario multiplicar cada elemento del
rengl6n i de A por el elemento correspondiente de Ia columna j de B y sumar
estos productos; observese que Ia multiplicaci6n de dos matrices esta
definida si y solo si el numero de columnas de A es igual al numero de
renglones de B, (n6tese que el numero de renglones de A y de columnas de
B puede ser cualquiera).
En general, si:
y :r b ! i rxr
10
su producto sera
ll
C = AB = '\' a ' b, c
"!; ==.
Ejemplo. 11-4
1 -1
AB = I~
l-2
4
2
-1
1
-2
;l
2
J
0
2
3 -2
2
1
Ejemplo. 11-5
a.-
A = J2 4 lJ
L3 1
B = l' 31 ·24 lJ AB = [ 10 20 j-
10 1o
BA = t12 14
l14 a
lJ
En Ia mayorfa de los casas en que A B esta definida. B A no lo esta.
b.-
Nota. 11-1
La division entre matrices no esta definida.
Definicion. 11-9
Se llama matriz cera o nula a aqueila matriz de cualquier arden cuyos
elementos son todos · cera' . Se representa par:
0 0 0
0 0 0
¢
0 0 0
11
Ejercicio il-1
Sean.
r~ - ~' 6 2
A=
j~
I 1
4 4:I
6 I'
~
B=
!5 4 3 01
7 2 5 1: C= :4
15 0 4 31,.,
1
L7 2 3i J ~
1 3 2 o) I
i? 5 ~I
i3 4 41
- :5 0 01
D=
I
! 1 6 7i E= io 6 o! F= l2 1l
I
j? 2 3J lo 0 1j l3 oj
L L
2-B ~ C porque b ~ c
= 38
r15
1
12 9 18 l
7.- kB =
l 21
3
6
9
15
6 ; J
8-A+D~[~
8
14
12
4
1!6 1J
12
:-2 4 4 l
9A-E~l ~; ~j
A - C no esta definida.
r104 43 43 l
10.- BC = l7929
I
24
4
40
19
J
BF no esta definida.
I
~~~ ~~ ~~21 25~g j'
11.- CB = 29 24
171 41 48 47
'-
Vemos que 8 C ~ CB
2.- 1 A = A
3.- A+ 0 =A
4.- A -A= 0
5.-A-rB=-B+A
6.- A ( B + C) = AB + AC
7.- (A T 8) c= AC + BC
8.- ( A + B ) + C = A + ( B + C)
13
!eyes que se satisfacen siempre y cuando las operaciones propuestas esten
bien definidas
Definicion 11-1 o
Dada Ia matriz A, Ia matriz AT, obtenida de A intercambiando los
renglones par las columnas en A, se llama Ia matriz transpuesta de A. Si
AT = :j aiJ'i!, el elemento aiJ. que aparece en el rengl6n i y Ia columna j de
Ar, es el elemento aji que aparece en el rengl6n j columna i de A.
Ejemplo 11-6
17 6 5 l
A= AT= 13o 2
7
1
2
j'
l4 1 3
(A + B )T = AT + BT
Definicion 11-11
La matriz adjunta es aquella matriz cuadrada formada par tps cofac-
tores de Ia transpuesta de Ia matriz cuadrada dada. Se escribe A
Ejemplo 11-7
AT = [ ~ ~ -~ 1
0 1 2 J
entonces:
14
,, 4 1 4 ;1
-'-
0
lo
I 7
2 I
. C· 2 I
:
13 ·- 2 -4 -2 1
Ar = 2 -2 ;2 3
I 1
:o Ar == !
-8 4 -2
I' 2 2 ! io I
I I
IL 12 -1() -3
3 -2 i2 -2 12 3 I
+ -'-
L iO 4 I
I !1 4 I i1 0
J
Definicion 11-12
La matriz identidad ( 1 ) . es una matriz diagonal cuyos elementos en Ia
d1agonai pr1nc1pal son todos unos.
r1
! I 0 0 0
:o 0 0
I= lo 0 0
I :
10 0 0
L
IA=A=AI
Definicion 11-13
Cuando una matriz Ia subdividimos en varias matrices mas pequenas,
a estas se les conoce con el nombre de submatrices.
Ejemplo 11-8
I
~a~
a.2 a,3 a.<l a,. A<
-i I
A=
i a2, a:>< a::3
a __
a2< l I ~- A"22
:a,.
I "
a3~ ~~ a~ -'I
15
donde
)a,.=a ..
Sea:
B= siendo: B2 =
AB
b, l
B= .
rb2
l...
bm J
16
Los vectores son importantes en Ia teoria de matrices, debido a que
cualquier matriz m x n puede ser subdividida ya sea en m vectores rengl6n
o en n vectores columna. con lo cual se pueden analizar importantes
propiedades de las matrices.
Si tenemos que:
donde las C son constantes, se dice que {3 es una combinaci6n lineal de las
x,. Tambien decimos que j3 es lineal mente dependiente de las x,, sip puede
ser expresada como una combinaci6n lineal de las x,.
~
las x;. Una relaci6n donde todas las C; '-" 0, se conoce como una relaci6n
lineal trivial; una relaci6n en Ia cual cuando menos uno de los coeficientes es
diferente de cera, se conoce como una realci6n lineal no trivial.
Definicion 11-14
Se dice que un conjunto de vectores son linealmente dependientes si
existe una relaci6n lineal no trrvial entre elias. En caso contrario, el conjunto
es linealmente independiente. Debe notarse que cualquier conjunto de
vectores que incluya el vector cera es linealmente independiente.
Si un conjunto {x;} de vectores es linealmente dependiente, entonces
rr
existe una relaci6n lineal no trivial de Ia forma 2:C; X; = 0 (no necesaria-
que es uno de los vectores del conjunto {xj, es una combinaci6n lineal de
los otros.
Si un conjunto {x;} de vectores se sabe que es linealmente inde-
pendiente y se obtiene una relaci6n lineal IC; X1 = 0, podemos entonces
concluir que todas las C . = o. ' =,
17
Es clare que el concepto de independencia lineal de un conjunto no
tendna sentido si un vector del conjunto pudiera incluirse un n(Jmero ar-
bitrario de veces en una posible relacion.
{x, x2 } en donde X1 = ~
Ejemplo. 11-1 0
Sea:
donde:
X,= (
~= 1,(
lS- ( 0,
1, 1,
0,
1'
0)
1)
1)
x 4 = ( 1, 1, 1)
por lo tanto:
X--
- x1+x2+x3
-----
4 2
Sin embargo, cualquier subconjunto de tres de los cuatro vectores es
linealmente independiente.
18
Ejemplo. 11-11
a -.- b -r 2c =0 .. (1)
a ~- b -r 2c +- d =
0 /0)
b + c =0 :(3)
a -r b +- 2c =0 .. (4)
de (3): b = -c
en (1): a + b - 2b = 0 de donde a = b
en (2) a - a - 2a + d = 0 de donde d = 2a
Ejemplo. 11-12
w, = (1.00,1)
w;;=(0,10,1)
w" = (0.0,1,1)
w4 ~c (1.1,1,1)
1Q
I 1 \1.~ ~8:C~S
a ~
d = 0 (1)
b - d =0 (2)
c + d =0 (3)
a - b- c + d = 0 .... (4)
de (1 ): a ~ --d
de (2): b = --d
de (3): c = --d
en (4)· -d-d-d.,.d=O = > d = 0
Definicion 11.15
El mayor numero de vectores linealmente independientes que puede
ser encontrado dentro de un con junto de vectores. se conoce como el rango
del conjunto.
Ejemplo. 11-13
Ejemplo. 11-14
Sea:
20
donde:
s.~\3.7.1)
s.=(6i4.2)
s. ·~ ( 2. 0, 3)
de (1) - 3(3):
-7c =0 ~ c = 0
en (1 ):
3a =6b =0 ~ b = -1/2a
si:
a = 1 ~ b = -1/2
par lo tanto:
s, - 1;2s2 = 0
a ( 3, 7, 1 ) ... b ( 2, 0, 3 ) ~c ( 0, 0, 0 )
de donde:
3a , 2 b = 0 .... (1)
?a =o (2)
a .,._ 3 b = 0 .... (3)
21
de (2):
a= 0
en (1):
b = 0
Teorema. 11-1
Cualesquiera n -1- 1 vectores, en un espacio vectorial de n dimen-
siones, son lineal mente dependientes.
Ejemplo. II~ 15
Definicion. 11-16
Se dice que un subconjunto de vectores linealmente independientes es
Ia base del conjunto de vectores a que el mismo pertenece, cuando todos y
cada uno de los vectores del conjunto original son una combinaci6n lineal
de los vectores en este subconjunto.
Una base no necesariamente es (mica
Ejemplo. 11~16
22
Ejemplo. II. 17
52 -- 2s.
53 -- 53
Teorema 11.2
Un subconjunto de r vectores linealmente independientes. tornados de
un conjunto de vectores, es una base del mismo si y solo si el conjunto tiene
rango r.
Demostraci6n:
Teorema. 11-3
Si un conjunto de n vectores tiene una base con un numero finito de
elementos. entonces todas las demas bases, en caso de existir, son finitas y
cuentan con ei mismo numero de elementos.
Teorema. 11-4
Un conjunto de n vectores en un espacio vectorial de n dimensiones,
es una base si y solo si es iinealmente independiente.
23
Definicion. 11-17
El rango de los rengiones de una matriz, es el rangu del conjunto de
sus vectores rengl6n. El rango de las columnas de una matriz es ei rango del
conjunto de sus vectores columna.
Ejemplo. 11-18
l
r1 0 p, = 3 (ver ejemplo 11-11)
I 1 -1 1 1
A= Pc 3 (a simple vista seve Ia depen-
12 2 1 2
=
/o
L
0 0
"
dencia entre ias columnas 1 y 4)
Teorema. 11.5
El rango de los renglones y el rango de las columnas de una matriz son
iguales. Es par esto que solamente se habla del rango p de una matriz.
Dada una matriz A, diferente de Ia matriz cera, (..cuando existe una
matriz A- 1 ,llamada el inverso de A, tal que:
A A- 1 = A- 1 A = I?
1
Si A no es una matriz cuadrada, entonces nunca existe tal matriz A- ,
de donde solo una matriz cuadrada puede tener inverso. Esto es debido a
que A- 1 deberia tener diferente numero de renglones en cada caso para que
los productos A A- 1 y A- 1 A estuviesen detinidos.
Definicion 11-18
Una matriz se llama no singular si su rango iguala tanto el numero de
sus renglones como el numero de sus columnas, en caso contrario se llama
singular. De aquf que solo las matrices cuadradas puedan ser no singulares.
Definicion 11-19
Una matriz m x m se llama no singular si sus vectores columna son
linealmente independientes.
24
Una matriz es no singular si y solo sf su determinants es diferente de
cera.
Teorema. 11-6
a) Si A es no singular, existe una matriz no singular A-', llamada inverso
de A, tal que A A-' = A 1 A = I
b) Si A es no singular y B es una matriz para Ia cual ya sea A B = I 6
B A = I, entonces B = A- 1 .
c) Solo las matrices no singulares tienen inverses.
Teorema. 11-7
Si A y B son no singulares, entonces:
a) A B es no singular.
b) (A B) -1 = B- 1 A- 1
c) A-' es no singular.
d) (A- 1)_, = A
e) Para a '# 0, (a A) es no singular y (a A)-' = a- 1 A- 1
Teorema 11-8
El rango de una matriz (no necesariamente cuadrada) nose altera si
esta se multilplica par una matriz no singular.
Ejemplo 11-19
Si: A=
r1
L..
2
3
2
3
4
~l
6 ' J
entonces A-t=
l ~ -2 -~
r-2 0
3
l
Comprobaci6n:
~ r
11 0 0 l 0 0
AA-1 = lo 1 oI A-1A = 16 1 0 I
IQ
L 0 1 I
J
lo
L..
0 1
J
Existen varies metodos para invertir una matriz:
25
1: ' \-1A ;RICES
Metodo 1.-
Ejemplo 11-20
En el ejemplo anterior:
I A! = 18 + 24 + 24-27-16-24 = -1
dado que:
entonces:
~ -1
-=-kll
-1
+ 11 21
12 3,
-1 J.
18-16 12-12 8-9
-1 -1
A_ 1 = _ 12-~\2 6-9 4-6
r -1 -1 -1
3-4
l
8-9 4-6
-1 -1 -1 J
26
-2 0 1
A~· = 0 3 -2
Ejemplo li-21
A= l
r 5 -4 l A ~ ~5 .,.. 4 "" -1
-1 J
-1
-1
-4
-1 I' -4 l
5 -5 I
-1 -1 L1
1
J
Metodo 2.
Ejemplo 11-22
.,
5 3 -1 I
I
A= 2 0 4 I
-2 3 I
J
5 2 -2 l I
AT= 3 0 3
I
-,• 4 1
J
-12 --6 12
AI = -10 3 -22
6 -21 -6
27
A ~ 0 - 24 - 6 - 0 - 60 - 6 co= - 96
N6tese que:
AI I ""A A~'
1 -1
l-12 -6 - -
12l 8 16 8
~221
A- - 1 5 -1 11
-96
-10 3 = 48 32
~~
48 I
I
Metodo 3.
[A I l = I; 2
3
3
4
1
0
0
1
l3 4 6 0 0
~ -2
l
2 3 1 0 0
I -1 -2 1 0
L o -2 -3 -3 0 1
c) Se busca dejar el segundo vector columna como {0, 1, 0} mediante
operaciones elementales.
28
1 0 -- ~r -3 2 0
0 1 2 2 -1 0
0 0 -2
0 0 -2 0 1
0 1 0 0 3 -2
0 0 -2
A_,=
l·-2
I 0
0 1 l
I -2
3 -21 J
L
Comprobaci6n:
A A_,
r~
13
L
2
3
4
:][-~ 0 1
3 -2
-2 1 !
=
l
j
rL
1
0
0
0
1
0
0 'I
0
J
r-~
0 1 l r, 2
~l
0 0 l
A-' A
-2
3 -2 I
J L3
12 3
4 6 J
=
[ 0
0
1
0
0
1 I
J
Teorema. 11-9
29
•1-2 ECUAClONES _:NEALES s:MUL 'ANEc".S
coi k. col k_
ro
0
0
0 0
a! '.k ~
0
--1 0
"
a''·· 2 - ·
a' 2.K . 1
2
~·
.l
I
I
.0 0 0 0 0 0 I
A' == I I'
I
lo 0 0 0 0 0 0 0 0 0 j
Ia cual satisface las siguientes condiciones:
a) Existe cuando menos un elemento diferente de cera en cada uno
de los primeros p renglo11es de A' y los elementos en los restantes
(m- p) renglones son cero.
b) El primer elemento diferente de cero que aparece en el rengl6n
i (i ~ p ) es un uno que aparece en Ia columna k,, donde
k, < k2 < ... < kP'
c) En Ia columna k el unico elemento :;eO es un uno en el rengl6n i
Teorema II -1 o
Una matriz A tiene rango p, si cuando menos uno de sus menores de
orden p es diferente de cera, mientras que todos los menores de arden
(p -,- 1 ), si existen, son cero. Una matriz nula se dice que tiene rango 0.
2x 1 + x 2 - X3 = 7 .... ( f3 )
El conjunto ordenado de valores x 1 "" 3, x 2 "" 2, x 3 = 1, se dice que
es una soluci6n de a {6 de j3), debido a que si se sustituyen estos valores en
Iugar de x 1, x2 , x3 en Ia primera ecuaci6n (6 en Ia segunda), se produce Ia
identidad 2 =
2 (6 7 =
7). La soluci6n (3, 2, 1.) se dice que satisface Ia
ecuaci6n a ( /3).
30
Supongamos el siguiente sistema de ecuactones lineales simultaneas
a 11 x 1 - a~ 2 x 2 + ... - a 11 x1 -
3 21 x 1 - a22 X2 - · .. · ,. a2 X -r
1 1
\'
L...a, x. = b para (i =- 1, . m)
,=:
6 en forma matricial:
A X= B
en donde:
Definicion 11·20
Una soiuci6n para ia i-ava ecuaci6n, es un conjunto ordenado de
nLJmeros (x' 1 . x' 2 , ... , x·nl· tal que:
31
Definicion 11-21
Un conjunto ordenado de numeros se dice que es una soluci6n de un
sistema de ecuaciones. sf y solo si es una soluci6n para todas y cada una de
las ecuaciones del sistema.
Definicion 11-22
Dadas las matrices A y B. Ia matriz aumentada A, B de un sistema de
ecuaciones lineales, se define como:
Definicion 11-23
Un sistema de ecuaciones que tenga soluci6n (unica o no) se llama
consistente. En caso contrario es inconsistente.
Definicion 11-24
Si en un sistema de ecuaciones una ecuaci6n es un combinaci6n lineal
de las otras, se dice que esta ultima es dependiente de las otras. La ecuaci6n
dependiente se llama redundante.
Definicion 11-25
Un sistema que no contiene redundancias se llama independiente.
Un sistema lineales claramente inconsistente, si es posible el encontrar
una combinaci6n lineal de las ecuaciones del sistema que que tenga Ia
siguiente forma:
0 X, + 0 >s + ~ 0 x" = d donde d ;;t:Q
32
Definicion 11-26
Un sistema can6nico, con un subconjunto ordenado de variables
llamadas bastcas, es un sistema tal, que para cada i, Ia i-ava vanable basica
tiene un coeficiente unitario en Ia i-ava ecuaci6n, y un coeficiente cera en
todas las demas.
Definicion 11-27
Dos sistemas se !Iaman equivalentes si un sistema puede ser derivado
del otro insertando o eliminando ecuaciones redundantes.
Teorema 11-1 0
Sistemas equivalentes tienen el mismo conjunto de soluciones.
Teorema 11-11
El sistema de ecuaciones lineales simultaneas A X =c B tiene soluci6n
(una o varias) si y solo si el rangode A es igual al rango de Ia matriz aumentada
A, B. Siempre que una soluci6n exista, todas las soluciones pueden ser
expresadas en terminos de v = n - p parametres independientes, donde
p es el rango de A.
J ustificaci6n:
33
soluci6n, que en este caso es unica
2.- Pix < n. Aquf existiran p ~~ vectores independientes y v = n-PrAJ
vectores dependientes (o parametres independientes), a los cuales
puede ser asignado cualquier valor arbitrario sin que el valor
escogido impida que exista una soluci6n para los vectores en !a
base, de aqui que en esta situaci6n existan una infinidad de
soluciones para el sistema de ecuaciones.
Esto es debido a que A' X = B' es tal que A esta en forma hermite.
Recordamos que xk, aparece solo en Ia i-ava ecuaci6n, ademas coeficientes
diferentes de cera aparecen solo en las primeras p ecuaciones. Dado que x,,
aparece en una sola ecuaci6n con coeficiente uno, a las restantes n ~ p
incognitas se les puede dar val ores arbitrarios para obtener los valores de xk,
correspond ientes.
Es importante el hacer notar que cuando p[A.BJ =piAl = r, donde
r < m, entonces tenemos que (m ~ r ) de nuestras ecuaciones deberan ser
combinaciones lineales de las otras r. De aqul que (m ~ r) ecuaciones
redundantes pueden ser eliminadas sin que esto afecte a Ia 61as soluciones.
Metodos de soluci6n.-
34
una de las n ecuaciones debe contener exactamente una de
estas variables ya que A A X = A-· B. de donde X = B.
La soluci6n deseada se lee directamente en las ecuaciones
resultantes
b)pl""' =P:"+1
En este caso. el sistema de ecuaciones no tine ninguna
soluci6n, debido a que el sistema es inconsistente. Este
caso no puede presentarse.
Ejemplo 11-23
2x 1 - 2x 2 j- 4x 3 - 22 . " ( 1)
--X~ ~ 3x - 9 .. (2)
2
4x 2 - x3 = 10 (3)
Eliminamos X1 de todas las ecuaciones menos de Ia primera.
35
unitario. dividiendo d1cha ecuaci6n entre (-5).
x, ~ 3x 3 = 21 ... (1 ")
x2 _;._ X2 = 10 (2 ')
X3 = 6 ... (3 ")
r -3 14 -12 l
- --
1-20 -20 --20 I
I -1 -2 -4
A-'=
1-20 -20 -20
I
1-4 -8 ___i_ I
_-20 -20 -20 J
36
r
i 3 -14 12 - 66 126 120
20
~--~
20
-
20
~-
Lx3 l6
Caso 2.- Cuando m7:n y/o A es singular
a) Pt,B! = PrAJ + 1
Ya v1mos que en este caso no hay soluci6n. Cuando se
aplica el metoda de Gauss-Jordan. se obtendra una
ecuaci6n en Ia cual el lado izquierdo se habra hecho cera,
mientras que ellado derecho sera diferente de cera, es decir,
nuestro sistema es inconsistente.
b) PrAs: = P,A 1 = n
Aquf. m > n (puesto que p = n)
En este caso tenemos (m- n) ecuaciones redundantes, las
cuales serim totalmente eliminadas par el metoda de
Gauss-Jordan, de forma tal que Ia soluci6n (mica sera iden-
tificada igual que antes (caso 1a)_
37
al sistema de ecuaciont:Os.
La transferencia de estas (n - r) variables al lado derecho.
dejara las r variables como una funci6n de los parametros
independientes.
Ejemplo 11-24
eliminando x 3 de (4"):
x, = 1 ................ (1 '")
= 1 ................ (2"')
x3 = 1 ................ (3'")
0 =-2 ............... (4"')
38
Ejemplo 11-25.
X2 - ~ x3 = 0 . (2")
X3 = 0 .............. (3")
Ejemplo 11-26
Caso 2 (c).
4x, + 4x2 ~ 4x 3 0
5x, -r 4X2 .,- 3x 3 -- X4 4
-2x, - 2~ - x 3 +2X 4 = -3
11 x, + ~ + 4x 3 + x. = 11
39
La matriz aumentada sera.
r
i 4 4 4 0 0
i 5 4 3 --1 4
1-2 2 -1 2 -3
l11 6 4 11
-~
11 0 -1 -1
I~ ~ ~ 2 -3 'J
lo o 3 6 -9
finalmente obtenemos:
~
0
1
0
0 -3
1
2 I l
0 1 2 -3
0 0 0
L
0
J
Como obtuvimos toda una hilera de ceres. Ia ecuaci6n (4) era redundante,
de donde el sistema de ecuaciones correspondiente es:
x, + X4 = 1
X;, - 3x 4 = 2
x 3 + 2x 4 = -3
o bien:
x. = 1 -- x4
X2 = 2 + 3x 4
x3 ~-3 - 2x 4
x4 = x4
donde x4 es el parametro independiente y
(x,. X;,. JS. X4 ) = (1, 2, -3, 0) + (-1, 3, --2. 1)x.
sera Ia soluci6n expresada en forma vectorial.
40
CAPITULO Ill
PROGRAMACION LINEAL
DEFINICIONES V TEOREMAS
Definicion .- 111.1.-
Uniendo dos puntas cualesquiera (x 1 ', x· 2 , ... , x'cl y (x" 1 , x" 2 , ... , x"n)
existe una colecci6n de puntas. conocidos con el nombre de 'segmento
lineal", tales que
(X 1 , X2 , .. , X")= (),x 1 '+(1-).)x,",llS'+(1-,.\)x2 ", ... ,,.\xn'+(1-))xn")
Sean:
(x'" x' 2 ) = ( 2 , 6 )
(x'' 1 , x'' 2 ) = ( 4, 3)
41
6 \"\ X l =-(2, 6)
~-'
4
3
V1=!4:11
2
2 3 4 5 6 7 8 X.
Definicion 111.2
Un conjunto convexo es una coiecc1on de puntas tales, que para
cualesquiera dos puntas en Ia colecci6n, el segmento lineal que los une esta
par completo tambien en Ia colecci6n.
Ejemplo 111-2
,_
l
I •: ,,,,
·ff--;;~
U
S CAE t ·,21
~EN'qO
S.J)
f-----'-----x.
' ,---x
42
Definicion 111-3
Un punta extrema de un conjunto convexo. es un punta en el conjunto.
el cual no cae en ningun segmento !meal que una cualesquiera otros dos
puntas en el conJunto
Geometricamente, un punta extrema es una esquina.
En Ia figura 111-3 los puntas extrema son: (0,0), (0.4), (3.4), (5.2), (5.0).
Con el objeto de poder llegar a una decision optima, todas las com-
binaciones posibles de distribuir los recursos a las actividades deben de ser
analizadas.
43
En te;minos generales, Ia programaci6n lineal puede ser utilizada para
problemas de optimizac16n, en los cuales las sigUientes condiciones se
satisfacen.
111-1.- MODELOS
Los modelos son representaciones de Ia realidad
Generalmente, se pueden construir modelos que son mucho mas
simples que Ia realidad, pero que aun se pueden utilizar para predecir y
explicar los fen6menos con alto grad a de exactitud. La raz6n estriba en que,
a pesar de que se suele necesitar de gran cantidad de variables para predecir
un ten6meno con exactitud pertecta, solo un pequeiio numero de estas son
las responsables de Ia mayor parte de Ia exactitud. La clave esta en encontrar
las variables correctas y Ia relaci6n entre elias.
44
Este tipo de modelos son particularmente utiies para representar
situaciones estaticas, o dinamicas en un momenta especifico de tiempo. pero
son generalmente diffciles de usar para representar situaciones dinamicas.
tales como Ia operaci6n de una fabrica.
A) Cuantitativos y cualitativos:
45
C) Estaticos y dinamicos:
D) Probabilisticos y deterministicos:
46
Es necesario aclarar que noes esta una diferencia basica que 1mplique
que el tratamiento dado a un problema en un caso. tenga un tratamiento
completamente diferente en otro. smo que por el comrario. el tratamiento es
el mismo (sa vo las diferencias debidas a Max. o Min ) Existe. ademas. una
1
Encontrar: x,, x 2,. ,xn (los niveiesde nactividades), tales que maximicen
Ia funci6n lineal:
a, x, ~
a.2 x2 -t· ~
a! X ~ +
a, .. X,. ~ b
z =ex
sujeta a Ia condicion de que .A.X ~ 8
47
Z se conoce como Ia funci6n objetivo y !as desigualdades lineales en
AX ~ B se !Iaman restricciones lineales.
a
De aqui en adelante se supondra que el modelo est en Ia forma generaL
a menos que se especitique de otra manera.
48
Vemos entonces. que las desiguaidades no rienen otro significado que
el de indicar matematicamente que Ia suma de las cantidades de recurso
escaso i utilizado en las n actividades. debera ser menor o igual que Ia
cantidad del mismo de que se dispone.
donde:
c1 ~- ganancia obtenida al producir una unidad del producto i
c 2 = ganancia obtenida al producir una unidad del producto 2
etc.
49
(j ~~ 1, 2, ..... n) ( i = 1, 2, ... m)
a,. x, _,_ s b,
a,2 x2 + + a, X + T
a,n X"
'
a,, x, + a22 x2 T + a2j XI .. T
a2n xn s b2
a,, X, + a,2 x2 -r ~
a,, XI + ...... I a rn x rr
s b,
Ejemplo 111-4
50
Ia cantidad x. que debera fabnca~se del producto 1 y ia cantidad '-c que debe
fabr;carse del producto 2. para que nuestra ganancia sea maxima
Soluci6n:
maximizar
SUjetO a.
X, $ 400
x2 $ 300
X. 4·
x2 $ 600
X. ::::: 0
x2 ::::: 0
Ejemplo 111-5
~.
Oo
horas:
Carttdad m 2 madera m 2 madera 1 Cons1dera- i·
hombre
que pro~ ArtiCUlO del H -
1 de1 # 2 CIOi'eS de Ganancia !I
consu- 'i
i duciremos consurn1da consurT11da ~ ·;enta ~!
m1das
i~~-~-~ 41
1
' x....... rnesas 5 2 3 2: 40 12 :i
t-------+- < ~~---4~
,. x5
:r-- --- ----,.-1
sillas
--t--
3 2 2:130
___'
. -----~
5 1'1
it xe eswtorios 9 4 5 2: 30 15 i!
ii-- -~---~-~~~---~- -~-- -+---~-~--11
,, ,,
iibreros 12 10 ~10
i(___x_-+ 10 "
·':' dispon;ble '
~-,-50_0_1
1000 800
+---~--j
I
===:h===-===~
Ejemplo 111-6.-
52
- :. :.... : .. ..:. '~;:
Solucion.-
Sea x. ~· r'-lo de cinturones tipo A. fabricados a! dia.
x2 ~c No. de cinturones tipo 8 faoncados al dia
I~
1
, LITRO DE i KILO DE DOCENA : REQUERIMIENTOS
TAM INA LECHE CARNE DE . DIARIOS MINIMOS
il
I
I
HUEVOS L ____ _
-A==--r 10 I 1 mg.
r-r
· - = = - - .:- --+----·----T·-----
c ! 100 10 ' 10 I so mg.
D i 10 100 -j- -10_ ___,_ _ _ _1O_m_g_ _
Soluci6n.-
Sea: x, = No. de litrcs de leche.
x c = No de kilos de carne.
x" = No. de docenas de huevos.
entonces tendremos:
53
Min. Z :_ 10x ~ 11Cx ~ SOx_
s.a.
X r XC - 1Oxh 2 1
100x, _,_ 1Ox, - 1Ox, 2 50
1Ox + 100x_ - 1Ox, 2 10
(x .. X 0 , xJ 2 0
Ejemplo 111-8.-
MAXIMA
WHISKY
CANT lOAD I COSTO POR
DISPONIBLE, I BOT ELLA
BOTELLAS/DIA I
I
A 2.000 i
$7.00
B 2.500 I $5.00
c 1.200 $4.00
54
Una cantidad ( x. ~ x. 2 ~ x. 3 ) de Blue Dot se produce y se vende
par $6.80 ( x .. - x, 2 -'- x 3
).
+ x~ 2 ) de 8 y ( x. ~· x" ) de C.
de A, ( X. 2 3 3
55
Aunque las anteriores condiciones son las basicas para poder es-
tablecer un modelo de programac1on lineal, en este curso estableceremos
otras dos que son necesarias para poder aplicar los metodos que veremos.
sin que esto signitique que no existan tecnicas que puedan resolver
problemas de programaci6n lineal que no satisfagan las condiciones
adicionales 4 y 5 (estas tecnicas solo seran mencionadas en nuestro curso).
Las condiciones adicionales son
4.- Divisibilidad.
5.- Determinismo.
111-3.1.- Proporcionalidad.
En el modele de programaci6n iineal se exige que tanto Ia funci6n
objetivo como Ia utilizaci6n de los recursos sean proporcionales al nivel de
Ia actividad Con esto queremos decir que si lo que deseamos es dL•plicar el
nivel de las actividades, bastara con duplicar las cantidades que intervienen
para el nivel unitario.
En el ejemplo 111-5, vemos que si producimos una mesa, consumiremos
3 horas hombre, mientras que si producimos 5 consumiremos 3 x 5 = 15
hrs-hm. (caso de proporcionalidad al nivel de Ia actividad en Ia utilizaci6n de
los recursos). Asimismo, ai producir una mesa ganamos $12.00. mientras
que al producir 5 ganaremos 5 x 12 = $60 00 (caso de proporcionalidad al
nivei de Ia actividad en Ia funci6n objetivo).
Se debe estar prevenido para aquellos casas en que el problema parece
tener todas las caracterfsticas de proporcionalidad, pero que no lo es. Par
ejemplo: el caso par todos conocido de un terreno que se desea cultivar; si
un campesino lo cultiva, se obtendra una cosecha de "x' toneladas, si lo
cosechan 10 campesinos se obtendra una cosecha de 1Ox" toneladas, pero
si Ia cult ivan 100 campesinos, Ia mas probable es que no se obtenga una
cosecha de "1 OOx" toneladas (rendimiento decreciente). La producci6n suele
tener una curva de Ia siguiente forma:
r CRNIIOAO
COSECHADA
I///
No. CR~PfSINOS
56
sin embargo, este problema ~imula a pnmera vtsta ser proporctonal. aun
cuando no lo es en todo su range.
Extsten ocasiones en que a proposito se asume que existe propor-
cionalidad, aunque Ia suposici6n no concuerde con Ia realidad. Estas
suposiciones son validas. siempre y cuando se calculen los efectos que dicha
suposici6n tendra en el resultado y se sepa que no lo afectara en forma
apreciable, sino que Ia aproximaci6n lograda servira a los prop6sitos del
programador .
En otras ocasiones. existen actividades que solo son lineales sabre un
:ntervalo determinado, p ej. el costa de puesta en marcha hace que el costa
de producci6n no sea linea: en cera para poder considerar estas actividades
en programact6n lineal, debemos estar seguros que el nivel de dicha ac-
tividad estara siempre dentro del intervale correcto.
c
( tI
0 si X=O
c
i
/
I
k + px si X>O I -~-
I
I
,/
/
p
• !/
111-3.2.- Aditividad.
La aditividad presupone que Ia medida total de efectividad y Ia
utilizaci6n de recursos resultantes de Ia operaci6n conjunta de las ac-
tividades, debe igualar las sumas respectivas de estas cantidades resultante
de Ia operaci6n individual de las actividades.
Vemos que s61o Ia presencia conjunta de Ia proporcionalidad y Ia
aditividad, pueden garantizar Ia linealidad de que hablamos ya con
anterioridad. La presencia de una sola de las dos condiciones, sin embargo,
no Ia garantiza.
La principal causa de Ia no-linealidad es Ia presencia de interacciones
entre las diferentes actividades, es decir, !a operaci6n de una actividad noes
independiente de Ia operaci6n de otra u otras actividades y al variar el (o los)
nivel (es) de esta (s) ultima (s), indirectamente modificamos el nivel de Ia
primera.
57
Supongamos que una campania precisa de una instalacion generadora
de vapor para poder etectuar su proceso productive. AI mismc rnpo utiliza
el vapor de salida para calentar cierto material. Suponiendo qL ' .::1 costa de
producir el vapor exclusivamente para el proceso productive fuese c,x, y que
el costa de calentar el material par un proceso independiente del vapor tuese
C2 X2 . el costa total serfa c. x, .,.. C2 X2 (existira linealidad), pero dado que el
material no es calentado por media de un proceso independiente, sino
usando el vapor de salida, el costa total de operar arnbas actividades sera:
111-3.3.- No negatividad.
M1entras que cuaiquier multiple positivo de una actividad es posible,
canudades negativas para una actividad no son r;osibles
Par ejemplo. supongamos el caso de una transportacion de cierto
producto. Aqui es imposible el transportar cantidades negativas de dicho
producto.
111-3.4.- Divisibilidad.
En intinidad de problemas, los cuales satisfacen las tres primeras
condiciones y son par lo tanto aptos de ser resueltos par programacion lineal,
nos encontramos con que las variables de decision tienen sentido
unicamente cuando taman valores enteros. Por ejemplo, asignar aviones a
diferentes rutas que deben ser cubiertas.
El procedimiento que trataremos mas adelante en el curso no conduce
(salvo en raras excepciones), a val ores enteros para las variables de decision,
raz6n por Ia cual, al aplicar este metoda de soluci6n, debemos permitir que
Ia soluci6n este expresada en valores fraccionarios.
Existen ocasiones en las cuales es absolutamente necesario el obtener
una respuesta en valores enteros (numero de edificios que debemos
construir}. Cuando este es el caso, exist en dos procedimientos para resolver
el problema:
a) Utilizando prograrnaci6n lineal de enteros. Este metoda garantiza
Ia obtenci6n del valor 6ptimo entero buscado, pero cuenta con el
inconveniente de que Ia obtenci6n de valores enteros para nuestra
soluci6n es una restricci6n dificil de manejar matematicamente.
58
el problema
X,
c
::::
4
1 - - - - - - - - - - , . ('i 3. 5)
3
Rec1o~ de /
2 so JC,8ne,.,
0
! e% CO"b
2 3
I
4
I
5
I
"'
I
7 k X.
metoda son obvias Las desventajas, sin embargo, !Iegan a ser en
extrema importantes en algunas ocasiones, ya que:
x2
6 '0 5\ :Jpt mo
- ' . entero Opt;rr:o
~ ·-- ~;; ·- ~-- (5, 4 9) no ente;o
~·~~~~~~~~-d
-~ ~-- :·:- ,. . .;, . . -------... . . -- . . . --- z
4 ······· ~-~~-~~~.":.-_,. __ .,. ::::::--::::z:
0 p 1 1rr:o
(5. 4) .-~---.z,
3 redondeado
2
0 1 2 3 4 5 6 7 8 X.
AI redondear (5, 3.5) a (5.3) 6 a (5,4) nos damos cuenta de que
ninguna de las dos es factible.
59
111-3.5.- Determinismo
Uno de los problemas mas comunes en Ia aplicaci6n practica de
programaci6n lineal es Ia dificultad de determinar el valor correcto de los
parametros del modelo (a,. b y c). los cuales dijimos con anterioridad que
deberfan ser constantes conocidas. Sin embargo, los valores que estos
parametros taman. a menudo estan afectados por eventos aleatorios im-
posibles de predecir. ya que estas coeficientes, en general. se utilizan en
madelos para seleccionar un curso de acci6n futura. raz6n par Ia cual los
mismos suelen ser predicci6n de condiciones futuras.
Es entonces necesario, para aplicar los metcxJos que veremos. el pcxJer
determinar, con razonable confiabilidad. los valores de dichas constantes.
Consideraciones.
La ecuaci6n y =ax (Ia cual es una relaci6n lineal), representa una linea
recta de pendiente ·a'.
Cualquier soluci6n (x, ax) para Ia ecuaci6n y = ax, sera un punto que
cae exactamente sabre dicha recta y cualquier pareja (x, t) que no sea una
soluci6n para esta ecuaci6n, sera un punta que no cae sabre Ia recta. Es
decir, cualquier soluci6n de Ia ecuaci6n y = ax. sera un par de numeros
(x, y), tales que para un valor dado de x (sea xc) .el valor dey (sea y 0 ) sera
a" veces mayor
60
y
v ~··· ;e ex'"
(x ax 1
(x ax)
ZONAl
~-------------------------x
Ejemplo 111·9.-
61
y
Las soluciones que satisfaces las desigualdades (1) y (2), seran las
situadas en Ia zona (a 2 ).
y
Las soluciones que satisfacen las desigualdades (1 ). (2) y (3). seran las
situadas en Ia zona (a 3 ).
y
6
0 B X
62
111-5.- METODO SIMPLEX.
Los problemas de programaci6n lineal a rnenudo cuentan con un gran
nl.lmero de vanab!es y restricciones. y ia importancia de un metoda etic1ente
para obtener una soiuci6n es de surna importancia.
Ei metoda simplex es el procedimiento utilizado para resolver
problemas de programaci6n lineal. Fue desarroliado por Jorge Dantztg en
194 7 y presentado en 1949. Este tue el primer metoda que permiti6 atacar,
ordenadamente y con un procedimiento rutinario, los problemas de progra-
maci6n lineal que antes habian sido muy dificiles de resolver.
Este metoda es un procedimiente algebraico que progresivamente se
acerca a Ia soiuc16n optima mediante un proceso iterative bien definido, hasta
que final mente esta se alcanza.
DebldO a las caracteristicas de este metoda, este es apropiado para su
utilizaci6n en las computadoras. aunque estas ultimas utilizan, por lo general,
lo que se conoce como Metoda Simplex Revisado, el cual veremos mas
adelante. El metoda Simplexes extremadamente simple y mediante este nos
es posible resolver problemas moderadamente complejos a mana. Sin
embargo, el fundamento que existe detras del metoda es mas complicado.
Definicion 111-5.-
Region de soluciones posibles es aquel conjunto de vectores que
satisfacen todas las restricciones.
Definicion 111-6.-
una soluci6n factible (o posible), es cualquier soluci6n (vector) que
satisface todas las restricciones. es decir, cualquier vector de Ia region de
soluciones posibles. Si existe una soluci6n factible para un problema, este
se dice que es factible.
63
Definicion 111-7.-
una solucion basica factible es aquella que corresponds con un punta
extrema de Ia region de soluciones posibles.
Una soluci6n basica no factible es aquella que corresponde con Ia
interseccion de dos o mas restricciones, fuera de Ia region de soluciones
factibles.
Una solucion oosica factible es el conjunto de ·m cantidades x (x ~ O)
~(X. = 0), que satisfacen el sistema de restricciones.
y de ··n-m · cantidades
b b d d
a) Soluciones factibles ( Nc oo )
b) Soluciones basicas factibles (5)
c) Soluciones no factibles ( No oo)
d) Soluciones basicas no factibles (5)
Definicion 111-8.-
Una solucion optima, es Ia mejor de todas las soluciones factibles, esto
es, aquella que maximiza Ia funcion objetivo.
Una soluci6n optima es aquel vector X0 , tal que CXc ~ CX para todas
las soluciones (vectores) X factibles.
64
1) El con junto de soluciones factibles es un conjunto convexo.
2) Si existe cuando menos una soluci6n factible. entonces tambien
existe cuando menos una soluc:6n bas1ca factible.
3) 1 .::; numero de soluciones basicas factibles < oo. Es decir. el numero
de soluciones basicas factibles es finito
4) Cuando menos una de las soiuciones 6ptimas (en caso de exist1r
varias), es una soluci6n basica factible
Debe quedar clara que una soluci6n optima no necesita ser una
solucion basica factible (es dectr un punta extrema). Esto sucede cuando
son varias las soluciones factibles que maximizan Ia funcion objetivo. o sea,
existe un segmento iineal de soluciones que Ia maximizan yen este segment a
lineal solo los puntas extremos son las soluciones basicas factibles, las
demas son soluciones factibles pero no basicas. La propiedad 4 dice que
cuando menos una de las soluciones optimas sera una soluci6n basica
factible, pero no dice nada sabre !as demas (en caso de que existan).
~
I n"-- =
65
La propiedad 3 surge del hecho de que tambien es finito el numero de
las restricciones y par lo tanto el numero de sus intersecciones o puntas
extremos.
La propiedad 4 provee Ia base fundamental del Metoda Simplex, ya que
nos indica que solo un numero finito de soluciones. las basicas factibles (de
entre el numero infinite de soluciones factibles) necesitan ser investigadas,
con objeto de llegar a localizar una solucion optima. De a qui que una solucion
optima pueda ser siempre encontrada examinando unicamente cada una de
las soluciones basicas factibles y eligiendo aquella que de un valor mayor a
Z.
A pesar de que las soluciones basi cas factibles exist en en numero fin ito,
en Ia mayorfa de los casos su numero es muy grande, razon porIa cual resulta
poco eficiente buscar entre todas elias hasta encontrar Ia soluci6n optima.
Es aqui donde entra en acci6n el Metoda Simplex, el cual, ademas de
que solo examina los puntas extremos de Ia region de soluciones factibles,
realiza su busqueda de una manera optima, ya que no exam ina a todos elias,
sino que parte de una soluci6n inicial, a partir de Ia cual busca otra mejor y
asi sucesivamente hasta encontrar el optima.
Ejemplo 111-10.-
66
lineales que constituyen ias restricciones resulta en el conjunto convexo de
puntas dado por el poiigono OABCD. Cualquier punta (x. y) dentro de este
poi1gono (region a) satisface el sistema de desigualdades ( ll AI polfgono
OABCD ie ilamaremos region de so!uciones factibles
X,>
X
60C
5CC 01
0
0
40C
X, = 300
3CO
200
100
n
X.
'-
'·" X1 + 2X 2 = Z
~-
67
El problema puede ser pensado como el de determinar aquella linea.
de entre Ia familia de rectas x. - 2x" = Z. que esta mas lejos del origen
pero que aun contiene cuando menos un punta dentro. o en Ia frontera del
polfgono OABCD
Como primer intento. hagamos Z = 100 y veamos si existen val ores
(x., X2 ) dentro de Ia region de so!uciones posibles que nos den ese valor de
Z Vemos que cualquier punta de Ia recta dentro de Ia region de soluciones
factibles nos da este valor de Z y que es aun posible alejarnos del orfgen sin
salirnos de nuestro polfgono.
Probamos con Z = 600, nuevamente un segmento de Ia linea x. _,. 2x 2
= 600 cae en Ia region a, par Ia que el valor maximo permisible para Z sera
cuando menos 600. AUn podemos alejarnos mas del origen.
Final mente, Ia recta 1000 = X1 + 2x2 sera Ia que satisfaga Ia condicion
de ser Ia mas alejada del origen, conteniendo cuando menos un punta en Ia
frontera de a. Este punta sera l;::t solucion optima buscada
(x* 1, x* 2 ) = (400, 300)
X,
600
500
68
2) Examina las aristas del vMice. para ver si al moverse por una de
elias. hasta el siguiente vertice adyacente. se aumenta e! valor de
Z. Si el recorrido sabre las aristas no aumenta el valor de Z. el vert1ce
en el cual estamos situados maximiza Z. Sial recorrer al menos una
arista se aumenta ei valor de Z. se pasa al paso tres.
3) Se escoge una de las aristas a lo largo de las cuales aumenta el
valor de Z y se sigue sabre ella hasta alcanzar el siguiente vertice
adyacente.
4) Se repiten los pasos 2 y 3 hasta que el valor de Z ya no pueda
aumentarse.
x, + ~
x ~ 0
= 600} ()
II
3
69
Similarmente:
X 2 :::; 300 se reemplaza por x2 .,_ x j = 300
X4 2: 0
y
x. + x3 = 600 }
X2 + x. =300
(a)
3x, + 4x 2 + x5 = 2400
X1 2:0; X2 2:0; X3 2:0; X4 2:0; X5 2:0
Definicion 111-9.-
La soluci6n obtenida al resolver el sistema de ecuaciones por m
variables en terminos de las restantes (n- m) variables, asignandoles un valor
de cera a estas ultimas, se conoce como una soluci6n basica. Si el vaior de
cada una de las m variables es :::::0, tendremos una soluci6n basica factible.
70
Si ei valor de cada una de las m variables es > 0, tendremos una soluci6n
basica factib1e no degenerada,
Las m variables escogidas se denominan como basicas o como las
variables en Ia base, Las (n m) variables restantes se conocen como
<
variables no basicas,
Ejemplo 111-11.-
La soluci6n:
x, = 600
x2 = 300
X5 =-DOO
es una soluci6n basica no factible para ei modelo de nuestro ejemplo anterior.
La soluci6n:
X, = 600
x2 = 150
X4 = 150
es una soluci6n basica factible no degenerada para el mismo ejemplo.
En forma semejante podemos obtener el valor de las restantes B
soiuciones basicas.
El Metoda Simplex no puede iniciarse mas que con un sistema cuyas
ecuaciones estan siempre en forma can6nica
Paso 1.-
J ustificaci6n.
71
variables para que formen Ia base. El inconveniente de este metoda. es el de
que entonces el sistema de ecuaciones ya no esta en forma canonica para
las variables en Ia base. razon par Ia cual debemos de efectuar dicha
transformacion para poder tener Ia solucion deseada (cosa que no tue
necesaria en el metoda anterior) exist1endo Ia posibiildad de que. despues
de haber ejecutado el trabajo de transformaci6n para tener Ia soluci6n
deseada. descubramos que dicha soluci6n no es factible (alguna b. < O) y
que sea necesario hacer un nuevo intento.
Las variables de holgura se usan para encontrar una solucion ba.sica
factible inicial rapidamente y por lo tanto no son indispensables en aquellos
casas en que tengamos ecuaciones o exista una soluci6n basica factible que
sea obvia para nuestras ecuaciones.
En general, Ia soluci6n basica factible inicial para el Metoda Simplex
sera siempre
x = 0 para j = 1 ,2, .. n- m
xn-m~· = b parai = 1,2, .... m
y Ia soluci6n oosica factible en Ia etapa testa relacionada con Ia soluci6n basica
factible en Ia etapa t + 1 de Ia siguiente manera: Una de las variables no oosicas
en Ia etapa t, toma un valor diferente en Ia etapa t + 1 y se denomina variable
de entrada' . Para compensar. una de las variables basicas en Ia etapa t se hace
cero en Ia etapa t -r- 1 y se denomina "variable de salida".
Las otras variables no basicas de valor cero se mantienen en cero; las
otras variables basicas no nulas, en general, siguen siendo diferentes de cero
(aunque su valor puede variar). Lo anterior es equivalente a pasar a un vertice
adyacente de nuestro conjunto de soluciones posibles. que incrementa el
valor de Z.
Paso 2.-
Detennine Ia nueva variable oosica entrante (xe ), Ia cual debera ser aquella
varwe no ba.sica, en Ia tunci6n objetivo actual, que satisfaga Ia condici6n:
C' e =max C' ! > 0
cuando dicha tunci6n se encuentra en Ia forma:
Z = C', x, + C' 2 x2 + ... -t- C' c xn ... (I}
Justifjcaci6n.
72
Es importante recordar siempre que x. es Ia variable nula, en Ia soluci6n
basica factible actual, que va a tamar un valor no nulo en Ia siguiente iteraci6n.
La regia de selecci6n para x. tambiem puede ser expresada asf. x. sera
aquella variable no oosica en Z que cuente con el mayor coeficiente entre
aquellos con signa positivo. estando ZenIa forma indicada en (1). Esto se debe
a que sic·. = max C' >. 0, al incrementar x., el valor de Z parece 1ncrementarse
mas rapidamente que incrementando cualquier otra variable.
El motivo de introducir x. a Ia base, es que si cuando menos existe un
coeticiente C' que sea positivo en Ia ecuaci6n (1), entonces es posible
(suponiendo que todas las b, > 0) el construir una nueva soluci6n basica
factible con un valor mayor para Z. Esta nueva soluci6n puede ser obtenida
incrementando el valor de una de las variables no basicas (x.) y ajustando
los valores de las variables basicas de acuerdo con este cambia.
Existen otros metodos para seleccionar x. (selecci6n arbitraria, proban-
do el efecto producido en el valor de Z par cada variable), pero Ia regia que
utiliza el Metoda Simplex, ha demostrado. en Ia practica, que casi siempre
requiere de menos iteraciones para llegar a Ia soluci6n optima.
Debe notarse que aunque asi lo parezca, Ia variable x. no necesaria-
mente es Ia que mas incrementa el valor de Z, debido a que las restricciones
pueden impedir que su valor aumente tanto como con variables con coeti-
cientes menores, como puede apreciarse en el siguiente ejemplo:
Ejemplo 111-12.-
Z = 3x 1 -r X2
x, :s 1
x2 :s 100
Aqui se selecciona x. = x, y sin embargo noes Ia variable que mas incrementa
el valor de Z.
Paso 3.-
73
Justificaci6n.
xbm = b''"-a''"'l.x.
Z = Z' 0 + c'. x.
Del sistema de ecuaciones (II), vemos que ellimite superior para el valor
que podra tamar x. en Ia ecuaci6n (i) sera:
<1> = 00
x.:-s <t>, b
{
a,.
sea:
<t> = min <t> I
74
PROGRAMACION L!NEAL
ape
=a:,.min>O '
aie
<!: 0
Paso 4.-
75
a) AI no contener Ia funckm objetivo una de ias variables no basicas.
no podemos juzgar el efecto que Ia misma tendrfa sabre ei valor de
Z al incrementar su valor.
b) Este hecho nos impide observar si alguna de las variables no
basicas pudiera incrementar el valor de Z, pues los valores C'.
cambian del ciclo tal t + 1, ademas de que algunos cambian de ser
positives a negatives y viceversa.
c) Este punta nos hace que no podamos determinar el efecto de
introducir una variable no basica a Ia base en el siguiente ciclo, ya
que este hecho puede afectar el valor de x. (en el ciclo t + 1). Ia cual
es ahara basica, con lo cual es imposible predecir e! efecto neto
sabre Z, ya que, en general, los valores de las variables basicas son
afectadas par los cambios en Ia base.
Paso 5.-
1a. lteraci6n.
Paso 1.-
76
pendientes (variables no basicas) segun las ecuaciones (II) tenemos
x, - 600-- X.
X~ 300 -- X,
x5 -- 2400 - 3x. - 4x 2
dado que a Ia variables no basicas les hemos as1gnado un valor cero.
obtenemos Ia siguiente soluci6n basica inicial:
variables basicas variables no basicas
--x --;; 600 - ---
X.= 0
-- ---
3
X4 = 300 X~= 0
X0 =2400
"
z 0
o en forma vectorial:
(x .. X 2 , X 3 , X4 . X 5 ) = (0, 0, 600. 300, 2400) Z=0
La anterior es una soluci6n basica tactible no degenerada que co-
rresponde al punta extrema (x. , X2 ) = (0 , 0)
Automaticamente hemos satisfecho el paso 4 y debido a que existen
c, > 0, el paso 5 nos indica seguir con el paso 2.
Paso 2.-
Paso 3.-
77
nunca se volvera negativa ( x, sigue siendo cera). (rp. = x)
De Ia ecuac1on (2). vemos que lo maximo que podemos incremerttar
x 2 es 300, antes de que X4 se vuelva negativa. (¢' 2 = 300)
La ecuaci6n (3) nos indica que lo maximo que podemos incrementar
X2 es 600, antes de que X5 se vuelva negativa. (¢ 3 = 600)
x2 = 300 x4 == 0
x 5 =1200
z :::: 600
78
ansta OA hasta el vert1ce A (O . 300)
A B
~
aqu esta:nos ahora
0 c
+----+--+--+---+--+--+---+-~ X.
Paso 5.-
2a. lteraci6n:
Paso 3.-
79
formamos dicho sistema de ecuaciones. pero en esta segunda iteraci6n Ia
haremos directamente de nuestras ecuaciones que es Ia forma en que
normalmente se hace.
De Ia ecuaci6n (1 '):
b1 600
, = -,-
a,. = - 1 = 600
rp' ~ x, $ 600
De Ia ecuaci6n (2 ')
a' 21 $ 0 ~ ¢ 2 = a: es decir: x, $ cx:
De Ia ecuaci6n (3 '):
b3 1200
¢3 = --;---- = - - = 400 ~ X
1
$ 400
a 31 3
basicas no basicas
Paso 4.-
Obtengamos ahara una forma can6nica para nuestra nueva base con
ayuda del metoda Gauss-Jordan:
80
Ordenando las variables
z ~
2 3x~ ~
1 '3X 5 1000 (0 ')
4/3X 4 - 1 ,3x5 ~
x2 = 300
X3 = 200
z =1000
Geometricamente, lo que hemos hecho es movernos a lo largo de Ia
ansta AS hasta el vertice B = (400.300).
X?
A B
Paso 5.-
81
Z' = 1000
x·. -= 400
x' 2 = 300
x" 3 = 200
valores que coinciden con los obtenidos por el metoda grafico.
Maximizar Z = CX
Sujeto a:
82
Esto se efectua naciendo usa de una tabla como !a que s1gue:
p P,,- p p p 0
p
-~~--~---~- ---~----~---~ -·
a .. a.2 a_ a. a a.
a __' a.
a21 a2? <'0 a:::'"' a2
"
a al2 a<3 a a a
a a,n a
M~ arT'\;;: arr-J I I
~- "~
No. X, ~ • ~ X4 ><:, P8 ¢
t I
v.b Ec Z P, P2 P3 P4 P5 , j
----+-------------------t--- ------+--~__,..._..~---;------~--------+--~--·------1
z 0 1 I -1 i -2 ' 0 ; 0 i
I
0 I
'
0 i
I
;
x3 0 1 0 0 ! 0 ! 600 ! ex
x4 2 0 0 0 1 i 0 300 300
x5 3 0 3 4 0 0 2400 600
83
Paso 4.- En una nueva tabla. cambtese x" par x 2 en !a columna v.b. y pongase
ia ecuaci6n i = p (ecuaci6n correspondiente a xJ dividida entre aPe (es decir
entre el coeficiente de nuestro pnmer pivotel en Ia nueva tabla En nuestro
caso aPe= 1. por lo que obtendremos
v.b No
Ec
z X. x2 X.,,.) x4 xs
p
.c
,..-----+----+--
z 0
x3 1 i
I
x2 I 2 0 0 0 0 300 I
'
xs 3
. v b I No.
1
__
z
_j -=.E:..:cc:_.,_'--+
z 0 1 i -1 0 0 2 nva(O) =antl0)+2ant(2)
1 0 0 i 0 0 600 : 600 . nva( 1 J =ant( 1)
2 0 0 1 0 1 I 0:300, rx 1
1.
v.b I No. I Z I
p :
x3 [ x4 i x5 i o :
1 1 • '
Ec. I x, x2
I --t-==+--+--+--i----r-:-.:-+---4----j
! z 1I o I 1
f-1
84
corno nuestra ecuacion (0) no tiene ya ningt.'m coeticiente negat1vo. nuestra
soluci6n sera 6ptin:a.
variables basicas variables no basicas
x:'~ · 4oo X:=
0
x: =<'
300 x;
= 0
x.~= 200
Z* =1000
Ejercicio 111-1.-
Sujeto a
X, :s 4
x< :s 3
X. +2 x_< :s 8
X, ~ 0
x2 ~ 0
v.b No. i Z : l i i :
. • Ec. i ! x, x2 I x3 : x4 x5 ; p o i ¢; i variables basicas no basicas
r------;- --;~,--+-,____,~-+---r-----t--~~
x = 4 - ·- ·x:~,;- 0
' Z 1 0 • 1 -2 -5 I 0 0 0 0 i 3
x3 0 1 0 : 1 0 0 4 rxl x. = 3
X
4
0 0 ( 1 I 0 1 0 3 I 3 i x5 = 8
0 ' ;;(i Q 0 4
I 1
•
:
I 8 1,
z= 0
. I I
· v.b i No.I Z
1Ec.' variables basicas no basicas
X = 3· -·
0 0 15 II X1 = 0
I~
1 -2 0 2
0 1 0 ! 1 0 4 I 4 I x3 = 4 x. =0
0 0 0 i 1 0 3 ex ! x5 == 2
0 i 0 0 I -2 2 2 I
z = 15
X 0 =X,
X s =X5
85
v.b
,,ariables basicas no basicas
z - x* =- 2- -- -- ~ --x* o·
X:=
L
3
2 0 0 3 x"' = 2
3 -2 2
Z* = 19
Ejercicio 111-2.-
Sujeto a:
- x, x2 ::; 2
x, + 2JS ::; 10
2x, - lS ::; iO
(x,, JS) ~ 0
1 0 0 0 1 i-1
/513'5 i6 ~= 2 x; = 0
2 0 0 0 I 2/5 j-1'5 I 2 x7" = 6
3 0 0 0 i 1/5 1 2;5 ! 6
z· = 22
86
z
2X - ~ = lB
' ' '
6, 2)
-.!<I + 2X,= 18
K
-2 -1 1 ~ 6 R g 1R
111-12.- COMPLICACIONES
En muchas ocasiones, nos encontramos que al construir nuestro
modelo matematico, este aparece con algunas variaciones respecto al
modelo general que hemos analizado hasta el momenta.
A continuaci6n veremos cuales son las variaciones mas comunmente
encontradas y las trataremos individualmente, sin que ello implique que no
se pueden presentar varias de elias simultaneamente, en cuyo caso se
aplican, en forma sucesiva. los metodos que se describiran para solucionar
estas variaciones.
Las variaciones mas ccmunes son:
1 - Minimizaci6n.
2.- Desigualdad con sentido invertido.
3.- Constantes negativas (b, < 0).
4.- lgualdades.
5.- Variables no restringidas en signa.
6.- Empate para entrar como variables basicas.
7.- Empate para dejar Ia base (degeneraci6n).
8.- Soluciones multiples
9.- Ausencia de soluciones posibies.
10.- Soluci6n optima sin limite.
87
En todos ios casas anteriores aun es posible utilizar eln1etodo simplex.
solo que sera necesario efectuar ligeros aJustes en nuestro rnodelo. Los
ajustes necesarios se trataran en cada caso particular.
111.-12.1.- Minimizaci6n
Ejemplo 111.-13.-
88
s.a
'X. < 8
3x, T
2x" :5 30
-X. T 3x~<' :5 12
(x. , x) ~ 0
para resolverlo, multipl!camos Ia func16n objetivo por (-1 ). con lo cual
obtenemos:
Max -~Z = x. ~ 2X 2
s.a
X, :5 8
3x. + 2x 2 :5 30
..._ 12
-X, 3xry :5
(x,, x 2 ) ~ 0
y lo resolvemos a.p:; . :<:mdo oi rr;;~todo general vista con anterioridad:
v b No :
. Ec -Z ' ~. v.b. vn...b~-
I z I o _, ' -2 : 0 c 0 ' 0 ' -z = o x. = 0
' X I 1 0 0 1 : 0 o 8 I o: x3 = 8 X2 = 0
: 3 I
x4 2 0 3 2 I 0 1 0 : 30 i 15 X4 = 30
xo 3 0 3 0 0 12 i 4 x5 = 12
I
: v b! No i .
' !
i
• I E , -Z 1 X, 1 x2 1 x5 P0 : rf;, vb. v.n.b.
I C' -~~- _-z-;;.:-8
z
1
o , -5/3 I o 0 i 0 2,'3 i 8 i X1 = 0
x3 = 8
1
1
3/111-2/1116 x, = 6
!
1/11 3/11 1 6 x2 = 6
89
.____ ____ ( 5' 5 J
__.------------;~--.
. ......
" ·, ~~-.
'"I ( s. 3)
X,
4
Ejemplo 111.-14.-
x1 + Xo< :s 9
(x,, x2) 2:: 0
Soluci6n:
'-,,X,+ 2~,= 12
'"- ..
X.
90
para poderio resolv£~r. multipl1carnos por (-/) Ia tuncion objetivo, y ob-
tenemos
Max Z -~ x. - Bx~
sa
-x ~ 4x~ ~ ~ 2
x. -+- 2x 2 ~ 12
x. ._ x2 ~ 9
(x. , x 2 ) 2:: 0
v.b No -Z
~----Ec • x. Xc X3 X4 Xc l pC ! rp vb. V.c.O~b.
j
r------ t- --··-:---r------t--;
z 0 -1 -8 0 0 0
1
0 I --Z = 0 X,= 0
x3 0 -1 4 j ' 0 0 ' 12 ' 3 x3 = 12 x2 = 0
I X 2 0 2 0 0 12 ' 6 X4 = 12
. 4 I
I x5 I 3 0 0 0 1 9 9
i I
x5 = 9
v.b -No I -Z
Ec x• x2 - x3 x4 · x5 pc cp v.b. v.n.b.
~-=-:;_+-----r-·-r- ·-·~---'·-~---- --· ::::z= 24 ·
Z o 1 -3 I o 2 · o o 24 X,= 0
o -1-4 I I 1•4: 0 0 3 x2 =3 x3 = 0
0 3/2 j 0 1 1 o 6 i 4 x. = 6
0 5l4: o -1/4: 0 6 ! 24!51 x5 = 6
j I
I v.b! No 7 ___yj), v.n.b.
i I r::c . x1 X. ;(3 x. ' x~, Po
x;
-L <
L----~- ___.. ___ - I ..J -Z*= 36 - 0
i z I
0 1 0 0 2 0 36
x2
I
I 1 0 0 1 1/6: 1/6 1 0
!
4
i x; =4 x*4 = 0
I i
x*
:
0 I:-1131 =4
i x, i
2 0 . I 213 'I 0 I
4
x*5 =
I x5 I 3 0 0 0 ! 1/6 i-516[ 1 I
91
',:-i2.- C8MPL.C~CIC~ES
Ejemplo 111-15.-
Ejemplo 111-16.-
7 > 5
-7 < -5
b, =-130
b2 "CO -5
Ejemplo 111-17.-
92
Max z ~~
X. - 2X 0
s.a
X. < 600
X < 300
3x. - 4x0 ~ :-2400
(x. x _ :l _> {j
X,.
0
0 t----------6~0f-O---- X.
-
x,
3x~
~
- 4x 2
-t- xc
~•~
x.
-t- xs
J(
=
~
600
300
-2400
0
} (I)
93
basica factibie iniciai obvia (pues x, < 0), par lo que deberemos de corregir
esta nueva complicaci6n segun el metoda que se explica mas adelante.
Debe notarse que el hecho de que hayamos puesto un ejemplo en el
cual al corregir Ia complicaci6n 2 caemos en Ia 3. no significa que este sea
el caso siempre, pues es posible encontrar problemas en donde no suceda
esto. Lo unico que nos enseria el ejemplo es que es ei caso mas facil de
encontrar.
Saul I. Gauss ha presentado un metoda para hacer minimo el numero
de variables negativas, en un conjunto de restricciones con varias desiguai-
dades del tipo 2:.
Veremos el metoda par media de un ejemplo
Ejemplo 111-18.-
Max ·z = 3x, + x2
s.a.
x. + lS 2: 4
2x, - x2 2: 6
X, 2: 4
4x 1 + 3JS ::s 24
(x, , JS) 2: 0
94
lntroduciendo las vanab!es de holgura para formar un sistema de
ecuaciones. ootenemos:
z - 3x, -
x2 0 ... .. (0)
X. - X. -
XJ 4 (1)
2x. -
x2 X~ 6 .. .(2)
X, - xs 4 .. .. ... (3)
4x, - 3x 2 -'-
xs 24 ... (4)
4x. +
3~ + x., = 24 ... (4)
95
basica no factible. con el consiguiente trabaJO extra).
b) lntroducir una variable artificial (una variable artificial es diferente a
una variable de holgura), !a cual entra restandose en aquellas
ecuaciones en donde existe una constante negativa y usar esta
variable artificial como Ia variable basica inicial para esta ecuaci6n.
Si tenemos:
a ,x, + a,2 x2 T- ~ anxn :S - b (b > 0)
introducimos primero una variable de holgura para eliminar Ia
desigualdad y a continuaci6n Ia variable artificial mencionada, con
Ia cual obtenemos:
a,x, -r a2x2 ~ + a,xc -r x,-xa ""-b,
cambiando signos:
-a,x,-a 2 ~- .. -a, X,, -x, -xa = +-b
0
Ejemplo 111-19.-
hacienda Xa = X6
- 3x,- 4~ + x5 - X6 = --2400
cambiando signos:
3x, + 4~ - x5 -r x6 = 2400
debido a que todas las demas ecuaciones se mantienen sin cambios, Ia
soluci6n basica factible inicial sera ahara:
basicas_~ _no_Msicas
Z = 0 X.= 0
X:J= 600 x2= 0
X4 = 300 x5 = 0
X5 = 2400
El problema que surge al aplicar este metoda, es el de haber modificado
nuestro problema original al agregar una variable de holgura y una variable
96
artificial con signos contraries. ya que x., - x_ puede tamar cualquier valor
arbitrario (desde- ,, hasta cc) con !o cual. el resultado neto es Ia eliminacion
operativa de Ia restricci6n origmal en donde surgio este cambia. Esto trae
consigo un agrandamiento del conjunto de soluciones posibles como es
clara en el siguiente ejemplo:
Ejemplo 111-20.-
Sea
3x, ,-4x 2 ~2400 ........ (1)
Ia cual se transforma a:
3x, + 4X 2 - x5 +- xe == 2400 ..... (2)
si hacemos:
3x, - 4x2 -= 1000
habremos ya roto Ia restriccion ( 1), sin embargo si hacemos
xs - x., = 1400
cumpliremos Ia ecuacion (2).
De lo anterior, es clara que es posible cumplir con Ia ecuaci6n
modificada sin satisfacer Ia restriccion original. Graficamente tendremos:
,
'
''
'
A 3
c~,----------------0~-----x. 1-----x.
97
Existe Ia posibilidad de que I ' solucion optima para ei problema
revisado no sea factible para el origmal {p.e el punta A con una funcion
objetivo como Ia recta II), en cuyo caso, al solucionar Ia complicacion 3.
habremos hecho que Ia solucion optima obtenida no corresponda a nuestro
problema original. Es por esto que el metoda simplex introduce un
procedimiento que tuerza el valor de las variables artificiales a cero, par lo
que nuevamente volvemos al problema original y Ia solucion optima que se
obtenga de esta forma. sf sera Ia que corresponde a este ultimo
El procedimiento que emplea el Metoda Simplex, es el asignar una
penalidad muy grande a las soiuciones factibles del problema revisado que
no coincidan con las del problema original (soluciones en Ia region
a 2 - a J La forma de asignar dicha penalidad es introduciendo -Mxa a Ia
funcion objetivo, estando esta en Ia forma Z ~ c,x . ..,.. ... + cnxn, de tal forma
que si xa ~ 0, el valor de Z disminuye enormemente ( M representa un numero
muy grande).
b) Alguna (s} variable (s) aparece (n) en Ia base optima del problema
modificado con valm mayor que cera. Cuando esta situacion se
presenta, el problema original no t1ene solucion (no tiene
soluciones tactibles}.
c) Alguna (s} variable (s) aparece (n) en Ia base optima del problema
modificado con el valor igual a cera (soluci6n degenerada). Esta
soluci6n tambien corresponde a Ia soluci6n optima del problema
original.
Ejemplo 111~21.~
98
z - X 2x Mx" --
0 (0)
X X, 600 ( 1)
~
x_ -- X" 300 ... (2)
"
3x. -r 4x~ -X, XC -~ 2400 (3)
i z 0 1 -(3NI~1\i-(4M.,-2! 0 : 0 ! M 0 -2400M
x3 1 0 1 0 0 0 0 600 rxi
x4 I 2 0 0 0 I I 0 0 300 3ool
x6 3 0 3 4 o / o /-1 I 1 2400 5oo I
v.b v.n.b.
-z -= -24ooM --x, = o
X3 = 600
300
X6 = 2400
99
v.o Ec
'No
z X. x2 I x, x, xs
p
rp v.b v.n.b.
~-+-- -
x;-:-~ 0
-----;--
: z
I
I 0 II 1 c 0 1 2 0 M I
1200 Z* = 1200
x5 I
I
1 0 0 0 3 4 ·-1 600 x* = 600
"
x: 0 =
x2 2
!
Ejemplo 111-22.-
Max Z = 3x, - x~
s.a.
2x, + ~
;:>: 2 ·--.t
x,:
-~-T
..
x, + 3~ ::::: 3
) ~·
L z' z'
x2 ~ 4
(x, x) ;:>: 0 2
.r--r \.--
X,
e 3 '4. 5
J x~
,
.
~
2X, X, 5
100
Quitanoo las desigualdades e introduciendo Ia variable artificial. ob-
tendremos·
s.a
2x. ~ x_< X.0 + X" 2
X. .,_ 3x2 ... x4 3
x< + xs 4
X ~ 0 u= 1,2,. 6)
En forma de tablero:
I·, b I Ec.
·No.
z x,
i x2 I x3 I x. xs x6 : Po i
---r----
z :
I
0 1 I ~3 1 ' 0
I
0 0 M 0
x I I 0 2 1 i -1 0 0 2
6
x4 2 0 3 0 0 0 I 3
xs ; 3 0 0 I 0 0 0 4 i
i I I !
1 v.b i Ec. Z . 1 p "'
. 'N ' x, x2 i x3 I. x4 xs i xo I o 't'l II
'-~' ...E:..L-+------+------+- --+1-------+1---4·---·
l i
I
0 1 0 5/2 i-31210;-0
j I
(M+3/2)
I' 3
1
1.,
I
x, I 1 0 ' 1 1/2 ,... 1/2 i 0 0 I 1/2 i 1 oc I
x4 i 2 0 0 5/2 f~1;2! 1 0 I -1/2 I 2 4
xs I 3 0 0 I 0 ' c I 0 I 4 I iX
ivbEc(z x, ; x2 i x3 x., xs x6 pa
1 No.I -i-----+- - ~ -+-----r- ......,....1
~ z· o i 1 0 1 10 0 3 I 0 M 9
X1 • 1 I 0 I 3 0 1 I 0 0 3
X •
3
2 0 0 5 1 2 0 -1 4
i x5 • 3 1 0 ; 0 0 0 0 4
101
111-12.4.- lgualdades en las Restricciones
102
cor' lo cual es clara que utilizando el metoda (c) tambien (al igual
que ei metoda (bl del caso 3) modificaremos el problema original.
por !o que es necesario volver a el. obligando a las variables
artificiales a tener un valor cero, para lo cua! se emplea el metoda
de 1a Gran M visto en ia compltcacion lli-12.3
Ejemplo 111.23.-
reg•or de
SOiUCiO'l8S
'-'----+-·~ 'actroles
Q ~--------------~~------x.
600
103
aplicando el metoda (C) tendremos
X,
z -- X. - 2x 2 0
x, X, 600 -C~1V
32S :inJ aue~.k:
X + X 300 ---
3x. + 4K. - X, - 2400 J!!•
J
z - x, - 2x 2 + Mx5 0
X, +- x3 600
x2 + x4 300
3x, - 4x 2 - )(5 2400
v.b v.n.b.
Z = -2400M x, =0
X3 = 600 X.. = 0
<
x. = 300
xs = 2400
104
\ b v no
Z c.= -600r2M-1) x. = 0
X2 = 600
X" =:: 0
X~e = 300
v b.,.___ _~v._n-=b
Z* = 1000 =0 x;
x; = 200 = o x;
x; = 300
x:" = 400
Ejempto 111-24.-
JS - X4 + x5 = 2 ........ (3)
4x. · 3JS + x 6 = 24 ........ (4)
·r;,·
x, - x2 2
105
formando ahara un s1stema can6n1co:
· v. b ~ Ec
_ _.,
, No.'
,_____,___
z
X, x2 x3 x4 x5 ' x5 I X' P, ¢
-+-~--~
z 0 1 ~(2M+3) (M-1) 0 M 0 0 !
0 -BM
' x3 1 0 ! 1 -2 -1 0 0 0 2 2
i i
X? I 2 0 2 -1 0 -1 0 0 6 3
I
I
xs 3 ' 0 -1 0 -1 1 0 0 2 2
x6 I
4 I 0 i 4 3 0 0 0 1 I 0 24 6
I v b Ec I
I
I
No. I
z x, )(2 ! x3 i x4 x5 )(6 x7 Po ¢
I
I
iI-2M+61'
~(M+4) I
i '--,---
II z 0 ! 1 0 0 -(M+3) I (2M+3) I 0 0
i x3 !0 0 ' -1 ! 1 i 0 i -1 ! 0 0 0 a
i 2 I 0 0 : 1 ! 0 1 -2
x7 \ :
:
! 0 2 2
)(1 3 I 0 i 1 -=-1 0 -1 0 0 2 c:
4 II 0
I i
! x6 0 7 0 4 --4 0 16 16(7
i
I v.b Ec. I i : i
No.
z I
! x, x2 ' x3 x4 I C/J, x5 x~ , po
z 0 1 --~--+--~-+-0--+---+-=-:--+1-~ I(M+--4-)t-~-~-4--tl-a---1
X"
~
1 0
~ ~ ~ =~ ~ ~
2 0
x2
x,
X6
3
4
0
0 o
o
o o -3 1o 1
I
I -7
!2
!
! , ;s
106
vo Ec
z X. X. X X
,--~------~~ --....-------t------.---~---T-------~--~·-~-
z 0 1 0 0 0 -1 2 0 1 '2 I M+ 1 2 15
X-: c (j 0 1 ·o 0 3 10 -; 1 10 n5 26
x2 2 0 C! 0 ;~ 5 0 1,'5 ~2/5 12:5 6
< 3 0 1 0 I] -3 10 0 1 10 3 10 21,5 I u
)C 4 0 0 0 0 -3. 10 1 10 ., 10 15 u
"
vb Ec ' p
·No ~ ' 0
f--------· ~~-----~__....--~-..,......._.__ _____,. ___ ~_ -·~---·--~-~
z· 0 1 0 5,'4 0 0 0 3;4 M 18
I
iI XJ * 0 0 -1:4 1 0 0 1 ·4 -1 2
X. .
x4 * 2
3
0
0
0 . 5,2
34
i 0
0 0
0
0
1/2
1.A
-1
0
6
6
x5 *' 4 0 0 3!4 0 I 0 1.4 -1 2
107
Ejemplo 111-25.-
problema que coincide ahara con nuestro modelo general, porno contar ya
con ninguna complicaci6n. por lo cual puede ser resuelto utilizando directa-
mente Ia rutina simplex, como lo hicimos con el problema original
Si:
Z = 5x, _._ SJS + 2x3
entonces x. y lS estaran empatadas para entrar en Ia base. pero si tenemos
que:
Z = Sx, + SJS - 7x 3
entonces x. = x 3 y no existira empate.
111-12.7 .· Degeneraci6n
Cuando al aplicar el paso 3 de ia rutina simpiex nos encontramos con
108
que 'P = rp, == 9 . -~ esto es. con que varias variables se hacen cera
simultaneamente cuando xe se incrementa. tenemos un empate para dejar
Ia base (degeneraci6n).
Cuando es este el case. se escoge arbitrariamente cualquiera de las
variables en empate para ser Ia variable de salida.
El efecto de tener un empate en ia variable de salida. es que todas las
variables en el empate que no sean escogidas para dejar Ia base quedaran
con un valor de cera en Ia siguiente iteraci6n, durante Ia cual tendremos que
¢, = 0 (para t = aquelllas i's correspondientes a las ecuaciones conteniendo
las variables en em pate); con Ia cual las variables en em pate continuaran
siendo candidatas para dejar Ia base en esta nueva iteraci6n, hacienda
ademas que Ia nueva variable entrante no pueda tamar un valor mayor que
cera, con Ia que el valor de Z permanecera sin cambia respecto a Ia iteraci6n
anterior.
De io anterior. cabe senalar que si existe degeneraci6n, es posible
tener una secuencia de iteraciones que no aumenten en nada el valor de Z,
y ha sido probado par Hoffman (1951) y par Beale (1955} que Ia rutina
simplex. en este caso. puede conducir a una recurrencia del mismo conjunto
de variables basicas despues de k iteraciones (y despues de 2k iteraciones.
etc. indefinidamente}. De aqui que sea evidente que ya noes posible el seguir
manteniendo que el metoda simplex necesariamente terminara en un
numero finito de iteraciones. A pesar de los ejemplos construidos par
Hoffman y por Beale, el caso de cfrculos viciosos mencionado no ocurre en
Ia practica. raz6n par Ia cual las reglas que existen para resolver el case de
degeneraci6n son casi siempre ignoradas. pues at:m en el caso de volver al
mismo conjunto de variables basicas, se pueden salir del circulo vicioso
cambiando Ia selecci6n de xs con respecto al ciclo anterior.
Ejemplo 111-27.-
z -x, - 2~ 0 rp
x, + x3 600 a
109
donde:
pero ahara:
¢ = rp 1 = ¢2 = 300
Escogemos x4 arbitrariamente
basicas
z ~x, + 2x 4 600 - - · no basicas
X2 =300 X.=O
~l
X, x3 600
X3 =600 X4 =0
x2 + x4 300
x-
5- 0
3x, - 4X 4 + x5 0
z =600
y obtenemos una soluci6n degenerada (pues X5 =0), con Ia que seve que
Ia variable en empate que no result6 escogida para dejar Ia base qued6 con
un valor de cera.
Como Ia soluci6n anterior no es aun optima, hacemos:
x. = x,
¢ = ¢3 = 0
y obtenemos:
basi(2as _ . _ng_b_?_~c_a~
Z 2/3x 4 + 1 /3x 5 = 600 x7= o x;=o
~
+ ~
+ x3 + 4/3x 4
+ X4
1/3x 5 =
=
600
300
x;= 300 x;=o
x;= 600
X1 -- 4/3X 4 + 1/3X5 = 0
Z*= 600
y los val ores de Ia soluci6n optima coinciden con la iteraci6n anterior, deoido
a que x, no pudo tamar un valor mayor que cera, con lo cual el valor de Z
tampoco pudo aumentar.
110
Ejempio 111-28.-
Max Z
s.a.
-
x. 4x= > 0
X. ~ 4
x" ~ 3
- 3x. - x2 ~ 0
(x, x2) ~ 0
I
v
z
x3
x4
b; Ec
'No.
0 II 1
2
z
0
0
--
x,
-1 '-4
x2
-4
0
x3
0
1
0
x4
0 ! 0
0
xs
0
0
I
x5
0
0
0
Po
0
0
4
¢,
a
u
I
X
5 3 0 0 1 0 0 1 0 3 3
x6 4 0 -3 ir 0 0 0 0 0
2 o ;' 1_ iI o o 1 o 1
o 4 4 1
I
1 v.b ! Ec. Z ,
I No. x• x2 I x3 i x4 . Xs x6
r----+.:..:..:c..;______-t----+---------+-------+--
z 0 1 0 I 0
1
0 l 0 i; 13/3 !--1
,
/3 iI 13 i
0 0 0 0 : 11 /3i 1/3 I 1 ~ ! 33
0 ! 0 i 0 0
,
i
r- I
1
-113 113 . 3 i 9 I
0 1 0 0 0 i 1.'3 11/3! 1
0 0 o o : 1 I o i 3 a'
111
J b Ec.
No
z X. X, X~ X~ x, p
~~
< 0
~ ........... ,----
z· 0 0 0 0 4 0 16
)(
3
. 0 0 0 ~1 4 0 8
xo* 2 0 0 0 I 0 3 --1 9
x,* 3 0 0 0 0 0 4
X/t: 4 0 0 0 0 0 3
Ejemplo 111-29.-
Sea:
X.
Max Z = 3x 1 -r 4X2
s.a. z,
''
x1 :5 600 '
~ s 300
3x. ~4~ :5 2400
(x1. x2) ~ 0
4. X.
112
De Ia representacion grafica es clara que cualquier punta sobre el
segmento lineal BC sera una soluci6n optima
Resolv1endo el modelo llegamos finaimente a:
x" ~
2400 x~=o----
x_ X~= 400
~
4/3X 4 1 '3x_0 200
x2 " X" 300
x*·= 300
2
x"'=O
5
x1 - 4;3X 4 _,__
1 3x 5 400 x;= 200
Z"'=2400
600
x:== 150
2*=2400
so!ucion que tambien es optima y que cuenta con el mismo valor para Z que
Ia obtenida en Ia iteraci6n anterior.
Si efectuamos una nueva iteracion, regresaremos a Ia primera solucion
optima obtenida (pues ~ahara es Ia que tiene coeficiente cera), por lo tanto
solo existen 2 soluciones bcisicas factibles optimas (puntas 8 y C en Ia
representacion grafica), mientras que existiran un numero infinito de
soluciones no basicas factibles 6ptimas (segmento BC en Ia representacion
grafica) que pueden expresarse como:
113
111-12.10.- Soluci6n Optima Sin Umite.
Cuando en ei momenta de obtener rp vemos que 'P. = o: para toda i,
esto indica que Ia variable x. puede aumentar su valor todo Ia que se quiera
sin hacer negativa ninguna de las variables basicas. De aqui que no exista
ninguna variable xs y que el valor de Z pueda aumentar sin limite. En aquellas
ocasiones en que nos encontramos con este caso en un problema real, es
buena el revisar que no hayamos omitido ninguna restricci6n y que nuestro
modele matematico este bien formulado.
114
Ia. X,+ a,j_ + + a .. .X.. r x-.~- =b 'i
=b~ i
1 - " '
~Z=O
=b b 2: 0
donde ·
b, 2: 0
~ 15
y
d = -(a. -t- a2 + .. +a~) (j = 1, 2. . n)
~w, =~(b. + b 2 + ... + b_J
siendo Ia solucion basica factible inicial
Xr+, = b }
(i = 1, 2, .. , m)
X = 0
D) Si Max (~W) = min W > o. esto implica que no existe una soluci6n
factible para el problema y el procedimiento termina. Si por el
contrario, Max (~W) = Min W = 0, iniciese Ia fase II, para lo cual:
donde:
d' ~ 0 y W' 0 0
116
si es que una solucion factible existe. Para factibilidad: W =0, por (3), lo que
signitica que cada x correspondiente a una d', 2: 0 debe ser cero: de aqui
que tales x puedan ser igualadas a cero y dejarse ya de considerar
Ejemplo 111-30.-
s.a.
x, T 2x2 T 3x 3 = 15
-2x, -
~
- 5x 3 = -20
x, + 2x2 -+- x3 -
x4 10
(x, , x2, x3, x4 ) 2: 0
Paso A-
s.a.
x, -t- 2x2 + 3'S 15
2x, T x2 + 5'S 20
x, ·t- 2~ + 'S + x4 10
(x, , ~. x3 , X4) 2: 0
Paso B.-
117
Max-W = 0
Z - x, - 2x 2 - 3x 3 ~ x4 = 0
X. ~ 2~ + 3~ - X< 15
2X 1 - x2 5x 3 - x5 20
en forma de tabla:
Ec l_w I
Ii v b !:No. I I 2 ,, x 1 i x2 x3 x4 x. xs x7 PO
~~~~--+~--~~--~~~~--~~-~--,--+~
i w I 0' I 1 i 0 0 i I 0 i () i () 1 1 0
i z i o !o 1 : -1 i -2 I -3 I 1 o o o o
I
x5 1 1 0 0 1 2 3 0 1 0 0 15
x6 I 2 0 0 2 1 5 0 0 1 0 20
x7 I 3 0 0 2 0 0 10
1 I I I I
v.b j Ec.l_w
, No.
Iz x, I x2 i x3 I x I! x5 I x6 I x7 ! PQ i: ¢
4
I
( :
I
i
w -s -s ~51
0' 1 0 -4 1
-1 0 0 0 I
I z ,
0 0 1 -1 1 0 0 0 i
-2 t -3 I
I xs 2 _;t
I x7 6
1
2
0
I
0 II 0
0
2 1 ! ,?
0
0
1
0 ,
0 0
0
15
20
5
4
I
X
3 I 0 I0 1
I
2 i 1 1
)
I 0 I 0 1 10 10:
vb 1:: J-wl z I i x, x2
I
I
I
I
' x3 ' x4
I
I
f xs xs I x7lpo!
I I
I
¢1 I
i
w 0' 1 0 0 -1 0 : 9/5
i
0 -9
z
xs I
,
0 0
0 0
0
0 0
1 0 13/5
1 i-3/5
0
0
12
3 15,7
x3 I 2 0 0 ~0 1 0 0 ; 1/5 0 4 20
x7 I 3 0 0 9 15 0 0 i-1/5 1 6 10/3
118
'/ ::, Ec
V·v' L
0
'c x, X ,, ~
\
;o
-~~
NO
-- -......--- -- ___. ____ ----+-------- . -
;.(:=:
---~------ -~--------l-
Vi rJ c
-6 7 0 0 :6 7 3J c c 7
z 0 0 c 0 0 ~
u -~ ~s
x2 0 c -1 7 C) 0 57 -3 7 •J 15 ( y
X 2 0 0 37 0 0 -1 7 2 ~ n
u 25 l
X-, 3 0 0 67 0 0 l 1 -9 7 4 i 157 - --'
~ 7
J.b -V
-W X. z )<2 f._ I(~ ,, X P,
w
f-~o
0
---.,.------------------- -~- ----~--
--' - - - - - - - - - -.t--
0 0 "v 0 8 0
z 0 0 -6.7' 0 0 c 16 7 -47 : 90/7
x2 0 0 1.7 i 1 0 0 57 -3,7 : 0 15/7
xJ 2 0 0 3:7 0 G -; 7 2'7 I
0 25!7
I I I
x, 3 0 ' 0 I 6!7 0 0 I -9/7 ! 4/7 I 1 15/7
Paso D-
Seve que Min W =0, par lo que el problema original tiene soluci6n y
podemos iniciar Ia fase II, para lo cual eliminamos Ia ec. (0'), asf como las
columnas correspondientes a x5 , X6 , X7 , y -W.
v.b. I Ec.:i Z !
rp
1
x, x2 "~ x, Po
'No,:
c----+=-r-- -----.-
zI 0
II
1 , --6; 7 I 0 0 0 I
90!7
I I
:v.b.! Ec. z x, p
)(2 )(3 )(4 0
1
I
No. 1
I Z* I 0 1 0 0 0 I 15
II X? * i, 0 0 1 0 i
1/6 512
I X-* I 2 0 0 0 1 l-1/2 35!14
I 3 !
I X~* I 3 0 0 0 ; 7!6 15/6
119
Ejemplo 111-31.-
Min Z = x. - X2
s a.
X,
-W + x5 -t- xs =0
z- X.
""' ~ =0
x1 - ~ - ~ + x5 = 1
x, - x. 'f' xs =2
en forma de tabla:
I I I I
I
z i x, I
X, x3 I x. II x5 II xe ! po
I
¢. :
'
0 i -2 1 1 0 0 i -3 -1
0 0 0 0 0
-
0 1 .. 1 -1 0 D 1
0 T c D -1 D 2 2
v b. 1Ec I
! No.
-W; z I
I
X;: x. ~5 xs Ip 0
rp
1
w 0' I c I
I
-1 2 0 -1
z 0 0 I
I
1 I
2 0 -1 0 -1
: x, 1 0 II 0 0 1 0 ex
I
x5 2 0 0 J -1 -1 1 II
!20
vb Ec
No -V'V z X'
--~ ----------+-------+--- -- ~--- _.., ____ ~---.----~ --+-- --~
WC 1 :JOOOD 8
z J 0 0 0 ' -1 2 -2 3
X, 1 J 0 0 0 ·1 0 2
x2 2 C 0 0 ' -1
J.b Ec
1\Jo
z x, X? x2 ' x.
c
1:
__.___.__ ________
__
p rp
~---,-----;
' ~_.
I
I Z . 0 1 . 0 0 -1 2 : -3
x, 1 0 I
0 0 -1 2 :X
'
x2 2 0 0 i 1 1 I 1
Max Z = x, + 2 X2
s.a.
x, ~ 600
~ ~ 300
3x, + 4><,. = 2400
(x,, ~.) 2! 0
121
obtenemos:
z- X, -2 ~ 0
X, ~
Xo 600
X..< + xt ~
300
3x, ~4x..
< - x. = 2400
(x~, x2 • x3 , X< xJ 2:: 0
Fase I
lntroducimos Ia ecuaci6n de no factibilidad (Max -W):
-W + xs 0
z- x, -2 ~ 0
X, +X:; 600
X., + x4 300
3x, +4x_, +Xs =2400
I l -WI z II x, II ~ I ~
v.b Ec.
No. I
I I ' x4 x51 Po
! ¢; I
!
w 0' 1 0 -3 -4 0 0 0 -2400 !
I
z 0 0 1 -1 -2 0 0 0 0 i
,.
I
I
x3 1 0 0 1 0 1 0 0 600 ex
:
x4 I2 0 0 0 0 1 0 300 1300
1 x5 i 3 0 0 3 I ~ 0 0 I
I 1 2400 1 600
I
1 v .b I1Ec. !' -W I z I x, I ~ I x3 i x4 I x5 i pc ~
[No.!
~-
1
I 0' I 1 I 0
1
I w I' -3 I
1
0 11 0 i' 4 I
I
0 I
1 200 I!
Iz I0 I0 1 I I -1 i 0 i 0 I 2 l 0 I 600
;
x3 ! 1 j 0 01110 1 1j0 OJ6001600
I; x2 J 2 I 0 O'Oi1f011IOI300IO::
I xs I 3 0 i 0 I 3 \1 0 0 . --4 1 ; 1200 i 400 I
122
v b Ec
'\JG --·--·-
Vv z '·' X,2 XC
'·
p
-__.__,-------+- --·~----~--- .. - - . - - -
\V ,~.
~ J 0 J J ~)
{_
r
D
,,,_ 2 3 1 < :ouo
x, "
~ :J () 0 4 .3 ~
' 3 }()'"'
4C..<..
X. 2 "
~ "u 0 0 1 300
X, 3 0 0 0 ~4 3 1 3 40C
1
v.b.; Ec. z : p
: No. , x, x2 : x3 x4 0
123
obtenemos
Fase I
;
v.b Ec.
, No.
,-W z i x. X2 )(3 xa X< X,
0
pc
t-- I
--+---------____,....._-----....,-
i w 0' 0 I 0 0 0 0 I
0 0
z 0 0 1 i -~ -1 0 0 0 0 0
X~ 1 0 0 3 I
2 0 0 0 20
" 2 0 ! 0 2 3 0 0 : 0 20
x4
,.
x6 3 0 0 2 0 0 .. 1
L 2
I i !
v.b Ec. W
I
1
I
1
I-
1
z i
x, i x2 x3 X )(5 '(6 Po ' rp I
I INo I i i
4
-+---1- ---+----- -··-+---·--+- --
w , 0' I 1 ' 0 : ~1 '
-.::.
i 0 0 1
I
0 -2 I
1
i I
\ z I0 I 0 I
i -1
I
I ~ 1 0 0 0 0 0
i x3 0 0 3 2 1 0 0 0 20 10
I x4 I 2
I 0 I 0 2 3 0 1 0 0 20 20'3 I
I xs I 3 0 0 2 0 0 ! -1 2 I 1
I
I
i v.b
•
I
1
Ec.
No,
1-w:
! I
Z x, : x2 x3 x4 x :<~ P
0
---·-----r--~--+------+--···-;-~---1
~~-
0 0 G 0 I 0 0 0 i
'z 0 0 1 :-1/2 i 0 0 0 !2 I 1 2 I
1 0 0 I 2 ' 0 1 0 < 18 I
2 0 o ; 112 ! o o 1 I 3;2 2 r<> 1' I
3 0 o 1:2! 1 o o ~1 '2, 1 2
124
' .b ~c
:0
?
Nc "
z - 0 - - 2
----~---
..-
~ y
3 -
2 J ::; - 2 1 7 ,• ... ,
"'
..
3 0 ~
"
::;
- 2
\'.D Ec
No
z X. X •.
<
X.
j I X:. ·:.._ :0
1
···-·-+~-
z 0 0 0 0 I -1
- : 2
1 0 J -4 0 3 14
:t2
i I
" '
x. 2 0 0 ·1 0 I 2 16 2
x, 3 0 0 0 -1 2 J.
I· v.b I Ec. 1: z I
I ·No.;
I X. x, x3 I x. ' xs PC <P
z 0 r~;3i 1/31 0
I I
II 0 I 1 0 I 20/3 I
I I I I
I
' x5 1 0 0 ~3! 1'31 0 I 14/3 I ·X I
ii x. 0 0 i1 513 i 2131 1
i 4 I
I
2
..
0
r
I 2013 i'
i x, 3 0 i 1 ! 2:3 i,3 l 0 0 I 20/3
I ~0
i 11 b i Ec I z
i
I
i x. x2 I
I
I
'-3
I
I x4 ) x5 PC
I l No: I I
zi i 115 i 1/5! 0
. I (
0 ; 1 I 0 0 8 I
I I
· • 1 I 0 0 0 ~1 1 5 4'5! 1 10 I
I xx51 1
I ?15; 3'5 fi 0 I
• 2 2 [ 0 0 j'-1 : I 4 !
I
I x," 1 3 I 0 ! 1 0 13/5 f-2J5! 0 4 I
111-14.- PROBLEMAS.-
A) USANDO El METODO GRAFICO RESUELVA LOS SIGUIENTES
PROBLEMAS
125
3x. - x2 ? 8
4x. - x2 s 19
X, + 3x 2 2: 7
x, x2 Disponibil1dad
Material 5 4 100 i
M. De Obra I 3
-+-- 8
"-----1
172 I
---j
l_ Egui~o 3 1._1 53
''
a) MinZ = 3x, + 4~
s.a. -2x. T
~ s 3
- x,
-
2x, -t-3~ 2: 6
'S 2: -3
(x,, ~) 2: 0
126
b) Max Z oc X. ~ X
s.a. -1i2X. - x ~ :
X, ~ 5
(x,, x) ~ 0
Maxz = 2x. ~·
x2
s.a x2 ~ 10
2x, ~ 5x2 ~ 60
X, ~
x2 ~ 18
3x, ~ ~~ 44
(x,, x2) ~ 0
1 .- MaxZ = X, + 4x2 + 6x 3
s.a. 2x, ~ 3~ + 2'S ~ 18
2x. + x2 + 4x 3 ~ 16
x, - 2~ +
x3 ~ 14
(x,. X2, X3) ~ 0
2- Max Z = 2x, + 3x 2 ·- x 3 + x.
s.a. 8x, ·~ 2x2 + 5x3 - x. ~ 30
3x. -r-4x 2 - :s + 2x 4 ~ 10
2x. + 'S + x3 - 3x 4 8
(x .. ~. JS, x 4 ) ~ 0
127
3- Max Z ~-
3x. ~
2X-
"
s.a. X. :s 4
x_< :s 6
3x, ~
2x.< :s: 8
(x .. >sJ ~ c
4.- Considere el problema:
MaxZ = 3x . ..,.. 2X,
s.a. 2x, ~ 4x 2 :s 22
- X~ ~ 4X, :':':10
2x. - X. ~ 7
x, - 3x 2 ~
(x,, X,) ~ 0
a) Resuelva este problema en forma gratica. lndique los puntos
extremos.
b) Resolverio mediante el metoda del Pivote.
2x, + 2~ + 3~ ~40
(x,, JS. x3 ) ~ 0
x, +JS :s 10
(x,. JS) ~ 0
128
Max Z ··~ 2x. T 5x 2
s.a. X, - 2x~ 2:10
"
4x. ~
5x 2 2: 20
X, :5 7
x2 ~ 4
,,
.::., - MaxZ = 3x. - 4x 2
2X, +3x2 2: 6
(x, x2 ) 2: 0
2~ :S 12
3x, +2~ 2: 18
(x~, ~) 2: 0
12.- Max Z = x, + 2x 2 - x3
s.a. 2x, + x< - 3x3 :S5
- 4X, -
~ + x3 :S4
x, + 3~ ~6
129
NOTA No hay restricci6n de no negativ1dad de ias variables.
a) Resuelvase por el metoda Simplex
b) Formulese este problema de tal forma que las variables solo puedan
tamar valores positives [(x,. x 2 , x 3 ) 2:: 0) y resuelva por ei metoda
simplex.
Min Z = 3x, - 4x 2 ~ x 3 ·- 2X 4
s.a. 2x, + x2 -r 2x 3 + X4 = 10
X3 - 2x4 =:; 10
x, - x2 + X4 2:: - 5
2x, + 3x 2 + X3 + X4 =:; 20
x, 2:: 0. x2 2:: 0. x3 2:: 0
130
debe fabricarse de cada modelo para maximizar las ganancias si
estas son $10 para A, $15 para B. $8 para C?
3.- Una empresa cuenta con 1000 toneladas de mineral b~, 2000
toneladas de mineral b, y 500 toneladas de mineral b3 . A partir de
dichos minerales pueden extraerse y fundirse los productos 1, 2 y
3. La empresa desea determinar Ia cantidad de cada producto que
debe fabricar, a partir de los minerales aprovechables, para ob-
tener el maximo provecho de Ia operaci6n.
131
!r MATERIA REQUERIMIETO DE MATERIA PRIMA EN C.A.DA , DISPONiBILIDAD DE ~
II PRIMA PRODUCTO i MATERIA PRIMA j 'I
Pan Pasteles Gall etas
11 Harinalkg) 12 16 8 200 ~
~
Leche(lts) 2 3 1 l 40
Huevos i 4 6 4 I 100
Levadura :, 2 2
!
30
UTILIDAD $ 3 5 6
El panadero desea saber que catidad de pan, pasteles y galletas
producir para obtener una ganancia maxima.
I 1 i :; i 3 i 4
Calorias l 2 I 0 i 1 i 3
Vitamin as --r---o !' 2 ! 1 ! 2
Carbohidratos 1I 1 T 0
I 1 1 2
Precio 4 I 10 I 12 I
15
132
figura. Los costas de envio par tonelada tambien se muestran en
Ia misma.
A.
T. A2 T;
40 Tons 30 Tons 50 Tons
2.- Supongase que lnglaterra, Francia y Espana produce todo E>l trigo,
cebada y avena del mundo. La demanda mundial para el trigo
requiere 125 millones de hectareas para poder ser satisfecha.
Slmilarmente, 60 millones de hectareas se requieren para el cultivo
de Ia cebada y 75 millones de hectareas para Ia avena. El total de
tierra disponible para estes prop6sitos, al igual que las horas de
trabajo requeridas para trabajar cada hectarea, asf como su costa
per unidad, sedan en las siguientes tablas:
133
lii-14.· PROBLEMAS
134
Personal
Hora del dia Carras Patrulla
Motorizado
12-16 300 600
16-20 400 800
20-24 500 1000
24-4 600 1200
4-8 400 800
8-12 200 400
PROOUCTO TR!MESTRE
1 2 3 4
Refrigeradores 2000 1500 3000 1000
Estufas 1500 1500 1000 1500
Lavadoras 1000 3000 1500 3000
135
dl-14_- "ROBLEMAS
136
CAPITULO IV
TRANSPORTE Y
ASIGNACION
IV-1.- PROBLEMA DE TRANSPORTE.-
Una importante clase de programas lineales, cuyo origen es ffsico y
econ6mico, puede ser formulado en terminos de una red compuesta par
puntas (nodos) conectados por rutas (arcos). sobre los cuales se efectua un
transporte (ftujo).
El clasico problema de transportaci6n surge cuando se debe deter-
minar un plan 6ptimo de embarque que:
a) Se origina en fuentes de suministro (almacenes), donde existe
disponible un determinado inventario fijo de un articulo.
b) Son mandados directamente a su destine final, en donde existe Ia
necesidad de una cantidad fija del articulo.
c) Agotan los inventarios y satisfacen las demandas, esto es, Ia deman-
da total iguala Ia oferta total.
d) Y final mente, el costo total debe satisfacer una funci6n obJetivo lineal,
es decir, el costo de cada em barque es proporcional a Ia cantidad embarcada
y el costo total es Ia suma de los costas individuales.
La tecnica de transportaci6n es aplicable a aquella suoclase de
137
problemas de programaci6n lineal. en los cuales los requerimientos y recur-
sos se expresan en terminos de solo un tipo de unidad.
La forma estandar de este problema fue formulada par primera vez par
FrankL Hitchcock (1941), en un articulo llamado The Distribution of a
Product From Several Sources to Numerous Localities.
Es importante hacer notar que muchos problemas que no tienen
ninguna relaci6n con Ia transportaci6n. se ajustan al modelo matematico del
problema de transporte.
Notaci6n:
:L a.= :L b1
I=~ j =i
138
En forma conc1sa tendremos:
Encontrar x (i = 1.2 .m j ~ 1.2 .... ,n), de forma tai que se minimice.
=
==.
CX ~ Z
m
\~
x = b U '-'. 1.2 .... n)
g
e
n m Cm1 Cm2 Cm3 Cmn am
Requerido b; b2 b3 bn
Tabla A..
Una importante caracteristica de este problema. es que el modelo
139
puede ser abreviado en Ia forma de una arreglo rectangular. como se muestra
a continuacion.
Destine
2
0 X11 X12
2 X21
g
e
n m Xm1
-----------~--------~----- - Xm2
---- ----~·
Xmn 3m
requerido, b2 bn
Tabla 8
Las tablas A y 8 se pueden combinar de tal forma que el rengl6n (i),
columna U), muestre los vaiores de x,1 y de C, 1 , con lo que tendremos:
Destine
2 n Ii d'1Sponluo8
'hl
, c1_1: c12
0 -+------- - - - ' + - - - - -
x2., ,,x22 'X2n :a2
2 I . . . ! ! I
[xmn lam
1
n m
crn1 cm2 : I cmn :
--+--- -t------+-------;-------+----
requerido~ !b
: n
140
Las ecuaciones rengl6n son "' x = a (i =- 1, 2, . . ml
Ejemplo IV-1.-
0 Destino
J K L M N
p 20 10 15 10 8
g
Q 15 30 5 12 14
e
n R i8 15 20 7 19
J ! 00 unidades
K 70 unidades
Total 580
L 140 unidades unidades.
M 1 80 unidades
P 200 unidades I
Total650
Q 300 unidades unidades
L -R
- - - - 150 unidades
~~~:__ _ _ _ _ _ _ __
141
Min Z 'c 20xG - 1Ox~, - 15x~c ~ 1OX 1 - 8x"" 0 ,_
Sujeto a:
X;p ~ XJ::l 1- X"'1 100
xKP -'- xK:J + x,<r=< 70
X~P -r XLO - X_q 40
1
_,_ x'A"
x.~F + x,~c: 180
xw + x,c: - x," 90
x p + x,P - x_, + x.,.p _, x,P
0
-s 200
XJ 0 T XKa ____._. xl 0 -+- x~ 0 - x~c s 300
xJR -~ xKR ' x._R - xf.•R -+- x\IP 5 150
x, 2 0 para todas i y i
Para formuiar este problema como uno de transporte, hacemos:
Origen i = Almacen i (i = RO,R)
Destine j = Cliente j U = J, K. L, M, N)
x 1 = Cantidad embarcada del almacen (i) ai cliente (j)
C. = Costas unitarios asociadas con x,,
a,' = Unidades disponibles en el almacen (i)
b = Unidades requeridas par el cliente (i)
' -
Nuestra tabla quedara ahara:
Destine
,---,---~--~----,-----------~ .........,
142
y el simbolo x. no aparece ni siquiera en Ia primera tabla. Aqui hemos
introducido el slmbolo como una representacion de que en su iugar debemos
poner el valor correspondiente.
Ejemplo IV-2.-
Monterrey 8 7 4 7
Guadalajara 10 3 6 8
Mexico, D.F 7 9 4 5
Puebla 8 6
143
Min Z = SxMP ~ 7xML -+- 4x~~R ~ 7x,,& • 1Cx_3 "' , 3x.:;L ~- 6x'-'F - 8x 83 •
sujeto a:
p L
Destine
!
R
I
-s Pisponible
I
---,
IXMP ;xML -!XMR ~ j20
M
Bl
'1~71
0 _l_ I
I
136
IXos I
D I
e
I
'
7' 91 ______41_;_____ 5 i ___ j
I
I;XPL I
I !16 I
n I p \xpp I!XPR lxPs I
! Mi M! 8! 6/ I
!
I I
requerido \c____
12 123 121 \29 I 85
En este caso conviene notar que como los requerimientos son iguales
a las disponibilidades, no hizo falta introducir ninguna variable figurada (ode
holgura).
Notamos tambiem que dos de lo caminos no 'pueden ser usados debido
a que ningun tecnico puede ser mandado de (P) a (P) o a (L).
E! problema de caminos inusables pude hacerse que ajuste con nuestro
problema general, considerandolos como caminos usables, pero
144
asignandoles un alto costa (Metoda de Ia Gran M) y reteniendolos en Ia
matnz. Si Ia soluci6n se intenta a mana, Ia Gran M puede no usarse, pero es
indispensable cuando Ia soluci6n utiliza otros medias.
En el problema de transporte, definimos las variables basicas factibles
igual que lo hicimos antes con el problema de programaci6n lineal general,
solo que aqui Ia soluci6n se efectua directamente en el arreglo rectangular
que hemos planteado.
Vemos que en un problema de transports tenemos (m + n) ecuaciones
(m orfgenes y n destines) Si empieamos el criteria desarrollado con
anterioridad al cubrir el problema de programaci6n lineal general,
deberfamos tener que nuestras soiuciones basicas factibles no degeneradas
deben contar con (m + n) variables mayores que cera (m + n entradas
diferentes a cera en nuestro arreglo rectangular). Esta consideraci6n resulta
de Ia omisi6n de una nueva restricci6n, con Ia cual contamos ahara, que es:
rn
2.: a, = 2: b1
·=1 !=1
y como:
n
·"'
L x IJ =a1
para ( i = 1, 2, ... , m )
j=1
m
)......., X lj =b1 para ( j = 1, 2, .... n )
:=1
tendremos:
m n -1 m n n -1 m
I
1=1
aI - .._,
)
I ::::1
bi 2: 2:
,::::, )::::1
Xlj - 2: 2:
j=1:==1
XIJ
m
2.: X para j = n
1=1
bn
145
Definici6n IV-1.-
146
ContinCJe er esta forma. paso a paso alejandose de Ia esqu,na
noroeste. hasta que f1nalmente obtenga un valor para Ia celda sureste
Esta regia selecciona tantas variables basicas como haya columnas
mas renglones me nos uno (m- n~ 1). debido a que en el paso final quedan
un rengl6n y una columna. los cuales son ambos eliminados.
Aplicando el metoda anterior ai ejemplo IV~ 1 •
\2) (5) (11) {;4)
j K L M N s
Paso 3.-
147
Eliminamos Ia columna 2 (5)
Reemplazamos a'. por a _ = (a' b) (100 - 70) --~ 30 (6)
Paso 4.-
Paso 7.-
Paso 8.-
148
=> X35 = 80 (19)
Eliminamos ia columna 5 (20)
Reemplazamos a 3 por a 3 = (a 3 - b' 5 ) = (150- 80) cc 70 (21)
Paso 9.-
149
se define como Ia diferencia aritmetica entre el element a mas pequeno y e1 siguiente
mas pequeno en ese rengl6n o columna Esta cantidad nos tndica Ia penalidad
unitaria rnintma en que se incurrina si fallamos en hacer Ia asignaci6n a Ia celda con
costa minima en ese rengl6n o columna. De aqui que este procedimiento. en forma
repetitiva, efectue Ia asignaci6n factible maxima en Ia celda de menor costa del
rengl6n o columna sobrante que tenga ia mayor diferencia.
Debido a que una asignaci6n ya afectuada no es cambiada, el
procedimiento termina tan pronto como Ia oferta y Ia demanda total se
agotan.
Los paso a seguir son los siguientes
1) Construya Ia tabla de costas y requerimientos y vaya al paso 3.
2) Use Ia tabla de costas y requerimientos para el problema que queda
despues de que las asignaciones previas (tanto ceros como
positivas) han sido efectuadas.
3) Calcule cada diferencia rengl6n y cada diferencia columna.
4) Seleccione el rengl6n o columna con Ia mayor diferencia (en los
em pates se escoge cualquiera en forma arbitraria).
5) Asigne tanto como sea posible a !a celda con menor costa en ese
rengl6n o columna.
6) Asigne ceros a todas las demas celdas de aquel rengl6n o columna
en Ia cual Ia oferta o Ia demanda se agate.
7) Efectue Ia (mica asignaci6n factible en todos aquellos renglones o
columnas que ya tengan solo una celda sin asignaci6n
8) Elimine todos aquellos renglones o columnas total mente asignados
y ya no los tome en cuenta en las iteraciones posteriores. Si ya no
existen renglones o columnas sin asignaciones, pare. pues ha
terminado; en caso contrario vaya al paso 2.
Aplicando el Metoda de Aproximaci6n de Voguel al Ejemplo IV-1.
tendrfamos:
u
I
K L M N s
0
p 20 10 15 10 8 0 200 i 8
'~0
Q 15 30 5 12 14 0 300 5
R 18 15 20 7 19 0 150 7
3 5 10 3 6 0
ft
150
--
'-'
--------------
K
--- -
M
--------
N
--~---------
s
p 20 10 10 8 OJ 200 8
0 15 30 12 14 07[ 160 12 <:=
R ! 18 15 7 19 0 150 7
L......_.._.._~-- - - ----------------~----~-
100 70 180 90 70
3 5 3 6 0
J K M N
p 20 10 10 8 200 2
0 15 30 12 14 90 2
'50 c
R 18 15 7 19 150 8<:=
100 70 180 90
----
3 5 3 6
J K M N
----------------------
70
p 20 10 10 8 200 2
0
0 15 30
L __ _ _ _ _ _ _ _ _ _ _ _ _ _ _ __ 12 14 90 2
100 70 30 90
----------~------
5 20 2 6
ft
J M N
90
---,
p 20 10 8 ! 130 2
!
0 ~--~~-----1_2___ .. 2_4__j 90 2
100 30 90
5 2 6
11
J M
I ·o
p ! 20 40 10 ¢:
:
90
Q . i5
L _ _ _ _ _ ____l
1/ 90 3
100 30
-------- . -----------r-- - -
5 2
151
Cabe aclarar que con practica todos los pasos se pueden hacer en Ia
misma tabla.
Los pasos anteriores nos dan Ia siguiente soluci6n
J ,
· - - - - .- - - - - , - - - - - +
K L
--+---
M N
-,-----~---
s
10 70 30 90 200
p 20i 10; 10 8.
15.-r-------"-+---------=----t--- 0
-- -+------------r--------~-;"- -- --r--·---
0
90 140 JO 300
-~----t- ~----~
15 30' 5 ' 12
~---------------J..-----t-
14 0. --
R .150 i 50
18 15 : 20 I 7 19 0
----,-------------t---------------L- ~ -----
i HJO 70 i 140 180 ,90 J 70 I
I
c -= 5020
sujeto a:
152
multiples de estas ecuaciones, conteniendo !a variable basica correspon-
diente con coeticiente unitario. se pueden sumar o restar a ia funci6n objet!vo
para eltminar las variables basicas de ella. Llamemos u (i = 1.2 ... m) al
multiple de Ia ecuaci6n rengl6n (i) y v U -~ 1.2. .n) al multiple de Ia ecuaci6n
columna U)
Tenemos entonces:
z "' "-
-
'·
,,
"'
- \ I X ~
<,
--::..=: u ~ ~
'\' \/
~
b
::=.< ··..:;
=
TABLA (A)
---~+- ~ ----;----
K
~·--·-~
M N s
I
100 l?O I ,2oo 1
p ' I
0
~
20 10 0
I
1
153
l:scogemos el rengl6n uno como ia ecuaci6n redundante y hacemos
u, cc 0, con lo cual i:1mediatamente obtenemos (recordando que
C s ~ u - vs) v, = 20. v 2 = 10 y v 3 = 15
Como C 23 c::), = 5 y vc = 15
"" = ~ U~ -10
"
Como C 24 = cc:M= 12 y u 2 = -10~v 4 = 22
Como C 25 = CO'< = 14 y u2 "" -10 ~ v5 24
Como C 35 = c"N = 19 y vs = 24 ~ u3 -5
Como C 35 ,= c"'s - 0 y u:.; - 5 V, ~
0
5
6
C., :::: U1 + V1 para toda (i, j )
X e =Xst
154
jonde x,. sera aquelia variable para Ia que
TABLP. (8:
K :_ cJ
--------r---~---~ --- - ----~------,
--i6, -5 200
p
_ 10_:____ s_,·~~~o-~~-oi
C.:
v: 30. 5 •300
0
I 15 30 ! 0 . 10 '
!-----+--------~~----- ·- ------t-----~----t~~-----+---------- i
3 10 ' 10 -1 0 I 15Q
R
18 ' 15 ::!0 7 i 5 1
-~~+-~ -----L~--+------__.._ _ _ _ __
100 ,70 140 180 .90 1 70
v: 2o i
~-_.]_-
1o :
__[__-~------I
1s 22 24
~~-~~-
I s !
155
restricciones. Para este fin. pondremos - !:fen algunas variables basicas y
-A en otras, y nada en las restantes.
Conviene hacer notar que cuando hacemos xe ·~ H. habremos intro-
ducido con esto una entrada no independiente, ya que ahara existira un
circuito cerrado que parte de xe y regresa a xe. Las entradas en las cuales
cambia de direcci6n dicho circuito,seran las que se veran afectadas por un
_._ fJ 6 un ·- A segun sea el caso. a fin de mantener Ia vaLdez de las
restricciones.
Aplicando Ia antenor a nuestro ejemplo obtendremos:
TABLA (C)
-·- --+---+---
J K L M N s
-l------,.cr----,----r------,
:100 i 70 --.-30- {)~ ! + fJ
1
: 200 I
p
I 20 i 10 15 10 I 8 I 0:
r-------r-1' :11 o..,... aC:J}: 180 :, 1 o- e(a) ;-------30-o- --l
l Q I 1 ,
i : 15 ! 30 i 5 ; 12 14 : 0:
t----+------+--·- ---- ~-------,---~ --i:-7-0--·1-1-15-0----!
I R I 18 • 15 i 20 I 7 I . 191 0I
L---~---~---~---T---~---~--~--~
Tambien aqui debe notarse que el trabajo efectuado sabre Ia tabla (C)
pudo haberse efectuado sabre Ia tabla (A). Nuevamente se separ6 para hacer
mas clara ei ejemplo.
Una vez que todas las fJ se han puesto. reemplazamos fJ par aquel valor
156
que haga cero primero algunas de las "anabies :Jasicas. Esto es. 6 toma el
valor X0 •. correspondiente ai vaior menor de las variables basicas a las cuales
se ies adjujic6 Ji de forma tal que xcr -8 se haga cera .A.si. tenemos que:
xo_ = x,.
Usando el valor de 13 asf determinado, todas 1as variables basicas se
recomputan para obtener Ia nueva soluci6n basica
2'! soiuci6n
TABLA (D)
), oo ..· ~ 70 . . 20 . : 10 / 1 : 200
~-P---+~·~·/_._.·_2_o4J~·-·--·_,o~11 /_,_./__./_,_s7j____1_o~l~·---~---o~l____ -4
157
~--~--~--~----TA_BLA (E)
J K L M N s u.
-+-~--
·-~---t-
J 20: 10 I 15 ! 22 I 8 -11!
Como C' 34 = -26 < C'" para toda (i,j), entonces:
x34 = xRM = x.
AI introducir & en Ia celda (3.4) y formar el circuito, tendremos que las/
entradas basicas a las cuales se les asigna- e son:
20- fJ celda (1 ,3)
180 - 8 celda (2,4)
80- fJ celda (3,5)
Vemos entonces que :
y que:
e = 20
Sustituyendo el valor de e. obtendremos Ia tabla (F).
3? soluci6n.
I +fi -211i
1 12 141 o 1~ 4 140 /~j 16Q -:Jfl -1Q 1. -5 !300 \
Q
fsl 1 30 , 5 1 . 1 .
----~----1-3~~----u-+l---2-o~:2-0--+-8V~60-----~-~./-7o--,-.~/!-15--o
1
1 1
R 18 I 15 ! 20 i / .· 7 / 19 ~- / 0I 11
100 170 ; 140 180 90 : 70
I . i 1,.
158
L = (100 X 20) • (70 X 10) ~ (30 X 8) ~ (140 X 5) ~ (160 X 12) -~ (20 X 7) +
(60 X 19)- (70 X 0) ~· 6840
(j = 60
X s =X 35
TABLA (G)
-··--~-----~---·
.J i K L M i N S 1 u,
I ---~~-5 i 70 ---.-- 10 '30 ~---T,()-·-1200 I
P . 20 1 1o 15 1o I .. .. 8 : · . o: oI
.~---t----+-- +---____,...---+-----~
:
1
100 I 20 140 2! 6 60 300 :
. Q 15 ·.' I •
~-+· 12J___..J_~ o• o
I! I '
I 30: 5 1
6i 8 18 150 I 14 3 150
1
Vemos que esta soluci6n es optima, debido a que C'. 1 > 0 para toda (i,j).
159
Nuestra soluci6n inic1al sera .
·-·-- -------~
J K L M
--r-----,---------.-~-------,----t--
N s ·-· ------
; 100 70 10 '180
p
20 10 15 10 8 0
r------- ---;.--~---- -----r---- ~------- ------ -~ ----------
z= 7220
Como no hemos modificado ios costas. nuestros C', deberan ser los
mismos que en el ejemplo anterior y Ia celda (1 5) corresponded! a ia variable
de entrada, con io que obtenemos
r-------------~--- --- T - - · - - - - - - - - - - - - - - -
1 J . K . L ,._/:
N S
r-----t-- T-~-·
1
I '100 70 180
p ..\:J-
2G 10 j•)
:·~
8 '
',
i--t-------~-.~----------------···-.{ ·-·----,-----···-- +---·
110-- 1;t 18C~ 1 0- A I .:300
0 c
15 30 14
·--~---· .. - T - - · · · ·-·.,.---. .
so ,?Q ~50
A
18 ' ~ 5
--- - - - - - - - - + - - --+------
20 ~-- ___19 i / 0
7Q 120 180 ,":.1".,_· 70
de donde vemos que x, 3 y x25 estan empatadas para dejar Ia base. Selec-
cionamos arbitrariamente una de elias ( ~ 5 en nuestro caso) y a Ia no
seleccionada (x, 3 ), para indicar que sigue siendo basica, le asignamos un
valor (t:) infinitesimal y esta variable se manipula como cualquier otra
asignaci6n, con Ia salvedad de que :
t + < {t si.tt=O
> o
c: Sl
160
obtenemos entonces·
J
--..---..,--·
K L M
_ __________.__ _________
N
-+-·------~~----'
s -------
70 r- fJ -1210 ~ 11 11'180
p
20. 10 10 8i o: 0.
·----+-------L.--~- '
5 30 120 ~fi 180 -e!--1~ 21:300
Q
15 30' 5 12( 14 0. 10 '
-----+--------~-- ~--.,---·
v .. 20 10• 15 a: -11
Z 0
= 7060
X c~ X
e 34
f) = €
,------,-- -,-- - - - - - - -
: J K L • M I N 1
s u
,-------+-P
[100;-Yho i 26t 14\10+&-/j 11:1so
; I 15 1 10 :
. .· . 20 I 10 . 8I 0 : 0 I
~~--~ o--2_,1-1,1~--4+!-,2-o--
.. ,-s-o---fr4 +-: 1_o+1. --_-s~!3-o-o-~
1_ _ _ _ _
z= 7060
&=80
161
~-- --~ ~-------- ------,-
,-------~
J K L M
---~--~----r- ------~-------~-,
N s u
8 '
15' 20 80-,lf I 21 70 - if i 150
R 20:
18: 15. 7! 191 0 -10
-r----
'100 70 ! 120 180 j90 JO
I
I
v.I 20: 10! 10' 17' 8 10i
______j
z = 5380
e = 20
J
;----_;__---+----,t------L----~ ---
K L M :-,: s uI
!
1
10 : 70 / 15 3 C}J 20 180 I
_P_+-__2_0+(_ _1_0-+1_ _1_5-+-r ~~~-- -----~- 0 0 i
i 100 15:120 ,80-1} : 1 A .:.:sj300 !
0 I // 15:
3o 1 - 5· 12 - 14 : o- 1 s !
~--i'------i'-----5+-i
8
-- 20 I 1 oo ~~A~----11T5o=-&11so-j
R
'-----r----~---+-1- - -
18 1s[ 20 !. 7;
--+-----r------t--
19 : . o oI
100 70 120 1 80 19() . 70
v 10 10l ol 7\ 81 0
z= 5180
Xe -== Xr,·~S
e = so
J K L M s
c·-· -- -- ---5 -io·-·--r----1 0 ~4-----~2~90----~20-=-~- -,-180 ---,
p
20 10 • 15 . 10 8 0 0
·------------~-----~-~-~.----cr---~----+----~---
I 8 10 : 20 150 16 5 150 ;
R 18 15 I 20 7I 19 0 . -5 !
--- ---+---L -r----------j__~--
100 70 120 . 180 ;90 ! 70
' I ;
v, I 15 10 5 ' 12 I 8' 0 :,
~--
z= 4930
e = 20
----------
J K i L M N s u
I
r - - _ T~· 7 i 70 /
I: 12 120 . / ,90 I
'
I p .I 20'/101 15 I /10 I ai
1oo l/ l8j12/o/. l~o ~~-
i0
L__
1
I
15 I 301/ 5 I/
i/ 121
I R 11. 8 ; 8I 20 i 150 .. I'
141
i
5 '150
I
i
I
181 15 20 I / / 7 19 I 0 il -3 II
---+---+--+--+---------<
[1oo !7o 180 'oo !7o
163
TRANSPORTE
Al X
A2
El cuadro es equivalente al r a <b (1\
modele matemat1co. Si
t b.<a(2)
MinZ = ,:S " c x
Jrengl6n i ( 1)
1=-l :==~ eiim1nar
lcolumna j (2)
sa. \' X ,, =a
..... ra b (1)
~=·
,.,.,
y actual1zar 1b a
12 )
)' X =b A3 kn Ia ult1ma celda
·=1 (xmn_l_:""~.:m ~ bn__
'..IETODO DE LA
ESQUif'</1 NOROES:r!O__
::uALJRO
DATOS ~
"JIC!Al
METODC DE INICIAl
APROX. DE IOG_E~ __ _
164
z z z
MIN
K-AVA
1~ ITERACION ~ ITERACION
ITERACION
C4A Haciendo
C3 Se determina x = min C' tJ < 0 I
. e X = + ()
C4 Se determma x5 e
Cerrando un "cir-
C5 Se recomputan las Se calcula el Valor de C48
cuito" y ajustando
otras Variables+--+ Z
sus entradas ba-
Basicas (V 8.) del para esta iteraci6n
sicas con -8y +Ben
circuJto.
sus esqumas
C6 Se repite el cicio +---"''ota: El valor de Z
hasta que no haya debera disminuir en
cad a iteraci6n o · C4C Reemplazando (:)
c 'i < 0 per el valor que
mantenerse = en
case de degenera- haga '0'' primero a
alguna(s) V.B. del
a6n.
"c1rcuito" se en-
cuentra X5
165
IV.-2.- PROBLEMA DE LA ASIGNACION.-
Exist en casas especiales del problema de transporte. los cuales presen-
tan ventajas de computaci6n adicionales sabre este. raz6n porIa cual se han
desarrollado metodos especiales de calculo para estos casas. Una de las
situaciones mas comunes, es aquella que se conoce como problema de Ia
Asignaci6n , el cual es el caso mas sencillo de todos los problemas de
programaci6n lineal.
El problema de Ia asignacion es un caso especial del modelo de
transporte, en el que existe una matriz cuadrada (m = n). con cada una de
las restricciones de disponibilidad y cada uno de los requerimientos iguales
a Ia unidad (a = 1 y bi = 1 para toda i y j ).
Un problema tipico es aquel de asignar personas a diferentes trabajos,
en donde tenemos tantos trabajos como personas y cada persona debe
asignarse a solo un trabajo. A una persona i se le asigna el trabajo j con un
costa C.'I
En el caso anterior:
(2) J: X.=
- 'I
1 (i = 1,2, ... , n)
166
E! objetivo de! problema de asignaci6n es e! de escoger las x quP
satistagan (1) (2) y (3), en tai forma que el costa total
Z = 2: ~ C.x
sea minimizado.
EJEMPLO IV-3.-
I
I
I
r-
T1 -I I
T2 T3 T4
! M1 10
I
I 16 12 8
!
M2 8 12 I
15 12
---+-
M3 15 !
I
13 13 11
I
M4 12 I 15
'
' _ _ _ ____l__
_L
10 7
___j
167
Vemos que si intentamos resolver el problema utilizando nuestro
modelo general de programaci6n lineal, obtendremos un modelocuya matriz
de coeficientes es enorme
=1
=1
=1
X41 ·~ ·:Gtl .x.,., =1
=1
=1
=1
·X44 =1
T1 T2 T3 T4
I
M1 I
121 81
r-
11 i
I I
M2 I
Bi 121 15:
I I
121
M3
I
:o 1,
I
I
11;
151 131 13l
I
M4
I ro !1 '
12i 151
!<
. 101 ?I
11 11 [1 11 4
I
Teorema IV-1.-
168
entonces una asignaci6n que hace mfnima Ia efectividad total en una matriz,
tambien minimizara Ia efectividad total de Ia otra matriz.
Demostraci6n intuitiva:
Si restamos A unidades de cada elemento del rengl6n (i), tendremos
que como cada soluci6n posible debe tener exactamente una asignaci6n en
el rengl6n (i), el costa total para ia nueva matriz sera siempre exactamente
A unidades menor que el que hubiesemos obtenido con Ia matriz original De
aquf que Ia soluci6n que minimice el costa total para una matriz, debera
tambien min1mizar el costa total para Ia otra.
Utilizando el teorema anterior, podemos transtormar nuestra matriz
original de costas en una que consista de elementos positives o nulos. Una
soluci6n optima consistira en hacer las asignaciones a los elementos con
valor de cero.
El procedimiento de soluci6n estriba, precisamente, en obtener una
matriz modificada como Ia mencionada.
Definici6n IV-2.-
A todos los elementos del rengl6n (i) se les resta el elemento mas
pequeno de dicho rengl6n.
En nuestro ejemplo, escogemos el numero mas pequeno del primer
rengl6n y lo restamos de cada uno de los elementos de dicho rengl6n. El
resultado es:
T1 T2 T3 T4 j
M1 2 8 4 0
J ''
M2 8 12 15 12
l
i
I M3 15 13 ! 13
M4 ~
12 15 10 7
169
Supongamos ahara que hemos asrgnado un trabajo a cada maquina,
cualquiera que sea Ia asignaci6n que hayamos hecho. el costa de Ia misma
con Ia nueva matriz, sera $ 8.00 menor que con Ia matriz anterior.
Procedemos ahara a restar el elemento minima en cada rengl6n res-
tante de todos los .elementos de su rengl6n, para obtener:
T1 T2 T3
----~---------------~------~-
2 8 4
0 4 7 I 4 .
---~-- r------
4 2 2 i 0 .
3 0
Paso 2.-
r-------·-,---------,---- ---,-----~
T1 T2 i T3 T4 i
! j .
170
( ) ya que ahi se efectuara una asignaci6n Marque con una X todos los otros
ceros en !a misma columna para tndtcar que no se pueden usar para hacer
otras asignaciones
B) Examinese a continuaci6n las columnas para encontrar ceros no
marcados (uno exactamente), los que se marcaran con un cuadrito ( ; y con
una X todos los demas ceros no marcados en el rengl6n correspond1ente.
Cl Rep1tase (A) y (B) sucesivamente hasta que una dP ::Jos situaciones
ocu:·ra
a) Ya no hay ceros sin marcar
b) Los ceros que quedan sin marcar son cuando menos dos en
cada rengl6n o columna.
Silo que resulta es el caso (a), tenemos una asignaci6n maxima (que
puede no ser completa).
Si resu1ta el caso (b), podemos usar el ingenio y/o tanteos para localizar
una asignact6n maxima.
T1 T2 T3 T4
2 6 2 0(1)
M1
(1) par (A)
M2 0 (2l 2 5· 4
(2) par (A)
M3 4 0 (3) ~ ~rJ
(3) par (B)
M4 5 6 (}01
[_
Vemos que tenemos Ia srtuaci6n (a), que cae dentro del caso li y lo
dejamos pendiente hasta que veamos el paso 4.
171
• -2.- ~P:OBL:::-Vl-"1.. :)E _;., ..:..S.G'-oAC ·;:::•.,
Ejemplo IV.- 4
7 0 5 4 3 7 0 5 4 3 7 0 5 4 3
0 0 3 8 6 0 Jl: 3 8 6 0 0 3 8 6
I
0 6 0 3 0 0 6 0 3 0 {} 6 0 3 0
0 0 7 0 i 0 0 7 0 .0 0 7 0
4 3 0 0 0 4 3 0 0 0 4 3 0 0_ _0:___j:
~'- - - - · · - - - - '
·-----
: ! l} ..
7 0 5 4 3 7 0 5 4 3 i 7 5 4 3
I
i 0 )I 3 8 6 0 ~ 3 8 6 0 .. a:' 3 8 6
:.tKJ
i
lu. 6 0 3 0 !Y 6 .0 3 {} 6 0 3 ff·
i
! {} 0 7 0
r
r {(' '•Lt 7 rrO Q 7 0 1
!
I
I
I
I
L4 3 fl 0 :G. __j
'I 4 3 s 0 4 3 9:, 0 (f
Nota.-
172
Paso 4.-
M4 ~-~----""--!~•• ..-(1)por(a) !s
L
6 Q-'i+-
t t
(2)por(b) par (e)
173
~ 2':::f {f)
·-------~ ~-- ------~~--------,
I 0'(41 ol3l
4- 0 -- 0- 1*- y obtene- 4 (3) por :A!
'
c· = s -r 8 -r 13 + 1o = 39
Ejemplo IV-5.-
A I!' 10 16 12
B 8 12 15
I
c ! 15 13
D I__ 12 15
Soluci6n.-
c 15 M 13 0 c r M 0
I Q __jI 4 3 ~vi
D 12 15 M D 0
--··-· -~------- ------------
174
TRONSPORTE V AS~G"'ACION
·----t --
M(1) M(2) M(1) M(2) M(3) M(F)
i
A 2 4 A 2 4 ,0
B .. 8 0 0 3
c 7 M 1 '' 0
!
I
~~-
I
c 6 M 0 0
D 4 3 M ~ 'JJ(
' I,.......____ ' D 3 2 M 0
~
II
' M(1) M(2) J(3) M~(F) M(1) M(2) M(3) M(F)
)<
A
B
c
I
2
j_-=o=J
6
I
l.._.____J
/~
M
4 ,if·
~
I
::s:
I
'Ol
~
1
i
~
.B IX A
c
:o
4
r--,
!
~I
2
0_ji
M iol
5
~--___.! X
3
Ejemplo IV-6.-
'~ •~...s ....,, C* =
<:
. • ~''!'Oi
.
10 + 12 + 13 + 0 = 35
, -<_.,:.
,, ; " J.
! ~
7 14 4 3 5 4 11 0 2
2 6 9 6 7 0 4 7-... 4 5
9
4
8
5
12
7
10
9
5
6
1
1
• 4
0
3 7
3
5
5
0
2
0 5 6 7 8
_,
0 5 6 7 8
1 1 1 1
II
• '
~
~
rt 1i
2
2
5
4
2
5
175
--------~
8 11 0 n 2
~ 1} 2 G
8 3 6 5 0
3 0 4
0 3 4
c· 4 6 T 5 + 5 T 0 = 20
ASIGNACION
DATOS:
MinZ =I ,=,
m
I c, xij
[a:l- . <
1.,
~osta aei e1emento
pequeno ds .:aoa rang16n
8 ·o"S elementos ae ese
rr"~8s
1=1
1 I 0 I-• rer.~to...,
0L.:_j
En cada colwmna aun S11" Os
~a as:Qf1aciO·.., a~~~
Ltb
CASO II - - - - - (C, - - - - CASO I
max 1ma .__,___. $upon1endo Que :a
maxtma 1caros s;marcar) as
para no H t .. ne "'"a as1gnac:O n
cornptela 1una por caae reng1on y
~:~~~:rr~;g~o~:,::~. :~:Cti'.IIOad una por cada columna) en 3A o JB
med1ante s~mato y restasJ ~bl Mercer columnas con
ceroa as1gnados o no er
~
on ,es poslcJones q~.oe e~
Matnz de Efect1..,1dad
e1ermmen Z rntn
J
<iN
176
Destino
, Dispon:-
2 'bil:dad
-+--~~-
c 6 4 2
L
A.
s 2 8 5 4
E
~---------+- -r---
Dernanda · 3 3
Ciudad
2 3 4 i GAS CI:>PONIBLE
p
E 4 7 9 10 8
F
~ 2 1 6 4 3 6 10
E I
~ 3 i 9 6 4 8 6
__
A --- -· -- I__
'lEOUERIMIENTOS l 5 3 8 4
9
~-----
16 4 '
''
. 19 8
2 c 12 10
8 \.} 8 8
3 12 3 23 6 30
----- ---~----
bj 5 4 9 8 22
4.- Use los elementos: { (1 ,4), (2,3), (2,4), (3,2), (3.5), (4, 1), (4,2). (4,3),}
a
como solucion inicial basica para = 0 y resuelva el problema
I
j•
2 3 4 5 ai
6 11 14 11 4 3 +a
2 7 17 9 10 5 15
3 9 5 7 18 0 14
4 9 13 13 15 7 10
bj 3 9+o 14 7 9
~Para que rango de valores de o Ia soluci6n basica optima obtenida
continua siendo optima?
5.- Considere el problema de transporte, teniendo Ia siguiente tabla de
requerimentos y costas.
I
Destino I
I
2 3 4 5 i Disponibilidad
I
8 6 3 7 5 ' 20
0
R
I 2 5 M 8 4 7 30
G
E 3 i 6 3 9 6 8 . 30
N
_4 _l__c>_ ____o__o_~o___o---i___2_o___
Demanda [ 25 25 20 10 20 j
178
b) Obtenga Ia solucion basica tactible :niciai siguiendo ei metoda
de Aproximaciones de Vogel.
·> Compare los valores de !a tuncion objetivo para esas
soiUCIOnes. EsCO)a Ia mejor y obtenga Ia SOIUCIOn optima.
2 3 8i
---------- __________ , __
5 14-2!. 7 15
2 6 7 8 20
~3 : 14 ---- ~. 24 3) 9 4
-.-~ ---+--~--·--· ----- ----~ ·- ----~~~~--
b 7 11 21
Mercado
I ; en Ia planta
--;--1-+-: I
2
9
3
5
4
6
5
7
(tons)
10
~ 2 I 2 8 10 13 25
N I
T 3 i 3 12 6 5 2 ~ 13
_A___
Requeri:-niento~
4 I 1 --~---~~~-3__ .~2_1___~----
i
en el mercado 1 30 22 17 19 12 :
I
/tons)
8.- Una firma que produce un solo producto, tiene tres plantas y cuatro
clientes. Las tres plantas produciran 3000, 5000 y 4000 unidades respectiva-
mente, durante el siguiente period a de tiempo. La firma tiene un contrato para
vender 4000 unidades al cliente 1, 300 unidades al cliente 2 y cuando menos
1000 unidades al cliente 3. Los clientes 3 y 4 quieren comprar Ia mayor
179
·_, 2 ~ .... ERC:<:;·c:; ;.;RCP 1....,t=:S~:JS
CLIENTE
2 3 4
---·--- - - - __ _ j
p
L
65 63 62 64
A
N 2 68 67 65 62
T
A 3 I 63 60 59 60
- - - · - ---·--
La directiva desea saber cuc'mtas unidades deben venderse a ios
clientes 3 y 4 y cuantas unidades deben enviarse de cada planta a cada
cliente, con el fin de maximizar el beneficia.
2 3 4 2 3 4 5
---,
A 7 8 3 6l A 2 4 6 8 5 !
I
I II
B 4 6 5 B 3 5 7 9
! I
c 9 9 7 4 i c 6 7 9 4 3
D I
2 5 6 7 D 3 5 4 7 4
E 4 3 4 3 4
·- -·-
2 3 4 5 2 3 4
A 6 5 3 7 /l, 4 3 5
B 6 2 4 3 5 B 3 9 3
c 7 2 3 c 7 3 2 2
D 7 5 4 6 6 D 6 4 7
E 9 3 5 E 5 4 3 3
180
TEORIA DUAL
zbr;
;?:.
,,.
' -·
f";-
..... , ., ..
·~ c 'c
CAPITULOV
TEORIA DUAL
.. ,li
V.1.-PROBLEMA DUAL
' I I
I 1 I ' i
Para cada problema de programaci6n lineal, existe otro problema de
programaci6n lineal muy estrechamente relacionado con el, al cual se le
denomina DUAL. AI problema originalle denominaremos de a qui en adelante
PRIMO.
Sea el siguiente caso:
Problema primo.-
(O) Max Zx
I 1 ,.
!
= C 1X1 + c~
I
I
+
t. 'I
s. a.
',01'""\ ............ t- .~d--_
181
El problema dual correspondiente, se obtiene asociando a cada
restricc16n una variable y a cada variable una restricci6n, trasponiendo los
renglones y las columnas de coeticientes en ias restricciones. intercambian-
do el papei de los coeficientes de Ia funci6n objetivo y las constantes,
invi:tiendo las desigualdades y minimizando en Iugar de maxim1zar.
Problema prima Problema dual
Restricci6n Variable
Variable Restricci6n
Es decir. existe una variable por cada una de las restricciones primas y
una restricci6n dual por cada una de las variables primas.
La relaci6n prima-dual podemos representarla par media de Ia
siguiente tabla:
PRIMO
;-~~r----------,----~-----.---~,
x,
D y,
u
A
L Y, a,,
Problema Dual.-
(0) Min Zy = b, Y, + b2 Y2 +
s. a.
(1) a,,Y, + a21Y2 + ... + a,, Y, + + am,Ym2!C,
(2) a,2Y, -t- a22Y2 T + a,2 Y. -t- ... + a""2Ym ~ c2
a,l Y, +
a2J Y2 + + a,! Y, + + a,.,, Y., 2: c,
_,..
(n) a,rY, +
a2"Y2 + ... + a,n Y, + ... arrnY"' 2: en
y, 2:0 (i = 1,2,3, ..... ,m)
182
Ejemplo V-1.-
Max Z x ~ 2x.
s. a
X. S 600
XC s 300
3x. + 4x 2 s 2400
X 2: 0 (j=1.2)
el problema dual correspondiente sera
s a
y 3y3 2: 1
y2 'j- 4y3 2: 2
y 2: 0 (i = 1' 2, 3)
PRIMO
183
Ejemplo V-2.-
Problema primo.-
Max Z - x. - ' x
s a
X, - 3 XC 4x_ .)
-- 5
X, - 2 X< s; 3
2 )(" X > 4
(x., x 2 ) ~ 0 x 3 no restringioa en s1gno.
Problema dual.-
184
Ejemplo V-3.-
Min Z, X. ·~ X X,
s.a
X. -3 x2 ~
4x 3 = 5
X, -2 x2 :::; "
0
2 X0 --
X ~ 4
(x. X,) - 0
s.a.
- x, ... 3 x 2 - 4x 3 = - 5
x. ·- 2 x2 :::; 3
- 2X + X3 :S- 4
obtenemos el dual:
Min-Z )
s.a.
-
Y, + y2 ~ -- 1
3y, -· 2y2 - 2y3 ~ - 1
- 4y, -r y3 = -1
(y2, y3) ~ 0;
Max -Z, = - x, - X2 - x3
s.a.
- x, ... 3 x2 - 4X 3 - 5
X, - 2 ~ :S 3
- 2x 2 + X3 :s -4
(x,, ~) ~ 0;
185
V-2-. COMPLEMENTALIDAD.-
Lema V-1.-
n n m m
Zx=I c X: :::; )
..... X ~a y = [pues c :<:: I a y en el dual]
1=: ::=. 1
m n
I Y, L a, x:::; [cambiando solo x, y y ]
m n
<::)yb=Z
...._. I y [pues L a, : x s ben el prima ]
1=1 :=1
186
De acue'OO con :a amenor notaci6n. nuestra funcion objetivo, cuando
hayamcs otitenido ia soluc16n CJOtima. puede representarse como
Lema V-2.-
· pues z'n., ,~No. de veces que Ia
ecuaci6n ( i ) se suma a Ia ec.(O).
Como Z, = 0 original mente. su valor Ia
a) z: = ~ z~ 1-1 b.
=·
obtiene de las sumas de las
ecuaciones ( i ) que se le hacen hasta
llegar a Ia soluci6n optima.
b) z.*I = \~
._.
z*c-h a' , par<:~ 0 ·.~ '.2 . n)
1=1
numero de veces que Ia ecuaci6n ( i )
se sumo a Ia ecuacion (0) par el coef.
de x: en las ecuaciones (1) a (m).
187
Tan p~·onto hayamas encantrado Ia .so!JC:~]rr uptima (usando Ia rutina
Simplex) para el problema prima entonces estamcs en posicion de obtener
Ia solucion optima del dual El valor de Ia variable dual ( i ), es igual al
caeficiente de Ia variable de halgura ( i ) del problema prima en Ia ecuaci6n
(0). Asi
y* = z·" para (i = 1.2 .... m)
Prueba.-
por 10 tanto
2, z~ + b ~ .L Y, b
( j = 1, 2 ..... n)
Prueba .-
a) Ellema V-1 dice que Z, :S ~· Ia cual implica que existe una solucion
optima finita para ambos problemas (pues max Z, :s: Zy y
minZ,:::: ZJ
b) Par el teorema V-2: ( y· = z· n ~, ) tenemos que:
'"
... ( {3)
de (a)y(f3):
z* = z*
' y
189
Teorema V-4.-
Prueba.-
m
z* = .2: z*o-h a ,, (j=1.2 . . n) [ por ellema V-2.b ]
Por definicion:
r'1
y~ +
'
= 2:=1 a , Y,* - c,
' '
({3) para (j=1, 2, ... ,n)
Sustituyendo ( a ) en ( f3 ):
v:-t, = (z"J -c:) L.Q.Q.D.
Corolario V-1.-
Prueba .-
Si tenemos que x; >0 (k = 1,2, ... ,n) esto implica que x; es una
variable basica y por lo tanto (z;-cJ = 0
Si tenemos que x: +. > 0 (i = 1,2, .... m) esto implica que x:+ es una
variable oosica, par lo que: z:-h = 0
!90
-~ '
t:' . .--., -
~ .
__ : ,__, -'.. ~
Teorema V-5.-
entonces:
z
X
= z y
es decir:
\~ CJ x, = 2: b, Y,
=·
Prueba.-
Usando Ia demostraci6n del lema V-2 y per el teorema V-2, pero sin
hacer aiusi6n a Ia optimaiidad, obtenemos que :
rn
=
z, =I 1=1
zr 11 b,
de donde:
z '
= z I
191
Ejemplo V-4.-
lterac16n
Sol bas:ca factlb!e so: basica :Jua! z
[Ecuac16n 10 I] complerr<e .taria.
--------------
iCOO,-'-;!
(0 .0 .600.300.2400\
[Z 600y.- 300y _ 0 0
[Z, -x. -2x 2 =OJ .:_240Cy =Q] ~
3
(0,300,600 0. 1200)
[Z-x. +2x 4 =600] (0.2.0.-1 .0) 600 600
---+--------------~-- ------- -----~----------1
!
z·. = 1000
= 1000] '
Z* X
192
V-3.- METODO SIMPLEX DUAL.-
Este metoda trata directamente con soiuc1ones basicas no factibles
mejores que el 6ptimo y trabaja para obtener soiuciones basicas factibles
Las soluciones complementarias. en el problema dual. son basicas tactibles,
pero no 6ptimas y se mueven hacia Ia optimaiidad.
Aplicamos este metoda cuando, para obtener una solucion basica
factible inicial. es necesario 1ntroducir muchas variables artificiales (cuyo
valor debera llegar a cera), razon par Ia cual puede ser mas conveniente el
comenzar con una soluci6n no factible · mejor que Ia optima' y aplicar el
Metoda Simplex Dual.
Otra aplicaci6n frecuente, es en el caso del analisis post-6ptimo,
cuando se introducen cambios en el modelo que hacen que Ia solucion
6pt1ma obtenida para el problema original ya no sea factible para el revisado,
siendo entonces mas rapido aplicar el Metoda Simplex Dual, partiendo de
esta ultima soluci6n, que utilizar el Metoda Simplex para resolver, desde el
principia. el problema revisado.
193
su coeficiente en Ia ecuaci6n (Ol entre el valor absoluto del coeficiente en esa
c c
ecuaci6n. As!. tenemos que: 1 --"-T = min T::-:"1 para a, < 0. Esto es
i a.l I a. I
equivalente a determinar Ia variable 'basica saliente en el problema dual. El
coeficiente en Ia ecuaci6n ( 0 ) que se hace cera primero. corresponde a Ia
variable en el problema dual cuyo valor se hace cera primero al incrementar
el valor de Ia variable entrante.
Paso 3.-
Paso 4.-
Ejemplo V-5.-
Max Z "' - 600x, - 300 x? - 2400x 3
sujeto a:
X, + 3x3 ~ 1
X
2
~ 4x 3 ~ 2
X; ~ 0 (i=1,2,3)
N6tese que el dual de este problema es nuestro ejemplo clasico, con
!a excepci6n de que aqul estamos Max. en Iugar de Min
La soluci6n sera:
194
Notese que Ia funci6n objetivo parece indicar optimalidad. pero nuestra
soluci6n noes factible. ya que
X, =0 X" =
x2 ~ 0 xs -2
x 2 =0 z = 0
vemos que:
y obtenemos
X, =0 x4 = - 1
X2 =2 x5 = 0
X3 =0 z =--600
en este paso:
pues (b = min b)
bi <0
Xe =X 3 pues (1200/3 < 600/1)
195
X~=O x:= 0
x; = 2/3 x: .-c 0
x; = 1,3 z; = --1 ooo
y; = 400 x:=o
y;=300 v;=o
v;=200 z;= 1ooo
solucion que es igual a Ia obtenida usando el metoda Simplex para el dual
de este problema.
Ejemplo V-6.-
x4 + x,., + xs =-3
--X, - x2 x:; x4 + x_0
+Xl -- 2
X, + x2 -2x 3 - 2x 4 - 3 x5 -x" 4
X =X
' 5
xe --x2
196
z - 5x. - 3x 4 - 7x 5 -- 2x_ -6
"
-1,2x_ ,- x2 - 1 2Xo --12x - 1 2x_ - 1 2x 6 - 3 2
-3/2X, - 1;2x_ - 1 '2x 4 -1 2X 5 - 1 '2x - x_ = -- 1 2
' "
3.2 X, -5/2 x 2 - 5.'2x" - 5/2x 5 - 1 t2X 5 -X = 512
"
x, -~x 7
Xe =X 3
z + 5X 1 ~3x
4 + 7x 5 -2x.c =-6 rp '
-- 2x, 'X_<' - xe + x7
3x. - x3 - x4 - xs - x6 -- 2x, oc
- 2x, - x2 x6 +x7
X. r 2x2 -t x3 -x4 - xs x6 3
- X. -r 5x 2 -5x 5 - 2x6 +Xa = 10
x.=x,
197
Xe =X 3
X2 =x,
z -r 5x. ~ 3x. 7x 5 -+ 2x 5 = -6
-2x. -r x2 - X-0 ~
x7
3x. ~x3 - x4 - x5 XC0 - 2X 7
Primo
Max Z,
s.a.
X, 2X 2 + 3X 3
-+- :52
3x. - ~ + 2x3 :59
-2x, + ~ -r 8x 3 s16
(x, x2 , x3 ) :::: 0
3x, + x2 +2X,
"
-2x, + x2 +8X3
Xe =X3
X5 = x 4
X =X..
e "
X5 =X 3
198
z X, x3 -L
2x" -
4
( Y:o I 1.2x, ~-x_
"
- 3i2x -1,2x,
( y2) 5 2x, - 1/2X. 0
- 1:2X~ - Xc -: 8
( y3 )· -- 5/2x. ~ 13/2x3 - 1 '2X, -x6 =15
Dual:
Min Z = 2y, -L
9y2 ~
16y3
s.a. '
y, + 3y2 - 2y3 ~1
2y, + y2 y3 ~
~4
3y, -,- 2y2 T 8y3 ~5
( y, y2' y3) ~0
Ys = Ye
Y. = Y,
z + 23/3y2 + 32/3y3 + 2/3y6 = -10/3
(X, )
- 7/3y2 + 14/3y3 .,.. y4 -1/3y6 = 2/3
( )(2) 1/3 y2 +13/3y3 +y5 -2/3y6 =- 2/3
( x4) y1 ~- 2/3y2 + 8/3y3 -1/3y6 = 5/3
Ys = Ys
Y. = Ys
z + 8y2 ~-
15y3 + Ys =-4
( x,) - 5/2y2 -,- 5/2y3 -+- y4 -1 /2y5
( )(3) - 1/2y2 - 13/2y3 -3/2y5 +Yo
I x4) Y, + 1,'2y2 .,.. 1/2y3 -1;2y5 2
199
V-4.-INTERPRETACION ECONOMICA.-
La reciente introducci6n de Ia programaci6n lineal en Ia econom1a
parece ser un anacronismo. Parecera 16gico que hubiese comenzado al-
rededor de 1758, cuando los economistas comenzaron a describir los
sistemas econ6micos en term1nos matematicos. Un crudo ejemplo de un
modele matematico de programaci6n lineal puede ser encontrado en el
Tableau Economique' de Quesnay. Sin embargo, no fue sino hasta los 1930's
en que se comenz6 Ia explotaci6n del modele de tipo lineal en economfa.
La mayor parte de los economistas matematicos se ocuparon con el
analisis de problemas te6ricos, asociadas con Ia posibilidad de equilibria
econ6mico y eficiencia distributiva, bajo condiciones competitivas y
monopolfsticas. Para tales estudios. encontraron el usa de funciones con-
vexas clasicas, con derivadas continuas. mas convenientes para las
demostraciones de condiciones de estabilidad que utilizando funciones
basadas en desigualdades lineales.
Hasta esos dias, el mundo econ6mico usaba el modele econ6mico
para describir, en una forma 'cualitativa en Iugar de cuantitativa", las
supuestas interelaciones dentro de un sistema.
La inspiraci6n del modele general de programaci6n lineal, fue com-
pietamente independiente a los desarrollos anteriores y tuvo un prop6sito
diferente. Surgi6 de las necesidades de programaci6n empiricas de Ia fuerza
aerea y de Ia posibilidad de generalizar Ia estructura practica simple del
modele de Leontief para este prop6sito.
La mayor contribuci6n de Leontief, fue Ia construcci6n de un modele
cuantitativo de Ia economia americana, con el prop6sito de trazar el impacto de
Ia poHtica gubemamental y de las tendencias del consumidor sabre un gran
numero de industrias, embebidas en una compleja serie de interrelaciones.
Cabe hacer notar que el obstaculo para el desarrollo de modelos
cuantitativos era Ia ausencia de computadores, lo cual limitaba el modele a
20 ecuaciones con 20 incognitas.
Para 1940, el analisis de regresi6n se empezaba a utilizar para medir
ten6menos econ6micos. Este heche di6 nuevo impulse a los modelos
cuantitativos.
El punta focal del anal isis de consume producci6n. es una disposici6n
de coeficientes llamados 'Ia matriz consumo-producci6n'' 6 'tabla
econ6mica·'. Una columna de esta matriz representa los requerimientos de
consume de varies artfculos para Ia producci6n de otro particular articulo
con valor de un peso.
200
Ex1ste exactamente una columna para cada articulo producido en Ia
economfa
Vemos entonces que Ia producci6n de un articulo. corresponde al
concepto de una actividad en el modele de programaci6n Lineal .Si los
factores de consume que aparecen en un rengl6n de Ia matriz se multiplican
par el total correspondiente a Ia producci6n de las industrias que compran.
los totales horizontales representan Ia distribuci6n del valor monetario de las
compras entre las industrias proveedoras. De aqui que ei modele haga
posible, no solo el determinar Ia relaci6n de producci6n para satisfacer Ia
demanda directa, sino que tam bien traza los efectos indirectos de los gastos
gubernamentales sabre cada industria.
En 194 7 T C. Hoopmans atrajo Ia atenci6n de los economistas hacia el
potencial de los modelos de programaci6n lineal, usando para ella los
modelos de transporte.
En 1951, Dorfman expres6. en terminos de programaci6n lineal. Ia
teoria econ6mica de una empresa bajo condiciones competitivas y
monopolisticas y compar6 Ia aplicabilidad de esta teoria con el tradicional
analisis marginaL
El numero de aplicaciones practicas continua creciendo y se usa Ia
programaci6n lineal para estudiar industrias especificas.
Mirando ahara en nuestro modele dual, podemos distinguir los siguien-
tes puntas
C esta dada en pesos par unidad de actividad ( j ), pues dijimos que
era el costa (o ganancia) de operar Ia actividad ( j ).
~~ esta dada en unidades de recurso ( i) pelf una unidad de actividad ( j ),
par lo tanto y debera estar dado en pesos par unidad de recurso ( i ), puesto
que:
I a,l Y, ;::: c,
.=1
Debido a las unidades que tiene y , se le conoce como 'el precio unitario
del recurso ( i )", esto es, Y, representa el valor IMPLICITO que el recurso tiene
para el que lo utiliza.
Asi, a'I y ,cuyas unidades son:
I
201
Tenemos asi que~ a. y representara el costa implfcito total de operar
!a actividad ( j ).
~ a. Y, ~ ci
=·
es decir, que cuando utilizemos nuestros recursos en Ia forma mas eficiente,
sera cuando como maximo podamos obtener c1 , pues si no utilizamos bien
el recurso, su costa implicito sera mayor que Ia ganacia o costa c
debe ser fijado de tal forma que se minim ice el costa implfcito total al usuario,
par lo tanto, y,* representa el valor implicito real por unidad del recurso
respect iva.
El valor implfcito del recurso ( i) es cera ( y,* = 0) cuando el suministro
de ese recurso no se agota debido a las actividades (esto es, xn~• >0)
(coroiario V-1).
Si consideramos las leyes de Ia oferta y Ia demanda en economfa,
vemos que cuando un articulo tiene una oferta excesiva, su precio debe ser
cera; en este caso tenemos lo que se conoce como artfculos libres (par
ejemplo el aire).
El corolario V-1, tam bien nos indica que el valor implfcito de los recursos
necesarios para una unidad de actividad ( j ), iguala a Ia ganancia unitaria,
siempre que Ia actividad ( j ) se opere a un nivel positive ( x.* > 0) .
m
202
si Ia cantidad de recurso ( i) fuese incrementada (disminuida) sabre un rango
dado (al rango de b sabre el cualla base optima original nose cambia y par
io tanto y* = zx -r nose cambia)
V-5.- CONCLUSIONES.·
De todo io vista con anterioridad. ha quedado clara que al resolver un
problema de programaci6n lineal. podemos escoger entre trabaJar con el
problema orig1nai (prima) o con su dual Debido a que. como regia general,
el numero de iteraciones requerido para resolver un problema de
programaci6n lineal es igual a una o una y media veces el numero de
restricciones. nosotros podemos, mediante una selecci6n adecuada, facilitar
Ia computaci6n, especial mente en aquellos casas en que existe una marcada
diferencia en el numero de renglones para cada uno de los dos problemas.
Ejemplo V-8.-
203
14.5
X •
.V:2~
(\ ~ HI
8.5 14.5 X 1
-2.5
Otro problema estriba en que Ia soluci6n redondeada puede estar muy
lejos de Ia solucion optima de enteros.
Ejemplo V-9.-
2 3
204
V-7.- EJERCICIOS PROPUESTOS.-
Problema 1.-
AT .1 0 .05 07 $0.21
BT .24 .12 .09 $0.48
CT .25 .21 .12 $0.60
Problema 2.-
205
$ 15 00 en A. S 30.00 en B y $ 20.00 en C El material y otros costas
directamente relacionados con !a produccion de estos modelos, son de
$100.00 para el estandar y S 150.00 para el modificado La compar1fa planea
vender el modelo estandar en $ 300.00 y el mod1ficado en S 400.00. El
departamento de lnvestigacion de Mercados de Ia companfa reporta que no
se debe esperar vender mas de 2000 unidades del estandar y 1500 unidades
del modificado.
Formule este problema como uno de Programacion Linea! si se desea
maximizar Ia utilidad. Formule el dual y de una interpretacion. De Ia
formulacion del prima y el dual en un solo tablero.
Problema 3.-
Una comparifa produce un articulo cuya demanda varfa mes con mes.
La materia prima y Ia disponibilidad de mana de obra presentan variaciones
estacionales. La informacion necesaria se presenta a continuacion.
I Coste de Mano de
Obra Dlsponibilidad de i Demanda Precio de
mes I ($/ton de art.) Mat Prima I del Articulo Materia Prima
(ton) I (ton) ($/ton)
I Tiempo Tiempo :. I
I
extra I
I normal
206
producir 12CJ y 600 tons. por mes durante t1empo normal y tiempo extra
respect1vamente En ios meses 4, 5. 6. 7. 8. 9. estas capacidades de trabajo
son 800 y 500 tons. respectivamente La producci6n de un mes puede ser
vendida en cuaiquier tiempo del siguiente mes o mas tarde. Los costas de
almacenaje son de$ 10 par tonelada de producto almacenado de un mesal
s1guiente La materia pnma no puede ser almacenada. par lo que esta debe
utilizarse en el mismo mes en el que se adquiere. Las operaciones empiezan
en el mes 1, con una existencia de 50 ton. del producto. AI final del mes 12,
Ia campania debe tener una existencia deal menos 50 toneladas del produc-
to. Determine una programaci6n optima de Ia producci6n, a traves de Ia
formulaci6n de un modelo de programaci6n lineal que minimice el costa total.
Formule el dual. De una interpretacion del problema dual.
Problema 4.-
i
Fabricante Fabricante Fabricante \ Fabricante ! Fabricante
1
1 i 2 ' 3 i 4 . 5
--------~--~ r
Planta 1 $0.5 i $07 $1.1 $1.5 j___
i
$1.6
Planta 2 $0.8 $0.6 $1.0 $1.2 $1.5
Planta 3 $1.0 $0.9 $0.9 $1.0 $1.6
207
Problema 5.-
Problema 6.-
208
Nota Se considera que el oenoao 1 sigue inmediatamente
desoues del pe r/odo 6
Cada persona trabaJa 8 horas consecutivas. Sea x1 ei numero de
personas que empiezan a trabaJar en el penodo t todos los d1as. El depar-
tamento de policia desea obtener una asignaci6n diaria que emplee el menor
numero posible de policias. tomando en cuenta que se cumplan todos los
rsquerimientos.
Formular un modelo de programaci6n lineal. de tal manera que se
obtenga una asignaci6n optima Exprese su dual y trate de darle una
interpretacion
Problema 7.-
Min Z = 2K -
3X~ + x3 - x4
<
s.a.
3x. -r
x2 - 2x 3 t- 3x 4 10
-t- 2x 2 -r 4X 3 - 2x 4 $ 12
2x, - 3x 2 + X~ ?:- 6
"
6x, + x, $ 15
X ?: 0 ( i = 1,2,3)
'
Problema 8.-
209
--, EJERCCiOS ?PC;:::l~...~t=S res
Problema 9.-
Considere el problema:
Max Z = 6x. + 8 ~
s.a.
5x, + 2x 2 :s: 20
x. + 2x2 :s: 10
x, :::: 0 (i = 1,2)
a) Obtener el problema dual.
b) Resuelva graficamente el problema prima y el problema dual.
ldentifique las esquinas que den soluciones basicas factibles y no
factibles para ambos problemas.
c) Use Ia informacion obtenida en el inciso {b) y construya una tabla
en Ia que se listen las soluciones basicas para estos problemas
(prima y dual).
d) Resuelva el problema prima par el metoda Simplex. Despues de
cada iteraci6n identifique Ia soluci6n basica factible para este
problema y Ia soluci6n basica complementaria para el problema
dual.
Problema 1o.-
210
a) Obtenga ei 6ptimo usando ei metodo Simplex
b) 2.euai es el modelo dual correspondienre"
c) 6eual es Ia solucion optima parae! problema dual? Use ei Metoda
Simplex Dua!
Problema 11.-
211
CAPITULO VI
METODO SIMPLEX
REVISADO
Vl-1.- PROCEDIMIENTO.-
El Metoda Simplex Revisado es una version especial del Metoda
Simplex y tambien fue desarrollado por G.B. Dantzig y sus colaboradores de
Ia Rand Corporation.
En tanto que cada iteraci6n del Metoda Simplex requiere que una nueva
tabla sea calculada y almacenada (proceso poco eficiente al ser implemen-
tado en una computadora digital, par Ia gran cantidad de memoria requerida
para almacenar cifras que no son utiles para decisiones posteriores), el
Metoda Simplex Revisado, diseriado para obtener exactamente las mis-
mas casas que el Simplex original, solo calcula y almacena las piezas de
informacion relevantes en cada iteracion, por lo que este procedimiento, al
ser mas conveniente, es el que se utiliza en las computadoras digitales. Con
este metoda se optimiza el tiempo de maquina y se permite, dada una
capacidad limitada de memoria, Ia solucion de problemas mas grandes que
con el uso del Simplex original.
La idea basica detras del Metoda Simplex Revisado, es el uso de un
con junto de numeros llamados "multiplicadores simplex" o 'precios ', junto
con Ia inversa de Ia base, utilizados para generar directamente, partiendo de
las ecuaciones originales, Ia informacion precisa requerida para las
decisiones de esa iteracion. Si consideramos que es comun que 90% o mas
213
de coeficientes cera existan en los datos originates. es evidente Ia bondad
del metoda al requerirse muchas menos multiplicaciones
La 1nformaci6n relevante en cada iteraci6n es Ia siguiente:
a) Los coeficientes de las variables no basicas en Ia funci6n objetivo
(ecuaci6n 0).
b) Los coeficientes que tenga Ia variable que entrara en Ia base (xe),
en todas las otras ecuaciones de las restricciones (ecs 1.2 ..... ,m)
c) Ellado derecho de las ecuaciones (b) para (i = 1 ,2, .... ,m).
A pesar de Ia eficiencia de este metoda, para calculos manuales se
recomienda el usa del Metoda Simplex Original.
El Metoda Simplex Revisado utiliza explfcitamente operaciones con
matrices, par lo que describiremos el problema usando notaci6n matricial
Existen dos versiones de este metoda. cuya (mica diferencia es Ia
manera en que Ia matriz inversa se almacena. Estas son:
1) La forma explfcita
2) La forma de producto
Recordamos que nuestro modelo matematico de programaci6n lineal.
una vez introducidas las variables de holgura, puede ser representado como:
(1) ... Max z = ex 6 ( Max z- ex = zo J
sujeto a:
(2) '\'
~
X
J
p )
= p 0
=1
p
I
= p 0
=
y: [a ,1 • b;. c 1
] son constantes.
El Metoda Simplex requiere que m de los vectores P sean inde- 1
pendientes.
214
Sea
iP p 1
, . s.-r, J
tai conjunto de vectores 1ndependientes eilos forman una base B del con-
junto de vectores (P . P2 . P_;
lB. B.•...
I .<'
B."'
I B2. B22 ... B2,..
(4) ..
II .
XBM J
debido a que todas las variables no basicas valen cera y las eliminamos de
consideraci6n.
En general, una forma can6nica se obtiene premultiplicando ambos
Iadas de las ecuaciones (2) por B-' . esto es:
+ (B- 1Pn)x • = B- 1p o
215
cc;ocEDif.l E~ ~c-
donde:
s~·p = P
(7) ... s~·p = P (j = 1, 2, .... n)
0 v
entonces:
x8 = s~,p = P
0 0
teniendo en cuenta que estos coeficientes se taman con el signa que tengan
dellado derecho de Ia ecuaci6n (0).
Si premultiplicamos ambos Iadas de Ia ecuaci6n (6) par y , ob-
tendremos:
donde:
( y P ) = cte.
216
en especial
de forma que Ia ecuaci6n (8) tiene los mismos coeficientes para las variables
bas1cas que Ia ecuaci6n (1) • de aquf que sumando Ia ecuaci6n (8) a Ia (1)
eliminamos !as vanab!es basicas, obteniendo:
donde•
c.. ==-c+ ;; P, == -c. + y(B- 1 P)
=-C + (I B- 1 ) P
=-C.+ .7 p I
y en particular:
( ver lema V-2b )
217
c' e = min c
e incrementando el valor de xe tanto como sea posible, hasta el valor
xe = x* e' en donde algun componente p del vector ( Pc - P-e xe ) cambie
de signa, mientras que los demas siguen no negat1vos.
manera:
I
b',
b' 2
l
raa ·1 2e
b' b'
p = p = x* = _::_e_=min '
l
0 e e a' pe a ,e >oa',.
b' n
a"'
8-' p
•
= p!
8-' P0 =P 0
p = 8-' p
entonces:
P=BP I
218
de donde:
'a .e
a'2e
(10) ... P.
I
I . i
La·,. J
.;. P Brn a ,....,e
donde:
,
{ a ' '• , a ,2e , . , a . ., • .f
'
es Ia representaci6n de P. en terminos de B.
Despejando PBp de Ia ecuaci6n (1 0):
a'
. - pe a' .-PB m ~
a
pe pe
es decir:
B*K
en donde:
k =
p a' pe
219
de tal forma que Ia relaci6n entre Ia vieja' y Ia nueva' base esta dada par:
B = [P,,Pa 2 , .,P 8 "'] = B*[U,U 2 , .. ,K, ... U"']
multiplicanda par s-~ a Ia derecha y par (8*)-' a Ia izquierda en ambos Iadas
de Ia ecuaci6n, obtenemos Ia relaci6n entre el inverso de Ia nueva base y el
inverso de Ia base previa, es decir :
en donde:
l km
Definici6n Vl-1.-
Teorema Vl-1.-
220
Partiendo de
p = B-' p
[U, U2 , ... , K, .. , Uml = [U 1 . U2 , ...• Up, ... , U"'] + [0, ... , K- Up, ... , OJ
y de Ia ecuaci6n (11):
(B*)- 1 = [U, U2 , ... , UP .... , U,]B- 1 + [0, 0 .... , K- Up, .... 0] B-'
par lo que:
donde:
K = K- Up= k~ k p -1
km
221
Llamando 8 a los renglones de s-·. es decir:
(B* )- 1 = B- + [0, 0,
1
.. , K , ... , 0] f3
p
L..
en donde se ve que (8*)·' difiere de B-' par una matriz K {JP, que es el
KfJ p =
222
Sea
B-' = [ p] I
entonces el elemento ( i, j) de K /3 0
es simplemente k f3P, De aqui que para
o=: g-1 p
I
+ ( K f3 p ) p I
= P; + Ka'P,
donde a'P: es el valor del componente pen Ia representaci6n de P, en
terminosde B"'_Asi, el nuevo Pt
ditiere del anterior Pen un vector pro-
_Q_orcional a K; el factor de proporcionalidad es el componente "p'' de
p
Las relaciones (11) y (16) son dos tormas de expresar Ia nueva inversa
en terminos de Ia vieja. La relaci6n (16) requiere, en general, de m 2 cambios
en los componentes deB-' , mientras que Ia (11) muestra que el proceso de
obtener (B*) 1 de B- 1 , multiplicandola porIa matriz elemental (12), requiere
s61o del conocimiento de los "m" componentes del vector K y su localizaci6n
en Ia columna 'p" de Ia matriz, con lo que se consigue un notable ahorro de
memoria.
Se ve que puede ser computacionalmente conveniente el representar
Ia inversa de Ia base como el producto de matrices elementales. Comenzan-
do con Ia base inicial, Ia cual es una matriz identidad, Ia inversa de Ia 2a. base
es tambien una matriz elemental, pero Ia inversa de Ia 3a.base y las sub-
secuentes, excepto en casas fortuitos, no son matrices elementales. Sin
embargo, podemos redefinir el problema despues de Ia 1a. iteraci6n, de tal
manera que Ia 2a. base sea tratada como una base original para el nuevo
problema. Con esto, Ia 3a. base y su inversa son matrices elementales, en
terminos de Ia 2a. En Ia misma forma, cada base es una matriz elemental en
terminos de Ia inmediata precedente.
223
Teorema Vl-2.-
MaxZ = x 1 + 2~
s.a.
x1 :S 600
~ :S 300
3x, + 4~ s 2400
(x,, x 2 ) ~ 0
aquf:
Y Po= [
2400
~~]
1a. lteraci6n:
Paso 1.-
=
100]
B
[0 1 0
0 0 1
224
Paso 2.-
Obtenemos ;' •
/, = r o. o. o 1
Paso 3.-
Obtenemos B-' :
Paso 4.-
y,B~' = (0,0,0]10
11 0
1
0] = [ 0, 0, 0 l
0
,o
L
0 1
Paso 5.-
P. P2 c, c2
Paso 6.-
Paso 7.-
Determinamos Pe :
225
Paso 8.-
Encontramos P0 :
I1 o o l r 6oo
lo 1 o i : 300
l 600
300 I
l
lo o 1 J l2400 J [ 2400 J
Paso 9.-
Paso 10.-
,------
I 8 _,
B. p = p
I ,,v {xg/, ( p 0 ),
1
----~-----
e 2 (xe ~ x2)
! x3 x4 xs
0 x3 600 0 0 0 0 0
0 0 p:vote (x, = x 4 )
II 00 )(4
)(5
300
2400
0
0 0
0
0 0 4
Paso 11.-
Actualizando, tenemos:
226
'.'::::.TC~O Sr".1P:._E;>., Rl? JrS;. :)C
-,
a', l a' 12 I
a·o: I - a'221 Ql
1 I
K1 =
ro1
a' pe = 1I
a'22 =
- a'3e J 2
4
- l4J
_ a> J
L a' pe a 22
Paso 12.-
2a. lteraci6n:
E, ~ r~ ~ n
Paso 1.-
Paso 2.-
f~l n1 j
Paso 3.-
cs;t,==B;'=E,I=[~
* o
~ ~1~~La o~ ~l=[~o ~ -4 1 1j -4 ~]
Paso 4.-
y2 E, == [
11
o, 2. o1 lo
0-4
0
1 j[
ol
~ ~ o, 2, o I
*
Paso 5.-
227
Ia~1 ~~l
c, c 2 c', c' 4
[ o.2,o J -[1,0]=[-1,2]
'3
L
Paso 6.-
Paso 7.-
P. = P 1 = 82 1 P 1 = E, P. =
*
Paso 8.-
Paso 9.-
*
[L! 0]0 [ 600l
1
300 = [ 6001
2400
300
J 1200 _,
*
~ Xs ( p = 3)
Paso 10.-
Z2 = y2 ( P0 )2 = [ 0, 2, 0] [ :g
1200
l = 600
Resumiendo:
l
r 1 - E1 82 * s2-1 - -
f2 Psf2 ( P0 )2 -K~~- ---- _ _ _ _ p =P (xe=x1)
e 1
I._______ 1 x3 x2 xs
I 0 x3 600 0 0 0 0 0 0
I 2 x2 300 0 0 0 0 0 0 0
II 0 1200 0 4 0 -4 3 pi vote
xs 0 -4'
(xs=xs)
Actualizando:
228
r
I, __
I
a·,.
a' pe
! ,, I l
1
I= I~!
3
I
a' 2e
~-~
1\= - a' pe
Paso 12.-
IL a·.,.
-
1
J
=
l -1
3
0
, -V3]
0
0 V3
3a. lteraci6n:
Paso 1.-
Paso 2.-
y3 = [ 0, 2, , 1
Paso 3.-
( 8 2* ). _, = B-,3 == E2 E, I = r~
* lO
Paso 4.-
229
Paso 5.-
P4 Ps c4 cs c'4 c's
Paso 6.-
Paso 7.-
No se necesita ya.
Paso 8.-
[~
0 0
(P,J,=B:P,=E,E,P,= 1 -113
0
l ['
J 00 1 0][600]
0 300
0 V3 -4 1 2400
=
[~
413 -1/3l
1
-4"3 v~ J
r 600
l300
2400
l l
=
r200
300
L400
Paso 9.-
No se necesita ya.
Paso 10.-
230
METOOO SIMPLEX REVISADO
Resumiendo:
Soluci6n Optima
8::13 l l l
' I
v.b.
__ ! I v.n.b 1
l i
z*= 1000! x*=o \
4 I
aquf:
[A,I]=[~ 2~ 0~ 0~
-1
g1] y p0 =
231
'/1-1. · PROCEDIMIENTO
1a. lteraci6n:
Paso 1.-
Paso 2.-
Y1 = ( 0, 0, 0]
Paso 3.-
B;· = [H ~]
Paso 4.-
Paso 5.-
P, P2 c, c2 c', c'2
Paso 8.-
Paso 7.-
232
Paso 8.-
= B:' po -= ~~0
(Po ) ,
[~
'-
Paso 9.-
(p = 3)
Paso 10.-
Z~ y,P,~[O,O,O) [~] =0
Resumiendo:
B, B-,1
i
Y, {xa}, (~), -·----- P=P (x 0 = x2 ) II
e 2
! x3 x4 X
,,~
r
i 0 x3 3 0 0 0 0
i
I
I 0 x4 4 0 0 0 0
,....,, I
i
I
I
0 xs 2 0 0 0 0 ~'2'
'
'"--.-'
'
__j
Paso 11.-
a' 1e
a' pe
1 r-·:,
a32 l
2.1
2
I
K1 =
a'2e
a' pe = 1-::: = 2
--,
1
a'pe J l a'32
J
l ij
233
Paso 12.-
2a. lteraci6n:
Paso 1.-
Paso 2.-
y2 = [ 0' 0' 3]
Paso 3.-
( B7 )- = 8; = E I = [~
1 1
1
* 0
Paso 4.-
Paso 5.-
P1 P5 c1 c5 c' 1 c' 5
234
Paso 7.·
pe = P = B-
1 2
*
1
P 1 "" E1 P1 = [~
0
0
1 -V2
0 V2
-v,lJ [ 1] = [ ~ 4
-1
912
-V2
l
Paso 8.·
:~ ][;] = [~]
0
1
0
Paso 9.-
X4 ( p = 2)
Paso 10.-
"'B-1
82
y1 {xs}2 ( J5;;)2 _ _._E1___ - - - 2- - P=P (x 0 = x2)
·----rc-'' , Xa x4 x2 • 2
Paso 11.-
Actualizando:
235
'JI· 1.· PROCEDIMIENTO
-,
a' 1e a' ~2
a'pe a'22
-1
3
1 1
~== a' = -
a' 22 == -2
9
l-::; j a' 32
a' 22
-1
9
Paso 12.-
~ [~
n
-V3
zig
E,
~
38. lteracion:
Paso 1.-
Paso 2.-
y3 = [ 0, 2, 3]
Paao 3.-
Paso 4.-
"' ~ y, B;' ~ y, E, E, ~ (
0
0, 2, 3] [~ -~ ~ [~ l
236
:1 -1 '3 -1/31
= [ 0. 2, 3] i0 2ig -tg! = (0. 7i9, 1CYgj
IO
L·
1/g 4'91
J
Paso 5.-
P. ps
1o ol
[0, 7/9, 10/9]l~ ~J [ 0, 0 l [ 7/g . 1(}9 J
Paso 6.-
Paso 7.-
No se necesita.
Paso 8.-
ollr01
0
2
-11 ]
01 -112 [31
4
1 0 0 1/2 2
= 1 -113 - 113]
0 zig -11g
[ 0 1/g o/'g
4
[3l = rlz;;
2 !
J
1
4"3
l
Paso 9.-
Ya no se necesita.
Paso 10.-
237
Resumiendo.
s;·- vb ·; n b
Ejemplo Vl-3.-
Max. Z = -3x, + ~ + 7 ~ - 3X 4 - x 5
s.a.
Sx, - 4 ~ + 13x3 - 2X 4 +- X5 = 20
x, - ~ + 5x3 - x4 + x5 8
(x,.~.~.x 4 ,x 5 ) 2: 0
238
con lo que el sistema queda:
Max. -W
Z ;.. 3x. - x 2 - 7x 3 c 3x" - x 5 -~ 0
5x - 4x 2 ~ 13x 3 - 2x, + X0 -"- x6 -~ 20
x. -- x2 ~ 5x 3 ·- X4 - x
5
- x, = 8
FASE 1:
1a. lteraci6n.-
Paso 1.-
B
l1
= lo
0
1
0!o
0 0 1
Paso 2.-
(, [0,-1.-11
Paso 3.-
Paso 4.-
Paso 5.-
[ 0, -1' -1 ] l5r3-1-7 3
-4 13 -2
1]
11 - [ 0, 0, 0, 0, 0] = [ --6, 5, -18, 3, --2 ]
1 -1 5 -1
Paso 6.-
239
Paso 7.-
Paso 8.-
( P,), = B;' P, = [~ !~ l Hl = Hl
Paso 9.-
min. ¢, = ¢ 2 = 20/13 ~X
s = X6 (p=2)
Paso 10.-
w, = y,
Resumiendo:
(P;J, = 1 o. -1,--1 I l
[2~ = -2a
8';1 I
(X 8 = ~) I
I
I o z 0 0 0 0 0 -7
I
0 0 0 ' '
,1~)
pi vote (x ~ x 6 ) i
0 0 0 5
• I
Paso 11.-
a'pe = a' 23 13
7
- ::1•1 pe 13
1 1
a'pe I= 13
a' 3e
a'pe
j - 1~ J
240
Paso 12.-
2a. lteraci6n:
Paso 1.-
j1 -7
r
ol
8=,013 o!
lo s 1J
I
Paso 2.-
y2 = [0,0,-1]
Paso 3.-
Paso 4.-
~] =[0,>'13,-1]
, 7/13
Jr 2 = y2 B;'= y 2 E,::[0,0,-1] 0 V13
[ 0 -5-13
Paso 5.-
l
3-1 3 1 0]
[ 0, 51'13, -1
[
5 --4 -2
1 -1 -1
1
1
1 - [ 0, 0, 0, 0, -1]
0
241
Paso 6.-
d'_0 -8/13 ~ x.
Paso 7.-
Paso 8.-
0 l i 0l
o ! 12o 1
1
j
I l. 8
j
I
Paso 9.-
Resumiendo:
y2 {xs}2 \P~)2
E
--K~
----
z
82
x3 x7
s--'
2
~
p =P (x.
e 5
~-
,x1
0 z 140/13 1 7/13'. 0 -7 0 7/13 0 20/13
I
I
0 x3 20/13 0 '1/13 0 0 13 0 0 1!13 0 1/13
-- I
I
/
pivote i
x7 4/13 0 t-5/13 0 5 0 -5/13 8/131
(xs=x7) I
Paso 11.-
a' pe 8/13
242
:
II a· .•
-- 2013 20
\ a oe &'~3 8
I
a'2e 1/13 1
K2 = = = --
a\ ae 8-'13 8
1 13 13
-
- 8
a' ce 8
Paso 12.-
0
1
-2Qi
-vs
s1
o 13/s J
Ja. lteraci6n:
Paso 1.-
[~ l lo s l B = I110 -7
13 1~
Paso 2.-
i'3 = [0,0,0]
Paso 4.-
2.43
Paso 5.-
3 -1 3
[ 0, 0. 0 J 5 -4 -2
0 0]
1 0 -[0,0,0,-1,-1] = [0.0,0, 1,1 J
[ 1 -1 -1 0 1
Paso 6.-
Paso 7.-
Paso 8.-
[~]
Paso 9.-
Paso 10.-
Resumiendo:
244
B:~
!! f3 {x3}3 ( f5;;)3 E2 83 j
K_c z x3 x5
FASE II:
4a. lteraci6n:
Paso 1.-
Paso 2.-
I4 [ 7 , --1 ]
Paso 3.-
Paso 4.-
n 4 -- 'Y4 B-1 - [ 7,
4 - -
1/a -V8]
1 ] [-Sis - [ 31:2, -~/2 ]
13/a -
Paso 5.-
245
.'' 2 METODO SrMPLEX REV SADO USANOO FASE 1 'f 0 ASE ,r_.
Paso 6.-
c' 2 -9/2
Paso 7.-
Paso 8.-
Paso 9.-
4/7
Paso 10.-
Resumiendo:
[ y4 {x } (x x } i
~~,----8--4------+-~~~-~-----------------~
3
:1
L_ 5
: 112 o 13/8 · 5 -5/8 13/8 7/s\
·, /
(xs = X5 ) 1
-
1
Paso 11.-
a' 1e r_~
r.,]
318
a' pe a' 22 7/a
K4 = = = =
7/8
lS'r
a' pe a' 22
246
Paso 12.-
Sa. lteraci6n:
Paso 1.-
8=113-4]
l 5 -1
Paso 2.-
i's = [ 7, 1 ]
Paso 3.-
Paso 4.-
Paso 5.-
Paso 6.-
247
Paso 7.-
Ya no se necesita.
Paso 8.-
Ya no se necesita.
Paso 10.-
Resumiendo:
r--
I
I Ys E4 B-15
I
{xs}s
(~)5 --- k4 i
~-
x3 x2
~
I 7 12;7 3;7 13 --4 -1;7 4!7
I x3
I x2 4;7 0 I
8;7 5 -1 --5;7 13;7
soluci6n optima
v.b. I v.n.b.
z· =88!7 I <=o
x;*=12;71\ x;=o
x =4!7
2
x;=o
248
C AP IT UL0 VII
249
'.)11.1. !NTHO~·_:c:::::-::::--,
250
introduciendo las variables de holgura:
Z- X, - 2x2 0 .. (0)
_,.
X, x3 600 ... (1)
3x, ~
4~ + X4 = 2400 ... (2)
v.b. v.n.b.
z -+-1/2x 4 = 1200 (0') Z*= 1200, x*, =0
Y,
El dual correspondiente sera:
'I +3Y.,=~,
introduciendo variables de holgura :
zv - Booy, - 2400y 2
y, _,. 3y2
4y2
y Ia soluci6n optima final seria (segun indican los coeficientes de Ia ec. (0)
final en el problema prima):
v.b v.n.b
251
Jli-2. CA\1810 c... C CwANDC x· ES \oC 3~~ C:.o
de (a):
252
(z',- c,) = 112
c, = 1
c. = 5/4
Debido a que :
x, Y,
v.b. v.n.b.
!!00
400
.______ _,_
~-= 1200 y, =0
.
~----
200
Y/= 1/2
.
Y. =0
y.
I
~y, .. ,S.,'4l
y3*=1/4 '
Ejemplo Vll-3.-
253
deseable separarse un poco de Ia soluc16n optima. con objeto de
obtener una soluci6n que requiera menos modificaciones esen-
ciales ante estos camb!os.
b) Los valores de las constantes ( a ), ( c ) yio ( b ), pueden ser
controlables hasta cierto punta, en cuyo caso deseamos conocer
los efectos que resultarian de cambios en sus valores.
c) Aunque los valores de los parametros del modelo no sean control-
abies, las estimaciones de sus valores pueden ser solamente
'aproximaciones", hacienda entonces importante el conocer
para que rangos de estos valores Ia soluci6n obtenida sera aun
optima.
d) Existen ocasiones en que se hace necesario introducir cambios en
el modelo original, ya sea porque se descubren errores u
omisiones, o bien porque disponemos de nueva informacion que
indica que las estimaciones de los valores de los parametres
deben ser revisadas.
Ur 9Z localizados los parametro crfticos, se pueden establecer
metodos :;Stadisticos para determinar el momenta en el que varian estos
parametres y poder hacer los ajustes necesarios en nuestro modelo.
Las tecnicas que veremos en este capitulo, hacen innecesaria Ia
resoluci6n del problema desde e! principia cad a vez que un pequeno cambia
se efectua en el modelo. En Iugar de esto, dada Ia solucion optima anterior
y su correspondiente sistema de ecuaciones. se determina si Ia misma base
sigue siendo optima o no. yen el caso de ya no serlo, usarla como punta de
partida para llegar rapidamente a Ia nueva solucion optima.
Ejemplo Vll-1.-
Max. Z = x, + 2~
s.a. 600
x, s600
400
3x, + 4~ s2400
(x,,x 2 ) so 200
250
"-NA_,sis DE SE~SIBIUCAD v "OST OPHJC
o:
A
u
ca-< -(z-c.)
J
'
I
• aki
y podemos notar que para a' I<! = 0, no afecta el coeficiente.
255
Por lo tanto:
-(z-c) -(z.-c)
max , :"5 il C8 :-:;min · , (para toda j fuera de Ia base)
a' kJ >0 a kJ k a' <0
~~
a kJ
Ejemplo Vll-4.-
soluci6n que noes 6ptima, por lo que volvemos a aplicar el Metoda Simplex:
x. = ~
200
256
ANALISIS DE SENSIBiL.;OAD v POST OPTiMC
x, + x3 = 600
:lS -- 3/4X3 + 1/4X 4 = 150
Ejemplo Vll-5.-
Determinar cuttles serlan los lfmites dentro de los cuales podrfan variar
c2 y c 3 para que Ia soluci6n obtenida en el ejemplo Vll-1 no deje de ser optima.
Para c2
par lo tanto:
Para c 3 :
por lo tanto :
-Vz~6.c 3 ~ oo
600
400
200
Z MODIFlCADA
dejando una forma can6nica para ~ y ~ :
Z + 1/3X4 = 800
x, + x3 = 600
3/4x, + ~ + 1/4X4 = 600
el mismo resultado hubiesemos obtenido si a Ia ec.. (0') final original le
257
Vll-3.- CAMSIO E ..... C .:uANDO x" ES SASICA.·
1
x2
600
Z MODIRCADA
200
0 . 1000 1200
Z MOOIRCADA~
X • =X 1
258
Vll-4.- CAMBIO EN b.-
Supongamos ahara que, una vez obtenida Ia solucion optima, nos
damos cuenta que el valor de b ha cambiado en ~b. hasta un nuevo valor
b. Dado que al cambiar b. por 5. no se afecta mas que ellado derecho de las
ecuaciones, los coeficientes (z'i - c1) seguiran siendo 2=0, sin embargo, el
vector Pc sf cambiara, pudiendose presentar 2 situaciones:
a) Todos los nuevas val ores de los componentes de Po siguen siendo
positives, en cuyo caso Ia soluci6n anterior sigue siendo optima.
b) Algun(os) nuevo(s) valor(es) de los componentes del vector Po es
(son) negativo(s), en cuyo caso Ia solucion optima anterior ya no
es factible y habra que aplicar el Metoda Simplex Dual para obtener
Ia nueva soluci6n optima.
De lo anterior, es clara que, cuando cambiemos el valor de b, a ~ , lo
unico que procede investigar es si Ia solucion modificada sigue siendo
factible.
Sea:
a\."., = coeficiente de x, ., en Ia ecuacion final k (k = 1,2, ... ,m)
b'. = Iado derecho de Ia ecuacion final k (k = 1,2, ... ,m)
Dado que xn +, , que es Ia variable de holgura correspondiente a Ia
ecuacion asociada con b,, aparece con coeficiente unitario en ella y cera en
todas las demas, es facil saber que multiples de b . han sido sumados, directa
o indirectamente, al lado derecho de las demas ·ecuaciones y este multiplo
es a·k , n -t-1. para k :;ei.
PorIa tanto tendremos:
m
b: = 2: a:.n+i b,
1=1
(k= 1, 2, ... , m)
259
b' -+ a·, r 6b ~ 0 (i = 1,2, .. ,m)
-b
6 b ~
I a*
' · 0 .,... '
yparaa· •. n •• < 0:
-b'
6b. s ·--·-
a*K, o +1
1
de donde:
-b', b',
max * s min ---
a*K, n +1
a *k, n +' >Q a K, n + I
a*K, n +1 <Q
Ejemplo Vll-6.-
max. [ ~0
0--,~ ~OJ s6bi s oo
- 2400 s 6 b, s co
Ejemplo Vll-7.-
260
por lo que ahara Ia nueva soluci6n es basica no tactible mejor que Ia optima
y debemos aplicar el Metoda Simplex Dual.
En este caso el problema original no tiene soluciones basicas tactibles,
par lo que en el Simplex Dual no hay vanable de entrada.
BASICA.-
Si ahara-
suponemos que tenemos un cambio en el valor de a.I] a
a 1 (~a"= a,- a 11 )U = 1,2, .... n) y que ya contabamoscon Ia soluci6n optima
para el caso original, sera evidente que este cambia solo puede atectar a los
coeficientes de x·. en el sistema de ecuaciones, permaneciendo inalterables,
respecto al caso original, todos los correspondientes a las demas variables.
Dado lo anterior y que x*, es no basica, Ia solucion optima original debera
seguir siendo factible y 16 unico que procede averiguar es si sigue siendo
optima o ya no, para lo cual solo revisamos el nuevo coeficiente de x·, en Ia
ec. (0) final. Si sigue siendo positivo, Ia solucion anterior sigue siendo optima,
pero si ahora ya es negativo, partimos de Ia solucion optima original y
escogemos a x como variable de entrada, aplicando el Metoda Simplex hasta
1
261
Ejemplo Vll-8.-
es decir:
600
400
200
262
como Ia nueva solucion ya no es optima, aplicamos nuevamente el Metoda
Simplex:
'·'· I mb
z· = , 500 x·3 =o
x· , =600 II x· 4 =O
x" 2 =450 ,
En caso de que despues del cambia Ia solucion anterior siga siendo
optima, ningun paso adicional debe llevarse a cabo.
263
siendo positive en todos sus elementos y el vector de coeficientes
correspondiente a Ia ec .. (0) modificada C* tambien es positive en
todos sus elementos; en este case Ia nueva solucion es todavfa
optima. Es decir:
si ~ ~ o y C* ~o :;.solucion optima
Ejemplo Vll-9.-
z- x 1
- ~ o
x, + lS 600
3x1 + ~ + X4 = 2400
264
Z --.- 1/2Xi + 112(- 2)x2 , 1/2x, = 1200
X, 600
3/4x, + [1 + 1/4(-2)]x2 = 600
es decir:
Z + 2x, X4 = 2400
x, ' x3 600
3/2x, + ~ + 1/2X4 = 1200
v.b. I v.n.b.
z• = 2400 I x·, = 0
x· 3 =500 I x·
' 4
=0
• I
X2=1200I
Graficamente tenemos:
_/RESTRICC!ON
'..--/ ORIGINAL
265
..'>~-6 ::.AM8\C E\ a ::u,~";::JO I. t:S SAS-:::;.
'
Ejempio.-VII-1 0.-
z x3 .._ 2x 4 = 14
-1- 1i2X 3 ..,_ 1/2x 4 4
X, -- 1/2x 3 ~ 1/2x 4 = 2
z'= 14 i
.
x· 3 =O
I x·4 =O
y:
LlaIJ ::~ - a 21 ) = ( 3 -
= ( a'J - a ) = ('L1 !j
1) =2
266
z 2(2)x, + x3 -+- 2x 4 = 14
1/2(2)x, + x~ + 1/2x3 .,.. 1 /2x 4 4
[1 + 1/2(2)]x~ - 1/2X3 + 1/2X4 ·""' 2
simpiificando:
z + 4x, x3 -1- 2x 4 = 14
x, + lS -+- 1/2~ + 1/2X4 = 4
2X, 1/2x3 + 1/2X4 = 2
z + 2x3 + X4 = 10 v. b. v. n.b
-X. +X,=2
2 4 6 8 10 12 14
267
Z 1 (-3/4)x 2 + >s -r 2x 4 = 14
[1 + 1/2(-3/4)]x2 + 1/2>s + 1/2x4 4
x, - 1/2(-3/4)x2 - 1/2>s + 1/2X4 = 2
simplificando:
Z + x 3 + 2X4
-3/4~ = 14
5/B>s + 1/2>s + 1/2x4 4
x, + 3/S>s - 1/2>s + 1/2x4 = 2
x, =-vs 1
x,
rX,+(1/4)X 2 =2
....-~___,~---+-+-~~x,
14
268
Vll-7.- ADICION DE UNA NUEVA RESTRICCION.-
Si por alguna omision, o bien porque las condiciones del sistema han
cambiado (p. ej. el gobierno impuso restricciones a Ia importacion de ciertos
componentes, o se presentaron cambios en las condiciones del mercado
que antes no existian. etc.). se hace necesario introd ucir una nueva
restriccion cuando el modelo ya ha sido resuelto. su efecto solamente puede
ser el de eliminar soluciones anteriormente factibles, incluyendo, tal vez, Ia
solucion optima original, sin que exista Ia posibilidad de anadir alguna nueva,
par lo que el valor de Z debe permanecer constante o bien disminuir.
Las 2 posibilidades que se pueden presentar en este caso son:
a) La solucion optima original es factible para Ia nueva restriccion, en
cuyo caso seguira siendo optima para el sistema modificado.
b) La solucion optima original no satistace Ia nueva restriccion. En este
caso. Ia nueva restriccion debe ser agregada al sistema de
ecuaciones final obtenido para el problema original, para lo cual y
en el case de una desigualdad. se ie agrega a esta Ia variable de
holgura correspondiente. misma que sera Ia variable basica para
esa ecuacion desde ese momenta y se eliminan de ella todas las
demas variables basicas, con objeto de tener un sistema
canonico.
Como Ia solucion original no satisface Ia nueva restricci6n, esto se
vera reflejado en que, al obtener el sistema can6nico, Ia variable de
holgura recientemente introducida (que actua como variable
basica), tendra un valor negative. Como Ia funcion objetivo sigue
indicando optimalidad y Ia soiucion es ahara no factible,
deberemos aplicar el Metoda Simplex Dual para obtener Ia nueva
solucion optima.
En caso de que Ia nueva restriccion amerite Ia introducci6n de una
variable artificial, esta debera tamar las funciones de variable basica
para esa ecuaci6n, pero debera penalizarse su existencia dentro
de Ia base mediante el Metoda de Ia Gran M, por lo que al formar
el sistema canonico deberemos eli.minarla de Ia funcion objetivo.
Ejemplo Vll-11.-
269
Tt!nfamos !liJ8 Ia soluci6n optima era
x· -= 0 ' x· 2 cc 600 ; x· 3 = 600 : x· 4 = 0 y z· = 1200
de donde es ciaro que esta no satisface Ia nueva restricci6n (pues x· 2 = 600),
por to que deberemos agregar at sistema final de ecuaciones Ia sigwente
ecuaci6n:
.(3)
donde:
X5 = variable de holgura.
obteniendose:
Z + 1/2x, + 1/2x4 1200 ... (0)
x, + x3 600 ... (1)
3/4x, + ~ ~ 1/4x4 600 ... (2)
x2 -r- x5 = 400 ... (3)
v. b. v. n.b
z + 1/2x, ~
1/2x4 1200
x, ·+- XJ 600 Z*= 1200 x, = 0
3/4X. + X? + 1/4x 4 600 X3 = 600 x4 = 0
- 3/4x, 1/4x 4 T
xs -200
X2 = 600
x5 = -200
Xe =X 1
v. b. I v. n b
Graficamente tendremos:
270
x.,
3X. +4X 2 =2400
271
Debido a que al introducir Ia nueva variable le hemos dado un
tratamiento de v.n.b. en Ia soluci6n optima originaL Ia obtenci6n de
sus coeticientes en el sistema final moditicado se hace mediante el
procedimiento indicado en los puntas 1 (cambios dec; cuando x·
es no basica) y 4 (cambios en a,: cuando x' 1 es no Mslca), en los
cuales los valores originates decAy a A (i = 1 ,2, .. ,m) son ceros.
1
L a~ +, ( aiA- acA.) = 2:
n a~ n +, a.,. (k= 1,2, ... ,m) [pues~ =OJ
;=1 i=1
Ejemplo Vll-12.-
272
ylasa·, .. son:
simplificando:
z + 1/2x, + 1/2x4 - 2Xs = 1200 ¢,
x, + x3 + x5 = 600 600
3/4x, + x2 ~ 1/4x4 + 2x5 = 600 300
x. = Xs
x, = x2
v. b. v. n.b
Z + 5/4x, + 'S + 3/4X 4 = 1800
z• = 1800 x. • = o
5/8X 1 - 1/2X 2 + X3 - 1/8X4 300
3/Sx, ~ 1/2'S + 118X4 + x 5 = 300
x.• = 300
273
CAPITULO VIII
TEORIA DE AEDES
275
VIII -1.- AEDES.-
Definici6n Vlll-1.-
Sea una red que eonsta de 4 nodes: (s,x,y,t) y 5 areas: (s,x), (s,y), (x,y),
(x,t), y (y,t). Oibujar Ia red.
Esta red se dice dirigida, debido a que cada area involucra una
orientaci6n especffica.
Vemos que en Ia anterior red dirigida, el arco (s, x) E A, mientras que
el (x, s) (£ A.
Es posible tambien Ia existencia de redes no dirigidas, en las euales el
conjunto A consiste de pares no ordenados de nodos. Asi mismo, se pueden
encontrar redes mixtas, en las cuales algunos areas estan dirigidos y otros
no.
En nuestro curos eliminaremos las siguientes posibilidades:
a) Arcos (x, x), que eonduzcan desde un nodo x hasta el mismo, raz6n
por Ia cual todos los areas que consideraremos se supondran de
Ia forma (x, y) con x ;t: y.
b) La existencia de areas multiples que unan x cony, por lo que todas
276
nuestras redes contend ran un area. como maximo. que vaya de un
nodo a otro.
En forma sirnplista podemos tener
Definici6n Vlll-2 .. -
Se considera como una red a una grafica lineal en Ia cual existe algun
tipo de tlujo en sus areas.
A continuaci6n se listan algunos ejemplos de sistemas que satisfacen
Ia definicion de red, mismos que podemos encontrar a nuestro alrededor·
Rutas de maneJO de
Centres de trabaJo Trabajos
matenales
Estaciones
Lineas Telef6nicas
telefonicas
Estaciones de
Tuneles Personas
Metro
Definci6n Vlll-4.-
En Ia figura del ejemplo Vlll-1, Ia cadena s. (s, x), x, (x, t), t, conduce
desde "s" hasta ''t''. Esta red no contiene ni un solo ciclo dirigido. Si Ia red
fuese no dirigida, entonces tendrfamos varios ciclos, entre elias:
277
-· =4EDt:S
Definicion Vlll-5.-
Sean x., X2 , .... x" una secuencia de nodos diferentes, ios cuales tienen
Ia propiedad de que ya sea (x, x .) 6 (x ... x) es un arco (i= 1, 2.... n-1).
Restringiendo para cada i solo una de las 2 posibilidades, llamamos a Ia
secuencia de nodos y arcos resultante, un 'camino" desde x 1 hasta xn.
Nuevamente si X. = xn tenemos un ciclo.
Vemos entonces que un camino difiere de una cadena, en que el
primero permite Ia posibilidad de recorrer un arco en un sentido opuesto a
su orientaci6n, al ir desde x1 hasta xn.
Definicion Vlll-6.-
s 0 0 0
X -1 0 0
y 0 -1 --1 0
0 0 0 -1 -1
Definicion Vlll-7 .-
La red del ejemplo Vlll-1, es una red conectada, pero dejarfa de serlo si
quitasemos los areas (s, x) y (s, y).
278
Definici6n VIII-B.-
Definici6n Vlll-9.-
Una red G = (N, A) con "n" nodos, es un arbol si tiene (n-1) areas y
ningun ciclo (es decir, es una red coneetada).
Efeetuando Ia demostraci6n por inducci6n, tenemos:
El teorema es evidentemente cierto para 2 nodos. Supongamoslo cierto
para 2, 3, ... , n-2, n-1, debemos demostrar que es cierto para "n'' nodos.
Lema Vlll-1.-
Baja Ia hip6tesis del teorema Vlll-1 , existe cuando menos un nod a que
es un final, es decir, un punta "p" con solo un area (p, q) que lo conecta al
resto de Ia red.
Demostracion del lema Vlll-1 :
Para eneontrar un nodo que sea un final, comience seleccionando
cualquier nodo, sea p,; este nodo esta unido cuando menos a otro nodo, sea
p 2 , mediante un solo area (de no ser este el caso, eliminando cualquier area
(p,, p) que una un par de los restantes n-1 nodos, obtendrfamos n-2 areas
sin dctos). Par nuestra suposici6n inductiva, esto formaria un arbol, par lo
que existirfa una cadena de areas uniendo P; con pi' Adjuntando el area (P;· pi)
a esta cadena, formarfamos un ciclo, lo cual serfa contrario a nuestra
suposicl6n. Debido a que existe un area entre p, y p 2 , muevase a p2 a Ia largo
del area (p,, p2 ). Pase a p~ a lo largo de otro area. (en caso de ser posible).
Dado que el numero de nodos es finito y no existen ciclos, procediendo de
esta manera llegaremos a un punta "p" que sea un final con solo un area (q, p)
uniendolo al resto de Ia red.
279
Ia suposicion inductiva. Si el punto p' eliminado y su arco (q. p) se insertan,
sera posible el conectar a "p' a cualquier otro nodo via el nodo q, con lo
cual probamos que una red con ' n·' nodes, (n-1) areas y sin ciclos esta
conectada y por lo tanto es un arbol.
De heche, cualquier pareja de las 3 condiciones:
a) G esta conectada.
b) G no tiene ciclos.
c) IA I = IN I~1
implica Ia tercera y caracteriza a G como un arbol.
Dada una red G = (N; A), sup6ngase que cada arco (x, y) E A tiene
asociado con 91 un numero real no negative c(x, y). Uamamos a c(x, y) Ia
"capacidad" del arco (x, y); intuitivamente, c(x, y) puede pensarse como Ia
representaci6n de Ia maxima cantidad de algun articulo que puede llegar a
"y", procedente de "x", por unidad de tiempo.
La capacidad de un arco no orientado es igual en ambas direcciones
[c(x, y)], mientras que en un arco orientado, es c(x, y) en Ia direcci6n de
orientaci6n y cero en Ia contraria.
lgnoraremos en nuestro curse lo posibilidad de nodes con capacidad
de flujo restring ida.
Definici6n Vlll-11.-
Si:
X EN
sea A(x), (antes x), el conjunto de todas las yEN tales que (y, x) E A:
A(x) = i y E N I (y, x) E A ~
similarmente, sea D(x), (despues de x), el conjunto de todas las yEN tales
que (x, y) E A:
D(x) = i yEN I (x, y) E A ~
280
Vlll-2.~ FLUJO MAXIMO.-
Considere una red que conecta 2 nodos (una fuente y un destine). par
media de varios nodos intermedios. Sup6ngase que el flujo que llega a un
nodo es igual al flujo que sale de el (conservaci6n del flujo). para todos
aquellos nodos que no sean Ia fuente y el destine. Supongamos que un flujo
f(x.y) par unidad de tiempo va del nodo 'x" al nodo y"; este flujo necesaria-
mente tendra que ser menar o igual a Ia capacidad c(x, y) del area. Es decir:
Sea Ia tuente (t) y el destine (d). Un tlujo estatico de valor F de "f' a "d",
satistace las desigualdades y ecuaciones lineales:
r F si x = fI
): f(x.y) - ) f(y,x;
y t::D(x) y EA(x)
=
l 0 si x #. f,dl ... (2)
-F si X = d
) f(x,y) - ) . f(y,x)
y Eo(x) y ~(x)
entonces, las ecuaciones (2) puede ser expresadas como: el flujo neto que
sale de Ia fuente es F y el flujo neto que ··sale" del destine es -F (o el flujo neto
que entra al destine es F), mientras que el flujo neto que sale de un nodo
intermedio es cera. Una ecuaci6n de este ultimo tipo, se llama ecuaci6n de
conservaci6n de flujo.
Dado un flujo F, llamamos f(x,y) al flujo en el arco (x,y). Cada flujo en el
area (x, y), ocurre precisamente en 2 ecuaciones de (2) y tiene un coeficiente
de 1 en Ia ecuaci6n correspondiente al nodo "x" y un coeficiente -1 en Ia
ecuaci6n correspondiente al nodo y. En otras palabras, Ia matriz de coefi-
cientes del sistema de ecuaciones (2), independientemente de Ia columna
de F, es Ia matriz de incidencia nodo-arco de Ia red.
Sabemos que:
281
Max F = \ .f(f,y)
, Eo(t)
s.a.
282
Definici6n Vlll-13.-
Sl (N;A) es una red y si F es un flujo entre (f) y (d) en (N; A), entonces
c(x, y) - f(x, y) es Ia capacidad residual" del arco (x, y) con respecto a F.
Teorema Vlll-2.-
283
Comience determinando todos aque!los nodos que pueden ser alcan-
zados desde Ia fuente, siguiendo un solo area con capacidad de flujo .> 0
que salga de (f) Para cada uno de estos nodos que fueron alcanzados,
determine todos los nuevas nodos (aquellos aun no alcanzados con
anterioridad), que puedan ser alcanzados desde ese nodo, a lo largo de un
area con capacidad de flujo > o Repita este procedimiento. sucesivamente.
con los nuevas nodos que vayan siendo alcanzados. El resultado sera Ia
formaci6n de un arbol de todos !os nodos que pueden ser alcanzados desde
Ia tuente, a lo largo de caminos con capacidad de tlujo > 0. Por lo tanto, este
procedimiento identificara un camino con capacidad de fluJO > 0 entre Ia
fuente y el destino, si es que existe cuando menos uno.
Ejemplo Vlll-5.-
Ejemplo Vlll-6.-
284
o bien:
o bien:
285
La busqued3 de nuevas caminos con capacidad de f!ujo ~' 0 que unan
(f) y (d). noes ya necesaria cuando iiegamos a un valor F que corresponda
a Ia minima capacidad de corte de todos aqueilos cortes que separan (f) de
(d), de acuerdo cone! teorema Vlll-2. ya que este nos garant1za haber liegado
a un flujo max1mo.
La 1a. parte del paso 3, tiene por objeto actualizar Ia capacidad residual
de los arcos considerados, para poder buscar nuevas caminos con
capacidad de flujo > 0, si aun existen. La 2a. parte, tiene par objeto poder
cancelar, total o parcialmente. algun flujo ya asignado a algun otro area, par
podersele asignar otra ruta que nos permita 1ncrementar mas el flujo total.
Es decir. puede suceder que un flujo asignado originalmente a traves de un
area, pudiera haber circulado a traves de otro area y que a! hacerlo par este
particular, impida que otro ftujo pueda hacerlo: sin embargo, de poderlo
cambiar al otro, abrimos nuevas caminos para incrementar nuestro flujo
total.
Ejemplo Vlll-7.-
286
de donde podemos observar que ya no quedan caminos con capacidad de
flujo > 0, par lo que nos conduciria a concluir que F"'"" = 2, lo cual es
totalmente falso, segun podemos observar a continuaci6n si efectuamos Ia
modificaci6n propuesta:
287
-
F=O
Paso 1.-
4 Fig. B
Paso 2.-
----
t, = 1
4 Fig. C
.... oo
Paso 3.-
Fig. D
es decir, c'(t, j) = c(i, j) ~f. y c'U, i) = cU. i) +- f, para todos los areas (i, j)
en ei camino escogido.
Pasos 1 y 2.- (C/ = 1)
0
'I ~1
~=1
-F=2
289
,;: - - -·~ -'-
0.,
290
Dado que no es posible llegar al destine. el proceso termina y hemos
encontrado un flujo maximo Frax = 3, que circula en Ia siguiente forma a
traves de Ia red:
F=3
-
F=3
coite 1
Ejemplo Vlll-9.-
292
• 1
0
VIII-3.-RUTA MAS CORTA.-
Metodo 1.-
Este metoda fue propuesto por G. J. Minty en 1957 yes aplicable para
el caso de redes no dirigidas. Simplemente construya un modele de hila de
Ia red, en el cuallas longitudes de los pedazos de hila sean proporcionales
a las longitudes de los areas. Tome Ia fuente con una mana yet destine con
Ia otra y jale, los hilos que queden tensos indicanin los arcos correspondien-
tes a Ia ruta mas corta
Metodo 2.·
?Q4
TEORIA DE qEDES
Ejemplo Vlll-10.-
Paso 1.-
Paao 2.-
295
factible ir de a a' g , a traves de 'd', a un menor, costa que directamente
Ademas, se puede ir a d directamente o a traves de c . En este caso de
empate, marque ambas rutas con lineas solidas.
G)
Paso 3,.-
Agregue todos los nodos que puedan ser alcanzados, partiendo de los
nodos considerados en el paso 2 y repita ese paso con respecto a elias.
lnserte las distancias:
296
~EOF<IA DE ?EDES
Paso 4.-
Metodo 3.-
2<37
Dada Ia informacion anterior. c:, de que manera identiticamos el nodo
que tiene Ia k-ava menor distancia al origen. a traves de su ruta mas corta ?.
es decir. 1.- cual de los nodos nuevas esta mas cercano a! origen ? . Para
calificar como candidate. un nodo nuevo debe estar conectado mediante
un solo area a alguno de los nodos en S (nodos originales). ya que de estar
conectado a traves de algun otro nodo nuevo, este ultimo estaria mas cerca
del origen que el primero. Ademas el nodo nuevo considerado. debera ser
el que se encuentra a menor distancia del nodo original de entre todos los
nodos (en caso de ser varios). que se encuentran unidos a el mediante un
solo area.
Debido a que el conjunto S contiene k nodos originales y solo con-
sideramos el nodo nuevo mas cercano a cada uno de elias, existen. cuando
mas, k candidatos para ser el nodo nuevo mas cercano al origen.
Para seleccionar el candidate ganador procedase asf:
1) Sea 'i un nodo enS
2) Sea o. su minima distancia al origen.
3) Sea j el nodo mas cercano a i y que no esta enS (si es que existe).
4) Sea d. su distancia desde i .
298
La anterior selecci6n. implica que el camino mas corto hasta is desde
el origen (con una distancia 6, , d), pasa par el nodo 's '.Para ver mas clara
esto, considerese cualquier otro camino entre is y el origen. Eventuaimente,
dicho camino debera llegar a algun nodo i en S. desde algun node i que
no esta enS (donde i puede ser is). Suponemos que las distancias a lo largo
del camino desde is hasta i son no negativas, de tal forma que Ia distancia
total hasta el origen, a lo largo del camino, no es menor que (6, +d) sin
embargo par (u):
Ejemplo Vlll-11.-
Destine
Elabore, para cada node, una lista de los areas que salen de ese node,
en arden ascendente de sus longitudes. Noes necesario ineluir los areas que
entran al origen o que salen del destine, si solo se desea calcular Ia ruta mas
299
corta entre estos dos, pero los segundos s1 deben incluirse. silo que se desea
es caicular Ia ruta mas corta entre ei origen y el resto de los nodos.
0 ;:. 8 c 0 ~ F ,--.
'..J
CG-3
k=O
El conjunto S consiste inicialmente solo del nodo 0 [S 0 = {0}].
k c~ 1
Debido al paso anterior, solo el nodo 0 es nodo original. por lo que
solamente Ia columna · 0" debe tomarse en cuenta (ya que indica las
distancias de todos los nodos que estan directamente conectados a el).
Rapidamente se descubre que el nodo A es e! mas cercano y par Ia tanto es
el candidate seleccionado para entrar a formar parte deS [S, = {O,A} ].
Para indicar que hemos encontrado Ia ruta mas corta al nodo A, anote
Ia distancia encontrada encima de Ia columna "A', encierre en un circulo Ia
entrada OA-1 de Ia columna ''0" y tache todos los areas que entran a A de
todas aquel!as columnas en las que aparezcan (BA-3, CA-3, DA-3).
1__1
0 A 8 c D E F G
OA--1~ :1
AB-3 BC-2 i CB-2 OC-2 i EF-1 GF-1
. 1
08-2 AC-3 '8~3- CD--2
I
! ~j-1! EC-3 GC-3
1! . i
AD--3 BG-4 CA-J-, DE-3
i ED--3 GB--4
I
I
CE-3
. CG-3
k=2
Los candidatos para el 2'2 nodo mas cercano ai origen, son aquellos
nuevas nodos que esten mas cercanos a los nodos originales {O,A}. E! nodo
mas cercano a 0 es 8 (parser el unico). mientras que el mas cercano a A es
tambien B (por ser el 1'2 no marcado de Ia columna 'A").
300
I
1 ! 2
2
0 A B c D E F G
X__?
I
OA-t_, -~3121
;!j DG-2
..
oB--2-. ~
. I
AG-3
--,! BG-2
~t1
GS':-2--j
CD--2
"•
OA:'-3~
I
:
EF-1
EG-3
GF-1
GC--3
I
! . 1
CA--3. 1i DE-3 I
I
' G8-'4 2
AD--3 BG-4 ED-3
I - I'
CE-3
!
CG-3
k=3
Buscamos ahara los nodos nuevas mas cercanos a los nodos
originales A y B (0 ya nose considera por tener una X su columna). Estes
los encontramos, mirando los arcos no marcados que se encuentran en
primer Iugar en las columnas correpondientes (es decir, en las columnas que
ya tienen distancias). Sus distancias son:
Nodo C via A - ( 1 + 3) = 4
Nodo 0 via A- ( 1 -r 3) = 4
Nodo C via B - (2 + 2) = 4
y debido al empate, selecionamos los nodos 0 y C (via A o B) como nuevos
miembros deS. Ahora S3 = {0, A, B, C, 0}.
Debido al empate, encerramos en un cfrculo las entradas AC--3 y A0-3
y Ia BC-2 de las columnas 'A' y "8" respectivamente, anotamos Ia distancia
de 4 encima de las columnas "C" y "'D", tachamos todas los demas arcos que
301
entren a C y a D r,CD- 2. OC -2. EC 3 ED-3. GC -31 'y ponemos una X debajo
de A para 'lid:car que ya no ia consideraremos mas
2 3 3
2 4 4
0 .A. 8 c D E "' G
2 3
X X
1 2, 3 ?
OA.-1 AB--3- BC--Z ca-:.2- :- DC -2 3:
.. EF-1 GF-1
2 3 3
OB-Z -. .<\C-~ 3
'· BA-3
1
CD--2
1 - Ec-..::3 3!
OA--3 --· GC--3
--· 31 1 3i 2
AD-3 ' BG-4 CA-3 _E0-:.3 -
DE-3 GB-4
CE-3
CG-3
k ~-4
lnvestigamos ahara cual es el nodo mas cercano a los nodos originales
en 8 3 (los nodos o· y A ya no se consideran par tener una X)
Tenemos:
Nodo E via C -- (4 + 3) = 7
Nodo E via 0- (4+3) =7
par !o que selecionamos al nodo G (via B) como nuevo nodo original. Ahara
8 4 = {0, A 8, C, 0, G}. Escribimos su distancia de 6 encima dela columna
de 'G' y encerremos en un circulo el arco BG-4, tachando todos los demas
areas no tachados que entran a G (CG-3); ponemos tambien una X en ia
columna 8":
1 2 3 3 4
1
·-- 2 4 4- 6---
0 A B c D E F G
X
2 X
3
--
X
4
302
k =5
Comparamos:
Nodo E via C - (4 ~ 3) -= 7
Nodo E vfa D- (4 ~ 3) = 7
Nodo F via G - (6 ~ 1) = 7
4 3 4 3 7 _5 4
6
0 A 8 C 0 E F G
X J X ._3 X 4 x 5 X _5. X 5 X
..A6-3~3i
4 GA:-3 ,I bE-3: 51 . ffi..3 ~ /Ga-4~-
. .
__ / 1 13G4.
-~ I -'L / l
! CE-3>51
{)0-341
•.
La ruta mas corta entre el origen '0' y ei destines 'f'·. puede ser
encontrada partiendo del nodo F hacia atras, siguiendo los areas encerados
en cfrculo. Asf encontramos:
F-G-B-0
OB ~ BG- GF
303
Destino
Ejemplo Vlll-12.-
Encontrar Ia ruta mas corta entre el nodo E y todos los demas en Ia red
del ejemplo Vll-11.
7 6 6 5 5 4 3 3
··- 3 3 o2 1 1 22
0 A B c D E F G
X
5 - X 6- X 6 X ·5- X 5 X 3 X -2 X -4
..·· 61
s:
OA,.:.1 --·r AO--t -J
jl
s0-2·.Pi
.. ·' I
.tB-s~ DC:.2c3i EF-r~!
1
FG-·.1 fj GF'-1· 1
~~
O&::z~J
~2~
.Aa::-3 4J SV2JI (DA~:1 1 EC~3~i
...
~3
3
• I '·._ - ~~ i
304
3
k=1
8, = [E,F]
k=2
82 = [E, F, G] node C via E = 0 + 3 = 3
node D vfa E = 0+3 =3
node G via F = 1 _._ 1 =2
k=3
83 = [E. F, G, C, D] node C vfa E = 0 +3 = 3
node D vfa E = 0 + 3 = 3
node C via G = 2 + 3 = 5
k=4
84 = (E, F, G. C, D, B] nod a B vfa C = 3 + 2 = 5
nodo A vfa D =3 +3 =6
node B vfa G = 2 + 4 = 6
k=5
85 = [E, F, G. C, D. B, A] node 0 vfa B = 5 ..,- 2 = 7
nodoAvfaC =3+3 =6
node A vfa D = 3 + 3 = 6
k=6
ss = [E, F, G, C, D, B, A, 0] nod a 0 vfa A = 6 + 1 =7
nod a 0 vfa B = 5 + 2 = 7
305
Vlll-4.- MINIMO ARBOL DE EXPANSION.-
Este tipo de problema es una var1aci6n del problema de Ia ruta mas
corta.
Dada una red conectada G. con n nodes. se pueden ir eliminando
arcos de G hasta que se obtiene un arbol: tal arbol se conoce como un arbol
de expansion de G. A nosotros nos interesa encontrar. de entre todos los
arboles de expansion aquel que tenga una longitud total minima.
Como antes, el conjunto de nodos de G y las distancias entre los pares
de nodos se conocen, pero ahara nose especifican los areas que concectan
dichos nodos. En este caso. en Iugar de encontrar Ia ruta mas corta a traves
de una red completamente definida. el problema consiste en encontrar
aquel!os areas para Ia red, que tengan Ia longitud total menor, al mismo
tiernpo que suministran una ruta entre eada pareja de nodes. Lo anterior se
logra eseogiendo los areas, de forma tal, que Ia red resultante forme un arbol
que conecte todos los nodes dados.
Si una red G y un arbol de extension T de G se especifiean, nos referimos
a los areas en T como areas del arbol' y a todos los demas como ·areas
fuera del arbol'.
Supongase que cada area (x.y) de una red conectada G. tiene
asociado con el un numero real a(x,y), que podernos pensar como Ia
'longitud' de dieho area . Supongase que T es un minima arbol de expansion
de G y sea: x,. x 2, ... , xk Ia cadena de areas del arbol T que une x. y x,. Aquf
(x,, x.J es un arco fuera del arbol T. Es clara entonces que:
ya que de ser menor que alguno de elias, el arbol serfa mas corte incluyendo
(x,, xk) en Iugar del otro.
Teorema Vlll-3.-
306
Vlll-4.1.- Metoda de Soluci6n para el Minima Arbol de
Expansion
1) Seleccione un nod a arbitrariamente y conectelo al nod a mas
cercano.
2) ldentifique aquel node no coneciado que se encuentre mas cercano
a un nodo conectado: conecte estes dos nodos. Repita este paso
hasta que todos los nodos hayan sido conectados.
Ejemplo Vlll-13.-
0 A B c D E F G
CG-3
k=1
El primer paso se inicia seleccionando arbitrariamente cualquier node.
Para facilitrar Ia comparaci6n con el problema de Ia Ruta mas Carta, es-
cogeremos el nodo "0".
El nodo "A" se identifica inmediatamente como el node mas cercano,
con solo mirar Ia parte superior de Ia columna "0", par Ia que el area OA se
agraga a Ia red. Para indicarlo, encierre las entradas OA-1 y A0-1 en un
circulo y ponga una X encima de los nodos "0" y "A", para indicar que son
nodes conectados.
307
X 1 X 1 X 2 x3 X 4 x7 X 6 x5
::; A 8 ~ D E F r;
OA-1 1 AO-t 1' B0-22 CB-2,3, DC-2 4: EF-1 7 ~""E-'
7 G=-1 6
OB-22 A£H.2 st~~ CD-24, D.A-3 4. E..C-3 7 t:G:_16 GC-3 5
sA~32' CA-3.3' QE-3-?' ED-3 7,
"
AC-33 I GB-4 5
A(}--3 4) f?G-45! CE--3 7,
CG-3.5!
k=2
Para identificar el nodo no conectado mas cercano a alguno de los ya
conectados (0 y A). solamente compare las entradas superiores no mar-
cadas para cada nodo conectado (08-2 y A8-3) Dado que 2 < 3, 08 es el
siguiente area que se agrega a Ia red. Para indicarlo, encierre 08-2 y 80-2
en un circulo y ponga una X encima del nodo B. Debido a que nunca
debemos agrager a Ia red un area que conecte 2 nodos ya conectados con
anterioridad, cualquier area potencial desde el nodo mas recientemente
conectado, que pudiera unirlo con los otros, debe ser eliminado de futuras
consideraciones. Par lo tanto tache AB-3 y BA-3.
k=3
El siguiente area se obtiene comparando AC-3 y BC--2, que nos indica
que 8C debe ser agragado a Ia red. Encierre 8C-2 Y C8-2 en un circulo y
ponga una X encima del nodo 'C". para indicar que ya esta conectado. Tache
AC-3 y CA-3.
k=4
Compare ahara AD-3. 8G-4 y CD-2. Se escoge CD como el nuevo
area de Ia red. Encerramos en un circulo CD-2 y DC-2. ponemos una X sabre
''D" y tachamos AD-3 y DA-3.
k=5
Compare BG-4, CG-3 y DE-3. Arbitrariamente, entre el empate, es-
cogemos CG como el nuevo area. Encierre CG-3 y GC-3 en un circulo y
ponga una cruz sabre G". Tache BG-4 y G8-4. N6tese que aqui solo se
escogi6 uno de los areas en empate. el otro no recibe ninguna marca.
308
k=6
Compare CE-3, DE-3 y GF-1. Seleccionamos GF como nuevo area
de Ia red. Encerramos GF-1 y FG-1 en un clrculo y ponemas una X sabre
F.
k=7
Compare CE-3, DE-3 y FE-1. Seleccione FE como nuevo area de Ia
red. Encerramos FE--1 y EF-1 en un circulo y ponemos una X encima de 'E"
Tachar CE-3, EC---3, DE-3 y ED--3.
Debido a que todos los nodos han sido ya conectados, hemos ter-
minado. La red resultante tiene una longitud total de sus areas igual a 12.
1 1 7
X 3- X- x- x- x-6
G
,GF:1' 6
te-3 '5
', _I
~.Ei
309
COn una /OflQitUd total c= 12.
::1 mismo resultada se hubiese abtenida. cualesquiera que hubiese sida
el noda ;nicial.
Ejemplo Vlll-14.-
7 5 6 4 3 2 9
X X X X X X X X X X
0 1 A 71 8 ,2 C s, D ~ ~
I, G .;2, H~~
'S F
L* = 2+1-2-r3-r-3-r1 +2-'-4-'-3 = 21
310
Vlll-5.- PERT
La tecnica de evaluaci6n y revision de programas, PERT (E_rogram
.§valuation and Review Iechnique), uno de los varios metodos de ruta critica
que existen en Ia actualldad, tuvo sus principios en Ia Grafica de Gant. PERT
se desarrollo pare el Proyecto Polaris en 1958, por Ia Oficina de Proyectos
Especiales de Ia Marina de EEUU, junto con Ia Lockheed Aircraft Corporation
yen colaboraci6n con Ia Booz, Alden and Hamilton (empresa consultora en
administraci6n) y a el se le atribuye el haber adelantado Ia terminaci6n del
proyecto en mas de 2 arias.
311
Paso 3.-
Paso 4.-
Paso 5.-
312
TEORIA DE REDES
;: ~,
requisites de
mane de obra
mano de obra
disponible
~
,---r-
:
.
I
313
Definicion Vlll-16.-
Una flecha que representa una mera dependencia de una actividad con
respecto a otra. se llama actividad artificial o actividad ficticia . General-
mente se representan mediante flechas punteadas y tienen asociadas a elias
una estimaci6n de tiempo de cera.
AI construir un diagrama de flechas, el planeador debe tener en cuenta
las actividades requeridas y sus respectivas reiaciones de tiempo, lo que
puede hacerse escribiendo una lista de las actividades del proyecto. En un
proyecto muy complicado, parece imposible anotar inicialmente todas sus
actividades: sin embargo, ias actividades adicionales aparecen a rnedida que
se desarrolla el diagrama de flechas. En seguida, ei planeador debe deter-
minar el arden 16gico de las actividades, osea Ia forma en que cada una de
elias se ajusta a las demas: '-hay alguna actividad que precede, o que sigue,
o que se desarrolle simultaneamente con otra actividad?. Final mente, es
necesario dibujar el diagrama de flechas para mostrar como se inter-
relacionan las actividades en el tiempo. El planeador debe vigilar las ac-
tividades que sean demasiado grandes o demasiado pequeflas. Es posible
que una actividad de gran tamafio pueda tratarse como mas de una, o que
muchas pequefias puedan combinarse en una sola.
Definicion VIII- 18.-
314
nodo. El nodo hacia el que todas las actividades se dirigen y ninguna sale.
es el evento correspondiente a Ia terminacion del proyecto.
Los eventos se numeran en serie. desde el principia hasta el fin de un
programa. La regia general para numerarlos. es que ningun evento puede
numerarse hasta que se hayan numerado todos los eventos precedentes, es
decir, no podemos numerar ningun evento hasta que hayamos numerado.
primeramente, Ia cola de cada flecha cuya punta senala el evento. El numero
de Ia cabeza de una flecha siempre es mayor que el de su cola.
Ejemplo Vlll-15.-
I
I
,::::: 6C'!'IVid8d
oF
f.ctiCia
I
I
I
F F=actlvldad
1
1 ct1c1a
1
I
~~
315
d) D depende de 8
E depende de 8
G depende de C.D y A
H depende de C, D y A
I depende de A
J depende de G y E
K depende de H. I y J
L depende de G y E
M depende de K
316
•
~B E
~ '\3;_}
-0
;o....-.;-C
317
S'l.t78
.o.--:.c-.-o
·~·
(__) B \_)
Vlll-5.3.- Estimaci6n de tiempos.-
La asignaci6n de tiempos de duraci6n a las actividades es indispen-
sable para completar 13. red PERT. .!.Debe hacerse esto basandose en ei costa
mas bajo posible, independientemente de Ia longitud de tiempo requerida, o
en ei tiempo mas corte posible. independientemente de los costas, o en algun
compromise entre los dos casas anteriores, o sabre cualquier otra base?
La version original de PERT, hacia Ia suposici6n realista de que el tiempo
requerido para ilevar a cabo cada actividad en el proyecto. es real mente una
variable aleatoria que cumple con alguna distribuci6n probabilfstica. Sin
embargo, debido a que esta suposici6n complica bastante el procedimiento,
una version simplificada. que considera estos tiempos como constantes
predecibles, se utiliza a menudo en Ia pnktica.
La estimaci6n dei tiempo requerido para una actividad. puede entonces
estar basada en un solo valor, que es el procedimiento empleado por el CPM,
o puede estar basado en un sistema de 3 estimaciones de tiempo, que es el
caso del PERT, tal y como veremos mas adelante.
La estadfstica nos muestra que Ia mayor parte de los grupos de datos
tienden a tamar una forma de campana cuando se trazan (Distribuci6n
Normal), mientras que otras !a hacen en forma asimetrica (Oistribuci6n Beta).
v.m.p. = valor mas probable
I~
j \~ i:l)
318
Con lo anterior en mente. es ahara evidente que, con objeto de obtener
los resultados deseados, cuando se usa PERT. es necesario estimar el valor
esperado (t.) y Ia variancia (d') de los tiempos requeridos para cada actividad.
Debido a que los conceptos de valor esperado y variancia pueden
resultar demasiado complejos para aquellos individuos calificados para
estimar los requerimientos de tiempo para las actividades, PERT usa un
procedimiento estimative simplificado, por media del cual cantidades in-
tuitivamente signiticativas son obtenidas, mismas que luego se convierten en
estimadores del valor esperado y de Ia variancia. Este procedimiento consiste
en Ia obtenci6n de 3 estimaciones del tiempo requerido para cada actividad:
a) El tiempo "mas probable", llamado "m" y que trata de representar
Ia estimaci6n mas realista del tiempo que Ia actividad puede con-
sumir.
b) El tiempo "optimista", llamado "a" y que representa el tiempo en el
cual Ia actividad puede ser completada si todo sale excepcional-
mente bien.
c) El tiempo "pesimista", llamado "b" y que representa el mayor tiempo
que pudiese necesitar Ia actividad, bajo las mas adversas cir-
cunstancias.
e m te
Como puede apreciarse al efectuar c6mputos en PERT, ia distribuci6n
del tiempo de duraci6n de una actividad es meramente hipotetica, debido a
que no es posible efectuar un muestreo estadlstico en ella. Una vez que Ia
actividad ha sido ejecutada, el tiempo real consumido por ella, seat*, puede
ser considerado como una muestra tomada de esta distribuci6n hipotetica.
Sin embargo, todos los c6mputos se efectuan con anterioridad a Ia ejecuci6n
de Ia actividad, por lo que como ya se dijo, Ia base de los c6mputos PERT
no envuelve un muestreo estadistico, sino que depende del juicio de Ia
persona encargada de Ia actividad de referencia, ya que se le pide que de
acuerdo a su experiencia, estime los 3 tiempos mostrados en Ia figura anterior
(a,m y b), los que a su vez seran usados para estimar el valor esperado y Ia
variancia de Ia distribuci6n hipotetica del tiempo de ejecuci6n de Ia actividad.
319
Con objeto de obtener valores confiables para los estimadores de Ia
duracion de ia actividad, conviene tener presente Ia siguiente:
a) Una de las suposiciones importantes en el Teorema del Umite
Central, es Ia independencia de las variables aleatorias con-
sideradas. Debido a que este Teorema es Ia base de los computes
de probabilidad en PEAT, los estimadores de a . m y b deben
obtenerse satisfaciendo esta condici6n de independencia. es decir.
deben hacerse independientemente de Ia que pueda ocurrir en
otras actividades del proyecto, que pueden a su vez afectar Ia
disponibilidad de mana de obra y equipo planeada para Ia activ1dad
de referenda.
b) Las estimaciones de "a",' m' y 'b' no deben ser influenciadas par el
tiempo disponible para completar el proyecto. es decir, es il6gico
el revisar todas nuestras estimaciones, reduciendolas, despues de
saber que Ia ruta crftica del proyecto es demasiado larga.
c) Con objeto de obtener una atmosfera que conduzca a Ia obtenci6n
de estimaciones no desviadas para "a', 'm' y' b , debe dejarse bien
clara que los mismos son estimaciones y no compromisos de
programa.
d) En general, las estimaciones de a , m y b" no deben incluir
toleracias para eventos que, par ocurrir tan poco frecuentemente,
nose puede pensar en elias como variables aleatorias (p.ej. aetas
de Ia naturaleza: fuegos, huracanes. inundaciones, etc.).
e) En general, las estimaciones de· a". ·m y b deben incluir toleran-
cias para eventos norma/mente considerados como variables
aleatorias (p.eJ. los efectos del clima).
320
Con objeto de obtener Ia estimacion del valor esperado, tiene que ser
evidente que e! tiempo mas probable (m) debe tener una mayor ponderacion
que el optimtsta (a) y el pesimista (b). Ciertamente, hay mas probabilidad de
que un programa se complete mas cerca del tiempo "m' que de los otros 2
tiempos extremes. La formula desarrollada para el tiempo esperado de una
actividad (t.) es :
a+ 4m +b
6
El tiempo mas tardio para un evento (TL), puede ser definido como Ia
fecha mas lejana en Ia cual el evento puede ocurrir, sin retrasar Ia terminaci6n
del proyecto mas alia de su tiempo mas proximo.
Ejemplo Vlll-16.-
1---3-·0
321
D son 6 y 9 respectivamente. Los tiempos mas tardios para los eventos A B.C
y D son 0.2 6 y 9 respectivamente.
Definicion VIII- 21.-
La holgura para los eventos del ejemplo Vlll-16 es de 0.2,0 y 0 para los
y
eventos A,B.C D respectivamente. La ruta crftica para esta red seria A -c
-+D.
AI efectuar los calculos de T~ y Tc. se siguen estas reglas:
d) El tiempo mas tardio para un evento del que salen varias acti-
vidades, es el menor de los tiempos mas tardios de iniciacion de
las actividades que salen del evento en cuesti6n.
322
a) a2 para el nodo inicial de Ia red se supone cera.
b) 2
a para el evento que sigue a Ia actividad en cuesti6n, se
obtiene sumando Ia variancia de Ia actividad a Ia variancia del
evento que le precede. Lo anterior para aquellos casas en que
solo una actividad conduce al evento.
c) Para eventos a los que !Iegan varias actividades, a 2 se calcula
a lo largo del mismo camino usado para calcular te, es decir, el
camino mas largo. En caso de empate, escoger el camino que
da Ia mayor varian cia.
Ejemplo Vlll-18.-
t,-.? =
1 +8 + 3
6
12
=-=2
6
~-.? = [(3~1)r 1
9
t2-:l =
1 +16 + 7 24
6
=-=4
6
~--3= r(7~1)r=1
t2_4 =
1 +8 + 3 12
=-=2 ~-4 = [(3~1)r - 9
t3~ =
6
1 +8 + 9
6
2 +16 + 6
6
18
=-=3
6
24
~~= [< 9 ~
1< 6
ur
-2)r
16
4
9
a;-e = u4~2)r
2+12+4 18 1
t5-e = =-=3
6 6 9
Actividad ..!L a2
1-2 2 1/9
2-3 4 1
2-4 2 1/9
3-5 3 16/9
4-5 4 4/9
5-6 3 1/9
323
yl!l-5.- PERT
t.=3
cr~1/s G)
-~·------~
TE TL
~·-·------
valor valor
EVENTO esperado d esperado d TL- TE
0 0 0 27/9 0
2 2 1/9 2 26/9 0
3 6 10/9 6 17/9 0
4 4 2/9 5 5/9
5 9 26/9 9 1/9 0
6 12 27/9 12 0 0
De los calculos anteriores, vemos que Ia ruta critica es: 1--. 2 .... 3 .-. 5
-+6 y que el camino 2 .... 4 -. 5 tiene una holgura de 1. Esta ruta crftica, se
obtuvo uniendo entre silos eventos con holgura cera.
324
-rt::ORtA DE AEDES
Ejemplo Vlll-19.-
Para Ia red del problema Vll-15 encontramos que Ia ruta crftica era
1-2-3-5-6. Sea r Ia suma de las variables aleatorias t. a lo largo de este
camino. Tenemos:
r = r,_ 2 + r 2_ 3 + r 3--5 -+- r5-€
Media de r = TE
TE = (t.k2 + (t.)2-3 + (t.)3--5 + (t.)5-{j
TE = 2 + 4 .,_ 3 .,_ 3 = 12
Variancia de T' = ~* = of* '-2
+ of*2-J + of*3--5 + of*
5-6
a= v'3 = 1.73
TE- 3a = 12-5.2 = 6.8
TE + 3a = 12 + 5.2 = 17.2
325
·Jdi-5. PERT
Ejemplo Vlll-20.-
t
\"""2
= 1 +8+3 =.!£.=2
6 6 a~~=[3~1r ~
t
l-:J
=2+12+4=_:1__!!=3
6 6 ~-:J = [4 -6 2] 2 91
t
~ - 4
= 2 + 16 + 6 = 24 = 4
6 6 a~-4 = [6 -6 2] 2 49
t
2
~
= 2 + 16 + 6
6
= 24 = 4
6 a~~= [6 ~ 2r :
t
3
~
= 1 + 16 + 7 = 24 = 4
6 6 ~~= [ 7 -6 1 J 2 99
t
3
= 1 + 12 + 5 = _:~__!! = 3 ~-e= [ 5-6 1 J 2 94
-e 6 6
2
t
4
=2 + 24 + 10 = 36 = 6
~~= [ 106- 2] 196
~ 6 6
t
4
-?
=3 + 16 + 11 = 30 = 5
6 6 a:~~ ["; 3]'- ~6
t = 3 + 16 + 5 = 24 = 4 = [ 5-6 3] = _91
r
5
-e 6 6 a;-e
t
e-e
=2+20+8=30=5
6 6 a~-e = [8 ~ 2 =~
t = 1 +8+3=_g=2 ~-e = [ 3 -6 1] 2 - 9
'7-e 6 6
t =2+12+4=_:1__!!=3 4 -6 2] 2 91
8 a~~ [
~ 6 6
326
TE TL
-------~-- ·-·- ----------
valor valor
EVEN TO esperado d- esperado d- TL- TE
0 0 0 22/9 0
2 2 1/9 6 6/9 4
3 3 1/9 6 14/9 3
4 4 4/9 4 18/9 0
5 10 20/9 10 2/9 0
6 6 5/9 9 10/9 3
7 9 20/9 12 2/9 3
8 14 21/9 14 1/9 0
9 17 22/9 17 0 0
La probabilidad de completar el proyecto en o antes de 19 semanas,
se obtiene entrando a las tablas con el valor de Z, cuyo valor es el siguiente:
19-17 2X3 6
z= ¥2279 = 722 = 4.69 = 1·28
327
Vlll-5. 7.- Consideraciones.-
Siempre existe una ruta crftica para el sistema y siempre puede ser
identificada trazandola a traves de los eventos con holgura cera. La ruta
crftica es obviamente significativa, debido a que si algun evento a lo largo de
este camino requiere de mas tiempo que su fecha esperada determinacion,
el evento terminal puede esperarse que se retrase en un lapso igual; par lo
tanto, a las actividades a lo largo de Ia mencionada ruta, debe asignarseles
maxima prioridad. Si el tiempo mas proximo esperado para el evento terminal
no es satisfactorio. entonces, mejoras deben preveerse para estas ac-
tividades. Por otro lado, los eventos que nose encuentran en Ia ruta crftica
tienen holguras positivas, par lo que pequefias desviaciones de su fecha
esperada de terminacion, probablemente no afectaran a Ia fecha de
terminacion del proyecto.
El planeador tiene a su disposicion varies metodos de ajuste, uno de
los cuales es el intercambio de hombres, maquinas y materiales (si es que
son comparables), de Ia ruta no crftica a Ia critica. Otro metoda de ajuste de
Ia red, consiste en reducir las especificaciones tecnicas del proyecto (p.ej.
reducir Ia cantidad de pruebas que requiere este). Silas actividades pueden
ordenarse de otro modo, puede ser posible acelerar Ia terminacion del
proyecto. El planeador puede aprovechar Ia sobreposicion de las ac-
tividades. Ademas, el empleo del tiempo extra, proporciona flexibilidad
adicional para el ajuste y Ia nueva planeacion de Ia red.
Es interesante hacer notar, que el analisis de Ia ruta crftica esta basado
solamente en valores esperados, par lo tanto, existe cuando menos una
pequeria probabilidad, de que otra ruta diferente de Ia considerada crftica,
pueda volverse Ia mas larga.
Programas muy complicados de computadora, han sido desarrollados
para complementar PERT; par media de ellos, uno puede determinar
rapidamente el efecto de un cambia propuesto, o el efecto de un retraso en
Ia programacion. Esta habilidad para medir y evaluar el estado actual del
proyecto, mejorando asr el control administrative, ha sido citada como uno
de los mejores valores de PERT
Numerosas tecnicas similares a PERT han sido propuestas en afios
recientes, bajo numerosos nombres diferentes. Algunos de estos entatizan
e1 costa y/o el nivel de desarrollo tecnico de las actividades, ademas de sus
tiempos transcurridos. El mas notable de estes, principalmente par que no
se deriv6 de PERT sino que se desarrollo en forma paralela con el, es CPM,
el cual se ocupa de Ia planeacion, programacion y control de costas de los
trabajos del proyecto.
328
Vlll-6.- SISTEMA DE REDES CON ACTIVIDADES
EN EL NODO.-
En este sistema, las actividades se representan graticamente por los
nodes, en Iugar de por las flechas. Aquf, las flee has se utili zan sola mente para
representar ia relacion de dependencia entre los nodos.
Ocasionalmente, en este tipo de red, se agrega un nodo fictfcio, que
tiene un tiempo de ejecucion cera, con objeto de concentrar varias ac-
tividades en alguna relacion logica (p. ej. lnicio del proyecto o terminaci6n
del proyecto).
Las principales ventajas de este sistema son:
a) Con el sistema de actividades en el nodo, las relaciones de
precedencia se muestran por media de lineas que conectan los
nodos de actividades, en Iugar de tener que traer ias cabezas de
las flechas de todas las actividades que anteceden, a un punta
comun (evento), asf como las ·colas de todas las actividades
sucesoras. Por esta raz6n, el sistema de actividades en ei nodo es
mas simple de construir y modificar.
b) Este tipo de red no reqwere ei uso de actividades fictfcias, lo cual
elimina problemas al trazar Ia red, especialmente para los prin-
cipiantes.
c) Permite el uso de redes de nodos preimpresas, a las cuales solo
se les agregan las flechas de precedencia.
d) La informaciOn es presentada en una forma mas compacta, hacienda
que las descripciones y tiempos computados sean mas faciles de
leer.
e) La numeracion de las actividades es mas simple, ya que solo
se usa un numero, con lo que Ia codificacion de las actividades,
por proyectos, habilidades de mana de obra, etc., se simplifica
grandemente. Cuando se usan las actividades en las flechas, una
actividad se identiflca por el par de numeros de los eventos en los
extremes de Ia flecha; esto crea un problema de codificacion en
los nodos de union cuando las actividades que alii se reunen
requieren diferente coditicacion. Esto solo puede resolverse
agregando actividades fictfcias, tal y como se puede ver en el
ejemplo Vlll-21 (b).
La principal desventaja de utilizar actividades en el nodo, como sistema
de construir una red, es principalmente que no es un metoda estandar. P ei.
existen muy pocos programas de computadora para este sistema, mientras
que existen muchlsimos para el sistema de actividades en las flechas.
329
Ejemplo Vlll-21.-
Actividad Depende de
8 A
c A
0 ByC
E c
G DyE
G)
~G..._..
330
ii) actividades en los nodos:
~~
~·····-p----~
331
AEDES
FlH!O ~AXJ_M(l HUTA MA§CORTA MINIMO J\R80_l_Di= EXPRI=SION iJf: R I
Se hene una red cuyos elementos son det Se ttene una red cuyos elementos son del Se tiene una red G cuyos elementos son Se t1ene IJna red cuyos elementos ~on del
sugui~nte tipo: siguiente tipo: del siguiente tipo: srgu1ente t1po:
Nodos A y 8 Nodos. 0 y A Nodos 0 y A Nouos· (evento~) 1.~ '/ 3
Arcus Onentados: (I· .8) o A
l t 2
2 1 {B,A) Arcos: con propiedad (2) Arcos: con prop1edad (21 Aicm.· (aciiVIdade:.) .A. yB
A. i ( o': Capacldad de los tliCOS I ' 2 (gralmente dtstancta) 2 (A , (gralmenle drstancia) I 8 18 es antftcial o f1ct t
I ,0 ,I A) '_0 l
\.- · onentados: C(a,B) = 2 ; c1a).
3
C(A,A) -1 Nota: Se va dlbujando de Nottt.. Se va d1bUJando de A Acl 1-2
NOTA: Se va dibUJ&ndo de 2a3 2a3 8 Act. £.1 3
1a 3
r,I,ETO[)()l()(31A I
METODOLOGIA METODOIOGIA MHOQlliOGI.4
0. Etectuar los cortes y determiner el que es - : Se consJdera que_:_ --.- -----=-- -~ -. d-,.-A~?1 ~ [;'1 1. Se e!abora l1sta de l~s ~~~o-S de cad_a nodo, b ttempo pe~1rn1;:;la
al Flujo MAX {T.VII-2) I Sk \ o) ~ /
en arden ascendenle de sus longitudes. j
. . -.. Tl · trernpo
Encontra1 un camtno
1.1 A partu de Ia Fuente se alcanzan los
I · - ----it(_~ 1 2. Se selecciona arbJtrariamente cuaJqu1eq
1. Se elabora hsiH de los areas de cada nodo nodo. j
m tiernpo
Cl_f.JfOXIHlddO
IJfOUable
lt tiempo rna& tardiO
nodos mtis cercanos a traves de un area en orden ascendente de sus longitudes 3. Se alcanza el nodo mits cercano a uno1 Se L..Onsidera· a t1empo optirnisld
con capactdad > 0 No tncluye arcos No 1nc/uye arcosj
de los nodos conectados fen caso de· Te tiempo esperado
1 2 A parttr de los nodes alcanzados. se, que entren al origen que entren a 0 sJ lol empate se el1ge uno de elias y lot. otros n nurnero de eventos
alcanzan los nodo& mas cercanos a: no recrben ninguna marca) y se tttchan: - valor esperedo
0 y que salgan del que se busca es lal
trav&s de un arco con capactdad > 0.
1.3 Se aplica 1,2 haste llegar al destine. I destine X si lo que
se busca es !a
R.M.C. entre 0 y c/ul
aqueHos arcos q:;e desde el nodo ma~~
recJertemente conectado pud1eran untrlo
de los restantes 1.-SeotJIIenenlos fe (a 1 4m b)ib ~1e
2. Astgn~r a/ cammo encontrado, un flujo --11 a otros que ya lo estuvreran. ESH:. PASO c/act
la.mf01macapac;Jdad de t~os sus &cos (C )
A.M. C. entre 0 y X nodos. i
SE REPITE HASTA AGO TAR LA TAll LA. 1
3.- Drsmtnuir en C Ia capactdad de c/arco del I2 Se alcanza el nodo mB.s cercano al origen y 4. Se ttene un &rbol T cuya longilud total es: 2. Se obt1enen las (:/t! !tl.J-1-i)/b]~J d t U:Hjd
camino se cancel an lodes los otros areas que entren Ia minima que puede tener cualquier acl
lncrementar en c•Ia capacidad de C/arco i al nodo alcanzado. Se anota su dlst. mas S.rbol que alcance a todos los nodi)S dei
del camino en senttdo contrario ·~ c.orta al origen Ia red G o sea: el MINIMO AHBOL DEl
REGRESAH AL PASO 1 HASTA QUE: ' 3 - Se alcanza el sigwente nodo mas cercano EJ(PANSION I
4. Nq. hay a cam! nos de Ia fuente al desttno las I al or~an {en caso de em pale se aJcanzan 1' rabla ACTIVIDAD, le, d'
(C O) se cterran todos los nodos que empatan) y se can- 3 · Se litH· Jan en Ia red ong1nal In~ '/dfore:::.
VERIFICAR: FACTIBIUOAD celan las •nlradas a e&e 6 esos nodo(sj de Te y cR obtentendose rle esl8, los
CONSERVACION DE FLUJO 1
1
alcanzado{s). ESl'E PASO SE REPifE datos para elnborar Ia 2~ tablti en el
ti.- Se venfica el PASO 0 · HASTA AGO TAR LA TABlA 4. Se obllenen de Ia red rnod1f1cad~n srguiente arden
6.- Se anoia el FLUJO MAXIMO 4 - Se leen en Ia tabla los areas que conecten el paso 3 a)La medra de T"' LT e 3' 1 Se l1stan lo& evenlos de 1 an
El ongen con el El origen con cada {RUTACRITiCA) , ,, TFid., 3.2Seanota>uv.o; acumulando
1L -- TL evento post - Te act. post destine nodo b) la vananc1a de f .£...02 (AU I A ·1 3.3 Se anotd sua- acumuJdndo
!1.- Se suman las distandas eleg1das y se aneta
~~~!~~~ ·~3 ~=(~:?,~eod:l?~~~~~:~s~· c~~~--~ vd~
b) vanos camtnos:
TE "" mayor TE de laa act , convergen al Ia R.MG (]:= d2 est J ¥:bla 4
evento Del origen al Del origen a c..ada 0.· Se dJbuja ladl::>trtbuci6n norma: n obten1do en el pdso 3.2 y re:.tan-
T,L menor T L de las act. salen del even to destine no do 6. Se supone (data) una fecha T 5 y se T do ha!:>lo. CE="RO para el e·1entu 1
C) TE:_ . TL para node terminal .. ' - - - - - - - - - - - - - - - - - ' calcula laprobablltdad de cumpJir en l Reg: (ver reg/as). .,
o antes . j J.!J Se procade de regre:.o ~:un I« (7
RtcGlAS~ , 7.-SecalculaZ (Ts lr*Ji(.J LJSde-. hauer1do su vafcH t::n n o ·v
a) cP. -· o pat a node 1niclal REG LAS TE y TL: del.,._ acurnulando hasta alcanzar er1 ei
b) Un csrntno C12- r.:R evto. ant. s2 act1v. ant ', evento J el vato1 obten1do ~n el ~c1::.o
c) varros camrnos:c1? del cam1no usado TE, al ~olo 1 camtno TA~A
~: g~l t:r~~r~~=~t~,f~r~~~~~~~~~~·l
1
s1 empate tamar ~ 1 grande TE - TE evento ant 1 react. a'lt. j b
Zen.