对于Python GEKKO中基于ODE或PDE的生态系统模型，最合适的求解方法是什么？
对于Python GEKKO中基于ODE或PDE的生态系统模型，最合适的求解方法是使用GEKKO自带的求解器。GEKKO是一个用于动态优化和非线性模型预测控制的Python库，它提供了一个高效的求解器来解决ODE和PDE问题。
GEKKO的求解器基于不同的算法和技术，包括非线性优化、动态优化、动态系统仿真等。它可以自动选择最合适的求解器来解决特定的问题，并提供了一些参数来调整求解器的性能和精度。
对于基于ODE的生态系统模型，GEKKO的求解器可以使用常见的数值方法，如欧拉法、龙格-库塔法等。这些方法可以在一定程度上近似解决ODE问题，并提供了一定的精度和稳定性。
对于基于PDE的生态系统模型，GEKKO的求解器可以使用有限元方法、有限差分方法等。这些方法可以将PDE问题离散化为一系列的代数方程，并通过求解这些方程来得到数值解。
总之，对于Python GEKKO中基于ODE或PDE的生态系统模型，最合适的求解方法是使用GEKKO自带的求解器。它提供了一系列的数值方法和技术，可以根据具体的问题选择最合适的求解器，并提供了一些参数来调整求解器的性能和精度。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云计算服务：https://cloud.tencent.com/product
- 腾讯云数据库：https://cloud.tencent.com/product/cdb
- 腾讯云服务器：https://cloud.tencent.com/product/cvm
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云物联网：https://cloud.tencent.com/product/iot
- 腾讯云移动开发：https://cloud.tencent.com/product/mobdev
- 腾讯云存储：https://cloud.tencent.com/product/cos
- 腾讯云区块链：https://cloud.tencent.com/product/baas
- 腾讯云元宇宙：https://cloud.tencent.com/product/mu
相关·内容
我已经找了一段时间了，但是在任何地方都找不到这个问题的答案，如果它是重复的，很抱歉！ 我已经开始构建一个基于xarray-simlab framework的python package，目标是提供一个模块化工具箱，用于构建可重现的、灵活的海洋生态系统模型。Xarray-simlab目前只支持显式步长来求解模型函数。为了更安全高效地解决复杂模型，我开始使用GEKKO作为求解器后端，因为模型语法似乎很适合。(注意:目前，我只需要随时间推移求解模型方程的功能，但我希望在以后的阶段利用GEKKO的优化功能将模型参数拟合到现场或实验室数据。) 包的当前原型创建了一个xsimlab流程类，该流程类将GE
浏览 18提问于2020-06-26得票数 2
回答已采纳
我正在尝试使用SymPy将一些枫树代码转换为python。下面是我要转换的枫树代码：
下面是我尝试用python代码来做同样的事情，但这是我所得到的。
from sympy import *
c, k, x, t, tau, Y, xi = symbols('c k x t tau Y xi')
u = Function('u')
F = Function('F')
xi = k*(x-c*t)
U = Function('U')
pde = diff(u(x,t), t) + 6*u(x,t)*diff(u(x,t)
浏览 5提问于2022-09-11得票数 1
回答已采纳
在使用GEKKO模拟进程时(例如，在示例15 here中)，如何设置状态变量的导数的初始值？我正在使用IMODE=4，但我也可以使用IMODE=7。 编辑我已经使用IMODE=5将ODE模型的参数与测量的输入和输出进行了拟合，我想预测超出测量时间点的模型输出。
浏览 13提问于2021-10-04得票数 1
为了解决一组复杂的微分方程，我需要在我的模型中包含时间。下面是一个简单的问题，演示了常量k=0.1和初始条件y(0)=10的问题。 ? 我在Python Gekko中尝试过，但不知道如何将时间作为变量。在Scipy ODEINT中，函数有时间和状态变量。在Gekko中，我将m.time定义为我希望看到解决方案的点，但在等式中使用m.time会出错。 import numpy as np
from gekko import GEKKO
import matplotlib.pyplot as plt
m = GEKKO() # create GEKKO model
m.time
浏览 24提问于2020-04-29得票数 2
1回答
我面对的是一个Hessenberg索引-2 DAE，我试图用python模块gekko来解决它。经过几天的尝试和错误之后，我认为我离一个能工作的代码不远了。但我刚刚意识到，也许gekko不能处理复数？
下面是一个最低限度的工作示例：
import numpy as np
from gekko import GEKKO
# Define the simulation and its parameters
g = GEKKO()
g.options.IMODE = 7
g.options.NODES = 1
# define the time array
n_steps = 100
Time
浏览 15提问于2022-05-09得票数 2
回答已采纳
在Python GEKKO中顺序求解模型(例如，使用IMODE >= 4)在使用GEKKO附带的max2和max3函数时失败。 这适用于np.maximum或标准max函数将GEKKO参数视为数组的用例，这并不总是预期的用法，或者在与整数进行比较时可能会产生错误。 最小代码示例： from gekko import GEKKO
import numpy as np
m = GEKKO()
m.time = np.arange(0,20)
y = m.Var(value=5)
forcing = m.Param(value=np.arange(-5,15))
m.Equation(
浏览 25提问于2020-12-02得票数 2
回答已采纳
1回答
我从Matlab中找到了一个调用python的函数。所以，我试着用这个函数从Matlab调用GEKKO包来求解ODE。我在将gekko选项(remote=true)定义为函数或变量时出错了。
我用Anaconda提示打开了Matlab。
m = py.gekko.GEKKO(remote==True);
m.time = py.numpy.linspace(0,20,100);
k = 10;
y = m.Var(5.0);
t = m.Param(m.time);
m.Equation(k*y.dt()==-t*y);
m.options.IMODE = 4;
m.solve(disp==T
浏览 0提问于2019-06-17得票数 4
回答已采纳
1回答
我正在使用GEKKO来解决Python语言中的MINLP问题，该问题涉及到与电力系统仿真软件PowerFactory的联合仿真。MINLP问题目标函数需要从内部调用PowerFactory，并在将值返回给MILP问题定义之前执行特定的任务。在等式约束定义中，我还需要使用两个Pandas数据框(具有10000x2单元格的表)来获取与约束表达式中决策变量的值相对应的特定值。但是，当我在目标/约束函数中编写涉及循环或Pandas 'loc‘或'iloc’函数的任何代码时，由于这些函数调用中GEKKO变量的一些问题，我会遇到错误。在这方面的任何帮助都将非常感谢。代码结构如下：
浏览 33提问于2020-08-01得票数 2
2回答
是否有任何免费的开源硬集成器Java库？
我目前正在使用来自库的Dormand 8(5,3)集成器。我正在为我的模型实现一个DAE近似，这样我就可以使用这个积分器，但是我想尝试一个在直线ODE系统上的刚性求解器来进行比较。
我到处找过，但似乎找不到任何直接的答案。类似于MATLAB中的ode15s求解器吗？
浏览 8提问于2015-05-11得票数 3
回答已采纳
1回答
神经网络MPC
、、、
因此，我想使用神经网络作为我在Python中学习的MPC控制的动态模型函数。我还没有找到像CASADI，GEKKO，do-mpc这样的开源优化包做这件事的例子/文档？有没有人有一些参考/建议来实现这一点？谢谢 编辑01 a)我已经尝试过CASADI + tensorflow模型CASADI有一个关于如何在CASADI中使用tensorflow模型的博客。我完全不确定我的实现是否正确，因为很明显，我没有得到预期的结果。b)在互联网上可以看到“mpc.pytorch”库，这是一个提供nn模型的mpc工具箱。不确定其性能C)基于CASADI的do-mpc正计划集成NN模型。d)正如@john ge
浏览 231提问于2021-01-14得票数 2
回答已采纳
我正在尝试Gekko的MINLP问题，然而，可用的求解器APOPT对解决我的问题没有帮助。可以在Gekko Python中使用SCIP解算器吗？
浏览 50提问于2020-07-06得票数 2
有没有人尝试用Modelica实现Navier Stokes偏微分方程(PDE)？我发现了空间基函数(SBF)的方法，它通过数值修改得到Dymola可以处理的常微分方程(ODE)。
致以敬意，
维克多
浏览 2提问于2015-06-17得票数 4
1回答
我正在学习如何使用 python包。作为第一步，我想模拟一个简单的向量ODE: dx/dt = A*x，其中A是矩阵，x是向量。我为GEKKO看到的所有ODE示例都是标量ODEs，在线中的数组示例没有说明在声明方程时如何合并.dt()方法。我知道，在声明方程时，可以使用列表，所以我认为这样的方法是可行的：
import numpy as np
from gekko import GEKKO
m=GEKKO()
m.time=np.linspace(0.,1.,10)
N=5
A=np.ones([N,N])
x=np.ones(N)
x=m.Var(value=x)
A=m.Param(val
浏览 1提问于2019-09-17得票数 1
回答已采纳
我试图在python2或3中运行gekko，我只需遵循一个修改参数的教程，并使用scipy.integrate来模拟我的MPC。它在我的一台计算机上工作，但不工作在我的NVIDIA TX2上。并在运行m.solve(disp=False)时得到"Exec格式错误“。
这两台计算机都有python2和python3。因为我也在运行ROS，所以我想使用python2来运行脚本。最初，我认为可能是使用python2解释导致了问题，所以我编写了另一个脚本，使用子进程将我的mpc和模拟器解释为python3。然而，这一问题仍然存在。问题似乎是在gekko包内。我确信这是一个环境错误，因为脚本在我
浏览 1提问于2019-09-17得票数 1
回答已采纳
1回答
我对编程很陌生，我的第一个Python项目是非线性编程。我正在使用的Gekko优化套件，我有一切正常运行，但需要指导如何准确地在本地运行它。下面是文档所提供的解释和代码，但是我需要一些帮助来帮助我自己如何准确地完成它，以及它到底意味着什么。请表现得像在向小孩子或金毛猎犬解释一样。
run目录m.path包含模型文件gk0_model.apm以及远程(默认)或本地(m=GEKKO(remote=False))运行优化问题所需的其他文件。使用m.open_folder()打开run目录。run目录还包含诊断文件，如如果解决程序找不到解决方案时生成的infeasibilities.txt。可以更改
浏览 7提问于2022-09-20得票数 2
我如何才能从GEKKO获得关于求解我的模型需要多长时间的输出？我知道基于，我可以让我的代码打印运行代码所用的总时间，但我不知道如何隔离求解器时间。
from gekko import GEKKO
import time
start = time.time()
m = GEKKO(remote=False)
x = m.Var(value=0)
y = m.Var(value=1)
m.Equations([x + 2*y==0, x**2+y**2==1])
s1 = time.time()
m.solve(disp=True)
e1 = time.time()
print([x.value[
浏览 19提问于2019-11-22得票数 4
回答已采纳
我认为这个问题有点长，所以我认为最好先考虑一下简化的版本：
有两个类A和B.B继承自A。B (add)中有一个成员函数需要使用A中的一个成员函数来运行。
class A;
typedef int(A::*operation)(int c);
typedef void (A::*myfunction)(operation, int);
class A
{
public:
int a;
int b;
int do_something(myfunction f, operation op)
{
(this->*f)(op, 1);
}
浏览 0提问于2021-06-13得票数 1
回答已采纳
如果我将一个布尔变量初始化为0，我会得到一个不正确的解(0)。如果我把它初始化为1，我得到了正确的解(1)。
# Squaring doesn't work
#######################################################
m = GEKKO(remote=False)
b = m.Var(lb=0,ub=1,integer=True, value=0)
m.Maximize(b**2)
m.options.SOLVER = 1
m.solve(debug=0, disp=True)
返回：
Successful solution
浏览 2提问于2021-02-21得票数 2
我试图复制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",
浏览 0提问于2019-07-08得票数 4
我正在用Python设置一个新的线性优化代码。不幸的是，我没有相同的结果与Pulp，Scipy和Gekko软件包。
我尝试用不同的包实现Python中的线性优化代码。
利用GEKKO进行优化
from gekko import GEKKO
import numpy as np
import matplotlib.pyplot as plt
m = GEKKO() # create GEKKO model
x = m.Var(value=0, lb=0, ub=400000) # define new variable, initial value=
浏览 0提问于2019-04-14得票数 4
回答已采纳
1回答
我正在运行长时间的模拟，这涉及到一个ode的集成。由于运行时间长，我决定产生后援，允许从检查点继续模拟。对于ode集成，我决定选择实例。这种方法适用于可重入的求解器“dopri5”和“vode 853”，但对于“vode”和“lsoda”则失败。对于我的特殊问题，'vode‘algo似乎是最快的，所以我感兴趣的是，是否可以绕过这个问题。我希望这是可能的，因为在重新启动时也只需要一个求解器实例。
运行以下两个代码片段将演示我的问题，请查看此
initialize.py
启动集成，并通过酸洗(使用dill，也可以对lambdas进行腌制)来保存ode解决程序实例。
import scipy
浏览 1提问于2018-01-04得票数 2
在Python Gekko中，有两种方法可以使用m.options和m.solver_options设置求解器选项。哪种方法优先，什么时候应该使用其中一种方法？
例如，我想为求解器设置 (OTOL)和 (RTOL)。Gekko使用哪一个(1e-7或1e-8)？
from gekko import GEKKO
m = GEKKO() # Initialize gekko
m.options.SOLVER=1 # APOPT is an MINLP solver
m.options.OTOL = 1.0e-8
m.options.RTOL = 1.0e-8
# solver settings
浏览 38提问于2019-11-26得票数 1
我试图用python解决一个优化问题，我开始研究GEKKO，因为它解决了非线性程序。我用术语编写了一个简单的模型来检查答案，这个模型与我从GEKKO (同一个模型)得到的答案没有相同的值。
Python代码：-
from gekko import GEKKO
# Initialize Model
smplmdl = GEKKO()
# Create Variables
x = smplmdl.Array(smplmdl.Var, 3, lb = 0)
a = smplmdl.Array(smplmdl.Var, 3, lb = 0)
Constant_Val = [10, 15, 2
浏览 1提问于2022-06-15得票数 2
我想我遇到了和这个堆栈溢出相同的问题。
昨天我能够解决我的Gekko模型，今天我甚至不能解决来自Apmonitor网站的例子。解决问题需要花费很长时间，我得到了以下错误：
ImportError: No solution or server unreachable.
Show errors with m.solve(disp=True).
Try local solve with m=GEKKO(remote=False).
我尝试在本地解决它，正如我在上面提到的堆栈溢出帖子中所建议的那样：
m = GEKKO(remote=False)
但是，我得到了以下错误：
Error: &#
浏览 0提问于2020-05-07得票数 3
回答已采纳
1回答
在利用GEKKO对多级蒸发站的过程数据建立ARX模型后，我想看看这个预测模型如何与LSTM预测控制相适应，但是我很难将LSTM模型(保存为.h5)连接到GEKKO上。
所以我要问两个问题
有可能和GEKKO一起使用LSTM吗？
是否存在这样做的文献或案例？
浏览 8提问于2020-09-30得票数 2
回答已采纳
我正在使用非matlab的ODE模拟软件来重现一个用matlab中的simbiology工具箱创建的模型。
其中一个问题是重复赋值的表示。有没有可能用一种只支持ODE系统的标准Runge Kutta (或其他迭代方法)来模拟重复赋值？或者，将具有重复赋值的模型重新表达为ODE系统是不可能的？
浏览 0提问于2014-04-17得票数 0
我正在尝试求解Gekko中的DAE，其中一些组件依赖于卷积的解。
这需要一个恒定的dt，但我确信这是在选项中的某个位置。因此，我想要做的是使用一个函数来记录数组中状态变量的当前值，并返回到该值的总和。下面是我尝试使用一个简单的ODE示例：
import numpy as np
from gekko import GEKKO
import matplotlib.pyplot as plt
m = GEKKO()
class adder:
def __init__(self,):
self.i=m.Array(m.Param, 10)
s
浏览 3提问于2022-04-06得票数 2
回答已采纳
我有一个Simulink R2015b模型包含SimElectronics和SimMechanics块。这种模式是可行的。我的一位同事要求将模型导出到R2015a。但是，当我执行这个模型时，我会得到以下错误消息：
组件: Simulink类:模型错误
*“型号名称”：没有为一个或多个Simulink转换器块提供足够的输入导数，以供所选择的求解器使用。隐式求解器(ode23t、ode15s和ode14x)通常比显式求解器需要更少的输入导数，而本地求解器不需要任何输入导数。您可以通过在Simulink-PS转换器块的输入处理选项卡上选择不同的选项(例如，打开输入过滤)来提供更多的输入衍生物：.‘
浏览 2提问于2016-03-11得票数 2
回答已采纳
我试图通过修改找到的代码来解决问题。
import sympy as sp
t = sp.symbols('t', real=True)
x = sp.Function('x', real=True)
diffeq = sp.Eq(x(t).diff(t), sp.sqrt(1 - x(t)**2))
res = sp.dsolve(diffeq, ics={x(0): 0})
math.se文章建议，分段函数满足方程，但也许我的代码没有添加所有必要的假设。我得到以下回溯：
Traceback (most recent call last):
File
浏览 11提问于2022-10-22得票数 0
回答已采纳
1回答
我正在使用Python解决一个大规模的优化问题。模型中的一些变量在优化过程中不需要改变。因此，我想知道在变量(m.Var)上使用m.Var方法(我希望是常量)与用m.Const声明的常量替换变量之间有什么区别。
我使用IMODE =3(默认)，因为我不做任何需要模拟或估计、导数、时间等的事情。
浏览 5提问于2022-08-09得票数 1
回答已采纳
1回答
我刚刚开始探索Gekko，并尝试模拟Lorenz ODE系统。不幸的是，我得到一个错误(“找不到解决方案”)，对于一个简单的std情况下，运行良好使用fine。
如果我只集成time=0.5而不是1.0，这个问题就能很好地解决
from gekko import GEKKO
import numpy as np
import matplotlib.pyplot as plt
m = GEKKO()
m.time = np.arange(0.0, 1.0, 0.01)
sigma = 10.; rho = 28.0; beta = 8./3.
x = m.Var(value=10);
浏览 10提问于2022-05-03得票数 2
回答已采纳
大家好，Gekko开发人员，
我想为我的热能存储系统建立以下方程模型，
它们既是线性的，又是微分的。一些变量也是时间相关的，作为外部参数出现。从Gekko的例子中，看起来我只能添加这种格式的方程，比如mdl.Equation(T_a.dt() == mt*c_p_w*(T_b-T_a)) / d_w*c_p_w，但是对于不涉及微分的方程，我想添加为普通方程，也可以添加mdl.Equation(m[t] == mc - m[t])。
附注:我如何添加这两种样式，以及哪种模式和求解器将允许这样做。？
谢谢,
诚挚的问候,
拉胡尔
浏览 4提问于2020-11-18得票数 3
我试着在Gekko中解一个有逻辑条件的微分方程。我知道Gekko不喜欢这些东西，但是我假设简单的if3()函数，为了在两个给定的表达式(R1，R2)之间切换，会得到一些解决方案。下面是一个失败的简单示例代码--解决方案未找到。
import numpy as np
from gekko import GEKKO
import matplotlib.pyplot as plt
A=1
B=1e-5
C=2
D=0.01
G=1
y0=120
m = GEKKO() # create GEKKO model
nt = 101
m.time = np.linspace(0,100,nt)
浏览 8提问于2021-03-17得票数 1
1回答
我一直在和Gekko(https://gekko.readthedocs.io/en/latest/)玩。为了测试它，我从这里的https://julia.quantecon.org/continuous_time/seir_model.html实现了Covid模型。但结果只是平面图。它看起来一点也不像链接的结果。有没有人看到我做错了什么？ import numpy as np
from gekko import GEKKO
import matplotlib.pyplot as plt
m = GEKKO(remote=False) # create GEKKO model
#
浏览 35提问于2021-02-08得票数 1
回答已采纳
我已经成功地使用GEKKO解决了各种问题，但有一件事我还不能确定，那就是目标函数或约束中的非线性项的限制是什么。我发现一些非线性的术语是可以容忍的，例如 m.Obj(sum(x * values) / sum(x)) 然而，似乎工作得很好。 m.Obj(sum(x * values + x ** 2 * other_values) / sum(x)) 似乎是无法解决的。有没有关于GEKKO处理某些形式的非线性的限制的文档，可以帮助我重新格式化我的问题？
浏览 23提问于2021-01-07得票数 4
我正在尝试使用GEKKO优化我的有限元分析桁架解决方案。所以我定义了所有的有限元解决方案作为一个函数，并尝试使力量作为variables.SO的目标是最小化力与一个约束(总和的所有力量=-100000).But后，我运行我的代码，我得到一个错误，我无法修复它谷歌也没有帮助(我是新手).I将不胜感激，如果有人能告诉我的问题is.Thanks提前。 try: from gekko import GEKKO除非：# pip install gekko import pip pip.main('install'，' GEKKO ') from gekko import
浏览 5提问于2020-01-23得票数 3
我用Matlab求解器解决了一个PDE问题。初始条件是ODE的解，我在一个不同的m.file中解决了这个问题。现在，我有一个大小为NxM的矩阵形式的ODE解。我怎么能用它做我的IC呢？这可能吗？当我使用for loop时，pdepe只将最后一次迭代作为初始条件。任何帮助都是非常感谢的。
浏览 0提问于2015-01-14得票数 0
回答已采纳
我有一个从S-函数调用热力学性质计算(CoolProp)库的模型.如果它的输入没有意义，这个库会产生一个错误并突然停止仿真，这有时发生在Simulink干扰状态时，特别是当模型中存在一个Step块时，即使它的输出是浮动的！我认为这是一种预测系统状态变化的尝试，我注意到减小步长会减少问题，但我无法在10 us运行20,000 s的模拟(是的，微秒！)时间步骤，当它已经花了15分钟与5s步长！
有没有办法捕捉这个错误(我想是通过调用error()引起的)，并告诉Simulink“嗯，您使用的步长太大了”？
编辑：我使用的是ode45求解器，但是除了模拟时间急剧增加之外，即使是ode23tb或ode
浏览 3提问于2015-03-26得票数 0
3回答
我对“物理系统建模”这门学科很陌生。我阅读了一些基本文献，并在Modelica和Simulink/Simscape中做了一些教程。我想问你，我是否正确地理解了以下内容：
符号操作是将一个微分-代数方程组(物理模型: DAE)转化为一个微分方程组(ODE)，可由标准求解者(Runge，Kutta，BDF，.)求解的过程。
也有可以直接解决DAE问题的求解器。但是Modelica (openModelica，Dymola)和Simscape将系统转换为ODE (为什么这种方法比直接DAE求解器更好？)
“平面Modelica代码”是转换的结果(= ODE)。
非常感谢你的回答。
浏览 4提问于2016-05-20得票数 7
我使用数值积分一个常微分方程组(ODEs)。请看我的笔记本电脑。
与API文档中的示例代码一样，函数ode_fn(t, y, theta)定义了待解决的ODE系统。我能够接受ode_fn wrt theta的梯度，并将ODEs与tfp.math.ode.BDF集成。
但是，当我尝试接受ODE解决方案的梯度结果wrt theta时，我会得到以下错误。当我用一组简单的ODEs替换ode_fn时，代码运行时不会出现任何问题。是否应该调整解决程序设置以避免此错误？
InvalidArgumentError Traceback (most recent call
浏览 1提问于2021-11-04得票数 0
回答已采纳
1回答
有时，一个求解器比另一个求解器能更好地解决我的Python Gekko应用程序。很难预测哪个求解器会表现得最好。因为Python Gekko支持使用m.GEKKO(server='http://{address}')在不同服务器上进行本地或远程求解，所以是否可以创建一个并行Gekko应用程序，该应用程序将在任意数量的计算机(包括本地计算机)上同时尝试所有求解器，然后在第一个进程成功返回时终止其他进程？我一直在研究Python的多线程和并行程序包。有没有能和Gekko一起很好地进行并行求解的工具？这是一个连续的原型：
from gekko import GEKKO
m = GEK
浏览 11提问于2019-11-16得票数 2
我在IPOPT (v3.12.10)中遇到本地不可行错误。
Converged to a point of local infeasibility. Problem may be infeasible.
我还使用m.options.SOLVER=1对APOPT求解器进行了尝试，并得到了类似的错误。
No feasible solution
这个问题有两个方程和三个变量的多个可行解。这个问题的解决方案是平面与球体相交的地方。
from gekko import GEKKO
m = GEKKO()
x=m.Var(); y=m.Var(); z=m.Var()
m.Equations([x**2+
浏览 31提问于2019-10-27得票数 1
我有一个gekko模型，它没有给出我期望的解决方案。为了调试目的，最好以可读的形式打印模型(就像从CPLEX打印LP格式一样)。
这个是可能的吗？
我的模型在Matlab中得到了很好的解决。
浏览 2提问于2022-04-06得票数 2
回答已采纳
1回答
我是gekko和APM的新手，我试着解决一个问题，但是在二元变量中的一个新方程似乎不可行，这个方程应该是可行的。以下是简化的APM模型：
Model
Variables
int_v1 = 0, >= 0
int_v2 = 0, <= 1, >= 0
v3 = 0
v4 = 0
v5 = 0
v6 = 0
End Variables
Equations
(0+int_v1)>=100
v3=((3.15)*(int_v1))
v4>=((int_v2)*(300))
v5=(0+((in
浏览 2提问于2020-04-24得票数 4
回答已采纳
我正在建立一个GEKKO模型，在那里我试图最小化一个函数。这是通过2D-fourier展开函数并显式计算所述泛函的Riemann和来完成的。在增加了fourier系数和网格点的数目之后，我遇到了一个问题：
我的代码看起来像这样(我忽略了一些求和和复杂部分)：
m = GEKKO(remote=false)
ur_x_t =[m.Var(value=param[i][2],lb=-0.01,ub=0.01) for i in range(nG)]
ur_y_t =[m.Var(value=param[i][3],lb=-0.01,ub=0.01) for i in range(nG)]
fo
浏览 1提问于2021-12-07得票数 2
1回答
我试图给gekko的一个变量赋值无穷大。我尝试过numpy的无限值和python自己的无限值，但由于gekko的识别问题，它仍然无法工作。
这个想法的主要目的是强迫变量严格等于0，至少在求解器的第一次迭代中是这样的。
from gekko import GEKKO
from numpy import Inf
model=GEKKO()
R=model.FV(value=Inf)
T=model.Array(model.Var,2)
Q=model.FV()
model.Equation(Q==(T[1]-T[0])/R)
model.solve()
我所犯的错误是：
Exception: @e
浏览 2提问于2022-08-24得票数 1
回答已采纳
假设我有一个pytorch模型，描述了某个多维系统基于自身状态x和外部执行器u的演化，所以x_(t+1) = f(x_t，u_t)，f是pytorch的人工神经网络。
现在我想解决一个动态优化问题，找到一个最优的u值序列，以最小化依赖于x的目标。类似于：
所有时间步长φ的最小和(X_t)
s.t.：x_(t+1) = f(x_t，u_t)
此外，我还对x中的一些变量有一些上界和下界。
有没有一种简单的方法可以使用像pyomo或gekko这样的动态优化工具箱做到这一点？
我已经编写了一些代码，将前馈神经网络转换为numpy函数，然后可以将其作为约束传递给pyomo。这种方法的问题是，每当神经网络
浏览 28提问于2021-05-26得票数 2
回答已采纳
1回答
我使用deSolve软件包来求解描述捕食-食饵动力学的微分方程。作为一个例子，下面是一个简单的L捕食者-食饵模型.我希望模型中的一些参数随时间而变化。我可以在event函数中使用ode参数来改变状态变量(例如猎物密度)，没有问题。
但我不能使用event参数来更改参数。
下面是简单的L模型，没有添加任何事件(工作正常)
# Lotka-Volterra Predator-Prey model from ?deSolve::ode
# define model
LVmod <- function(Time, State, Pars) {
with(as.list(c(State, P
浏览 6提问于2022-11-10得票数 1
回答已采纳
我用Gekko编写了一个模型，但遇到了一个错误。所以，我开始逐行检查，发现了一个我不能理解的问题。可以执行下面的代码(显示了部分模型代码，输入文件)，并且结果也是正确的。但是，当我取消对第31行的注释时，出现了一个错误。会不会是因为我按顺序使用了逻辑条件(ab3、max3)？有没有办法解决这个问题？
import numpy as np
from gekko import GEKKO
# Read weather data
weatherData = np.load("Atlanta_TMY3_climate.npy")
# create GEKKO model
m = G
浏览 21提问于2021-02-16得票数 1
1回答
这个问题在一定程度上侧重于经济优化，在某种程度上关注python的实现，但也许社区中的一些人能够提供帮助。我试图在Python的GEKKO平台上实现一个标准的连续时间宏观经济储蓄模型，但一直未能解决这个问题。我已经以GEKKO的中提供的经济例子为例，并且适应了基本的储蓄决策模型，但是事情还没有完全解决。该模型最大限度地利用消费效用之和，其中消费+投资=产出。例如，最大积分(U(y-i))。输出y= k^ALPHA投资= dk/dt+delta*k。
有人能说出为什么我的代码不能被破解吗？这个平台有能力解决这样一个模型吗？我没有见过很多经济学家使用这个平台来解决模型的例子，但我不确定这是因为这个
浏览 2提问于2020-03-23得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
云直播
对象存储活动推荐
腾讯云开发者
