martes, 21 de septiembre de 2010

PROGRAMACIÓN LINEAL - Método Simplex



El tema fué extraido de la página web q nos proporcionó el docente encargado del curso. http://www.investigacion-operaciones.com


Método Simplex.

·El método del simplex fue creado en 1947 por el matemático George Dantzig .
El método del simplex se utiliza, sobre todo, para resolver problemas de programación lineal en los que intervienen tres o más variables.
El álgebra matricial y el proceso de eliminación de Gauss-Jordan para resolver un sistema de ecuaciones lineales constituyen la base del método simplex.

Es un procedimiento iterativo que permite ir mejorando la solución a cada paso. El proceso concluye cuando no es posible seguir mejorando más dicha solución.
Partiendo del valor de la función objetivo en un vértice cualquiera, el método consiste en buscar sucesivamente otro vértice que mejore al anterior. La búsqueda se hace siempre a través de los lados del polígono (o de las aristas del poliedro, si el número de variables es mayor). Cómo el número de vértices (y de aristas) es finito, siempre se podrá encontrar la solución.

El método del simplex se basa en la siguiente propiedad: si la función objetivo, f, no toma su valor máximo en el vértice A, entonces hay una arista que parte de A, a lo largo de la cual f aumenta.

Con miras a conocer la metodología que se aplica en el Método SIMPLEX, vamos a resolver el siguiente problema: 



MaximizarZ= f(x,y)= 3x + 2y
sujeto a:2x + y 18
2x + 3y  42
3x + y 24
x0 , y 0
Se consideran las siguientes fases:

1. Convertir las desigualdades en igualdades
Se introduce una variable de holgura por cada una de las restricciones, para convertirlas en igualdades, resultando el sistema de ecuaciones lineales: 



2x + y + h = 18
2x + 3y + s = 42
3x +y + d = 24

2. Igualar la función objetivo a cero
- 3x - 2y + Z = 0

3. Escribir la tabla inicial simplex
En las columnas aparecerán todas las variables del problema y, en las filas, los coeficientes de las igualdades obtenidas, una fila para cada restricción y la última fila con los coeficientes de la función objetivo: 



Tabla I . Iteración nº 1 
BaseVariable de decisiónVariable de holguraValores solución
xyhsd
h2110018
s2301042
d3100124
Z-3-20000

4. Encontrar la variable de decisión que entra en la base y la variable de holgura que sale de la base
  1. Para escoger la variable de decisión que entra en la base, nos fijamos en la última fila, la de los coeficientes de la función objetivo y escogemos la variable con el coeficiente negativo mayor (en valor absoluto).
    En nuestro caso, la variable x de coeficiente - 3.
    Si existiesen dos o más coeficientes iguales que cumplan la condición anterior, entonces se elige uno cualquiera de ellos.
    Si en la última fila no existiese ningún coeficiente negativo, significa que se ha alcanzado la solución óptima. Por tanto, lo que va a determinar el final del proceso de aplicación del método del simplex, es que en la última fila no haya elementos negativos.
    La columna de la variable que entra en la base se llama columna pivote (En color azulado).
     
  2. Para encontrar la variable de holgura que tiene que salir de la base, se divide cada término de la última columna (valores solución) por el término correspondiente de la columna pivote, siempre que estos últimos sean mayores que cero. En nuestro caso:
          18/2 [=9] , 42/2 [=21] y 24/3 [=8]
    Si hubiese algún elemento menor o igual que cero no se hace dicho cociente. En el caso de que todos los elementos fuesen menores o iguales a cero, entonces tendríamos una solución no acotada y no se puede seguir.
    El término de la columna pivote que en la división anterior dé lugar al menor cociente positivo, el 3, ya 8 es el menor, indica la fila de la variable de holgura que sale de la base, d. Esta fila se llama fila pivote (En color azulado).
    Si al calcular los cocientes, dos o más son iguales, indica que cualquiera de las variables correspondientes pueden salir de la base.  
     
  3. En la intersección de la fila pivote y columna pivote tenemos el elemento pivote operacional, 3.
