用scipy solve_ivp求解具有初始条件的方程

文章/答案/技术大牛

发布

1回答

python、scipy、equation、differential-equations

我很难适应scipy的solve_ivp。假设我们有一个二阶线性微分方程，例如弹簧(y'‘= -k**2*y)。条件是弹簧在位置0(时间0)速度为v0。我如何使用初始条件来解决它？y: function1(t, y, k = 10) #function with only t and y solution = solve_ivp

浏览 13提问于2019-01-13得票数 5

回答已采纳

1回答

solve_ivp (ODE求解器)的一个实现

python、scipy

这是我关于堆栈溢出的第一篇文章。因此，我在SciPy的文档中看到了这个solve_ivp求解器的例子。正在解决的问题如下：>>> hit_ground.direction = -1 >>> sol = solve_ivp(upward_cannon,

浏览 0提问于2020-01-23得票数 1

1回答

带实数的数字:函数调用与编译

python、scipy、numeric、numerical-integration

我是新来的。特别是，我要综合一个(复杂)微分方程组。我注意到scipy.integrate.solve_ivp可以做这项工作。由于我的方程依赖于某些数值参数，所以我编写了一个生成器类来建立方程和初始条件。from scipy.integrate import solve_ivp def __init__(self, parameters):numerical_function d

浏览 2提问于2019-10-31得票数 0

1回答

如何求解非线性方程的根？我要使用SciPy吗？

python、scipy

我试着找出给定时间的t。我在上面的方程式中有所有必要的变量，除了Ts。我意识到我可以对t进行集成，这使得事情变得更容易一些。从那里开始，我该如何解决Ts呢？

浏览 5提问于2019-11-02得票数 0

1回答

是否可以使用solve_ivp在每个时间步长访问解决方案数组？

python、scipy

我正在使用scipy.integrate中的solve_ivp来求解一个方程组， Xdot = A*X + B*U，其中X是大小为4*1的数组。def func(t,X): return Xdot 并使用求解器作为 solv = solve_ivp(convf, tspan,Z0,method='BDF')solution_array=s

浏览 43提问于2020-04-26得票数 0

1回答

如何使用Python (SciPy)解出一定数量的时间步长的ODE？

python、scipy、ode、differential-equations

我正在尝试使用SciPy来解决一组ODE。分配给我的任务要求我在500个时间步长内求解微分方程。如何使用SciPy实现这一点？到目前为止，我已经尝试使用scipy.integrate.solve_ivp，它给了我一个正确的解决方案，但我无法控制它运行的时间步数。t_span参数允许我配置t的初始值和终值，但我实际上对此并不感兴趣--我只对我集成了多少次感兴趣。(例如，当我使用t_span = (0, 500)运行我的

浏览 9提问于2021-05-13得票数 1

1回答

戴莫斯:如何显式调用solve_ivp方法来求解动力学方程？

scipy、openmdao

为了验证所使用的分析模型，我希望计算一些初始条件下的轨道。由于动力学方程有点复杂，我想知道是否还有其他方法来解决这些<em

浏览 1提问于2021-02-25得票数 0

回答已采纳

1回答

如何根据Python中最适合于数据集的参数来确定微分方程的未知参数？

python、scipy、curve-fitting、ode、differential-equations

我试图用python将不同的微分方程拟合到给定的数据集中。出于这个原因，我使用了use包，分别是solve_ivp函数。只要我对微分方程中包含的参数(b= 0.005)有粗略的估计，例如：from scipy.integrate import solve_ivpdydx= [-b[0] * y[0]] xspan= n

浏览 4提问于2021-02-03得票数 1

回答已采纳

1回答

odeint或solve_ivp中的初始条件

matplotlib、odeint

我是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

1回答

用scipy.odeint求解预定义函数的非线性常微分方程

python、scipy、ode

我想解一个形式为的非线性常微分方程 Theta2 = (C + j(Theta2))**-1 * (f(t) - g(Theta1) -h(Theta0)) 其中f()、g()、h()和j()是已经定义的函数Theta2和Theta1是Theta0对时间t的二阶和一阶导数。我一直在使用SciPy.odeint函数求解不带j(Theta2)项的方程，代码如下： from scipy.integrate import odeint de

浏览 17提问于2019-09-12得票数 0

1回答

如何求解不同初始条件下的微分方程

python、numpy、matplotlib、scipy、odeint

我只想通过改变初始条件来连续求解一个微分方程。我试了很多次，但没有找到合适的方法来做好这件事。有谁能分享一下这方面的情况吗。出于您的考虑，我在下面给出了可以求解微分方程的代码：import numpy as np plt.plot(x, sol[:, 0], '

浏览 0提问于2019-07-12得票数 4

回答已采纳

1回答

在solve_ivp函数的文档中，字母k是什么意思？

python、numpy、scipy、vectorization、scientific-computing

Solve_ivp是一个初值问题的求解函数。三言两句求解一个系统的初值问题。该函数对给定初值的常微分方程组进行数值积分。在solve_ivp函数文档(Sicie参考指

浏览 2提问于2020-02-18得票数 3

2回答

如何在Python语言中更新参数并传递给solve_ivp

python、scipy、ode

在我的常微分方程函数中，我需要迭代地求解一个参数的方程，直到在每个时间步收敛。我希望传递最新的参数值作为下一个时间步的初始值，这样当函数对参数进行迭代更新时，花费的时间会更少。但我想不出该怎么做。ODE函数的代码结构如下： if t==0: awhi

浏览 2提问于2020-08-08得票数 0

1回答

使用scipy的solve_ivp在时间步之间运行代码

python、scipy

我正在将我的代码从使用scipy的odeint转换为使用scipy的solve_ivp。solver.successful() : # do other operations 这个方法允许我在每个时间步之后存储依赖于解决方案的值我现在切换到使用solve_ivp，但不确定如何使用solve_ivp求解器来实现此功能。有没有人<

浏览 3提问于2019-10-18得票数 1

1回答

用Python耦合微分方程

python、physics、differential-equations、astronomy

我正试图用python来求解大地测量轨道方程系统。它们是耦合的普通方程。我尝试过不同的方法，但它们都给了我一个错误的形状(在绘制r和phi时，形状应该是一些周期函数)。知道怎么做吗？这是我的常量c = 0.0020053761Lz= 0.000024 #Angular m

浏览 0提问于2020-03-21得票数 5

回答已采纳

1回答

当first_step参数足够大时，solve_ivp中的参数不能工作

python、numpy、scipy、ode

我遇到了一个很奇怪的例子，我在first_step中测试了ode求解器solve_ivp中的scipy参数，当它跨越一个成为2(即first_step=2)时，它只是被求解者忽略了。所以我测试的代码是当我用first_step=1.9调用求解

浏览 4提问于2022-08-28得票数 0

回答已采纳

1回答

TypeError: fun()缺少1个必需的位置参数：'a‘

python、numpy、scipy、solver

我知道有很多关于这个错误的主题，我已经在其中的许多主题上试图理解这个简单的系统是怎么回事。这是我的代码，求解一个非常简单的方程来测试solve_ivp和odeint的效率。import numpy as np import matplotlib.pyplot as plt from scipy.integrate import solve_ivpTraceb

浏览 31提问于2020-04-25得票数 0

回答已采纳

1回答

向量微分方程

python、matrix、vector、ode、differential-equations

我试着用Python解决一个向量微分方程，但我总是得到一个错误，我不是很理解。下面是我的代码：from scipy.integrate import solve_ivp sol=solve_

浏览 4提问于2020-05-15得票数 1

回答已采纳

1回答

在Matlab中用事件函数求解常微分方程

matlab、ode

我现在正在做一个项目，我已经对多体机械系统的运动进行了建模。这是一个复杂的、非光滑的系统，系统的运动并不总是由同一组方程描述，这意味着如果我用数值方法解决这个问题，我必须有时停止积分(并用新的初始条件和方程重新启动它)。尽管如此，我还是写了一些代码，现在我只是停留在如何将所有这些放在一起:我希望常微分方程求解器(例如ode45)求解我的方程，然后

浏览 9提问于2019-02-21得票数 0

回答已采纳

1回答

解渐近微分方程组的语法

sympy

我浏览了堆栈交换中关于使用渐近符号求解具有初始条件的微分方程系统的文档和问题。具有初始条件V(0) = V0和I(0) = I0 我在堆栈交换中浏览了

浏览 1提问于2014-10-03得票数 6

点击加载更多