INVESTIGAÇÃO OPERACIONAL
 8ª Aula

Método Simplex aplicado ao Problema de Transportes
Como os Problemas de Transportes são um dos tipos de problemas de programação linear, por isso, é
possível resolvê-los pelo método Simplex dado nas 2 aulas anteriores. No entanto, dada a especificidade destes
problemas, o método simplex pode ser simplificado – Método Simplex dos Transportes.

Preparação do Método
Após construir o quadro dos coeficientes das restrições para o método simplex, converter a função objectivo para
a forma de maximização e introduzir as variáveis artificiais z1, z2, ..., zm+n, obter-se-á o seguinte quadro simplex:
Variáveis
Equação
Básicas
Z
(0)
Coeficientes
Z
...
xij
...
zi
-1
cij
M
0
1
1
0
1
...
zm+j
M
...
Lado
Direito
0
(1)
...
zi
(i)
si
...
zm+j
(m+j)
1
dj
...
(m+n)
Cecília Rocha # 1
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Preparação do Método
Falta agora realizar algumas operações algébricas antes da 1ª iteração para eliminar os coeficientes das
variáveis (artificiais) básicas iniciais da linha (0) que sejam diferentes de zero.
Após essas operações, a nova linha (0) terá a seguinte forma:
Variáveis
Básicas
Equação
Z
(0)
Coeficientes
Z
-1
...
xij
Cij – ui - vj
...
zi
M - ui
...
zm+j
M - vj
...
Lado Direito
m
n
i 1
j 1
 si ui   dj v j
Onde:
 ui – múltiplo da linha original (i) que tem de ser subtraído (directa ou indirectamente) à linha original (0)
no método simplex, durante todas as operações que levam ao quadro actual
 vj – múltiplo da linha original (m+j) que tem de ser subtraído (directa ou indirectamente) à linha original
(0) no método simplex, durante todas as operações que levam ao quadro actual
 Se xij é uma variável não básica, então cij – ui – vj é interpretada como a taxa a que Z se irá alterar à
medida que xij aumenta
Cecília Rocha # 2
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Preparação do Método


Em primeiro lugar, não são necessárias variáveis artificiais porque se pode obter uma solução básica
inicial com métodos auxiliares simples
A linha (0) pode ser obtida sem utilizar qualquer outra linha, calculando os valores actuais de ui e vj
directamente. Dado que cada variável básica tem de ter coeficiente zero na linha (0), os valores de ui e
vj podem ser obtidos pela resolução de um conjunto de equações:
cij – ui – vj = 0



Cecília Rocha # 3
para cada i e j em que xij é variável básica
A variável básica de saída pode ser identificada facilmente sem utilizar os coeficientes das variáveis
básicas de entrada, assim como, a nova SBA pode ser detectada imediatamente sem se realizarem
nenhumas operações algébricas.
Deste modo podemos prescindir de quase todo o quadro do método simplex.
Além dos dados de base (parâmetros cij, oferta si e procura dj), o método simplex dos transportes só
precisa da SBA inicial, dos valores actuais de ui e vj e dos valores resultantes da operação cij – ui – vj
para as variáveis não básicas xij. Estes dados podem ser organizados num quadro denominado –
Quadro Simplex dos Transportes.
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Preparação do Método
Formato do quadro simplex dos transportes
Iteração
?
1
2
1
2
Destino
...
...
n
Origem
c11
c12
c1n
c21
c22
c2n
cm1
cm2
cmn
Oferta
s1
s2
...
m
Procura
d1
d2
sm
dn
Vj = cij - ui
cij
xij
Cecília Rocha # 4
ui = cij - vj
Variável Básica
cij
cij – ui - vj
Z=
Variável Não Básica
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Inicialização do Método dos Transportes e de Distribuição
O objectivo da inicialização é obter uma Solução Básica Admissível (SBA) inicial.
Antes de começar este processo de inicialização, temos de garantir que:
Número de Variáveis Básicas = m + n - 1
A razão para esta situação é que as restrições têm a forma de igualdades e o conjunto das m+n restrições
tem uma que é dispensável (por exemplo, uma das restrições de procura é igual à soma das restrições de
oferta menos as outras restrições de procura). Assim, qualquer solução básica aparece no Quadro dos
Transportes com m+n-1 variáveis rodeadas com um círculo, em que a soma em linha e coluna, corresponde
à oferta e procura, respectivamente.