5. Encontrar los coeficientes de la nueva tabla.
Los nuevos coeficientes de x se obtienen dividiendo todos los coeficientes de la fila d por el pivote operacional, 3, que es el que hay que convertir en 1.
A continuación mediante la reducción gaussiana hacemos ceros los restantes términos de su columna, con lo que obtenemos los nuevos coeficientes de las otras filas incluyendo los de la función objetivo Z
También se puede hacer utilizando el siguiente esquema:




Fila del pivote:
Nueva fila del pivote= (Vieja fila del pivote) : (Pivote)
Resto de las filas:
Nueva fila= (Vieja fila) - (Coeficiente de la vieja fila en la columna de la variable entrante) X (Nueva fila del pivote)
Veámoslo con un ejemplo una vez calculada la fila del pivote (fila de x en la Tabla II):  



Vieja fila de s2301042
------
Coeficiente222222
xxxxxx
Nueva fila pivote11/3001/38
======
Nueva fila de s07/301-2/326
 



Tabla II . Iteración nº 2
BaseVariable de decisiónVariable de holguraValores solución
xyhsd
h01/310-2/32
s07/301-2/326
x11/3001/38
Z0-100124

Como en los elementos de la última fila hay uno negativo, -1, significa que no hemos llegado todavía a la solución óptima. Hay que repetir el proceso:
  1. La variable que entra en la base es y, por ser la variable que corresponde al coeficiente -1
  2. Para calcular la variable que sale, dividimos los términos de la última columna entre los términos correspondientes de la nueva columna pivote:
    2:1/3 [=6] , 26:7/3 [=78/7] y 8:1/3 [=8]
    y como el menor cociente positivo es 6, tenemos que la variable de holgura que sale es h.
  3. El elemento pivote, que ahora hay que hacer 1, es 1/3.
Operando de forma análoga a la anterior obtenemos la tabla: 




Tabla III . Iteración nº 3
BaseVariable de decisiónVariable de holguraValores solución
xyhsd
y0130-26
s00-70412
x10-1016
Z0030-130

Como en los elementos de la última fila hay uno negativo, -1, significa que no hemos llegado todavía a la solución óptima. Hay que repetir el proceso:
  1. La variable que entra en la base es d, por ser la variable que corresponde al coeficiente -1
  2. Para calcular la variable que sale, dividimos los términos de la última columna entre los términos correspondientes de la nueva columna pivote:
    6/(-2) [=-3] , 12/4 [=3], y 6:1 [=6]
    y como el menor cociente positivo es 3, tenemos que la variable de holgura que sale es s.
  3. El elemento pivote, que ahora hay que hacer 1, es 4.
Obtenemos la tabla: 



Tabla IV . Final del proceso
BaseVariable de decisiónVariable de holguraValores solución
xyhsd
y01-1/20012
d00-7/4013
x10-3/4003
Z005/40033

Como todos los coeficientes de la fila de la función objetivo son positivos, hemos llegado a la solución óptima.
Los solución óptima viene dada por el valor de Z en la columna de los valores solución, en nuestro caso: 33. En la misma columna se puede observar el vértice donde se alcanza, observando las filas correspondientes a las variables de decisión que han entrado en la base: D(3,12) 


* Si en el problema de maximizar apareciesen como restricciones inecuaciones de la forma: ax + by  c; multiplicándolas por - 1 se transforman en inecuaciones de la forma - ax - by  - c y estamos en el caso anterior
* Si en lugar de maximizar se trata de un problema de minimizar se sigue el mismo proceso, pero cambiando el sentido del criterio, es decir, para entrar en la base se elige la variable cuyo valor, en la fila de la función objetivo, sea el mayor de los positivos y se finalizan las iteraciones cuando todos los coeficientes de la fila de la función objetivo son negativos 

No hay comentarios:

Publicar un comentario