GEKKO中使用BPOPT求解器时的奇异矩阵问题及求解器选型咨询

一、指定特定求解器的方法 #

当你设置solver=0时GEKKO会自动选择求解器，但要指定MINOS这类特定求解器，有两种简单直接的方式：

• 直接设置求解器编号：GEKKO给每个内置求解器分配了固定编号，MINOS对应的是编号1。你只需要在创建GEKKO模型后，直接设置solver参数即可：

  m = GEKKO()
  m.solver = 1  # 明确指定使用MINOS求解器
  

  顺便给你整理了常用求解器的编号对应表，方便你快速切换：

  • 0: 自动选择（默认行为）
  • 1: MINOS
  • 2: BPOPT
  • 3: IPOPT
  • 4: APOPT
  • 5: SNOPT
  • 6: CLP

• 添加求解器自定义选项：如果需要给MINOS设置额外参数（比如最大迭代次数、收敛精度），可以用solver_options属性来配置：

  m.solver_options = ['max_iter 1500', 'feasibility_tol 1e-7']
  

二、BPOPT奇异矩阵警告的处理建议 #

你遇到的*** WARNING MESSAGE FROM SUBROUTINE MA27BD *** INFO(1) = 3 MATRIX IS SINGULAR. RANK= 581错误，本质是求解器在矩阵分解时发现了线性相关的行/列，导致无法正常计算。这里给你几个排查和解决的方向：

• 先检查你的约束条件，看看有没有重复定义的约束、或者完全可以被其他约束推导出来的冗余等式约束，这些是导致矩阵奇异的常见原因。
• 尝试调整变量的初始值，不合适的初始点可能会让求解器进入奇异的数值区域，换个更贴合问题实际的初始值说不定能解决问题。
• 既然你已经发现GEKKO有多个求解器可选，直接切换到MINOS、IPOPT这类求解器也是个可行的方案——不同求解器的数值处理策略不同，可能能避开这个奇异问题。

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