Pyomo“数据=无”输出错误

Question

非常感谢您调查我的问题！我是 pyomo 新手，现在尝试用它来解决具体模型。

这是我的部分代码：

——（已编辑）

def objective_rule(model):
        ans = sum(model.DAEB[t] * model.DAEP[t] for t in model.t)
        ans -= sum(model.DARUP[t] * model.RU[t] + model.DARDP[t] * model.RD[t] for t in model.t)
        ans += sum(0.5 * (sum(model.penalty[t, w]) + sum(-model.RTEP[t] * (model.DAEB[t] - ((model.RTRD[t, w] * model.RU[t]) +
               sum(model.veh_pwer_dem[t, v, w] for v in model.v))))) for t in model.t for w in model.w)
        ans += 0.95 * (model.epslon + 1 / (0.5) * sum(0.01 * model.miu[w] for w in model.w))
        return ans

• t 定义为一组 [0:23]，一天中的时间小时；
• w被定义为[0,1]，是不同的场景；
• v 是 [0:29]，30 个不同的样本。

在定义目标函数之前，我还定义了一堆与上述模型变量相关的约束。

除了最后一步，一切都很顺利，pyomo给出：

——（已编辑）

Traceback (most recent call last):
  line 61, in <module>
    myresult = result.solve(project, pricefile, reg_dispatch, SOC=0, SOC_margin=0.05)

  line 56, in solve
    opti_model, result = self.Fr_optimal_bidding_optimization(self.vehicles, pricefile, reg_dispatch, SOC)

  line 347, in Fr_optimal_bidding_optimization
    model.objective = Objective(rule=objective_rule, sense=minimize, doc='minimize the total cost')

  line 483, in __setattr__
    self.add_component(name, val)

  line 849, in add_component
    val.construct(data)

  line 307, in construct
    tmp = _init_rule(_self_parent)

  line 337, in objective_rule
    sum(model.veh_pwer_dem[t, v, w] for v in model.v))))) for t in model.t for w in model.w)

ERROR: Rule failed when generating expression for objective objective:
  line 337, in <genexpr>
    TypeError: '_GeneralVarData' object is not iterable
    sum(model.veh_pwer_dem[t, v, w] for v in model.v))))) for t in model.t for w in model.w)

ERROR: Constructing component 'objective' from data=None failed:
TypeError: '_GeneralVarData' object is not iterable
    TypeError: '_GeneralVarData' object is not iterable

首先，我对_GeneralVarData 对象感到非常困惑，不知道它指的是什么，其次，我的目标函数看不到问题，所以如果有人可以提供任何类型的建议将不胜感激！

（我也是“stackoverflow”的新手，如果我的问题没有明确说明，请在前面道歉！）

腾

Answer 1

异常是由倒数第二行中的以下代码触发的：

sum(model.veh_pwer_dem[t, v, w])

这试图将（标量）model.veh_pwer_dem[t, v, w] 视为可迭代对象并对其键（索引）求和。您的表达式的逻辑结构似乎还存在其他问题。我强烈建议将表达式分成更小的块，这将有助于您保持清晰。例如：

def objective_rule(model):
    ans = sum(model.DAEB[t] * model.DAEP[t] for t in model.t)
    ans -= sum(model.DARUP[t] * model.RU[t] + model.DARDP[t] * model.RD[t] for t in model.t)
    ans += sum(0.5 * ( sum(model.penalty[t, w] for t in model.t for w in model.w) + # ...
    #...
    return ans

关于_GeneralVarData 对象是什么。 _GeneralVarData 对象是一个内部 Pyomo 结构，用于表示作为索引 Var 组件的一部分的实际优化变量（您可能偶尔还会看到 _SimpleVar，这是用于非索引 Var 组件的内部实现）。作为 Pyomo（以及一般 Python）中的一项规则，以下划线 (_) 开头的方法、属性和类对于实现来说是“私有的”。