如何修复Python GEKKO最优控制Cod中的“未找到解决方案”错误

Question

我试图重现K.Renee Fister和Jennifer Hughes Donnelly 2005年发表的论文“免疫疗法：最佳控制理论方法”图1中的结果。为此，我使用Python的GEKKO包编写了一个数值最优控制解算器。我使用了与论文中相同的初始条件、控制边界、参数值和模型方程。但是，当我运行代码时，我得到以下错误：

Traceback (most recent call last):
  File "xxxxxx", line 45, in <module>
    m.solve(disp=False) #solve
  File "xxxx", line 783, in solve
    raise Exception(response)
Exception:  @error: Solution Not Found

我希望程序的输出提供两个图：一个是ODE动态特性图，另一个是最优控制解决方案图。在

我尝试过用多种方法修改代码：修改目标函数、时间步数和改变最优控制模式，但是每次都会得到相同的错误。下面是我使用的代码：

^{pr2}$

此代码是通过修改this Youtube video中提供的示例GEKKO代码而派生的。任何帮助解决这个问题将不胜感激！在