JuMP固定变量后获取@NLexpression值报Nothing has no field status如何解决

错误原因 #

你调用JuMP.value(k1)报错的核心原因是未触发求解器计算流程：Juniper求解器实例在你调用optimize!()前没有生成任何求解结果，相关状态字段为空，JuMP.value()尝试从空结果中读取数值就会抛出type Nothing has no field status错误。即便所有变量都已固定为定值，JuMP也不会主动触发求解器计算，需要手动调用求解命令。

解决方案 #

根据你的使用场景可选择两种方案：

• 方案1：直接原生计算（无需调用求解器，推荐）
  你所有变量和参数都是固定值，完全不需要走优化求解流程，直接用Julia原生语法计算表达式即可：

  # 直接代入数值计算，不需要调用JuMP接口
  hp_val = xhat[3]
  Lp_val = xhat[1]
  np_val = xhat[2]
  k1_val = hp_val / (1 - f(Lp_val - sig*sqrt(np_val)) + f(-Lp_val - sig*sqrt(np_val)))
  

  该方法运行效率最高，不会触发求解器相关的任何报错。

• 方案2：先调用求解器再取数
  如果你需要保留JuMP的接口调用逻辑，在定义完k1之后补充optimize!()触发求解，再调用JuMP.value()即可：

  # 原有定义k1的代码后补充：
  optimize!(ALT)
  # 可选：先校验求解状态，避免其他异常
  if termination_status(ALT) in [MOI.OPTIMAL, MOI.LOCALLY_SOLVED]
      k1_val = JuMP.value(k1)
      println("k1取值为：", k1_val)
  else
      println("求解异常，终止状态为：", termination_status(ALT))
  end
  

补充说明 #

你用@variable(ALT, hp == xhat[3])的方式固定变量是可行的，等价于给变量添加了等式约束；如果后续需要动态调整变量固定值，也可以用fix(hp, 新数值; force=true)函数实现，更灵活。

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