GEKKO与模式4和模式7的结果不一致

Question

用2个节点解决问题。

给定问题

    m = GEKKO()
    m.time = np.linspace(0, 5, 2)
    x = m.Var(1)
    m.Equation(x.dt() == 1)

如果使用模式7求解

    m.options.IMODE = 7
    m.options.SOLVER = 2
    m.solve(disp=False)

    x = [1.0, 1.0]

如果使用模式4求解

    m.options.IMODE = 4
    m.options.SOLVER = 2
    m.solve(disp=False)
    
    x = [1.0, 6.0]

Answer 1

当IMODE=7只有一个时间步长时，Gekko报告解决方案的方式似乎存在一个错误。切换到两个时间步长为IMODE=7和IMODE=4提供了一致和正确的解决方案。

from gekko import GEKKO
import numpy as np
m = GEKKO()
m.time = np.linspace(0, 5, 3)
x = m.Var(1)
m.Equation(x.dt() == 1)
m.options.IMODE = 7
m.options.SOLVER = 2
m.solve(disp=False)
print(x.value)

当为IMODE=4时，解决方案也是[1.0, 3.5, 6.0]。