如何在Gekko中为变量分配无穷大值？求解报错求助

Gekko不支持将Inf（无论numpy还是Python原生）作为变量初始值，因为求解器仅能处理有限数值输入。针对你的核心需求，提供以下三种可行方案：

1. 用极大有限值替代无穷大 #

选择一个足够大的有限值（如1e10）近似无穷大，既能让Q=(T[1]-T[0])/R在首次迭代时趋近于0，也能被求解器正常识别：

from gekko import GEKKO
model = GEKKO()
# 用极大值模拟无穷大效果
R = model.FV(value=1e10)
T = model.Array(model.Var, 2)
Q = model.FV()
model.Equation(Q == (T[1]-T[0])/R)
model.solve()

2. 直接强制Q初始值为0 #

如果核心目标是首次迭代时Q严格为0，无需通过R间接实现，可直接固定Q的初始值：

from gekko import GEKKO
model = GEKKO()
R = model.FV(value=1) # 后续测试场景可灵活调整此值
T = model.Array(model.Var, 2)
# fixed_initial=True确保求解器首次迭代不修改Q的初始值
Q = model.FV(value=0, fixed_initial=True)
model.Equation(Q == (T[1]-T[0])/R)
model.solve()

后续测试其他场景时，只需将fixed_initial设为False或调整R的数值即可。

3. 重写方程避免分母 #

将分式方程转化为乘法形式，彻底规避分母相关的数值问题，同时提升求解器的计算稳定性：

from gekko import GEKKO
model = GEKKO()
R = model.FV(value=1e10)
T = model.Array(model.Var, 2)
Q = model.FV()
# 用乘法形式替代分式，避免无穷大相关报错
model.Equation(Q * R == T[1] - T[0])
model.solve()

这种写法更适合后续测试不同R值的场景，不会出现分式带来的数值波动。

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