12.07.2015 Views

Programaç˜ao Linear - Notas de aula - CEUNES

Programaç˜ao Linear - Notas de aula - CEUNES

Programaç˜ao Linear - Notas de aula - CEUNES

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

CAPÍTULO 7. O MÉTODO SIMPLEX 587.4 Melhorando uma solução básica viávelComo dito anteriormente, o método Simplex busca por soluções básicas viáveis ótimas.Assim, nosso interesse é obter uma solução básica viável melhor que uma corrente. Como cadasolução básica viável correspon<strong>de</strong> a uma base B, <strong>de</strong>vemos obter uma outra base B ′ melhorque B. Ora, as matrizes básicas são formadas por colunas <strong>de</strong> A. Deste modo passaremos <strong>de</strong>uma base B para outra B ′ trocando uma <strong>de</strong> suas colunas por uma coluna não básica, <strong>de</strong> N.Devemos estabelecer então critérios para entrada e saída <strong>de</strong> variáveis da base.Suponhamos que temos uma solução básica [ B −1 b 0 ] t≥ 0 e consi<strong>de</strong>remos o mo<strong>de</strong>lo PLanterior, escrito nas VNB’s. Como sabemos, a FO só po<strong>de</strong> ser melhorada se existir j tal quez j − c j > 0. A fim <strong>de</strong> diminuir ao máximo a FO, escolhemos o índice k com maior z k − c k . Esseserá o critério para entrada <strong>de</strong> uma variável na base, ou seja, escolhemos a variável não básicax k tal quez k − c k = max {z j − c j },j∈Re a aumentamos. As VNB’s restantes continuam iguais zero, e então o novo valor da FO seráTambém, <strong>de</strong> x B = b − ∑ j∈R y jx j = b − y k x k segue que⎡⎢⎣x B1x B2.x Br.x Bmz = z 0 − (z k − c k )x k . (7.2)⎤⎡=⎥⎦ ⎢⎣⎤b 1b 2.b r⎥. ⎦b m⎡−⎢⎣⎤y 1ky 2k.y rk⎥. ⎦y mkx k , (7.3)on<strong>de</strong> a componente i <strong>de</strong> y k é <strong>de</strong>notada por y ik . Vamos analisar o quanto po<strong>de</strong>mos aumentar x kmantendo viabilida<strong>de</strong>, ou seja, x B ≥ 0. Se y ik ≤ 0 então x k ≥ 0 po<strong>de</strong> aumentar arbitrariamente.Agora, se y ik > 0, o valor máximo para x k <strong>de</strong> modo a manter x Bi ≥ 0 é b i /y ik (<strong>de</strong> fato,0 ≤ x Bi = b i − y ik x k ⇒ x k ≤ b i /y ik ). Assim, aumentamos x k a seu valor máximo permitido,x k = b { }rbi= min ; y ik > 0 . (7.4)y rk 1≤i≤m y ikObserve que x k ≥ 0, e na ausência <strong>de</strong> <strong>de</strong>generação x k > 0, pois b r > 0 (lembre-se <strong>de</strong> que naausência <strong>de</strong> <strong>de</strong>generação, b = B −1 b = x B > 0). Substituindo x k = b r /y rk na equação (7.3),obtemosx Bi = b i − y iky rkb r , ∀i = 1, . . . , m. (7.5)Assim x Br = 0 (verifique!). Escolhemos então a variável x Br para sair da base. Assim, escolhidox k para entrar, o critério para saída da base é a variável x Br tal que{ }b rbi= min ; y ik > 0y rk 1≤i≤m y ikse houver um y ik > 0 (o caso em que não houver será estudado adiante, e indicará que PL éilimitado).O processo acima passa da matriz básicaB = [ a B1 . . . a Br−1 a Br a Br+1 . . . a Bm]

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!