使用ODE实现动态模型中的共享参数
ODE是Ordinary Differential Equation（常微分方程）的缩写，它用于描述动态系统中的变化规律。在动态模型中，共享参数是指多个子模型之间共同使用的参数。
使用ODE实现动态模型中的共享参数的过程如下：
- 定义动态模型：首先，根据具体的问题和需求，定义动态模型的数学方程。这些方程描述了系统中各个变量之间的关系和变化规律。
- 确定共享参数：在动态模型中，确定需要共享的参数。这些参数可能是影响系统行为的重要因素，需要在多个子模型中共同使用。
- 设定初始条件：为动态模型中的变量设定初始条件。初始条件是模型开始运行时各个变量的初始取值。
- 选择ODE求解方法：根据模型的特点和求解需求，选择合适的ODE求解方法。常用的求解方法包括欧拉法、龙格-库塔法等。
- 实现共享参数：在求解ODE的过程中，将共享参数作为输入，确保多个子模型共同使用这些参数。可以通过函数参数传递、全局变量等方式实现参数共享。
- 求解ODE：使用选择的ODE求解方法，对动态模型进行求解。根据设定的初始条件和共享参数，求解模型在一定时间范围内的变化过程。
- 分析和优化：根据求解结果，分析模型的行为和性能。如果需要，可以对模型进行优化，调整共享参数的取值，以达到更好的效果。
ODE在动态模型中的应用非常广泛，包括物理模拟、生物学建模、经济学模型等领域。通过ODE求解动态模型，可以更好地理解和预测系统的行为，为决策和优化提供支持。
腾讯云提供了一系列与云计算相关的产品，包括云服务器、云数据库、云存储等。这些产品可以帮助用户快速搭建和部署云计算环境，提供稳定可靠的基础设施支持。具体产品信息和介绍可以参考腾讯云官方网站：https://cloud.tencent.com/
相关·内容
我想在数值格式的每一步提取状态变量的值。
deSolve R包中的ode()函数，使用实现的ODE求解器之一对一个常微分方程组进行数值求解。为了做到这一点，它使用了一个动态调整的集成步骤，该步骤基于每个集成末尾的本地截断错误的值。用户基本上指定了所需的输出时间，它定义了时间步长网格，它可以有一个等于或小于输出时间的步骤。
如果我们以Lotka Volterra捕食者-猎物模型(含逻辑食饵)为例：
LVmod <- function(Time, State, Pars) {
with(as.list(c(State, Pars)), {
Ingestion <- r
浏览 2提问于2019-03-31得票数 2
回答已采纳
我正在使用非matlab的ODE模拟软件来重现一个用matlab中的simbiology工具箱创建的模型。
其中一个问题是重复赋值的表示。有没有可能用一种只支持ODE系统的标准Runge Kutta (或其他迭代方法)来模拟重复赋值？或者，将具有重复赋值的模型重新表达为ODE系统是不可能的？
浏览 0提问于2014-04-17得票数 0
我认为这个问题有点长，所以我认为最好先考虑一下简化的版本：
有两个类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
回答已采纳
1回答
我是PyMC ()的初学者，正在尝试构建一个具有动态组件的模型，本质上是在每次调用模型时求解一个小的常微分方程(ODE)系统。
我尽我所能搜索了谷歌和(废弃的) PyMC邮件列表，得出了模型结构(pk_model.py)和调用者(pk_fit.py)：
运行pk_fit起初似乎是可行的，但随后会转储许多与ODE解算器相关的错误消息。ODEs和求解器的独立版本(未与PyMC代码集成)可以很好地工作。
在不同的打印语句中，看起来代码一开始是工作的，但是在某个时刻，未知的参数参数(kcp，kpc，ke)从标量变为numpy.ndarrays。这似乎是问题的一部分。
作为PyMC的新手，我认为我在做
浏览 3提问于2013-07-28得票数 0
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
回答已采纳
我对共同模拟这个话题很陌生。我熟悉以下定义(基于Trcka“建立ANDHVAC/R系统模拟的COMPARISON比较")：
准动态耦合，也称为松散耦合，或称orping Pong耦合，其中分布式模型按顺序运行，其中一个模型使用耦合模型的已知输出值(基于前面时间步骤的值)。
全动态耦合，也称强耦合，大葱耦合，其中分布式模型在每一时间步长内迭代，直到误差估计达到预定的容限。
我的问题是:FMI/协同模拟是一种松散耦合的方法吗？什么是FMI/模型交换？据我理解，这不是一个强耦合的方法。我的理解是，在模型交换中，导入FMU的工具是收集所有的ODE和代数方程，工具用一个求解器来求
浏览 4提问于2017-05-04得票数 4
我试图使用solve_ivp()函数来解决带有状态相关args的ODE问题。其思想是在给定轨道当前状态的情况下，每次更新ODE。
给定ODE模型如下：
def model(t, y, arg1, arg2):
'''
Some Dynamics model that includes arg1 and arg 2
'''
return dy_dt
arg1和arg2是使用其他函数计算的：
def gen_args(y):
'''
Calculate arg1 and arg2
浏览 1提问于2019-04-26得票数 1
回答已采纳
2回答
我有一个通用函数，它定义了我计划使用集成的一种ODE形式，例如：
def my_ode(K, tau, y, u):
return K*u/tau - y/tau # dydt
我的代码中有几个对象，它们都具有在my_ode中定义的表单的动态性，但是具有唯一的参数K和tau。我希望能够向my_ode传递一个惟一的句柄，在初始化对象时已经设置了这些参数，所以当我更新对象时，我所要做的就是在一些模拟时间( t )中做一些类似soln = odeint(my_ode, t, y, u)的事情。
例如，如果我定义了一个类：
class MyThing:
def __init__(self,
浏览 3提问于2016-01-14得票数 7
回答已采纳
我有两个类，其中一个(类ode)有一个私有成员，它是一个函数(指向函数的指针)：
class ode
{
private:
double T;
double (*sfn)(double,double,double*);//I mean this
int Ne;
public:
ode(double t,double (*f)(double,double,double*),int N)
{
Ne=N;
T=t;
sfn=f;
} //.....
}
现在我想构造一个ode对象，在它的构造中，函数(
浏览 0提问于2015-09-19得票数 0
我正在为硕士学生编写一个练习，我希望他们能够键入一个ODE系统，并使用deSolve进行自动评估。
我现在的问题是根据用户输入定义在R中动态评估的函数。假设我想从文本文件ODE.txt中读取，其中包含的ODE系统的语法与我使用的直接写入它的语法相同：
##function (differential form)
ODEfun <- function(time=seq(1:times_range), state=state, parms=parms) {
with(as.list(c(state, parms)), {
source(&#
浏览 11提问于2022-01-28得票数 1
1回答
双主机SI模型与R中的变参数
、、、、
我目前正在研究一种两宿主流行病学SI模型。也就是说，一个没有恢复室的隔间模型。
我对R还是比较陌生的，但是在大部分使用MATLAB之后，我正在开发一个不错的理解。然而，在寻找有用的资源方面，我遇到的问题是如何改变两个不同的输入参数，这样我就可以对它们进行检查，甚至可以对这些变量进行三维绘图，或者绘制相位图，以确定人口是否死亡。
所以，更具体的说，我想在0到1之间改变mu，0到1之间变alpha，我可以“即插即用”，但是我希望能够显示一个更动态的结果，并认为在我的车轮房里有一个工具是很方便的。
总之，这是我到目前为止的代码：
# Here we will load the required pa
浏览 8提问于2022-09-19得票数 1
在我的模型中，每个智能体在每个节拍上求解一个ODE系统。我采用了欧拉法(类似于NetLogo中的系统动力学模型)来求解这些一阶常微分方程。然而，为了获得稳定的解决方案，我被迫使用非常小的时间步长(dt)，这意味着使用这种方法进行模拟的速度非常慢。我很好奇，有没有人有建议可以更快地解决这些常微分方程？我正在考虑实现Runge-Kutta (具有更大的时间步长？)正如这里所做的()。我也会考虑使用R扩展和在R中使用ODE求解器，但同样，ODE是由每个代理求解的，所以我不知道这是不是一种有效的方法。
我希望有人对这些方法的性能有所了解，并能提供一些建议。如果没有，我会尝试分享我的发现。
浏览 0提问于2015-04-10得票数 1
2回答
我需要用Simulink来解决这段代码，我不知道怎么做。我只知道如何使用ODE求解器来实现它。
y'' - y' - 2y = e^(3x)
y(0)=1, y'(0)=2.
我重写了方程，得到了一首颂歌：
y' = f(x,y)
y(x0) = y0
y'1 = y2
y2= e^(3*x) + y' + 2y
使用ODE解算器。
如果有人能用Simulink模型帮我解决这个问题，我会很感激的。
我知道如何在Matlab中使用像ode23和ode23s那样的ODE求解器，但是我不知道如何使用Simulink模型来解决这个问题。
浏览 0提问于2015-06-01得票数 2
回答已采纳
我用Dymola建立了一个简单的模型，在不同的阶段使用了不同的方程，但令我惊讶的是Dymola无法处理这个简单的模型。经过指数降阶后，系统似乎是奇异的。
我的问题是：
在使用Modelica时是一个常见的问题吗？如果是这样的话，我应该如何修改这个模型，我知道当我用der(x)=0替换x=100方程时，模型工作得很好，但是我想知道对于更复杂的情况是否有一个通用的规则。如果能提供详细的解释，我将不胜感激。
下面是模型的代码：
model ErrorWhenUsingIf
Real x(start=100);
equation
if time<=0.5
浏览 5提问于2020-12-08得票数 3
回答已采纳
当我不断地得到形状误差时，我试图为我的最少的模型和方程绘制数据。我只是想猜测符合数据的参数，但如果看不见图表，我就不能这么做。
以下是我到目前为止所拥有的：
from scipy.integrate import odeint
Import numpy as np
Import matplotlib.pyplot as plt
#data
didata1 = np.loadtxt("diauxic1.txt")
time = didata1[:,0]
pop2 = didata1[:,1]
# model equations
def diauxic_ode(x,t,para
浏览 3提问于2020-01-09得票数 0
回答已采纳
我有以下颂歌：
x_dot = 3*x.^0.5-2*x.^1.5 % (Equation 1)
我正在使用ode45来解决它。我的解是dim(Kx1)的向量(通常k= 41，由tspan给出)。
另一方面，我从(1)建立了一个近似于模型的模型，但是为了比较第二个模型的精度，我想用ode45来求解它(求解第二个ODE)。我的问题是，第二个ode是离散的：
x_dot = f(x) % (Equation 2)
f是离散的，不是像(1)中那样的连续函数。我为f提供的值是：
0.5644
0.6473
0.7258
0.7999
0.8697
0.9353
0.9967
1.0540
1.107
浏览 3提问于2013-11-01得票数 1
回答已采纳
2回答
Modelica是否将DAE转换为状态空间形式(ODE系统)？
如果我想要为一个系统建模:例如，一个带有动态管道和水槽的简单模型，来源:中等简单的不可压缩液态水。状态变量是什么(压力，温度，..?)？
谢谢你的帮助!
浏览 6提问于2015-12-16得票数 1
2回答
我正在致力于将一个模型从MATLAB翻译成Python。该模型的关键在于MATLAB的ode15s。在MATLAB执行中，ode15s具有标准选项：
options = odeset()
[t P] = ode15s(@MODELfun, tspan, y0, options, params)
作为参考，和MODELfun一样，y0是一个向量(大小为98)。
我的Python尝试的等价物如下所示：
ode15s = scipy.integrate.ode(Model.fun)
ode15s.set_integrator('vode', method = 'bdf'
浏览 3提问于2014-06-27得票数 3
我正试图用求解器ode23s (完全相同的算法在Matlab中完美地运行)来解一些八度的乐谱，但是我收到了一个“x”未定义的错误，而当我使用求解器ode23或ode15s时，这并没有发生，尽管这是我要做的唯一改变。
下面可以看到上述代码和错误：
model = [Fg1 Fg2 ui mi ... gi]; %model is a vector of (1,77)
options = odeset('RelTol',1e-4,'AbsTol',5e-6*ones(16,1));
[T,Y] = ode23s(@odesys_test,0
浏览 1提问于2022-01-16得票数 0
2回答
我正试图解决以下几个问题：
function [eta, sol] = compressible_similarity_wo
global Gamm Ma Pr omega;
Gamm = 1.4;
Ma = 2;
Pr = 0.7;
omega=0.76;
global eta_max_ode;
eta_max_ode = 20;
opt = optimset('Display','off','TolFun',1E-20);
F = fsolve(@(F) eva
浏览 4提问于2020-02-06得票数 1
回答已采纳
首先，我对BPEL非常陌生。我有一个托管在Tomcat7.0中的webservice和一个部署在Apache ODE中的BPEL。我已经使用Eclipse JUNO完成了所有这些工作。现在，当我无法从eclipse调用web服务时，因为Apache ODE无法在Tomcat实例运行时启动。此外，我不能在Apache ODE中托管web服务，也不能在Tomcat中托管BPEL流程。现在如何调用Web服务？
浏览 6提问于2012-10-03得票数 1
回答已采纳
我试图在每个空间网格单元上应用一个常微分方程组(ODE)。因此，每个景观单元都有一个相关的ODE模型。易感蚊子(Sv)、暴露蚊子(Se)和感染蚊子(St)的数量在每一时间步骤上都被更新，ODE模型与一个基于离散agent的动物运动模型相耦合。下面是一个运行ODE模型的示例：
library(deSolve)
mod1 <- function(out_tab, time_step, var){
Sv <- out_tab[time_step,c("Sv")]
Ev <- out_tab[time_step,c("Ev")]
Iv &l
浏览 0提问于2018-04-17得票数 6
回答已采纳
我有一个长度为N的向量x，我想用它的值来解微分方程：dy/dt = x - 4*y。对于每个步骤，我希望ode求解器函数使用向量的一个值，然后在下一步使用矩阵的下一个值。
为此，我尝试将向量声明为全局变量，并在ode求解器中这样使用它：
global x
tspan = 0:0.01:10;
[t,filter_coef] = ode45(@ode_filter,tspan,0);
这样的解题功能：
function dx = ode_filter(t,fil)
global x
dx = x - 4*fil(1);
end
产生跟随错误
ODE_FILTER retur
浏览 0提问于2019-05-19得票数 0
回答已采纳
2回答
是否有任何免费的开源硬集成器Java库？
我目前正在使用来自库的Dormand 8(5,3)集成器。我正在为我的模型实现一个DAE近似，这样我就可以使用这个积分器，但是我想尝试一个在直线ODE系统上的刚性求解器来进行比较。
我到处找过，但似乎找不到任何直接的答案。类似于MATLAB中的ode15s求解器吗？
浏览 8提问于2015-05-11得票数 3
回答已采纳
我在用MVVM和Prism和Silverlight。我对一种模式有多种不同的看法。在我编写更多视图时，他们的ViewModels似乎重复了许多与处理这个模型相关的通用代码。与其在所有VM中重复相同的通用代码，我更倾向于将其推回模型中(这可能会混淆过多的关注点)。或者进入某个普通的ViewModel基类？或者我的VM需要在它们和模型之间建立第二层的“共享VM”？这个单一的共享实例，二级VM将合并多个常规VM共享的行为和状态。
对这些问题和可能采取的办法有何评论？
谢谢各位的评论。我可能应该告诉您更多关于具体的“共享”VM代码的信息。
我可以看到在VM基类中放置了一些未来的代码，但是我所看到的特
浏览 4提问于2009-11-15得票数 1
我正在尝试编写一个ODE求解器(类似于MATLAB中的ODE45 )。我希望这个求解器是“独立的”，即我可以传递给它我想要它求解的颂歌的名称，然后它运行那个特定微分方程的代码。
到目前为止，我有我的主类(程序的主要部分)和一个名为ODEsolver的类( ODE求解器所在的位置)。我希望能够从我的主类调用ODEsolver，并让ODE查看由字符串参数引用的方法(我想要求解的ODE的名称)
我设想的方式是：
public class Main {
double y[];
double x0,xf,y0;
x0 = 0;
xf = 10;
y0 = 1;
ODEsolve
浏览 1提问于2011-10-08得票数 0
回答已采纳
1回答
我有一个衬套的数学(参数)模型。模型的输入是位移(正弦波)，以力作为输出(在相应的时间步长)。
模型参数采用遗传算法辨识，力方程采用ode15s求解。
当假设一组参数(不使用GA)时，求解器给出输出，但是当从GA调用时，我得到以下错误n no.of时间：
警告:矩阵是奇异的，接近奇异或严重缩放。结果可能不准确。RCOND = NaN.
最后
警告: t=1.570796e+000失败了。在t时刻，如果不将步长减小到最小允许值(3.552714e-015)以下，则无法满足积分公差。
我有点困惑，因为代码可以很好地工作而不使用GA，而不使用GA。我真的很感谢你的投入。
提前谢谢。
浏览 1提问于2012-06-30得票数 0
1回答
我试图用以下方式绘制范德波尔方程的三维解：
t作为x轴
mu为y轴；
y(t)为z轴。
我正在用来解这个方程。
我的2D代码是：
tspan = [0, 6000];
y0 = [1; 1];
Mu = 1000;
ode = @(t,y) vanderpoldemo(t,y,Mu);
[t,y] = ode23s(ode, tspan, y0);
plot(t,y(:,1))
xlabel('t')
ylabel('y(t)')
我想要这样的3D情节：
tspan = [0, 6000];
y0 = [1; 1];
Mu = 1:1000
浏览 2提问于2017-05-06得票数 0
我正试图用JiTCODE来解决一个动态的食物网。网络的一个方面是，接受阈值的人群被设置为零。所以我得到了一个不可微的方程。有办法在JiTCODE中实现这一点吗？另一个类似的问题是网络的Heaviside依赖性。
示例代码：
import numpy as np
from jitcode import jitcode, y, t
def f():
for i in range(N):
if i <5:
#if y(N-1) > y(N-2): #Heavyside, how to make the if-statement
浏览 4提问于2021-08-07得票数 0
我正在使用FME包来将捕食者模型与我的数据相匹配。目前，我只是在学习和测试代码，还没有适应任何东西，遵循这个例子：
当执行全局敏感度分析时，我得到一个错误，并且我不知道它来自哪里。
当我在R会话中复制粘贴上面的示例时，它可以工作。但是我不能让我自己的代码工作。我基于FME上的另一个文档编写函数，因此语法略有不同。
# model function
RMmodel <- function(t, y, parms) {
#y <- c(N = 30, P = 3)
derivs <- function(t, y, parms) {
with(as.list(c
浏览 0提问于2019-05-21得票数 0
1回答
我目前正在尝试寻找一种机器学习算法，可以预测数学模型(MM)中使用的大约5- 15个参数。MM有4个不同的常微分方程(ODE)，并且将添加更多的常微分方程，因此需要更多的参数。大多数参数是可以测量的，但其他参数需要猜测。我们知道所有的15个参数，但我们希望计算机猜测5个甚至10个。为了测试参数是否正确，我们在MM中填写参数，然后用数值方法计算ODE。随后，我们用我们知道(并想要猜测)的参数计算模型的计算值与我们猜测参数的MM的计算值之间的误差。计算模型ODEs的值被多次执行，ODEs实时表示一分钟，我们计算24小时，因此计算1440次。
目前我们正在使用粒子过滤器来猜测变量，这工作得很好，但我
浏览 7提问于2019-12-26得票数 0
我知道如何在MATLAB中使用ode15s或其他ode求解器，但我不确定的是，这段代码(来自)对我来说似乎很模糊：
[VOI, STATES] = ode15s(@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS), tspan, INIT_STATES, options);
更具体地说，以下(？)是什么意思：
@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS)
该函数的头文件"computeRates“如下：
function [RATES, ALGEBRAIC] = computeRa
浏览 1提问于2011-03-29得票数 3
回答已采纳
2回答
它是否有可能采取一个步骤与枕叶代码积分器？我知道step参数是scipy.integrate.ode.integrate的，但我无法确定如何设置t参数，以保证它能够完全完成一步。显而易见的选择是t=inf，但这会导致ODE失败：
from scipy.integrate import ode
obj = ode(lambda t, y: -y)
obj.set_initial_value(1)
y_new = obj.integrate(inf, step=True)
-
DVODE-- At T (=R1) and step size H (=R2), the
co
浏览 0提问于2016-04-24得票数 0
回答已采纳
我的问题与有关，它解释了如何用VexCL实现boost::odeint，以实现并发(完整的代码可以找到)。
下图显示了我如何看待ODEINT的迭代：
现在我扪心自问，到底是什么/或者它的哪一部分是VexCL中的并行化的？
我的印象是，ODE部分是一个单一的任务，因为在给定的示例中，所有ODE方程都在一个块内。也许集成部分运行在三个并行任务中。这导致了四个任务，其中(我认为) ODE任务是一个瓶颈(因为方程可能变得非常大)。
如果这是正确的，我想知道，如何改进这个并发。我认为横向结合ODE和INT是有意义的。这就产生了3个任务，每个任务都不能在这个级别上进一步减少。
浏览 2提问于2020-02-11得票数 0
回答已采纳
1回答
我需要一个关于如何处理不同的几何形状(静态的，动态的)在opengl的小指南。
我试图开发一个图形引擎(作为学习)，我对opengl有一些想法(我做了一些事情)，我需要问一些google没有回答的问题：
在vbo中存储所有几何数据的正确方法是什么，然后绘制它。我是如何将这些数据组织成vbo的？我想这个：
全部分开
->静态和动态数据
2-记住，在世界上有不同类型的三维几何静态，动态等，这将是最理想的方式来处理它们？动态元素在更新每个帧时都应该有自己的vbo和ibo，而静态元素应该被处理吗？
假设游戏世界由3个静态立方体(具有相同纹理)、5个静态锥(每一个纹理不同)、5个动态特征(每个帧中
浏览 3提问于2013-11-15得票数 2
回答已采纳
我想在现代Fortran的两个类之间传递一个过程指针。这个过程指针应该
从第二个对象中调用，以访问第一个对象的组件，而不使用它作为虚拟参数。
这里有一个明确的例子，想象一下对ODE解决程序执行一个面向对象的包装：
module test_funptr
implicit none
public
type, public :: ode_solver
integer :: NEQ = 0
procedure(ode_api), pointer, nopass :: f => null()
contains
pro
浏览 6提问于2022-04-01得票数 3
回答已采纳
免责声明:在上交叉发布
目标：I试图在R中数值求解Lotka码，在中使用de sde.sim()函数。我想使用sde.sim()函数，以便最终将这个系统转换为SDE。所以最初，我开始使用一个简单的ODE系统(Lotka模型)，没有噪声项。
ODE系统：
的初始值为x= 10和y= 10。
α、β、delta和γ的参数值分别为1.1、0.4、0.1和0.4 (模拟此)。
试图解决问题：
library(sde)
d <- expression((1.1 * x[0] - 0.4 * x[0] * x[1]), (0.1 * x[0] * x[1] - 0.4 * x[1]))
浏览 0提问于2019-02-12得票数 1
我正在尝试将我的实验数据拟合到一个包含两个耦合的ODEs的模型中。我已经创建了一个单独的函数来定义两个ODE，其中a是包含我想要匹配的2个常量的数组。
function dy = CalcCoupODEs(t,y,a)
%dy= matrix with all 2 ODEs
dy = zeros(2,1); % variable definition
dy(1) = 1.41.*a(1) .* ((y(2)*y(1)^2)^2)*y(1) - a(2) .* y(1)*t ; %ODE for R
dy(2) = -a(1)*(y(2)*y(1)^2)^2; %ODE for n
end
然
浏览 1提问于2018-06-18得票数 2
回答已采纳
我正在使用ODE求解器来解决我的ODE系统。现在，我还想通过改变参数来计算各种参数的灵敏度。我的基本ODE解决程序是这样的:函数：
function F = ODE_site(t,y)
%%Parameter block
k1f=1e8;
k1b=1e5; %%and so on
%%Elementary rate block
r1=k1f*y(1)-k1b*P*y(5);
and so on
%%Mass balance block
F=[ r1-r5-r78+r86 ;
and so on ];
%% End of function
主要的ODE解决程序如下所示：
optod
浏览 7提问于2015-03-23得票数 0
回答已采纳
我正在尝试使用SIR模型对流感流行进行建模：。
它基本上需要我同时解3个微分方程。
我的任务是将图表与一组数据进行拟合，我只是想知道我应该如何进行这项工作。
我已经成功地建立了一个系统来解决如下的ODE：
function dydt = sir_ode(t,y,p)
B = p(1);
r = p(2);
S = y(1);
I = y(2);
R = y(3);
dydt = [-B*I*S; B*I*S - r*I; r*I];
以及调用它的m文件：
time = 1:24;
y0 = [400 1 0];
tspan = time;
p0 = [.01 .5];
[t,y] = ode4
浏览 4提问于2013-05-23得票数 1
回答已采纳
1回答
引用类的成员函数的指针。
、、、
本质上，我想要的是将成员函数的引用传递给另一个类的另一个函数。示例代码如下所示：
#include<iostream>
using namespace std;
class algorithm
{
void (*ODE)(double * timeDeri, double * var);
void simulation()
{
//use ODE to do simulation
}
parent(void (*someODE)(double * timeDeri, double * var))
{
浏览 1提问于2014-01-16得票数 0
回答已采纳
我们有一个用于Revit的设计自动化应用程序，可以处理BIM360中的模型，我们希望使用自己的流程打开模型，以便可以修改打开的工作集，等等。我们在找到了该链接。但是，根据这一点，您必须将rvt文件的名称硬编码为特定名称才能打开它。但是，我们有一些地方在导出的数据中使用了模型的名称，因此我们希望保留最初显示的名称。
是否可以从应用程序中引用命令参数？现在，我们已经为活动上的模型定义了如下参数：
"rvtFile": {
"zip": false,
"ondemand": false,
"
浏览 7提问于2021-07-06得票数 2
我正在尝试做与本文中描述的相同的事情：https://science.sciencemag.org/content/sci/early/2020/04/07/science.abb4557.full.pdf 我使用GNU Octave。 以下是ODE系统的函数文件： function f = COVID19ODE(t,y0)
alpha = 3.07*0.38;
beta = 0.38;
kappa = 0.5;
kappa0 = 0.5;
ydot =@(t,y,a) ([-alpha*y(1)*y(2) - kappa*y(1);
浏览 26提问于2020-04-12得票数 0
我在matlab中实现了一种滤波算法，在参数传递方面遇到了一些困难。
我的动力是由功能控制的：
function zdot = VortexDynamics(t, z, vorticity)
for ii=1:length(vorticity),
zdot(ii,1) = 1i/(2*pi) * sum( vorticity([1:ii-1,ii+1:end]) .* ( z(ii) - z([1:ii-1,ii+1:end]) ) ./ (abs(z(ii) - z([1:ii-1,ii+1:end]) ).^2 ) ) ;
end
end
我最初用ode45计算轨迹，没
浏览 3提问于2014-06-13得票数 0
回答已采纳
1回答
我试图用dopri5积分器在scipy.integrate.ode中解决一个简单的例子。如文件所述
这是一个显式的龙格-库塔方法的(4)5阶由于休眠和普林斯(步进大小控制和密集输出)。
这应该能行。下面是我的例子：
import numpy as np
from scipy.integrate import ode
import matplotlib.pyplot as plt
def MassSpring_with_force(t, state):
""" Simple 1DOF dynamics model: m ddx(t) + k x(t)
浏览 4提问于2015-10-03得票数 2
1回答
我试图为一个质量可以忽略不计质量的物体建立一个基本的引力模型。我遵循了ODE套件文档中提供的示例，但是我得到的结果显然是荒谬的。
下面是我用ode45调用的函数
function dy = rigid(t,y)
dy = zeros(4,1);
%Position
xx=y(1);
yy=y(2);
%Radius
r=(xx.^2+yy.^2).^0.5;
%Constants
M=10^30;
G=6.67*10^-11;
%dX/dt
dy(1)=y(3); %vx
dy(3)=-M.*G.*xx.*r.^-3; %ax
%dY/dt
dy(2)=y(4); %vy
dy(4)
浏览 3提问于2015-02-14得票数 1
回答已采纳
我有一个函数(对于SIR模型)，然后有一个脚本来解决这个函数，并将它与我试图适应这个模型的数据进行比较。因此，我试图运行一个for循环来更改函数中的一个参数，以便优化fit。我想知道如何在for循环中更改我的(r)和(a)参数，而不必手动更改它们：
function ydot=epidemic(t,y)
r=0.000001;
a=1/3;
ydot=zeros(3,1);
ydot(1)=-r*y(1)*y(2);
ydot(2)=r*y(1)*y(2)-a*y(2);
ydot(3)=a*y(2);
end
和
[t,y]=ode45('epidemic',[0:222],
浏览 7提问于2016-03-21得票数 1
回答已采纳
1回答
我正试图解一个二阶微分方程。我已经设法解决了一个角度(函数中有一个因素取决于一个角度)使用ode45在matlab中。在我的程序脚本中，也有与ode45的输入参数相同的角度。
我想要的是求解不同角度的方程。我用一个for -循环进行了测试，我只需在ode45中运行给定的角度(作为输入参数)，这样就可以绘制出不同的曲线。但是，我不太确定它是否正确，因为我没有改变ode45解决的函数的角度。
所以问题是，对于每个循环，我是如何改变函数内部的角度的。
浏览 3提问于2014-04-18得票数 0
回答已采纳
我能在Scilab中找到ode解算器(ode23、ode45和ode113)吗？我在MATLAB中使用了这些求解器，但是我不知道Scilab中是否有相同的选项。提前谢谢。
浏览 17提问于2022-04-07得票数 1
1回答
我们可以将R中的deSolve包用于常微分方程( ODE )，但是，我无法找到一种方法来求解两个嵌套的ODE方程。
b'(t) = beta - k*b(t);
a'(t) = alpha -b(t)*gamma;
其中'的意思是分化。那么，我们如何解决a和b呢？由于a'是b的一个函数，我们必须同时求解a和b。
我有个错误：
lsoda中的错误(y，time，func，parms，.)：使用的求解器组合不能嵌套。
当我试图将b的解决方案添加到a的ode解决方案中时。
浏览 8提问于2011-04-12得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
