CPLEX官员委任优化模型报‘Not an array type’错误求助

我构建了一个调用Excel表格数据的CPLEX模型，用于在资源限制条件下优化委任官员数量。目前已解决其他所有错误，但仍存在一处「Not an array type」错误，位置在代码的Progression Dynamics部分（已在代码中标记）。用尽入门知识和AI建议都未解决，以下是完整代码：

// 集合定义
{string} MSL = ...;
{string} Type = ...;
{string} FY = ...;

// 参数
float initPop[MSL][Type] = ...;              // 各层级、类型的初始人数
float gain[MSL][Type] = ...;        // 每年各层级、类型的外部新增人数
float progressRate[MSL][Type] = ...;      // 晋升到下一层级的比例
float retainMSL[MSL][Type] = ...;  // 未晋升但留任的比例
float CmsnRate[MSL][Type] = ...; // 委任比例
float costPerYear[Type] = ...;   // 每类人员的年度成本

float budgetLimit = ...;
int capacityLimit = ...;
int mission = ...;

// 决策变量
dvar int+ cadet[MSL][Type][FY];       // 每年各层级、类型的受训人员数量
dvar int+ commissions[Type][FY];             // 每年各类型的委任官员数量

// 目标函数
maximize sum(t in Type, y in FY) commissions[t][y];

// 约束条件
subject to {

   // 人员晋升动态约束
  forall(y in FY: y != "2026") {
    forall(l in MSL: l != "MSL1") {
      forall(t in Type) {
        // 从上层级晋升而来的人员
        cadet[l][t][y] ==   // 此处为「Not an array type」错误位置
            progressRate[prev(MSL,l)][t] * cadet[prev(MSL,l)][t][prev(FY, y)]
          + retainMSL[l][t] * cadet[l][t][prev(FY, y)]
          + gain[l][t][y];
      }
    }
    // MSL1层级仅包含留任人员和外部新增
    forall(t in Type) {
      cadet["MSL1"][t][y] ==
          retainMSL["MSL1"][t] * cadet["MSL1"][t][prev(FY, y)]
        + gain["MSL1"][t];
    }
  }

  // 定义委任官员为从MSL4和MSL5层级晋升的人员
  forall(t in Type, y in FY)
    commissions[t][y] == CmsnRate["MSL4"][t] * cadet["MSL4"][t][y] + CmsnRate["MSL5"][t] * cadet["MSL5"][t][y];

  // 年度预算约束
  forall(y in FY)
    sum(l in MSL, t in Type)
      costPerYear[t] * cadet[l][t][y] <= budgetLimit;

  // 年度容量约束
  forall(y in FY)
    sum(l in MSL, t in Type)
      cadet[l][t][y] <= capacityLimit;

  // 年度最低委任人数约束
  forall(y in FY)
    sum(t in Type) commissions[t][y] >= mission;

}

错误原因分析 #

1. 集合未声明为有序类型：OPL的prev()函数仅对有序集合生效，你定义的MSL和FY是普通字符串集合，没有明确顺序，导致prev(MSL,l)和prev(FY,y)无法返回合法的集合元素，进而让数组索引类型非法，触发错误。
2. 参数维度不匹配：gain参数定义为二维数组gain[MSL][Type]，但约束中使用了三维索引gain[l][t][y]，维度不匹配也会引发类型错误。
3. 初始值约束缺失：没有处理第一年（FY的首个元素）的cadet初始值，导致循环中访问prev(FY,y)时可能出现越界问题。

解决方案 #

1. 将集合声明为有序集合 #

修改集合定义，添加ordered关键字明确顺序：

{string} MSL = ordered {"MSL1", "MSL2", "MSL3", "MSL4", "MSL5"};
{string} FY = ordered {"2022", "2023", "2024", "2025", "2026"};

如果从Excel读取集合，需确保Excel中数据按顺序排列，再用ordered声明。

2. 修正参数维度 #

将gain参数修改为三维数组，匹配约束中的索引使用：

float gain[MSL][Type][FY] = ...;

3. 补充初始年份约束 #

添加初始年份的人员数量约束，覆盖循环未处理的首个FY元素：

// 初始年份人员数量约束
forall(l in MSL, t in Type)
  cadet[l][t][first(FY)] == initPop[l][t];

内容的提问来源于stack exchange，提问作者Brent Huber