Сумма переменных не совпадает с проблемой смешивания - PullRequest
Новая
       46

Сумма переменных не совпадает с проблемой смешивания

0 голосов
/ 31 января 2019

Я пытаюсь воспроизвести это смешение пример , но используя меньшее количество переменных, эта часть работает просто отлично: 

import pulp
from pulp import *

# Creates a list of the Ingredients
Ingredients = ['CHICKEN', 'BEEF', 'MUTTON', 'RICE']

# A dictionary of the costs of each of the Ingredients is created
costs = {'CHICKEN': 15, 
         'BEEF': 12, 
         'MUTTON': 17, 
         'RICE': 12
        }

# A dictionary of the protein percent in each of the Ingredients is created
proteinPercent = {'CHICKEN': 17, 
                  'BEEF': 2, 
                  'MUTTON': 16, 
                  'RICE': 8
                 }

# A dictionary of the fat percent in each of the Ingredients is created
fatPercent = {'CHICKEN': 10, 
              'BEEF': 14, 
              'MUTTON': 13, 
              'RICE': 16, 
              }

# Create the 'prob' variable to contain the problem data
prob = LpProblem("The Whiskas Problem", LpMinimize)

# A dictionary called 'ingredient_vars' is created to contain the referenced Variables
ingredient_vars = LpVariable.dicts("Ingr",Ingredients,0)

# The objective function is added to 'prob' first
prob += lpSum([costs[i]*ingredient_vars[i] for i in Ingredients]), "Total Cost of Ingredients per can"

# The  constraints are added to 'prob'
prob += lpSum([proteinPercent[i] * ingredient_vars[i] for i in Ingredients]) >= 15.5, "ProteinRequirement"
prob += lpSum([fatPercent[i] * ingredient_vars[i] for i in Ingredients]) >= 12.3, "FatRequirement"


prob.writeLP("WhiskasModel.lp")
prob.solve()
# The status of the solution is printed to the screen
print ("Status:", LpStatus[prob.status])

# Each of the variables is printed with it's resolved optimum value
for v in prob.variables():
    print (v.name, "=", v.varValue)

# The optimised objective function value is printed to the screen
print ("Total Cost of Ingredients per can = ", value(prob.objective))

Она рассчитывает оптимальное количество, требуемое для каждого ингредиента: 

Status: Optimal
Ingr_BEEF = 0.0
Ingr_CHICKEN = 0.77916667
Ingr_MUTTON = 0.0
Ingr_RICE = 0.28177083
Total Cost of Ingredients per can =  15.068750009999999

Но это не добавляет до 100%, когда я добавляю для этого ограничение к коду: 

prob += lpSum([ingredient_vars[i] for i in Ingredients]) == 100, "PercentagesSum"

Я получаю такой результат: 

Status: Optimal
Ingr_BEEF = 100.0
Ingr_CHICKEN = 0.0
Ingr_MUTTON = 0.0
Ingr_RICE = 0.0
Total Cost of Ingredients per can =  1200.0

Что неверно, поскольку оно не удовлетворяет другим противоречиям.

РЕДАКТИРОВАТЬ

Кажется, я неправильно истолковал это, я думал: если яхочу произвести 3 единицы, сумма входных данных должна быть 3.

Я думаю, что это так: 

# The constraints are added to 'prob'
prob += lpSum([ingredient_vars[i] for i in Ingredients]) == 3, "PercentagesSum"
prob += lpSum(ingredient_vars["CHICKEN"]) <= 2, "CHICKEN"
prob += lpSum(ingredient_vars["BEEF"]) <= 1, "BEEF"
prob += lpSum(ingredient_vars["MUTTON"]) <= 1, "MUTTON"
prob += lpSum(ingredient_vars["RICE"]) <= 1, "RICE"

Где 2,1,1,1 это количестводоступно для каждого сырья.

1 Ответ

0 голосов
/ 31 января 2019

Что не так, поскольку оно не удовлетворяет другим противоречиям. [Sic]

Какое ограничение нарушается?Если вы посмотрите, как вы определили ограничения, то обнаружите, что все они удовлетворены. 

prob += lpSum([proteinPercent[i] * ingredient_vars[i] for i in Ingredients]) >= 15.5, "ProteinRequirement"

То, что при решении 100 единиц BEEF означает, что у вас 2 * 100 = 200 единиц белка- намного больше требуемого 15,5. 

prob += lpSum([fatPercent[i] * ingredient_vars[i] for i in Ingredients]) >= 12.3, "FatRequirement"

Что при решении 100 единиц BEEF означает, что у вас 14 * 100 = 1400 единиц белка - намного больше требуемого 12,3.

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

...