Ampl Cplex QP报错：Hessian非半正定且无基，求排查建议

首先，这个错误本质是CPLEX的QP求解器要求目标函数的Hessian矩阵（二次项的二阶导数矩阵）必须半正定（对应凸优化问题），如果矩阵不定，求解器无法找到合法的基来启动求解。结合你的产销运利润优化模型，给你一套逐步排查的方案：

1. 先锁定目标函数的二次项问题 #

很多时候这个报错源于目标函数的二次项设计不合理：

• 先把目标函数里的所有二次项（比如产量的平方项、变量相乘项）暂时注释掉，换成线性目标再尝试求解。如果能正常运行，说明问题肯定出在目标的二次项上。
• 检查二次项的系数：比如生产的边际成本递增（二次项系数为正）会让Hessian正定，这是合法的；但如果是需求曲线带来的销售额（比如price[i] * sales[i]，而price[i] = a - b*sales[i]，最终销售额是a*sales[i] - b*sales[i]^2），这里的二次项系数是负的，会让Hessian不定（目标函数是凹的）。这种情况下，你需要告诉CPLEX允许求解非凸QP，在AMPL里添加参数设置：

  option cplex_options 'allow_nonconvex_qp 1';
  option cplex_options 'qpmethod 2';  # 切换到内点法，对非凸问题兼容性更好
  

2. 逐行排查约束中的隐性非线性 #

你提到怀疑非线性约束，但可能有些非线性是你没注意到的：

• 用AMPL命令show constraints;把所有约束的完整表达式打印出来，逐行检查：
  • 有没有变量相乘的项？比如transport[i,j] * cost_per_unit[i,j]如果cost_per_unit是变量（而不是参数），这就是二次项，会把约束变成二次约束，导致问题成为QCQP（二次约束二次规划），如果约束是非凸的，也会触发Hessian相关错误。
  • 有没有不小心用了^、log、exp这类非线性函数？即使是变量的平方项（比如x[i]^2 <= ...）也是二次约束，需要确认这类约束是否是凸的（比如平方项系数为正的话是凸约束）。
• 用AMPL的check;命令做自动诊断，它会提示约束的凸性问题、变量上下界冲突等，很多隐性问题都能在这里找到线索。

3. 逐步简化模型定位冲突点 #

如果上面的方法没找到问题，就用“减法排查法”：

• 先把模型简化到最小可行版本：比如只保留一个产地、一种产品（沙发或椅子），只保留核心的生产、运输、销售约束，去掉所有非必要的约束和变量。如果这个简化模型能正常求解，再逐步添加变量组、约束组，每次添加后运行求解，直到报错出现——这时候刚添加的内容就是冲突来源。
• 另外，你之前尝试移除等号，不妨试试完全注释掉某一类约束（比如所有销售约束，或者所有运输约束），而不只是改等号，看是否能解决问题，这样能快速缩小范围。

4. 检查变量定义与求解器参数 #

• 确认所有变量的类型正确：比如连续变量有没有被误设为整数？变量的上下界是否合理（比如有没有设x[i] >= 1000但实际生产能力只有500，导致可行域为空？不过可行域为空的报错通常不一样，但也可能间接影响Hessian的构建）。
• 除了前面提到的allow_nonconvex_qp，还可以尝试调整CPLEX的其他参数，比如preind 0关闭预处理器，看是否能绕过某些预处理导致的误判。

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