嘿，我有一个模型，可以给出微分返回值。, z 我的求解器可以求解这个方程，并在相应的时间值给出x和y。t = np.linspace(0, 10, 100)result = solve_ivp(model, [np.min(t), np.max(t)], f0, t_eval=t) 但是现在我想要我的z的解决方案，它不应该被求解器积分。有没有可能只对
浏览 19提问于2021-01-26得票数 0
我想用数值方法解决一个复杂的常微分方程系统。最初，我使用NDSolve，但是因为我最终需要进行大量的参数搜索，所以我决定改用Python，特别是SciPy的solve_ivp函数。然而，我很快就遇到了以下问题：
起初，当
浏览 3提问于2020-07-12得票数 1
这是我关于堆栈溢出的第一篇文章。因此，我在SciPy的文档中看到了这个solve_ivp求解器的例子。正在解决的问题如下： 1.11088891e-01 1.11098890e+00 1.1109989
浏览 0提问于2020-01-23得票数 1
回答这一问题的过程已经在下面的链接问题中开始，但这一专题具体是关于整合一个功能，并得到了答复。所以我增加了一个新问题。
from numpy import linspaceax2.plot(s.t[::2], s.y[1][::2], label=
浏览 3提问于2020-11-26得票数 0
回答已采纳
分配给我的任务要求我在500个时间步长内求解微分方程。如何使用SciPy实现这一点？t_span参数允许我配置t的初始值和终值，但我实际上对此并不感兴趣--我只对我集成了多少次感兴趣。(例如，当我使用t_span = (0, 500)运行我的方程时，求解器积分907次。)下面是我的代码的简化示例：
def simple_d
浏览 9提问于2021-05-13得票数 1
1回答
Solve_ivp是一个初值问题的求解函数。t_span，y0，method=‘RK45 45’，t_eval=None，dense_output=False，events=None，vectorized=False，args=None，**选项)参数有趣，系统右侧可调用.
浏览 2提问于2020-02-18得票数 3
1回答
我一直在尝试使用odeint和solve_ivp来求解以下三个耦合微分方程系统。(这不是确切的系统。我对它进行了修改，只是为了提高问题的可读性。)但现在我想用下面的方式修改它:当u，v或w中的任何一个低于阈值时，它被重置为零，然后让系统进化。每当这三个参数中的任何一个自动超出阈值时，都会发生这种情况。进化系统的“规则”保持不变。稍后，我将把它推广到由数百个耦合方程组成的系统。
浏览 2提问于2020-09-28得票数 0
1回答
我想知道，对于一个微分方程，ODEINT和solve_ivp之间有什么区别？它们之间有什么优点和缺点呢？f1 = solve_ivp(f, [0,1], y0) #y0 is the initial point
f2 = odeint(f, y0, [0, 1], args=(a, b)) # a and
浏览 2提问于2019-03-13得票数 5
回答已采纳
1回答
我想集成一个控制问题，即形式为dx/dt = A.f(t)的常微分方程，其中x(t)是R^3中的函数，f(t)是R^4中的函数，A是矩阵3x4。在我的特殊情况下，f(t) = F'(t)，即函数F的时间导数。此外，F是1周期的。因此，在区间0，1上积分ODE应该再次得到起始位置。然而，像scipy.integrate的solve_ivp这样的方法根本不考虑这种周
浏览 19提问于2021-01-22得票数 0
回答已采纳
1回答
抛射运动的奇异值
、、
在给定初始条件的情况下，我绘制了抛射运动的曲线图。
我坚持的是让tmax依赖于初始条件，所以如果你改变初始条件，tmax就会变成一个合理的值。其次，我如何找到撞击的距离、飞行时间和撞击时的速度。这是我的代码。谢谢。
浏览 24提问于2019-10-28得票数 0
2回答
空间离散化是在方程部分进行硬编码的，时间演化用一个由der(phi)组成的项来实现。
当使用变步长算法时，Dymola的时间积分总是稳定的吗？如果没有，我能做点什么吗？Dymola的Euler积分算法是显式的还是隐式的？
浏览 3提问于2014-03-24得票数 2
回答已采纳
1回答
具有可变边界的超数组集成
、、、
使用scipy.integrate.romberg的速度要快得多，几乎是瞬间的，而四边形是慢的，因为你必须循环它或矢量化它。我的函数很复杂，但是为了证明目的，假设我想从a到b集成x^2，但是x是一个计算x的标量数组。# Fails since the function doesn't seem to我认为唯一的方法是在numpy中重新实现算法本身，然后使用
浏览 7提问于2013-12-05得票数 1
1回答
我正在尝试在simulink中以矩阵形式求解一组微分方程。我的问题是，我的两个积分器的输出应该是每个积分器的"x“和"y”，然后我需要将这些x和y值放入一个矩阵中以求解它。所以我在想，如果我尝试做的事情是不可能的，或者我有错误的方法来求解方程。我想解决的是：
[dx
浏览 3提问于2015-01-29得票数 0
1回答
我试图在Python中实现Galerkin不连续方法来求解偏微分方程。这种方法涉及到大量的常微分方程组的形式。
其中$i=-1，1美元。为了进行这些积分，我定义了函数DotProduct(G，H)如下
from scipy.integrate
浏览 4提问于2022-02-10得票数 0
我对matlab ODE函数的作用感到困惑。如果我有系统动力学，为什么我不能自己手动集成它并得到结果。ode45函数的确切用途是什么？
浏览 6提问于2020-02-16得票数 0
我正在尝试使用Python的solve_ivp来解决一首颂歌。但是，我希望根据当前解决方案和以前的解决方案之间的比较，动态地更改我的ODE的右侧。这背后的想法是，我的右手边是一个向量场，我想通过根据前一个解决方案的方向反转右手边来确保向量场的方向性。
实现方法如下:我想检查右手边函数定义中前一个解决方案和向量场之间的点积。因此，我需要访问ODE求解器</e
浏览 44提问于2019-03-28得票数 0
我在matlab中用ode45和ode15s写了一个解决系统的代码。我想知道我是否可以在我的脚本中使用多核(或并行代码)来提高代码的速度。谢谢
浏览 3提问于2013-05-14得票数 5
回答已采纳
1回答
编辑t_list = np.linspace(0, 30, 100)[1. 1.端编辑yi = [1.0, 1.0, 1.0]X = ret[:, 0]它印的是nan nan nan n
浏览 10提问于2021-09-19得票数 1
回答已采纳
我有一个Simulink xPC目标应用程序，它包含几个不同采样率下具有离散状态的块，以及一些使用连续状态的部分。我保持连续状态的目的是为了更好的数值积分。产生问题的原因:一个块正在以非常快的速率(500 hz)读取设备。应用程序的其余部分可以而且应该以较慢的速率(例如，25或50 Hz)运行，因为以最高速率运行它将是过度杀伤力，而且因为处理器根本无法将整个应用程序周期压缩到较快速率的.002秒中。
浏览 4提问于2008-09-04得票数 2
我想在任意大的稀疏2D网格上使用不同的局部规则来处理离散的和连续的细胞系统。在这个框架下，简单的像素图形编辑器(在Windows 95中认为mspaint )与交互活动相结合也是可能的。以下是我的另一个梦想规范：
计算规则可以用C++ (或某些着色语言、UCDA、OpenCL等)编写，支持机器整数和浮动，并在任何可用处理器(GPU、CPU、FPGA、计算机集群)上执行--如果一切都在本地运行如果规则可以用某种声明性/脚本语言
浏览 0提问于2020-12-16得票数 -1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
云直播
对象存储
实时音视频活动推荐
腾讯云开发者
