Python Pulp - 无法对非方阵建模答案

【问题标题】：Python Pulp - Unable to Model Non-Square MatrixPython Pulp - 无法对非方阵建模
【发布时间】：2019-11-02 08:14:29
【问题描述】：

我在使用 excel 中的输入数组和 PULP 中的变量数组设置约束时遇到问题。

看来该模型仅适用于方阵，而我的最终代码有一个 365x24 的矩阵。下面的代码有一个 5x6 的矩阵，并且在运行说 index out of range 时会引发错误。

我已经提供了我目前正在使用的代码。

Sheet11 矩阵

91  37  36  38  33  16
1   36  59  29  23  4
25  74  72  39  69  1
22  68  48  70  12  41
98  86  75  16  99  12

Sheet12 矩阵

59  63  66  57  4   15
26  33  75  71  21  2
37  88  89  1   90  3
91  48  27  24  23  14
68  13  61  37  77  20

from pulp import *
from pandas import *
import pandas as pd
import numpy as np
import xlrd

model = pulp.LpProblem("Basic Model", pulp.LpMinimize)

YPER = 5
HE = 6

yearlyhours = []
yearlyhours = [(i,j) for i in range(YPER) for j in range(HE)]

book = xlrd.open_workbook('Stack.xlsx')
sheet11 = book.sheet_by_name('Sheet11')
sheet12 = book.sheet_by_name('Sheet12')
sheet13 = book.sheet_by_name('Sheet13')

TEST = [[sheet11.cell_value(i,j) for i in range(YPER)] for j in range(HE)]
YAPR = [[sheet12.cell_value(i,j) for i in range(YPER)] for j in range(HE)]
MAPR = [[sheet13.cell_value(i,j) for i in range(YPER)] for j in range(HE)]

YAHL = pulp.LpVariable.dicts("YAHL", (range(YPER), range(HE)), lowBound=0, cat='Continuous')
MAHL = pulp.LpVariable.dicts("MAHL", (range(YPER), range(HE)), lowBound=0, cat='Continuous')

##OBJECTIVE##
model += pulp.lpSum([YAPR[i][j] * YAHL[i][j] + MAPR[i][j] * MAHL[i][j] for i in range(YPER) for j in range(HE)]), 'Sum_of_Value'

for i,j in yearlyhours:
    model += pulp.lpSum([YAHL[i][j] + MAHL[i][j]]) == ([TEST[i][j]])

LpSolverDefault.msg = 1
model.writeLP('Opt.lp')
model.solve()
pulp.LpStatus[model.status]
print("Status:", LpStatus[model.status])
obj = value(model.objective)
print("Total Cost: ${}".format(obj,2))
print('\n')

当更改为非方阵时，会给出以下错误消息：

IndexError: 列表索引超出范围

感谢任何帮助，因为我仍在学习。

【问题讨论】：

如果您在这里没有得到答案，您可以考虑在Operations Research（或.stackexchange.com）上提问。
从错误消息的声音中，您尝试引用列表末尾之外的元素。找出哪一行引发了该错误，然后打印出列表的大小以及您引用的最大索引。请记住，python 是零索引的，因此 5 元素列表的最后一个元素是元素 [4]。

标签： python optimization pulp

【解决方案1】：

您要读取数据的代码产生 6 行和 5 列。变量字典结构是 5 个键，然后内部 dict 有六个键。因此，当您尝试同时迭代两者时，会出现索引错误。方阵之所以有效，是因为 5 x 5 与 5 x 5 相同，所以倒序无关紧要。

>>>YAHL
{
 0: {
   0: YAHL_0_0,                                                                                                         
   1: YAHL_0_1,                                                                                                            
   2: YAHL_0_2,                                                                                                            
   3: YAHL_0_3,                                                                                                            
   4: YAHL_0_4,                                                                                                            
   5: YAHL_0_5
 },                                                                                                          ...                                                                                                          4: {
  0: YAHL_4_0,                                                                                                         
  1: YAHL_4_1,                                                                                                            
  2: YAHL_4_2,                                                                                                            
  3: YAHL_4_3,                                                                                                            
  4: YAHL_4_4,                                                                                                            
  5: YAHL_4_5}} 
 }
}
>>>TEST
[[91.0, 1.0, 25.0, 22.0, 98.0],
 [37.0, 36.0, 74.0, 68.0, 86.0],
 [36.0, 59.0, 72.0, 48.0, 75.0],
 [38.0, 29.0, 39.0, 70.0, 16.0],
 [33.0, 23.0, 69.0, 12.0, 99.0],
 [16.0, 4.0, 1.0, 41.0, 12.0]]

>>>for i in range(YPER):
       for j in range(HE):
           print(TEST[i][j])
           print(YAHL[i][j])
91.0
YAHL_0_0
1.0
YAHL_0_1
25.0
YAHL_0_2
22.0
YAHL_0_3
98.0
YAHL_0_4
---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
<ipython-input-14-a9538e188b1e> in <module>
      1 for i in range(YPER):
      2     for j in range(HE):
----> 3         print(TEST[i][j])
      4         print(YAHL[i][j])
      5

IndexError: list index out of range

【讨论】：