Добавьте количество в качестве ограничения в задаче оптимизации смешивания - PullRequest
Новая
       8

Добавьте количество в качестве ограничения в задаче оптимизации смешивания

0 голосов
/ 01 февраля 2019

Я повторяю этот пример проблемы смешивания: https://www.coin -or.org / PuLP / CaseStudies / a_blending_problem.html

Со следующими данными: 

import pulp
from pulp import *
import pandas as pd

food = ["f1","f2","f3","f4"]
KG = [10,20,50,80]
Protein =       [18,12,16,18]
Grass = [13,14,13,16]
price_per_kg =  [15,11,10,12]

##            protein,carbohydrates,kg

df = pd.DataFrame({"tkid":food,"KG":KG,"Protein":Protein,"Grass":Grass,"value":price_per_kg})

И это код: 

deposit =  df["tkid"].values.tolist()

factor_volumen = 1



costs =  dict((k,v) for k,v in zip(df["tkid"],df["value"]))
Protein =  dict((k,v) for k,v in zip(df["tkid"],df["Protein"]))
Grass =  dict((k,v) for k,v in zip(df["tkid"],df["Grass"]))
KG =  dict((k,v) for k,v in zip(df["tkid"],df["KG"]))

prob = LpProblem("The Whiskas Problem", LpMinimize)
deposit_vars = LpVariable.dicts("Ingr",deposit,0)
prob += lpSum([costs[i]*deposit_vars[i] for i in deposit]), "Total Cost of Ingredients per can"



prob += lpSum([deposit_vars[i] for i in deposit]) == 1.0, "PercentagesSum"
prob += lpSum([Protein[i] * deposit_vars[i] for i in deposit]) >= 17.2, "ProteinRequirement"
prob += lpSum([Grass[i] * deposit_vars[i] for i in deposit]) >= 11.8, "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))

Эта часть просыпается, но мне нужно добавить еще одно противоречие - сколько килограммов я хочу произвести.

Я пытался сделать это с двумя ограничениями: 

## total KG produced == 14
prob += lpSum([KG[i] * deposit_vars[i] for i in deposit]) == 14, "KGRequirement"
### Can´t not use more that 8KG from deposit 1
prob += lpSum([KG[i] * deposit_vars[i] for i in deposit[0:1]]) <= 8, "KGRequirement1"

Я получаю эту ошибку: 

Status: Infeasible
Ingr_f1 = 0.83636364
Ingr_f2 = 0.11818182
Ingr_f3 = 0.045454545
Ingr_f4 = 0.0
Total Cost of Ingredients per can =  14.30000007

Но для этого можно использовать депозит 4,поэтому я думаю, что ограничение не является правильным.

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

prob += lpSum([KG[i] * deposit_vars[i] for i in deposit]) == 14, "KGRequirement"

И что средневзвешенное значение ингредиентов обеспечивает безопасностьа также

prob += lpSum([Protein[i] *KG[i] * deposit_vars[i] for i in deposit]) >= 17.2*14, "ProteinRequirement"

Это правое слово сейчас: 

prob += lpSum([Protein[i] *KG[i] * deposit_vars[i] for i in deposit]) >= 17.2*14, "ProteinRequirement"
prob += lpSum([Grass[i] *KG[i] * deposit_vars[i] for i in deposit]) >= 11.8*14, "FatRequirement"
prob += lpSum([KG[i] * deposit_vars[i] for i in deposit]) == 14, "KGRequirement"
prob += lpSum([KG[i] * deposit_vars[i] for i in deposit[0:1]]) <= 8, "KGRequirement1"

1 Ответ

0 голосов
/ 02 февраля 2019

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

prob += lpSum([KG[i] * deposit_vars[i] for i in deposit]) == 14, "KGRequirement"

И что средневзвешенное значение ингредиентов также обеспечивает безопасность этого условия.

prob += lpSum([Protein[i] *KG[i] * deposit_vars[i] for i in deposit]) >= 17.2*14, "ProteinRequirement"

Это правдоподобие: 

prob += lpSum([Protein[i] *KG[i] * deposit_vars[i] for i in deposit]) >= 17.2*14, "ProteinRequirement"
prob += lpSum([Grass[i] *KG[i] * deposit_vars[i] for i in deposit]) >= 11.8*14, "FatRequirement"
prob += lpSum([KG[i] * deposit_vars[i] for i in deposit]) == 14, "KGRequirement"
prob += lpSum([KG[i] * deposit_vars[i] for i in deposit[0:1]]) <= 8, "KGRequirement1"
...