对5个常微分方程进行建模并绘制模型以显示这5个方程
常微分方程 (Ordinary Differential Equations, ODEs) 是描述物理、工程、经济等领域中变量关系的数学方程。建模和绘制模型可以帮助我们理解和预测实际问题的行为。下面是对5个常微分方程进行建模并绘制模型的完善答案：
常微分方程1：一阶线性微分方程 一阶线性微分方程的一般形式为 dy/dx + P(x)y = Q(x)，其中P(x)和Q(x)是已知函数。它可以用来描述一些物理系统中的衰减或增长过程。 应用场景：电路分析、放射性衰变等。 推荐腾讯云产品：云服务器（https://cloud.tencent.com/product/cvm）
常微分方程2：二阶线性非齐次微分方程 二阶线性非齐次微分方程的一般形式为 d^2y/dx^2 + p(x)dy/dx + q(x)y = r(x)，其中p(x)、q(x)和r(x)是已知函数。它可以用来描述振动系统、电感电容电路等。 应用场景：机械振动、电路振荡等。 推荐腾讯云产品：云数据库 MySQL 版（https://cloud.tencent.com/product/cdb）
常微分方程3：二阶非线性微分方程 二阶非线性微分方程的一般形式为 d^2y/dx^2 = f(x, y, dy/dx)，其中f(x, y, dy/dx)是已知函数。它可以用来描述自由落体运动、混沌系统等。 应用场景：物体自由落体、混沌系统模拟等。 推荐腾讯云产品：云函数（https://cloud.tencent.com/product/scf）
常微分方程4：一阶非线性微分方程组 一阶非线性微分方程组的一般形式为 dx/dt = f(x, y)、dy/dt = g(x, y)，其中f(x, y)和g(x, y)是已知函数。它可以用来描述生物种群动力学、化学反应动力学等。 应用场景：生物种群模型、化学反应动力学模拟等。 推荐腾讯云产品：云原生容器服务（https://cloud.tencent.com/product/tke）
常微分方程5：高阶微分方程 高阶微分方程是指阶数大于二的微分方程，例如三阶、四阶等。它们可以用来描述弹性体力学、电磁场分布等复杂问题。 应用场景：弹性体力学模拟、电磁场分布模拟等。 推荐腾讯云产品：云硬盘（https://cloud.tencent.com/product/cbs）
以上是对常微分方程进行建模并绘制模型的完善答案。请注意，答案中的腾讯云产品仅作示例，实际选择产品应根据具体需求进行决策。
相关·内容
我是netlogo的新用户。我有一个反应系统(转化为常微分方程)，可以用Matlab求解。我想在netlogo中开发相同的模型(用于与matlab结果的比较)。我对time/tick感到困惑，因为netlogo使用"ticks“表示时间增量，而Matlab使用秒为单位的时间。如何将我的matlab秒转换为刻度数？有人能帮我写代码吗？模型是：
A + B ---> C (with rate constant k1 = 1e-6)
2A+ C ---> D (with rate constant k2 = 3e-7)
A + E ---> F (with rate co
浏览 2提问于2014-12-30得票数 1
有没有人尝试用Modelica实现Navier Stokes偏微分方程(PDE)？我发现了空间基函数(SBF)的方法，它通过数值修改得到Dymola可以处理的常微分方程(ODE)。
致以敬意，
维克多
浏览 2提问于2015-06-17得票数 4
谁能告诉我如何选择ODE求解器来设计汽车模型？在simulink中有各种常微分方程求解器。我可以知道每个ode求解器的区别以及这些求解器应该用于哪种类型的应用程序吗？
如果我得到一些关于这方面的指导，那将会有很大的帮助。
在此，我要对所有的祝福表示感谢。
浏览 0提问于2014-12-12得票数 1
我在NetLogo中为一个模型编写了一个简单的代码。同时，通过文献中的常微分方程对模型进行了深入的研究。现在我想比较一下NetLogo和Matlab (用于求解微分方程)得到的一些模型图。我使用了"ticks“命令来增加NetLogo中的时间，就像Matlab使用时间(以秒为单位)一样。为了比较NetLogo和Matlab得到的图形，我应该记住哪些预防措施(变化)。
浏览 2提问于2014-09-02得票数 2
回答已采纳
1回答
我有一个矩阵，其中包含通过常微分方程计算的行星运动的数据，然后绘制，我需要单独绘制每个行星，是否可以通过仅使用一个矩阵的集合函数来实时完成这一点，或者我是否必须相应地分解矩阵？(太阳，水星，金星，地球，火星，木星，土星，天王星，海王星，冥王星的x，y，z位置和速度有60列)这个矩阵中我需要绘制的列是每六组中的前三列
谢谢
克里斯
浏览 0提问于2013-02-12得票数 0
2回答
public static void sortByNumber(Course[] list) {
Course temp = new Course();
boolean fixed = false;
while(fixed == false) {
fixed = true;
for (int i = 0; i<list.length-1; i++) {
if (list[i].getNum() > list[i+1].getNum()) {
temp.setNum(list[i].getN
浏览 4提问于2015-12-02得票数 0
回答已采纳
我是Python3.4的新手，正在考虑使用Python来解决细胞生理学中的大型反应扩散隔室模型，并且想知道使用Pythons odeint例程可以求解多少联立的一阶常微分方程。我的问题需要同时解决至少300个一阶常微分方程。odeint能做到这一点吗？
浏览 0提问于2015-10-26得票数 1
引言
利用ODE45对二阶常微分方程组进行数值求解，用Matlab对一些动态系统进行了数值模拟。我从Mathworks那里找到了一个很好的教程(最后是教程链接)，介绍了如何做到这一点。
在本教程中，方程组在x和y中是显式的，如下所示：
x''=-D(y) * x' * sqrt(x'^2 + y'^2)
y''=-D(y) * y' * sqrt(x'^2 + y'^2) + g(y)
上述两个方程均为y'‘= f(x，x'，y，y')
问题
然而，我遇到的方程组中，变量不能显式地求解，如
浏览 4提问于2017-01-24得票数 1
在我的模型中，每个智能体在每个节拍上求解一个ODE系统。我采用了欧拉法(类似于NetLogo中的系统动力学模型)来求解这些一阶常微分方程。然而，为了获得稳定的解决方案，我被迫使用非常小的时间步长(dt)，这意味着使用这种方法进行模拟的速度非常慢。我很好奇，有没有人有建议可以更快地解决这些常微分方程？我正在考虑实现Runge-Kutta (具有更大的时间步长？)正如这里所做的()。我也会考虑使用R扩展和在R中使用ODE求解器，但同样，ODE是由每个代理求解的，所以我不知道这是不是一种有效的方法。
我希望有人对这些方法的性能有所了解，并能提供一些建议。如果没有，我会尝试分享我的发现。
浏览 0提问于2015-04-10得票数 1
我有一个代数方程的模型，偏微分方程和常微分方程。我用有限差分法做了偏微分方程的离散化，就像Tiller的例子6.9，用Modelica介绍PhysicaL建模一书。
然后，我的模型的初始条件是在时间零点离散化的所有不同点的数组。在我的模型中有一个输入：
问题是我需要在模型外部设置这个数组，并出现以下错误：'FMUException: valueref和value的长度不一致。‘
有谁知道是什么导致了这个错误，以及如何解决它？谁能给我提供一个例子或一段代码，说明如何在JModelica中实现这一点？
谢谢你，卡琳
浏览 5提问于2020-09-21得票数 1
1回答
Julia是否支持边值微分代数方程？我有一个隐式的常微分方程，它的可变质量矩阵有时是奇异的，所以我必须使用DAEProblem。我的问题是x1(t)和x2(t)的两个耦合的二阶常微分方程，通过设置x1'(t) = y1(t)和x2'(t)=y2(t)，我已经将它们转换为四个一阶方程。我在域的开头和结尾有x1和x2的值，但没有y1或y2 anywhere的值，因此我需要DAE和BVP值。 This github的帖子建议这是可能的，但我担心我对机制的理解还不够深入，无法理解如何将DAEProblem与BVProblem结合起来。 我已经成功地按照数字食谱编写了多个拍摄代码来解决
浏览 30提问于2021-07-08得票数 2
1回答
我正在尝试绘制一个常微分方程组的解。代码是：
tspan = [0 10];
z0 = [0.01 0.01 0.01 0.01];
[t,z] = ode45(@(t,z) odefun3(t,z), tspan, z0);
plot(z(:,3))
为什么输出绘制在区间0,60上，而不是代码中的0,10？
浏览 3提问于2018-06-10得票数 0
1回答
我试图将三个不同的函数绘制成一个图，一个常微分方程的近似解，一个常微分方程的精确解和相对误差。我的函数只输出这些图形中的一个。我不确定如何一次绘制所有这三个参数，以及如何更改k的参数，以便在.02、.1和1.5时绘制k。
我在一个单独的文件中有这个函数。这是近似的解决方案：
function dydt = prblm1_fun(t,y)
k = 0.02
if t < 1
dydt = y(1)*k + sin(pi*t);
elseif t < 2
dydt = y(1)*k + 0.5;
else
dydt=y(1)*k + exp(t-2)/2;
end
此函数尝试求解近似值
浏览 5提问于2017-03-15得票数 0
1回答
我有一个关于解一阶复微分方程的问题。我使用Runge-Kutta，但答案似乎不正确。
这是我的方程式：
y'=exp(-2*t)-i*y
常微分方程的结果是可以的，但对于复杂的方程，它似乎不是正确的。
我用gfortran编译了我的代码，并用gnuplot绘制了没有i的常微分方程的图表。这张图对于方程的解是正确的。但我在方程式中加入了i，图中只有一条直线。另外，gnuplot只绘制真实的部分。
编译代码时没有任何编译时错误。
这是我的代码：
program kutta
implicit none
real:: a=0.0, b=8.0,y=0.1,h,t
integer::n=40
c
浏览 1提问于2016-01-20得票数 2
Python中有哪些库程序包可用于动态系统的模拟和分叉研究？ 我想研究常微分方程中的分叉现象。我知道Fortran中有AUTO。我还在皮顿找到了PyDSTool。PyDSTool的问题是它与64位系统和Python3(以及numpy等)有一些不兼容。我想知道Python中是否有用于此目的的最先进的库。谢谢。
浏览 43提问于2019-05-24得票数 2
回答已采纳
4回答
我正在用c++编写程序来数值求解微分方程(常微分方程和偏微分方程组)。原因很简单--对于大规模计算，c++比matlab快得多。然而，当涉及到绘图时，matlab显然是赢家。有没有什么方法可以让matlab像在c++中绘图一样？
浏览 1提问于2014-05-26得票数 4
我有一个关于在Julia中如何训练神经ODE网络的通用问题。是否从定义了标准常微分方程的tspan中采样数据点，并在这些tspan上计算参数更新？换句话说，在训练过程中是否发生了一些混洗和批处理，或者是否计算了tspan中所有数据点的损失？
浏览 46提问于2020-02-19得票数 1
3回答
我正在尝试解决仅限于正解的常微分方程，即：
dx/dt=f(x)
和x>=0在一起。
在MATLAB中，这是非常容易实现的。R是否有什么变通方法或包来将解空间限制为仅为正值？
这在我的情况下是非常关键的，不幸的是没有其他选择。我现在找了一段时间，但没有任何成功。:-(
浏览 0提问于2017-11-21得票数 3
我是matlab的新手，想要解决常微分方程的
x'(t) = a - bx(t) + x(t)^2 * y(t) - x(t)
y'(t) = bx(t) - x(t)^2 * y(t)
其中我知道x(0) = y(0) = 0.9，并且a= 1，b=为1.9；
理想情况下，我希望最终绘制出解决方案。我在网上看了看，但我想不出该怎么做。任何帮助都是最好的。谢谢。
浏览 1提问于2015-04-30得票数 0
4回答
从一条曲线的作图出发，可以得到该曲线的参数方程？
特别地，假设x={1 2 3 4 5 6}是x轴，y= {a b c d e f...}是相应的y轴。我有图(x，y)。
现在，我如何才能得到描述绘制曲线的方程？可以显示从样条插值开始的参数方程？
谢谢
浏览 0提问于2010-09-30得票数 2
1回答
我目前正在尝试寻找一种机器学习算法，可以预测数学模型(MM)中使用的大约5- 15个参数。MM有4个不同的常微分方程(ODE)，并且将添加更多的常微分方程，因此需要更多的参数。大多数参数是可以测量的，但其他参数需要猜测。我们知道所有的15个参数，但我们希望计算机猜测5个甚至10个。为了测试参数是否正确，我们在MM中填写参数，然后用数值方法计算ODE。随后，我们用我们知道(并想要猜测)的参数计算模型的计算值与我们猜测参数的MM的计算值之间的误差。计算模型ODEs的值被多次执行，ODEs实时表示一分钟，我们计算24小时，因此计算1440次。
目前我们正在使用粒子过滤器来猜测变量，这工作得很好，但我
浏览 7提问于2019-12-26得票数 0
我是matplotlib/python的新手。我想用odeint或solve_ivp在某个区间内解决一些简单的常微分方程，例如segment0，10，其中“初始条件”是t_0= 5，而不是t_0 = 0。在我看来，odeint和solve_ivp都将初始条件作为值y_0，也就是说，y_0是函数在指定时间数组t开始时的值。
在我的积分间隔中间有一个条件，我如何解决一个常微分方程？
浏览 2提问于2020-05-17得票数 1
2回答
我是Python初学者。我正在尝试切换matlab中的一些程序。我需要解一个刚性常微分方程，它的输入都是矩阵。在matlab中，我使用
[ttT,uT] = ode23s('SST',t,fT);
浏览 0提问于2012-01-05得票数 4
回答已采纳
1回答
有人问我这个问题，即使经过大量的研究，我也找不到答案。
当我使用实体关系(ER)设计数据库，然后将其重新绘制为UML类模型(用UML的数据建模概要文件进行扩展)时，这两个模型有什么不同，这两个模型如何(以它们的方式)对db开发人员有用？关于多重性，外键。
浏览 0提问于2014-01-01得票数 3
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
2回答
为了说明具有无穷解族的一阶常微分方程的柯西问题，我将绘制参数化的解，并有可能通过滑块控制参数值。
完整的柯西问题是y'=sqrt(|y|)，y(0)=0，参数解是y_c(x):={0，if c=>x；(x-c)^2/4，if x=>c}。
所以我会得到一个y=y_c(x)的曲线图，用一个滑块来控制c的值。
谢谢。
浏览 1提问于2012-04-15得票数 4
回答已采纳
1回答
正如标题所说，我想知道如何从Matlab常微分方程中检索二阶导数。我试图模拟的系统用12个一阶微分方程来描述。
让我们来解释一下：
...
[T,Y] = ode113(@sixdofsolver,time,Y0,options,settings);
T= 1700x1向量
Y= 1700x12矩阵
现在，如果我这样做了：
[dY] = sixdofsolver(T,Y,settings)
dY = 12x1向量
我原本期望有一个和Y一样大小的矩阵。
我做错了什么？
浏览 0提问于2013-10-05得票数 0
我应该写一个MATLAB函数，它取一阶常微分方程，形式为y‘(T)= a*y(t) +b，初始点y(t0)=y0作为输入，并计算解的前15个点。还会绘制前15个点的解决方案曲线。我们要解的方程是:y‘(T)= 4*y(t)+1，初始点y(0)=0。
对于这个函数，我写了下面的代码，但这给了我一个关于y的错误。我应该如何正确地实现euler函数？而且我也不能确定如何绘制解决方案曲线。
function E=euler(f,y)
%Input - f is the function entered as a string 'f'
% - a and b are the left a
浏览 0提问于2012-12-17得票数 0
回答已采纳
嗨，我正在尝试实现一个简单的塑性有限元代码。我不清楚如何设置通过拉格朗日乘数应用的位移控制。在线性问题的情况下，我执行了以下操作来应用显示控制，并且它起作用了。如果能量平衡是线性的，那么将得到以下方程：
然而，在非线性问题中，K是K(u)并且是非线性的。因此能量平衡方程及其导数如下所示
感觉不对劲。我该如何评估它？如何在非线性情况下修改方程的矩阵形式？那么矩阵方程是如何修改的呢？我使用添加了用latex输入的图像。如果不清楚，请让我知道。
浏览 9提问于2021-06-24得票数 0
我搜索了JAGS的手册和一篇关于常微分方程( ode )的帖子()。我认为JAGS可以，因为它类似于WinBUGS (通过WBdiff接口完成)。然而，如果我让JAGS在我的ode代码中阅读，它甚至无法识别D(y[...], t)表达式。
JAGS能处理ode吗？也许我错过了JAGS的插件，比如WBdiff。
浏览 10提问于2022-10-15得票数 1
回答已采纳
我对具有不同初始条件的常微分方程的多个解决方案进行平均，因此所有解决方案在相同的时间具有值是很重要的；例如，增量为0.01。
我一直在使用数值食谱3 (nr3)中的ODE例程。它们执行自适应大小步长，并使用计算的值进行相同顺序的插值。我不能使用它们，因为它们与boost冲突。有没有其他类似的例程？
我看过GSL，它很不错，但它没有内置插值。我可以这样做的一种方法是求解具有自适应大小的常微分方程，然后运行Akima插值。但是看起来nr3解决方案会更快更准确。
浏览 0提问于2011-10-05得票数 2
回答已采纳
我想画出三维常微分方程组的数值方法的结果。我的输出是这样的(假设我们计算了三个步骤)：
import numpy as np
v= np.array([[1,2,3], [4,5,6], [7,8,9]])
其中，每个三元组中的第一个值是x坐标，第二个值是y坐标，第三个值是z坐标。
我想用最简单有效的方式在3D网格上绘制这些点。问题似乎在于数据应该像np.array([[1,4,7], [2,5,8], [3,6,9]])一样格式化。
浏览 4提问于2016-11-01得票数 6
回答已采纳
所以我们试着解决简单的常微分方程。我们使用文本。这是本书(第18页)中的示例代码：
> with(DEtools); with(plots);
> eq1 := diff(x(t), t) = y(t); eq2 := diff(y(t), t) = -sin(x(t));
> init := [[x(0) = 1, y(0) = 0], [x(0) = 0, y(0) = 2]];
> DEplot([eq1, eq2], [x, y], t = 0 .. 110, x = -2*Pi .. 2*Pi, init, scene = [x, y], stepsize
浏览 1提问于2011-05-13得票数 0
1回答
我有一个目标函数，它根据一些参数来计算非线性常微分方程系统的实际值与我的预测值之间的差值。
residuals_vector = actual_values-predictions
P参数N残差的[residuals_vector]=objfun(actual_values,time,parameters)，在Matlab语言中有没有一种简便的方法计算该系统的灵敏度矩阵P x N
浏览 0提问于2019-03-31得票数 0
1回答
我有一个常微分方程组，它有两个吸引子，一个在(1，0)，另一个在(-1，0)。我想在笛卡尔坐标中绘制一个吸引盆地，其中有两种颜色，显示随着时间趋于正无穷大，坐标的每个点上的吸引子将结束。但是，我不知道如何使用matplotlib绘制这样的图形。这是我现在所得到的：
from scipy.integrate import ode
import numpy as np
import matplotlib.pyplot as plt
from numpy.linalg import norm
"""
The system of ODE:
x' = y
y'
浏览 3提问于2018-10-18得票数 0
2回答
SciPy可以用scipy.integrate.odeint或其他软件包求解常微分方程，但要在函数完全解完后才能给出结果。但是，如果ode函数非常复杂，程序将需要很长时间(一到两天)才能给出完整的结果。那么，我如何控制它求解方程的步骤(当方程还没有完全求解时，打印出结果)？
浏览 0提问于2019-11-26得票数 2
2回答
我目前正在评估Dymola 2020建立太阳能发电厂模型，但我很难在使用自定义连接器的模型中计算出方程的数量。
以下是代码：
connector Port
flow Real Q;
Real P;
Real T;
end Port;
model Inlet
parameter Real Q = 1;
parameter Real P = 2;
parameter Real T = 3;
Port a;
equation
a.Q = Q;
a.P = P;
a.T = T;
end Inlet;
我在模型中设置了3个方程，但是有一个误差显示有4个方程。方
浏览 3提问于2019-11-11得票数 2
回答已采纳
1回答
我有一个形式为$\frac{dy}{dx}=f(x，y)$的常微分方程，我设法在Matlab中绘制了它的解。我想在曲线上添加“前进方向箭头”。下面所示的曲线是x(t)和y(t)的参数图，每个曲线都有1500个步长。
左边的黑色箭头是使用图中窗口中的"Insert Arrow“选项插入的。我想添加箭头，如在黑白图片所附。
有什么办法可以做到吗？任何帮助都将不胜感激。
浏览 0提问于2016-05-19得票数 0
我想在COMSOL中求解具有特定边界条件的Schrödinger方程。作为一个常微分方程，schrödinger方程是这样写的( 1D)：
af''(x) + b(x)f(x) = Ef(x)，
其中E是将由边界条件确定的未知常数。
我不习惯使用COMSOL，所以我不知道是否有可能解决这个问题。到目前为止，所有用于求解微分方程的模板都包含一些通用形式，其中您必须在每个项之前指定常量的值。这不适用于上面的特征值问题，其中E是未知的。有谁知道如何将微分方程指定为特征值方程，其中E是未知的？
浏览 40提问于2016-04-04得票数 2
2回答
我目前正在用deSolve来解一个常微分方程组，我想知道是否有任何方法可以防止微分变量值降到零以下。我看过一些关于在向量、数据框架等中将负值设置为零的帖子，但由于这是一个生物学模型(T细胞计数变为负数是没有意义的)，我需要从一开始就阻止它发生，这样这些值就不会扭曲结果，而不仅仅是替换最终输出中的负数。
浏览 5提问于2017-01-14得票数 8
回答已采纳
1回答
在libgdx 2d中，我们通常不渲染不在屏幕内的精灵，然而在3d中，从我读到的内容中，你绘制了整个场景作为关卡，我的场景相当大，我想知道是否应该处理屏幕外的模型？又是如何做到的？
浏览 31提问于2019-09-25得票数 1
回答已采纳
1回答
我正在尝试在Matlab中模拟方程y“+ 4y = 0，初始条件y(0) =1和y'(0) =0。下面的代码片段显示了上述常微分方程的位置作为时间函数的工作近似值：
clear
syms y(t)
%Differential Equation Conversion
P = odeToVectorField(diff(y, 2) == -4*y);
M = matlabFunction(P,'vars',{'t','Y'});
%Position Approximation
pos_solution = ode45(M,[0
浏览 3提问于2016-05-24得票数 0
1回答
我正在寻找一个(可以是付费版本)软件，允许我：
绘制给定三维模型的二维截面，在任何特定平面上裁剪
包含足够信息的2D横截面允许我进行线跟踪。
软件应该有一个API，可以让我直接从C++/C#编程图形用户界面(如WinForms)绘制横截面图。
为了详细说明，我实际上正在从事一个涉及超声波检测的项目。本质上，我的程序接受输入参数，如厚度，半径，长度和两个管道交汇的角度。因此，我将用管道的参数方程来表示(而不是点云格式)。
我要找的是一种在三维模型中可视化超声波光束路径的方法。
浏览 0提问于2015-02-18得票数 1
1回答
我目前正在解决一个涉及解决一些二维常微分方程的问题。我的问题在于计算与弹簧挂钩的测试质量的位置，到二维电势的测试位置。
我想用鼠标在一些矩形上自定义支持vía绘图的位置，支持将运行的路径，最好是在潜在的等高线地图上绘制。
我是python的新手，不知道最好的方法是什么。我正在为这个项目使用numpy和matplotlib，以防它可以用这些包完成。
非常感谢。
浏览 0提问于2016-03-20得票数 0
是否有可能使用scipy.integrate.odeint对任何常微分方程进行时间倒序积分？如果可能的话，谁能告诉我odeint中的争论'time‘应该是什么？
浏览 0提问于2012-11-05得票数 2
回答已采纳
在我的常微分方程函数中，我需要迭代地求解一个参数的方程，直到在每个时间步收敛。我希望传递最新的参数值作为下一个时间步的初始值，这样当函数对参数进行迭代更新时，花费的时间会更少。但我想不出该怎么做。ODE函数的代码结构如下：
from scipy.integrate import solve_ivp
def run(t, y):
if t==0:
a = 1e-8
nn = 0
while nn<=100:
nn = nn +1
#update a until convergence
浏览 2提问于2020-08-08得票数 0
2回答
对框2D中的一袋液体进行建模
、、、、
我将如何在框2D中对一袋膨胀的流体进行建模，即哪些物理方程在建模中将是有用的？当触摸时，装水的袋子可以移动，但这是唯一的互动。任何地方，方程式或模型都会非常感谢！
浏览 2提问于2017-06-06得票数 2
首先，我对java非常陌生，抱歉。我试着模拟地球绕着太阳转。在绘制了结果之后，地球在一次旋转后就明显地离开了轨道！
我已经检查和双重检查了常数，如太阳和地球的质量，以及初始速度和位置。我不知道我做错了什么，这些方程也得到了同事和讲师的证实。
代码由4个类组成：
发现。
浏览 2提问于2015-12-27得票数 0
回答已采纳
1回答
我必须在matlab中积分一个8变量的常微分方程。我的模拟时间是5e9，时间步长是0.1。但它显示内存错误。我使用的是i7核心,2.6Ghz中央处理器和8 8GB内存。如何模拟大量时间样本的ODEs？
浏览 1提问于2017-08-26得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