Procedimento para obter uma Solução Básica Inicial



Das linhas e colunas em consideração seleccionar a próxima VB, de acordo com algum critério
Atribuir à VB um valor que corresponda à oferta ou procura remanescente, a que for menor
Eliminar essa linha ou coluna dos cálculos seguintes (se a linha e colunas se anularem em simultâneo, escolha
a linha para ser eliminada, posteriormente, a coluna servirá para atribuir o valor zero a uma variável degenerada)

Cecília Rocha # 5
Se só resta uma linha ou coluna, então o processo termina com a consideração de todas as variáveis
ainda sem valor atribuído.
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Exercício Exemplo (recordar)
Suponha que Inglaterra, França e Espanha produzem todo o trigo, cevada e aveia disponível no mundo. A procura mundial
de trigo corresponde à produção de 125 milhões de acres de solo. Com o mesmo objectivo são necessários 60 milhões de acres para
cevada e 75 milhões de acres para aveia. O total de solo agrícola disponível para este propósito, em Inglaterra, França e Espanha é de,
respectivamente, 70 milhões, 110 milhões e 80 milhões de acres. O número de horas de trabalho necessárias para produzir 1 acre de
trigo é de 18h em Inglaterra, 13 em França e 16 em Espanha. No caso do cevada são necessárias 15h em Inglaterra e 12h em França
e em Espanha. Para o aveia são precisas 12h em Inglaterra, 10 em França e 16 em Espanha. O custo da hora de trabalho para
produção de trigo é de 3 u.m., 2.4 u.m. e 3.3 u.m., respectivamente em Inglaterra, França e Espanha. Para a produção de cevada o
custo da hora de trabalho será de 2.7 u.m., 3.0 u.m. e 2.8 u.m. em Inglaterra, França e Espanha. No caso da aveia haverá um custo da
hora de trabalho de 2.3 u.m. em Inglaterra, 2.5 u.m. em França e 2.1 u.m. em Espanha. O problema é definir a melhor distribuição da
produção em cada país, de forma a satisfazer as necessidades mundiais de trigo, cevada e aveia mas minimizando o custo de
produção total.
a) Formular este problema como um Problema de Transportes, construindo o quadro de custos e requisitos;
b) Utilize uma rotina automática do SOLVER para encontrar um solução óptima para o problema;
c) Utilize o método de Vogel e o método do Custo Mínimo para determinar uma solução básica admissível inicial;
d) Resolva pelo Método dos Transportes.
Cecília Rocha # 6
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Alternativas para escolher uma Solução Básica Inicial

Método do Canto Noroeste
Destinos
2
1
1
Origens
2
3
Procura
vj





Cecília Rocha # 7
54
Oferta
3
40.5
27.6
36
25
70
70
31.2
55
52.8
110
55
33.6
33.6
5
125
60
ui
75
75
80
Z = 10 164
Começar pela variável x11
Se houver ainda oferta disponível, passar para a variável xi+1, j
Se só houver procura disponível, passar para a variável xi, j+1
Prosseguir até obter todas as variáveis básicas (as que têm um círculo) e todas as outras variáveis (não básicas) serão
zero.
O valor da Função Objectivo Z = 54*70 + 31.2*55 + 36*55 + 33.6*5 + 33.6*75 = 10 164 u.m.
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Alternativas para escolher uma Solução Básica Inicial

Método de Vogel





Origens
Calcular as diferenças, em linha e coluna, entre os 2 valores menores
Seleccionar a maior diferença global (considerando as linhas e colunas em conjunto)
Na coluna da maior diferença, escolher o menor custo e atribuir à variável correspondente o menor valor entre a
oferta e a procura
Eliminar a linha ou coluna respectiva
Repetir o procedimento
1
2
3
Procura
Diferença entre colunas
1
Destinos
2
3
54
40.5
27.6
31.2
36
25
52.8
33.6
33.6
125
60
75
Seleccionar x21 = 110
52.8 – 31.2 = 21.6
36 – 33.6 = 2.4
27.6 – 25 = 2.6
Eliminar a linha (2)
1
Origens
1
2
3
Procura
Diferença entre colunas
Cecília Rocha # 8
Destinos
2
3
Oferta
70
110
80
Diferença entre
linhas
40.5 – 27.6 = 12.9
25 – 31.2 = 6.2
33.6 – 33.6 = 0
Oferta
Diferença entre
linhas
54
40.5
27.6
70
40.5 – 27.6 = 12.9
52.8
33.6
33.6
80
33.6 – 33.6 = 0
15
60
75
Seleccionar x13 = 70
54 – 52.8 = 1.2
40.5 – 33.6 = 6.9
33.6 – 27.6 = 6.0
Eliminar a linha (1)
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Alternativas para escolher uma Solução Básica Inicial

