CPLEX OPL代码问题咨询：不等式约束触发内部错误排查

问题代码 #

int nm=4;
int r=5;
range num=1..nm;
range hr=1..r;
int scale=100;
float m[hr]=[1026.8,1061.65,1096.5,1068.35, 1040.2];
float n[hr]=[736,713,690,667,644];
// Declaration of Decision variables
dvar float+ e[num][hr];
dvar float+ f[num][hr];
dvar float c[hr];
dvar float d[hr];
using CP;
  dexpr float b[i in hr]=(((m[i]+c[i]-d[i])^2)+n[i]^2);
  dexpr float nb[i in hr]=sqrt(b[i]);
  dexpr float amax=max(i in hr)(b[i]);
  dexpr float nf=1/sum(i in hr)nb[i];
            
// Objective
minimize  amax*nf;
subject to {
 forall(i in hr){
   c[i]==sum(k in num)e[k][i];
   d[i]==sum(k in num)f[k][i];
 } 
 forall (i in num){
   forall(j in hr){
     0<=e[i][j]<=1;
     0<=f[i][j]<=1;
     e[i][j]+f[i][j]<=2;
   }
 }
}

问题现象 #

将约束e[i][j]+f[i][j]==2改为e[i][j]+f[i][j]<=2时触发内部错误，但等式约束下可正常求解。

原因分析 #

1. 求解器选型不匹配：代码中使用using CP;调用CP Optimizer求解器，该求解器主打离散/组合优化场景，对连续型非凸优化问题的处理能力有限。
2. 约束放松后的问题特性变化：
   • 当约束为e[i][j]+f[i][j]==2时，结合0<=e[i][j]<=1和0<=f[i][j]<=1的上下界，e[i][j]和f[i][j]只能取1，变量取值被固定，目标函数的非线性复杂度大幅降低，求解器可正常处理。
   • 当约束改为<=2时，变量拥有更多取值空间，目标函数包含平方、平方根、倒数、最大值等非线性操作，组合形成非凸优化问题，CP Optimizer处理这类问题时会出现数值不稳定或内部逻辑错误，最终触发报错。

解决方法 #

• 更换求解器：删除代码中的using CP;语句，改用CPLEX默认的连续优化求解器（支持QCP/非凸问题求解），该求解器对连续非线性问题的适配性更强。
• 优化目标函数建模：引入辅助变量重新表达目标函数（比如用辅助变量替换amax，将最大值约束转化为一组不等式），降低求解器处理非线性项的复杂度，提升求解稳定性。

内容的提问来源于stack exchange，提问作者SUBHADARSHINI PANDA