1回答
我有一个非线性常微分方程系统，它有3个参数作为输入:µ、target和omega。我尝试通过最小化观测值与ODE系统(x5)输出之间的最小二乘误差来估计这些参数。我尝试了不同的目标函数以及不同的IMODES，但解决方案从未收敛。此外，由于我总是找不到解决方案的错误，我看不到最接近的最佳参数是什么。你对我的优化问题有什么建议吗？谢谢!这是我的yData 我的代码是： import n
浏览 26提问于2021-01-18得票数 3
1回答
我正面临着一个问题，在这个问题中，我必须解决一个方程系统：\部分B\δy= h(A) k(B)你还有其他更有效的想法吗？或者可能有更优雅的方式，比如组合自变量x和y，积分一个向量函数，然后撤消更改
浏览 3提问于2014-10-03得票数 0
1回答
我正在使用GEKKO来估计微分方程的参数，并且我将其中一个变量限制在0到1之间。但是，当我求解常微分方程时，我得到的值超出了这个变量的界限，所以我想知道是否有人知道GEKKO是如何找到解的，因为这可能会帮助我解决问题。下面是我用来拟合数据的代码。然而，之后，我尝试用scipy.integrate.solve_ivp求解这个常微分方程，用我得到的u的</em
浏览 7提问于2020-11-04得票数 2
引言x''=-D(y) * x' * sqrt(x'^2 + y'^2)
y''=-D(y) * y' * sqrt(x'^2 + y'^2) + g(y)上述两
浏览 4提问于2017-01-24得票数 1
我正在尝试最小化函数与数据点在不同时间点上的差异。因此，目标函数是模型(我的函数)和数据点在不同时间的差值的平方和。 我的模型有解析的一阶和二阶导数。如何向Gekko Python提供这些派生工具？
浏览 22提问于2021-02-04得票数 2
回答已采纳
1回答
我想对模型中的4个参数执行动态参数估计，如这个链接的示例4所示：。我与我的应用程序的区别是，有4个微分方程，其中包含4个参数。其中只有一个方程有实验数据。我已经调整了示例中给出的用于此应用程序的代码，但它对每个参数返回0。from gekko import GEKKO
import matplot
浏览 3提问于2022-07-24得票数 2
2回答
我有9个方程，随时间变化的系数为gfunction dy =tarak(t,y)g = 0.1*exp(-((t-200)/90).^2);其中系数G = 3.16和g = 0.1*exp(-((t-200)/90).^2)是随时间变化的系数和时间我得到了错误的输出y(1)，y(2)的负值。有没有人可以试着用ode45解
浏览 2提问于2015-04-01得票数 2
1回答
我对GEKKO不熟悉，所以我的问题可能很愚蠢。我想估计一下我的模型的一些参数，这是一个SIR模型。我看过很多次关于约翰·赫登格伦教授的文件和视频，但是盖科对我来说还是很难理解。这是我的密码：tab_data = [275.5,317,457.34,646,888.67,1236.67,1619.34,2077.34]
total_active_data59999385.33333334,59999158
浏览 0提问于2020-11-13得票数 2
最近，我一直在对一些动力学系统进行模拟，其中所有的动力学变量都是相互依赖的。因此，为了模拟动态，我在小的时间步长dt<<1上执行了循环，并在每次迭代中改变了数量。我得到了很好的结果，但由于迭代过程缓慢，模拟可能需要相当长的时间。一般来说，我听说应该避免像我使用过的for循环，因为它们会大大减慢模拟速度。然而，另一方面，我不知道如何在没有迭代的情况下以小的时间步长进行模拟。因此，我问你:对于一个动态系统，其中每个量都必须在超
浏览 1提问于2014-06-14得票数 1
1回答
我试图使GEKKO最大化的功能“利润”，但我无法使它工作。我总是以下面的错误告终。我需要GEKKO解决选择最佳的价值(Weight1，Weight2，Weight3.因此，函数“Weight8”返回的“结果”是尽可能高的。每个权重可以是0 (0%)到15 (15%)之间的任意值，所有权重的总和必须是100 (100%)。
所有的DataFrames都需要计算每个挑选(Pick1，Pick2，Pick3)的最终利润。最后，“结果”是每个
浏览 11提问于2022-05-23得票数 2
大家好，Gekko开发人员，
它们既是线性的，又是微分的。一些变量也是时间相关的，作为外部参数出现。从Gekko的例子中，看起来我只能添加这种格式的方程，比如mdl.Equation(T_a.dt() == mt*c_p_w*(T_b-T_a)) / d_w*c_p_w，但是对于不涉及微分的方程，我想添
浏览 4提问于2020-11-18得票数 3
我用四阶龙格-库塔法求解带有NumPy阵列的duffing振子的微分方程，但是我收到了一个错误。RuntimeWarning: overflow encountered in double_scalars 有谁知道可能导致溢出错误的错误来源，以及我可以尝试解决它的方法吗？n-1] + h/6 * (k1 + 2* k2 + 2 * k3+ k4)
xn[n] = xn[n-1] + h/6 * (l1 + 2* l2 + 2* l3 + l4) 这是我的</em
浏览 41提问于2021-04-27得票数 0
我的模型是常微分方程系统，系数是具有参数的函数。因此，如果你知道如何做到这一点，就使用数据。
浏览 4提问于2021-05-17得票数 0
我想用CSTR的稳态浓度数据拟合反应常数(k0和EoverR)与反应器温度的关系。我只是用一个简单的线性方程来生成操作数据来拟合。(Ca_data =-1.5*T_反应堆/100+ 4.2)import matplotlib.pyplot as plt
from gekko import <e
浏览 3提问于2019-11-01得票数 3
回答已采纳
我现在正在做一个项目，我已经对多体机械系统的运动进行了建模。这是一个复杂的、非光滑的系统，系统的运动并不总是由同一组方程描述，这意味着如果我用数值方法解决这个问题，我必须有时停止积分(并用新的初始条件和方程重新启动它)。尽管如此，我还是写了一些代码，现在我只是停留在如何将所有这些放在一起:我希望常微分方程求解器(例如ode45)求解我的方程，然后
浏览 9提问于2019-02-21得票数 0
回答已采纳
使用Gekko来拟合数据的数值ODE解。import numpy as npCa_data = [2.0081, 1.5512, 1.1903, 0.7160, 0.2562,
浏览 3提问于2020-03-09得票数 4
回答已采纳
1回答
如何使用Sympy解矩阵微分方程？
我有一个方程y'(t) = A*y(t) + B，其中A是3x3矩阵，y(t)是1x3向量，B是1x3向量。更具体地说，我正在研究一个使用微分方程在3D中移动点的计算机图形学问题。我在3D空间y(t)中有一个点，一个3x3旋转矩阵和一个平移向量。y'=Ay+B方程是我正在研究的实际问题的简化，但我甚至无法使用Sympy来求解y'=Ay+B。我在Sympy中寻找闭合形式的解(而不是数值解)
浏览 3提问于2018-08-25得票数 4
我正在学习如何使用GEKKO进行基于实验室批式反应器数据的动力学参数估计，该数据实质上包括A、C和P三种物种的浓度分布，就我的问题而言，我使用的是我以前在中介绍的一个模型，该模型与从单个数据集估计参数有关我的最终目标是能够使用多个实验运行进行参数估计，利用在不同温度、物种浓度等情况下可能收集到的数据。由于单个间歇反应器实验的独立性，每个数据集都具有在不同时间点采集的</em
浏览 3提问于2021-01-13得票数 2
回答已采纳
2回答
我正在尝试切换matlab中的一些程序。我需要解一个刚性常微分方程，它的输入都是矩阵。在matlab中，我使用
[ttT,uT] = ode23s('SST',t,fT);
浏览 0提问于2012-01-05得票数 4
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
