VinBR

PO - 27/11/2012

May 28th, 2013
58
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.73 KB | None | 0 0
  1. # # # # # # # # # # # # #
  2. # Pesquisa Operacional #
  3. # Prof. Alexandre Cunha #
  4. # notasdogarcia.blog.com #
  5. # # # # # # # # # # # # #
  6.  
  7. Vamos aprender a trabalhar problemas de programação Linear Inteira
  8.  
  9. Existem vários outras modelagens que podemos fazer.
  10.  
  11. Z = max cTx
  12.  
  13. x € X contido em Z^k (Inteiros)
  14.  
  15. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  16.  
  17. Seja um Prob. De programação Inteira.
  18.  
  19. Com heurísticas podemos encontrar limites superiores.
  20. Ou limites inferiores.
  21.  
  22. Esses limites garantem que o valor ótimo é menor/maior ou igual à eles.
  23.  
  24. Dessa forma se encontrarmos um limite superior igual ao inferior, teremos certeza de que obtivemos a sol. ótima.
  25.  
  26. O limite inferior conseguimos com o primal, e o superior com o dual por exemplo.
  27.  
  28. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  29.  
  30. A dificuldade dos PPLs Inteiros é sua natureza NP-Completo.
  31.  
  32. Para contornar esse problema, realizamos a Relaxação do problema.
  33. de forma que conseguimos um novo PPL que não seja inteiro por exemplo.
  34.  
  35. Nossas garantias quanto a esse novo problema são:
  36.  
  37. 1) O dominio do novo prob. contem o dominio original. (senão podemos perder a sol. ótima)
  38. 2) Substituimos a func. objetivo max ou min por outra que assuma um valor não menor/maior
  39. que a original em qualquer ponto do problema original.
  40.  
  41. 3) Obtemos a sol. ótima desse novo problema.
  42.  
  43. Essa nova sol. ótima nos dá um limite superior/inferior para o nosso problema.
  44.  
  45. // Nota: O que significa seu novo prob. ser inviável? logicamente o prob. original também é.
  46.  
  47. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  48.  
  49. Seja o Prob. inteiro:
  50.  
  51. max 4x1 -x2
  52. 7x1 -2x2 <= 14
  53. x2 <= 3
  54. 2x1 - 2x2 <= 3
  55.  
  56. x € Z² positivo
  57.  
  58. O vetor (2,1) é uma sol. viável, logo z>= 7.
  59.  
  60. O valor ótimo da relaxação é 59/7 para x = (20/7,3)
  61.  
  62. Como queremos uma sol. inteira temos que:
  63.  
  64. z < 59/7 => z = 8
  65.  
  66. - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  67.  
  68. Formulação 1:
  69.  
  70. max cTx
  71.  
  72. x € P1 intercessão Z^n positivo
  73.  
  74. Formulação 2:
  75.  
  76. max cTx
  77.  
  78. x € P2 intercessão Z^n positivo
  79.  
  80. Então estamos tentando solucionar um mesmo problema de duas formas diferentes.
  81.  
  82. A relaxação do 1 é:
  83.  
  84. max cTx
  85. x € P1
  86.  
  87. Relaxação 2:
  88.  
  89. max cTx
  90. x € P2
  91.  
  92. Possibilidades:
  93.  
  94. 1) P1 contido em P2
  95. 2) P2 contido em P1
  96. 3) Existe um x em P1 tal que x não pertença à p2 e o contrário.
  97.  
  98. Em casos genéricos (considere diversas instancia diferentes do problema)
  99. iremos ter as seguinte relações de dominancia:
  100.  
  101. Em (1) o limite de P1 <= limite de P2
  102. Em (2) o limite de P1 >= limite de P2
  103. Em (3) os limites podem ser diferentes para cada caso. -> não há dominancia.
  104.  
  105. Agora vejam a sacada.
  106. Supondo o caso (3). P1 e P2 contém a sol. ótima do original em seu dominio.
  107. E ambos discartam partes do dominio do outro => essas partes não contém sol. do prob. original.
  108.  
  109. Sabendo isso podemos juntar os dois problemas, obtendo um dominio que é a intercessão dos dominios de P1 e P2.
  110.  
  111. E obteremos um novo dominio que é mais próximo das soluções inteiras do prob. original.
  112.  
  113. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  114.  
  115. Problema de programação combinatória:
  116.  
  117. max 7x1 + 4x2 + 5x3 + 2x4
  118. 3x1 + 3x2 + 4x3 + 2x4 <= 6
  119. x € B⁴ (binário em 4 dimensões)
  120.  
  121. Uma solução é: x = (1,1,0,0)
  122.  
  123. Note que isso é um problema da mochila.
  124.  
  125. Um jeito de de fazer a relaxação desse problema é assumir o seguinte:
  126. Primeiro ordenamos as variáveis em ordem de eficiencia.
  127.  
  128. Onde eficiencia é o tanto que é melhorado no objetivo (coeficiente) sobre tanto de espaço que ocupa coeficiente em A.
  129.  
  130. Vamos preenchendo a mochila com os itens de melhor eficiencia até que um deles não vai caber completamente na mochila.
  131.  
  132. Agora temos uma solução que não é a ótima do original mas é uma sol. da relaxação.
  133.  
  134. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  135.  
  136. // Fim da aula!
Advertisement
Add Comment
Please, Sign In to add comment