如何通过for循环在gekko中制作消耗品？

问题描述

根据a related question中的提示，将字典用作gekko的源输入，我做了以下工作；

import pandas as pd
from gekko import GEKKO
my_vars = ['x1','x2','x3']

# stored as dictionaries

Cost = {'x1':100,'x2':125,'x3':80}
Min = {'x1':0,'x2':0,'x3':0}
Max = {'x1':70,'x2':40,'x3':15}

LP = GEKKO(remote=False)
va = LP.Array(LP.Var,(len(my_vars)))  # array
vd = {}                                # dictionary
for i,xi in enumerate(my_vars):
    vd[xi] = va[i]
    vd[xi].lower = Min[xi]
    vd[xi].upper = Max[xi]

# Cost function
LP.Minimize(LP.sum([Cost[xi]*vd[xi] for xi in my_vars])) 

# This also works as a dictionary
LP.Equation(LP.sum([vd[xi] for xi in my_vars])==100)

LP.solve(disp=True)

for xi in my_vars:
    print(xi,vd[xi].value[0])
print ('Cost: ' + str(LP.options.OBJFCNVAL))

成功了。

但是假设我的数据帧格式的变量系数为LHS，标准数据为RHS，

例如；

Coef = [['x1',33,8.8,0.2,0.1],['x2',24,4.5,2,0.5],['x3',82,'','']]
LHS = pd.DataFrame(Coef,columns=['Name','p1','p2','p3','p4'])
LHS.set_index('Name')
and
goals = [['p1',30],['p2',5],['p3',0.7],['p4',0.2]]
RHS = pd.DataFrame(goals,columns=['property','goal'])

我想在for循环中使用LHS和RHS数据动态构造约束。 我做了如下，

for i in range(len(RHS)):
    LP.Equation(LP.sum([LHS.iloc[list(my_vars).index(xi),i]*vd[xi] for xi in my_vars]) >= RHS.iloc[i,1])

但是有错误。

---

APMonitor，版本0.9.2
APMonitor优化套件

---

 --------- APM Model Size ------------  

每个时间步包含
对象：6
常数：0
变量：28
中间体：0
连接数：24
式：22
残渣：22
@error：模型表达式***错误 函数字符串的语法：无效元素：x1

位置：7 v9-（（（（x1）*（v1））） ？

回溯（最近通话最近一次）：

文件“”，第1行，在 LP.solve（disp = True）

文件 “ /opt/anaconda3/envs/py37/lib/python3.7/site-packages/gekko/gekko.py”， 2130行，在解决 引发异常（apm_error）

异常：@error：模型表达式***函数语法错误 字符串：无效元素：x1

位置：7 v9-（（（（x1）*（v1））） ？

我的代码有什么问题？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）