Método de Vogel
(repetido)
Origens
1
1
2
3
Procura
Diferença entre colunas
Origens
1
2
3
Procura
Destinos
2
3
Oferta
Diferença entre
linhas
54
40.5
27.6
70
40.5 – 27.6 = 12.9
52.8
33.6
33.6
80
33.6 – 33.6 = 0
15
60
75
Seleccionar x13 = 70
54 – 52.8 = 1.2
40.5 – 33.6 = 6.9
33.6 – 27.6 = 6.0
Eliminar a linha (1)
1
Destinos
2
3
52.8
33.6
33.6
15
60
5
Oferta
Diferença entre
linhas
80
33.6 – 33.6 = 0
Seleccionar x31 = 15
Seleccionar x32 = 60
Diferença entre colunas
Seleccionar x33 = 5
Z = 8 340
Cecília Rocha # 9
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Alternativas para escolher uma Solução Básica Inicial

Método do Custo Mínimo


Seleccionar o menor custo de todo o Quadro de Custos (atribuir à variável x23 = 75)
Procurar sequencialmente os menores valores possíveis

(x13 não pode ser, dado que já se satisfez toda a procura)

(x21 = 35, valor excedente da oferta)

(x32 = 60, satisfaz toda a procura)

( x12 e x22, não podem ser, dado que já se satisfez toda a procura)

(x31 = 20, valor excedente da oferta)

(x11 = 70, valor da oferta)
1
1
2
3
Origens
Procura
Cecília Rocha # 10
125
Destinos
2
3
54
70
40.5
27.6
31.2
35
36
25
52.8
20
33.6
90
70
60
60
33.6
75
Oferta
75
70
110 35
80 20
Z = 9 819
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Exercício Exemplo (recordar)
Suponha que Inglaterra, França e Espanha produzem todo o trigo, cevada e aveia disponível no mundo. A procura mundial
de trigo corresponde à produção de 125 milhões de acres de solo. Com o mesmo objectivo são necessários 60 milhões de acres para
cevada e 75 milhões de acres para aveia. O total de solo agrícola disponível para este propósito, em Inglaterra, França e Espanha é de,
respectivamente, 70 milhões, 110 milhões e 80 milhões de acres. O número de horas de trabalho necessárias para produzir 1 acre de
trigo é de 18h em Inglaterra, 13 em França e 16 em Espanha. No caso do cevada são necessárias 15h em Inglaterra e 12h em França
e em Espanha. Para o aveia são precisas 12h em Inglaterra, 10 em França e 16 em Espanha. O custo da hora de trabalho para
produção de trigo é de 3 u.m., 2.4 u.m. e 3.3 u.m., respectivamente em Inglaterra, França e Espanha. Para a produção de cevada o
custo da hora de trabalho será de 2.7 u.m., 3.0 u.m. e 2.8 u.m. em Inglaterra, França e Espanha. No caso da aveia haverá um custo da
hora de trabalho de 2.3 u.m. em Inglaterra, 2.5 u.m. em França e 2.1 u.m. em Espanha. O problema é definir a melhor distribuição da
produção em cada país, de forma a satisfazer as necessidades mundiais de trigo, cevada e aveia mas minimizando o custo de
produção total.
a) Formular este problema como um Problema de Transportes, construindo o quadro de custos e requisitos;
b) Utilize uma rotina automática do SOLVER para encontrar um solução óptima para o problema;
c) Utilize o método de Vogel e o método do Custo Mínimo para determinar uma solução básica admissível inicial;
d) Resolva pelo Método dos Transportes.
Cecília Rocha # 11
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Teste de Optimização
Uma Solução Básica Admissível é óptima se e só se
cij – ui – vj  0,
para todo (i, j) em que xij é variável não básica.


