Где я иду не так в следующем коде LP? - PullRequest
Новая
       14

Где я иду не так в следующем коде LP?

1 голос
/ 13 мая 2019

Я пытаюсь решить задачу LP с двумя переменными с двумя ограничениями, где одно - это неравенство, а другое - ограничение равенства в Scipy.Чтобы преобразовать неравенство в ограничении, я добавил в него еще одну переменную с именем A. 

Min(z) = 80x + 60y

Ограничения: 

0.2x + 0.32y <= 0.25
x + y = 1
x, y <= 0

Я изменил ограничения неравенства с помощью следующих уравнений:добавление дополнительной переменной A

0.2x + 0.32y + A = 0.25
Min(z) = 80x + 60y + 0A
X+ Y + 0A = 1

from scipy.optimize import linprog
import numpy as np

z = np.array([80, 60, 0])
C = np.array([
[0.2, 0.32, 1],
[1, 1, 0]
])
b = np.array([0.25, 1])
x1 = (0, None)
x2 = (0, None)
sol = linprog(-z, A_eq = C, b_eq = b, bounds = (x1, x2), method='simplex')

Однако я получаю сообщение об ошибке

Неверный ввод для linprog с method = 'simplex'.Длина границ несовместима с длиной c

Как это исправить?

1 Ответ

2 голосов
/ 13 мая 2019

Проблема в том, что вы не предоставляете границы для A.Например, если вы запустите 

linprog(-z, A_eq = C, b_eq = b, bounds = (x1, x2, (0, None)), method='simplex')

, вы получите: 

     con: array([0., 0.])
     fun: -80.0
 message: 'Optimization terminated successfully.'
     nit: 3
   slack: array([], dtype=float64)
  status: 0
 success: True
       x: array([1.  , 0.  , 0.05])

Как видите, ограничения выполнены: 

0.2 * 1 + 0.32 * 0.0 + 0.05 = 0.25  # (0.2x + 0.32y + A = 0.25)

, а также 

1 + 0 + 0 = 1  # (X + Y + 0A = 1)
...