08.02.2015 Views

Programação Linear (e rudimentos de otimização não-linear)

Programação Linear (e rudimentos de otimização não-linear)

Programação Linear (e rudimentos de otimização não-linear)

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.

notas <strong>de</strong> aula – versão 64 - Jerônimo C. Pellegrini<br />

94 CAPÍTULO 5. ANÁLISE DE SENSIBILIDADE<br />

Assim, a condição <strong>de</strong> otimalida<strong>de</strong> é<br />

c q − c T B Ba′ q ≤0<br />

c q − c T B B(a q + ∆e p ≤0<br />

( )<br />

c q − z q − ∆ c T B Be p ≤0<br />

( ) ∑<br />

c q − z q − ∆ c i B ip ≤0<br />

i<br />

( ) ∑<br />

−c q + z q + ∆ c i B ip ≥0<br />

i<br />

Isolamos ∆, levando em consi<strong>de</strong>ração dois casos:<br />

5.4 Nova variável<br />

∆ ≥ (z q − c q )<br />

∑<br />

i c ,<br />

iB ip<br />

∆ ≤ (z q − c q )<br />

∑<br />

i c ,<br />

iB ip<br />

se ∑ i<br />

se ∑ i<br />

c i B ip > 0<br />

c i B ip < 0.<br />

Se uma nova variável x n+1 é adicionada ao problema, sem mudanças nos<br />

coeficientes já existentes em A, b e c, teremos uma nova coluna a n+1 em<br />

A e um novo elemento c n+1 em c.<br />

Po<strong>de</strong>mos tomar o tableau que usamos para obter x ∗ e adicionar a nova<br />

coluna com a n+1 e c n+1 . Teremos também que calcular c n+1 − z n+1 . Isso<br />

já nos dará a informação que queremos: a solução x ∗ continuará sendo<br />

ótima somente se c n+1 −z n+1 ≤ 0. Caso não seja, po<strong>de</strong>mos imediatamente<br />

incluir a n+1 na base e usar o algoritmo Simplex para obter uma nova solução<br />

ótima.<br />

5.5 Nova restrição<br />

Suponha agora que uma nova restrição tenha sido adicionada ao problema.<br />

Se a solução ótima x ∗ satisfizer a restrição, ela claramente continuará sendo<br />

ótima. Trataremos então do que acontece quando nossa antiga solução<br />

não obe<strong>de</strong>ce a nova restrição (não é mais viável).<br />

Representamos a nova restrição por<br />

∑<br />

a m+1 x j = b m+1 ,<br />

Versão Preliminar<br />

j≤n

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

Saved successfully!

Ooh no, something went wrong!