Assim, só é necessário calcular todos os valores de cij – ui – vj
Deve-se iniciar o cálculo pela linha ou coluna que tenha mais variáveis básicas, como forma de facilitar os
cálculos
Destinos
2
1
1
Origens
2
3
Procura
vj
Cecília Rocha # 12
54
40.5
27.6
7.2
12.9
0
36
25
24
13
31.2
110
0
52.8
15
0
33.6
60
0
Oferta
3
70
70
110
33.6
5
0
125
60
75
V1 = 52.8
V2 = 33.6
V3 = 33.6
80
ui
U1 = C13-V3
= -6
U2 = C21-V1
= -21.6
U3 = 0

Esta já é a Solução
Óptima, no caso de
utilizarmos como Solução
Básica Inicial a obtida pelo
Método de Vogel
Z = 8 340
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Processo Iterativo

Teste de Optimização
Destinos
2
1
1
Origens
2
3
Procura
vj


70
0
31.2
35
0
52.8
20
0
125
V1 = 0
40.5
27.6
5.7
36
-7.2
25
24
0
33.6
60
0
60
75
33.6
0
75
V2 = 33.6 – 52.8 = -19.2 V3 = 33.6 – 52.8 = -19.2
Oferta
ui
70
U1 = 54
110
U2 = 31.2
80
U3 = 52.8
Z = 9 819
Ainda não é a solução óptima !
Escolher a VB Entrada


Cecília Rocha # 13
54
3
Como cij – ui – vj representa a taxa a que a função objectivo irá evoluir à medida que a variável não básica xij
aumenta, a VBE deverá ter um coeficiente cij – ui – vj negativo para diminuir o custo total.
A VBE será a que tem coeficiente cij – ui – vj mais negativo, ou seja, x13
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Processo Iterativo

Escolher a VB Saída





Aumentar a VBE irá provocar uma reacção em cadeia, na cadeia por nós definida.
Assim, passaremos a ter células receptoras e células fornecedoras, representadas no Quadro dos Transportes
pelos sinais + e –
Neste caso iremos utilizar a cadeia marcada a vermelho
O valor a transferir das células fornecedoras para as receptoras é dado pelo mínimo (x11 = 70, x23 = 75), ou seja, 70
VBE = x11
Destinos
2
1
1
Origens
2
3
Procura
vj
Cecília Rocha # 14
54
70
0
31.2
35
0
52.8
20
0
+
3
40.5
27.6
5.7
36
-7.2
25
24
0
33.6
60
0
+
75
33.6
0
125
60
75
V1 = 0
V2 = 33.6 – 52.8 = -19.2
V3 = 33.6 – 52.8 = -19.2
-
Oferta
ui
70
U1 = 54
110
U2 = 31.2
80
U3 = 52.8
Z = 9 819
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)

Processo Iterativo

Identificar a nova SB Admissível
Destinos
2
1
1
Origens
2
3
54
40.5
27.6
20.2
25.9
0
36
25
24
0
31.2
0
52.8
Origens
2
3
Procura
vj
Cecília Rocha # 15
20
0
Procura
vj
1
105
+
-
33.6
60
0
Oferta
3
5
33.6
125
60
75
V1 = 0
V2 = 33.6 – 52.8 = -19.2
V3 = 25 – 31.2 = -6.2
1
Destinos
2
3
40.5
27.6
7.2
12.9
0
36
25
24
13
31.2
110
0
52.8
15
0
33.6
60
0
+
-13
54
70
70
5
-13
U1 = C13-V3
= 33.8
110
U2 = 31.2
80
U3 = 52.8
Z = 8 405
Oferta
70
33.6
ui
125
60
75
V1 = 0
V2 = 33.6 – 52.8 = -19.2
V3 = 33.6 – 52.8 = -19.2
70
ui
U1 = C13-V3
= 46.8
110
U2 = 31.2
80
U3 = 52.8
Z = 8 340
2001/2002
INVESTIGAÇÃO OPERACIONAL

8ª Aula (cont.)
Cecília Rocha # 16
2001/2002
Download

8ª aula da Engª Cecília Rocha