estudos, Goiânia, v. 35, n. 2, p. 225-246, mar./abr. 2008.algoritmo simplex tem complexidade exponencial (pior caso) parauma certa escolha na entrada da base, apesar de funcionar bem naprática. A partir daí, surgiram outras classes de algoritmos para PLcom complexidade inferior à do simplex. Um deles, comcomplexidade polinomial e que funciona bem na prática, é ométodo de pontos interiores, publicado por Karmarkar [12] em1984.O que faremos aqui será reescrever o que já existe naliteratura sobre PL com uma perspectiva para discutir acomplexidade, em número de iterações, em algoritmos de pontointerior-inviável.Dividiremos este trabalho da seguinte forma: naseção 2, apresentaremos os problemas de PL primal, dual e primaldual;na seção 3, apresentaremos os problemas teste queresolveremos com nossas implementações; na seção 4,apresentaremos idéias básicas dos algoritmos simplex primal, duale primal-dual, dos algoritmos de pontos interiores (primais-duais)de passos curtos e preditor-corretor, e dos algoritmos de pontointerior-inviável(primais-duais) preditor-corretor e homogêneo eauto-dual para PL, a página do LabPL referente ao algoritmosimplex tabular primal fases 1 e 2, e o enunciado do algoritmo deponto-interior-inviável homogêneo e auto-dual; na seção 5,apresentaremos algumas implementações e, finalmente, faremosalgumas análises comparativas em torno destas implementações.OS PROBLEMAS <strong>DE</strong> PROGRAMAÇÃO LINEARO PPL primal, no formato padrão, é o seguinte problemade Otimização:T(P) minimizar c xsujeito a :Ax=bx ≥ 0,m n mnonde são dados A∈ R× , b∈R e c ∈ R . Supomos, semperda de generalidade, que posto(A) = m e 0 < m < n.226
Associado a todo PPL existe um outro PPL chamado dual.E o dual de (P) é definido por:onde(D)maximizarsujeito a :ATbyTy+s ≥ 0,my ∈ R é o vetor de variáveis duais e incluímosnexplicitamente um vetor com componentes não negativas s ∈ Rdenominado folga dual.Relacionando o par de problemas primal (P) e dual (D),existe um outro PPL chamado primal-dual, que consiste emencontrar, se existir, uma solução para o seguinte sistema deequações e inequações:s=cestudos, Goiânia, v. 35, n. 2, p. 225-246, mar./abr. 2008.(PD)ATyAx = b+ s = cxs = 0x, s ≥ 0,onde xs = 0 significa x = 0,j = 1,…,n.j s jAs condições de otimalidade para o par de problemasprimal e dual coincidem com as condições de Karush-Kuhn-Tucker, a saber: x é uma solução ótima de (P) se, e somente se,existe um par (y,s) tal que o sistema (PD) é satisfeito.Consideremos os problemas (P), (D) e (PD). O conjuntoviável e o conjunto de pontos interiores viáveis do problemaprimal (P) são, respectivamente,Xn= {x ∈ R ; Ax =b, x≥0 }eX0= {x ∈ X; x > 0 }.O conjunto viável e o conjunto de pontos interioresviáveis do problema dual (D) são, respectivamente,227
- Page 1: IMPLEMENTAÇÃO DE ALGORITMOSSIMPLE
- Page 5: (KM)Pminimizarsujeito a :xxε xε x
- Page 8 and 9: estudos, Goiânia, v. 35, n. 2, p.
- Page 10 and 11: 0( x,s,µ ) ∈ F × R+֏ δ(x,s,µ
- Page 12 and 13: de pontos interiores, para se obter
- Page 14 and 15: ( yestudos, Goiânia, v. 35, n. 2,
- Page 16 and 17: gerados pelo algoritmo. A última s
- Page 18 and 19: interior-inviável para a resoluç
- Page 20 and 21: ReferênciasBEALE, E. M. L. Cycling
- Page 22: YE, Y.; TODD, M. J.; MIZUNO, S. An