Drake中调用Gurobi求解MIQP不可行时的约束定位与问题调试方法

一、Drake对不可行约束定位的支持现状 #

• 目前Drake已实现Gurobi求解状态的基础透传，你现在在gurobi.cc适配层能获取的GRB_INFEASIBLE、GRB_UNBOUNDED、GRB_INF_OR_UNBD三类核心返回状态，就来自这层原生映射逻辑。
• Drake暂未封装开箱即用的“精准定位不可行约束”功能——也就是说求解完成后不会自动在返回结果中标记导致不可行的具体约束项，也不会自动生成可读的不可行原因报告。
• 但Drake开放了两个关键能力，可对接Gurobi原生的不可行诊断逻辑，完全不需要脱离原有框架的优化搭建流程：
  • 支持透传所有Gurobi原生求解参数，可在求解触发前开启不可行极小冲突子集（IIS）计算逻辑，不需要改动原有问题的代码结构
  • 保持Drake侧约束、变量对象与Gurobi原生模型对象的索引一一映射，求解完成后可直接读取Gurobi模型的IIS计算结果，通过索引关联回Drake侧的约束句柄，实现不可行约束的定位

二、大型MIQP不可行问题的高效调试技巧 #

（以下方法均不需要依托搭建小型测试用例，适用于百/千级变量、万级约束的大规模模型调试）

• 分层预校验大幅缩小问题范围
  不要直接等完整MIQP求解跑完再判断不可行：在搭建完所有变量和约束后，临时将所有整数/二元变量松弛为连续变量，把目标函数设为常数0，先求解纯连续模型的可行性。这一步可快速区分不可行原因是来自连续域的约束冲突，还是来自整数/二元变量的离散约束矛盾，能将后续定位工作量减少至少一半。预校验阶段可开启Gurobi的数值精度参数NumericFocus=2，排查是否因模型系数量级差过大、缩放失衡引发的数值层面假不可行问题。
• 松弛变量法快速锁定冲突约束簇
  不需要逐个删除约束排查：对所有等式、不等式约束统一添加非负松弛变量，将原有约束形式从Ax ≤ b调整为Ax - s ≤ b, s ≥ 0，把目标函数修改为最小化所有松弛变量的求和值。求解该松弛模型后，松弛变量取值大于0的约束集合，就是存在冲突的约束范围。调试时可对不同约束的松弛变量配置差异化惩罚权重：对自定义物理约束、参数配置类约束设置高惩罚权重，对通用的边界、常规等式约束设置低惩罚权重，可进一步收窄需要排查的约束范围。
• 调用Gurobi原生IIS能力精准定位
  通过Drake的参数透传通道设置gurobi.ComputeIIS=1，求解触发时Gurobi会自动计算最小不可行约束子集。针对超大规模模型可调整参数IISMethod=0，使用启发式逻辑快速计算IIS，不需要追求数学意义上的最小不可行集，只要能定位核心冲突簇即可。拿到IIS输出的约束索引后，可直接通过Drake的接口取回对应的约束对象，查看约束的具体表达式、绑定的变量与参数。
• 排查常见的隐式不可行诱因
  在通过循环批量生成变量、约束的场景下，优先检查三类容易被忽略的矛盾点：
  • 变量的上下限设置是否存在逻辑颠倒，比如连续变量的下限值大于上限值、二元变量的边界设置与整数类型规则冲突
  • 约束中传入的参数值是否存在量级错位，比如米级单位的约束被传入毫米级数值、弧度制的角度约束被传入角度制数值，这类错位会直接导致可行域与预期完全不符
  • 整数变量的取值范围是否超出了离散空间的可覆盖范围，比如要求二元变量同时满足两个互斥的逻辑约束

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