如何查看GEKKO的变量、参数及中间值以用于调试？

1. 求解后查看变量与中间值 #

求解完成后，所有变量（m.Var）和中间变量（m.Intermediate）都会生成计算值，直接访问其.value属性即可查看：

• 单个变量/中间值：print(ro.value)、print(phi.value)
• 变量数组：如果ro是用m.Array(m.Var, n)创建的数组，可通过ro[i].value查看单个元素，或ro.value查看整个数组的所有值
• 示例代码：

  m.solve(disp=False)
  print('Phase shift phi values:', phi.value)
  print('ro array values:', ro.value)
  

2. 求解前验证参数与模型定义 #

求解前变量还未得到最优解，但可以验证初始值和常数，以及检查中间变量的表达式是否正确：

• 常数（m.Const）：直接访问.value即可查看，比如print(c.value)
• 变量初始值：变量默认初始值为0，若你设置了初始值（如ro = m.Array(m.Var, 5, value=1.0)），可通过ro.value查看初始值
• 中间变量表达式：因为求解前中间变量还未计算数值，可通过.expr属性查看其定义的表达式，确认逻辑是否正确，比如print(phi.expr)会输出c * ro的具体表达式结构

3. 保存/查看完整模型结构 #

如果需要全面检查模型的所有元素（变量、常数、中间变量、约束等），可以导出模型文件：

• 打开模型临时文件夹：调用m.open_folder()，会自动打开GEKKO生成临时文件的目录，其中的.gk0文件包含了完整的APMonitor格式模型定义，能看到所有模型元素的细节
• 保存模型到本地：使用m.save('my_model.apm')将模型保存为本地文件，直接打开该文件即可查看所有变量、常数、中间值的定义和约束关系

4. 逐行验证模型构建过程 #

在构建模型的代码中插入打印语句，逐段验证参数和变量定义：

• 定义完常数或变量后立即打印其值，比如：

  c = m.Const(2.5)
  print('Constant c:', c.value)
  ro = m.Array(m.Var, 3, value=0.1)
  print('ro initial values:', ro.value)
  

• 对于复杂表达式，拆分步骤后检查每一步的表达式，避免逻辑错误

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