如何在CPLEX中编写指定约束（已知前两个决策变量结果）

需求说明 #

已知CPLEX中两个决策变量的输出结果：

• first_beam_in_each_row[rows] = [1,5]（共2行，每行的首个波束编号为1和5）
• beam_nomchannel[beams] = [2,2,2,2,2,2,2,2]（8个波束每个占用2个信道）

需要编写约束使beam_firstchannel满足：

• 若波束i是某一行的首个波束（即i属于first_beam_in_each_row的取值），则beam_firstchannel[i] = 1
• 否则，beam_firstchannel[i] = beam_firstchannel[i-1] + beam_nomchannel[i-1]（当前波束的起始信道等于前一波束的起始信道加上前一波束占用的信道数）

修正后的完整代码 #

int Nbeams=8;
int Nchannels=16;
int Nrows=2;          

range beams=1..Nbeams;
range channels=1.. Nchannels; 
range rows=1..Nrows;

dvar int first_beam_in_each_row [rows]      in (1..Nbeams);       
dvar int beam_nomchannel[beams]             in   (channels);
dvar int beam_firstchannel[beams]           in   (channels);

// 假设已得到的变量输出值（可用于测试）
// first_beam_in_each_row = [1,5];
// beam_nomchannel = [2,2,2,2,2,2,2,2];

Subject to 
{
    // 对每个波束i设置约束
    forall(i in beams) {
        // 判断i是否是某一行的首个波束
        if (exists(r in rows) first_beam_in_each_row[r] == i) {
            beam_firstchannel[i] == 1;
        } 
        // 非首波束的情况，避免i=1时出现越界
        else if (i > 1) {
            beam_firstchannel[i] == beam_firstchannel[i-1] + beam_nomchannel[i-1];
        }
    }
}

关键说明 #

1. 首波束判断：用exists(r in rows) first_beam_in_each_row[r] == i检查当前波束i是否属于任意一行的首个波束，逻辑清晰且符合OPL语法
2. 约束逻辑拆分：分两种情况明确约束规则，同时增加i>1的判断，避免i=1时出现i-1=0的数组越界问题
3. 语法合规性：使用==定义约束等式，符合CPLEX的OPL建模规范；通过forall循环遍历所有波束，确保每个波束的约束都被覆盖

内容的提问来源于stack exchange，提问作者MO MO MA