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)
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 />
108 CAPÍTULO 6. OUTROS MÉTODOS<br />
A nova solução é viável: além <strong>de</strong> ser positiva,<br />
E seu valor é<br />
c T x ′ = 0.975 + 0.0125 + 1 = 1.19875,<br />
maior que o que tínhamos antes (1.75, conforme a equação 6.1).<br />
Continuamos até que a norma <strong>de</strong> c p seja suficientemente pequena (por<br />
exemplo, menor que 0.001).<br />
◭<br />
Descrevemos a seguir o algoritmo em pseudocódigo.<br />
repita:<br />
D ← diag(x)<br />
P ← I − DA T [AD 2 A T ] −1 AD 2 c<br />
se ||c p || ∞ < ε<br />
PARE -- solu ção suficientemente boa<br />
d ← D.c p<br />
se d i < 0 para todo i<br />
PARE -- problema ilimitado<br />
{<br />
}<br />
δ ← ε δ min j − x j<br />
d j<br />
: d j < 0<br />
x ← x + δd<br />
6.2.2 Métodos <strong>de</strong> barreira<br />
Notas<br />
Descrições do algoritmo do elipsoi<strong>de</strong> po<strong>de</strong>m ser encontradas nos livros <strong>de</strong><br />
Griva, Nash e Sofer [GNS09] e <strong>de</strong> Papadimitriou e Steiglitz [PS98], que discutem<br />
inclusive a representação computacional (o algoritmo calcula raízes<br />
quadradas, sendo necessário representar irracionais).<br />
O algoritmo usando affine scaling apresentado neste texto é uma simplificação<br />
do algoritmo <strong>de</strong> Karmarkar, elaborado por Narendra Karmarkar<br />
[Kar84].<br />
Exercícios<br />
Versão Preliminar<br />
Ex. 44 — Mostre como obter uma solução viável inicial para o algoritmo<br />
affine scaling.<br />
Ex. 45 — O algoritmo affine scaling po<strong>de</strong>ria ser usado em um problema<br />
<strong>de</strong> otimização com função objetivo não-<strong>linear</strong>, já que só precisa do vetor