我是一名工程专业的学生，我正在试图弄清楚如何使用scipy.integrate模块中的odeint函数(我只在MATLAB中使用过ode45 )。我在尝试数值求解一个简单的二阶质量，弹簧，仪表盘系统。下面是我写的代码(特别是我使用Jupyter Notebook并运行最新版本的Python 3)： import numpy as np
from scipy.integrat
浏览 35提问于2020-09-12得票数 0
回答已采纳
在某些库中，像"scipy.integrate“这样的库调用像"odeint”这样的函数(集成函数)必须表示为"odeint(func，y0，T，...)“其中"func“是必须用两个参数(y，t0)预定义的函数的名称，y向量和t0是标量。问题是，如果已经定义的函数"func“是用两个参数指定的，但顺序与"
浏览 1提问于2017-04-25得票数 1
1回答
编辑:这个问题是由于一个错误造成的，修正在中。
在开发和测试代码时，我可能会犯一个简单的错误，比如一个NameError。当我使用scipy.integrate.odeint时，odeint会打印错误消息，但是无论我请求多少时间步骤，都要继续集成，所以我得到了许多相同的错误消息。我假设它具有这种行为，以便在算术错误发生时继续进行(例如，除以零)，
浏览 1提问于2013-05-07得票数 5
回答已采纳
3回答
import numpy as npimport mathB=np.dot(b,b.T) return math.exp(np.dot(y,B))y= odeint(g, y0, t, args=(B,))
我试着用下面的代码解决y的动态系统。当它
浏览 0提问于2018-04-05得票数 2
我研究如何用odeint求解微分方程，为了便于测试，我尝试求解以下ODE (下面的示例来自)：import numpy asnp def model(y,t):y0 = 5.0
t = np.linspace(0
浏览 9提问于2019-09-22得票数 0
odeint与以下设置工作良好；import matplotlib.pyplot as pltplt.plot(time, x[:,0], time, x[:,1])
但是，当我想定义一个参数化的模型，然后将参数传递给model()函数时，当odeint调用它时，我会遇到这样的错误
浏览 3提问于2016-03-29得票数 0
回答已采纳
1回答
我正在使用Python包来求解ODE方程。但是，我需要知道在python中的中间状态，或者换句话说，ode求解器的轨迹。from scipy.integrate import odeint在这里，输出只是给我最后的状态，而不是中间的步骤，我如何才能得到中间的步骤？谢谢 solution = ode
浏览 7提问于2022-11-15得票数 1
回答已采纳
我有一首二阶颂歌：其中a = 5、b = 10和u''(t)是一个函数： if t < 5:30*np.sin(t*np.pi/5) -20*np.sin(np.pi*(t-5)/6)t = np.linspace(0, 11x, v = S -
浏览 12提问于2022-03-10得票数 0
这是一个相当普遍的问题，其中一部分可能适用于任何耦合常微分方程的数值模拟，一部分可能仅适用于scipy.integrate Python库中的odeint方法。首先，odeint如何使用手动输入的雅可比矩阵( Dfunc参数)，为什么它会使大型ODE系统的速度如此之快？其次，更确切地说，对于我的具体问题，如果雅可比函数略有错误，odeint会产生
浏览 31提问于2020-03-11得票数 1
我试图解n(多元指数)的任意值的Lane-Emden方程。为了使用SciPy，我将二阶ODE表示为由两个耦合的一阶ODE组成的集合。我有以下代码：import numpy as npapp.launch_new_ins
浏览 1提问于2017-02-22得票数 2
回答已采纳
1回答
这个问题是关于scipy.integrate模块中的odeint函数的。所有这些都与odeint (scipy.integrate模块)集成在一起，它工作得很好。我
浏览 1提问于2012-05-17得票数 2
我试图用scipy.odeint()方法来求解二阶偏导数函数。但是我想对k的许多值尝试这个解。为此，我在列表k中包含了我想要的值，并通过一个循环将这些值作为参数插入到最终解决方案中。但是，我收到了一个错误import numpy as np
from scipy.int
浏览 3提问于2020-10-17得票数 1
回答已采纳
1回答
lambda x,y: (x)/((x**2) + (y**2))**(3.0/2.0),0,1) listx.append(f(x))问题是，要解决这个方程500次需要一个明显的长时间加上，它没有给我足够数量的绘图点。所以，如果我把它提升到1000，它需要的时间太长，而且在某些情况下，我不得不杀死蟒蛇，因为它没
浏览 4提问于2014-02-21得票数 1
回答已采纳
我正在尝试使用scipy中的odeint来集成代码中定义为f的函数。函数f以θ，t和K为参数，它们定义在函数f的下面。y是结果，我得到了错误。错误的原因是2维的theta。我无法执行集成。import numpy as npfrom scipy.integrate import odeint
f是要集成的</
浏览 0提问于2012-10-26得票数 1
回答已采纳
我对python和scipy比较陌生，因为我是MATLAB的转换者。我在对scipy.integrate中的odeint函数进行快速测试时，遇到了这个潜在的bug。考虑以下代码片段：from scipy.integrate import odeint
from scipy.interpolate import interp1dbounds_error=False, fill_value=0)
浏览 0提问于2011-04-13得票数 4
回答已采纳
我试图使用Kronecker函数来绘制来自ODE的输出，该函数只能在特定的时间= t1时变为“active”。这应该会给出锯齿状的响应，其中初始值呈指数衰减，直到t=t1再次上升，然后再次下降。然而，当我绘制这个图时，它看起来就像求解器将Kronecker函数看作是所有时间t的零。在Python中是否有这样做呢？from scipy import KroneckerDelta
import
浏览 3提问于2019-10-10得票数 1
回答已采纳
1回答
我正在努力理解scipy.odeint是如何工作的，但我有一些问题。例如：from scipy.integrate import odeint该代码应计算如下形式的微分方程组：
哪里
xdot_0 = -(a_0
浏览 0提问于2018-09-11得票数 0
回答已采纳
1回答
我正在尝试用scipy.integrate.odeint从最终时间的边界条件求解一个常微分方程系统，并一直工作到初始时间(如这里所描述的：Backward integration in time using然而，odeint迭代为负的时间值-超出了我正在寻找的实际解决方案的范围-这会导致错误，因为我的ODE依赖于时间的平方根，并且我的函数返回复数值而不是实数。这里
浏览 58提问于2019-06-18得票数 1
2回答
我写了一个用odeint解微分方程的程序。但它有一个问题。当我将Cosmopara设置为np.array([70.0,0.3,0,-1.0,0])时，它给出了一个警告：invalid value encountered in sqrt和invalid value但是我检查了这一行，没有发现任何错误。如果Cosmopara = np.array([70.0,0.3,0.0,-1.0,0.0])，Y不应该改变，但它改变了。此外，如果我选择Cosmopara = np.array([70.0,0.3,0.1,-1.0,0.1])
浏览 1提问于2015-04-25得票数 1
1回答
我试着把a=1的F-L方程解成a= 0。然而，当我试图用odeint求解它时，当a和t的值很小时，我得到的导数的平方根是负值，弄得一团糟。是我使用了错误的方法来解决它，还是我的方程式中遗漏了什么？import numpy as npimport matplotlib.pyplot as plt
浏览 6提问于2021-06-15得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
