我想用Python模拟带电粒子通过磁场时的运动。我尝试使用scipy.integrate中的odeint函数，但它似乎没有像我预期的那样工作。以下是我在给定初始条件的情况下所期望的结果：但这是我通过模拟得到的结果：问题来自于我对odeint函数的实现。任何帮助都会被删除。下面是我的代码：
from scipy.integrate import odei
浏览 21提问于2018-03-03得票数 1
1回答
当我使用scipy.integrate.odeint时，odeint会打印错误消息，但是无论我请求多少时间步骤，都要继续集成，所以我得到了许多相同的错误消息。
浏览 1提问于2013-05-07得票数 5
回答已采纳
2回答
我是这里的新手，也是Python的新手。我正在为常微分方程组的数值解编写一个程序。我的问题是:我需要在函数中实现一个循环。之后，我需要集成这个函数(odeint)，但是它不能工作。下面是不带循环的(工作)代码：from pylab import * # for plotting commands a = 2.0
return array([ y*a*(1-y/10)]
浏览 2提问于2014-03-02得票数 1
回答已采纳
我目前正在学习一些拉格朗日力学，这在很大程度上进行得很好，但我决定尝试一个计算机练习，我不知道如何让Python求解这些方程。求解一个普通的一阶方程很容易，因为你只需创建一个函数，将它设置为某个值，然后使用odeint。没问题。但是现在我得到了一个二阶方程，我不知道如何告诉odeint或函数，我要找的是两次微分的东西。我尝试使用不同数量的积分和微分函数(符号函数和其他函数)，并说了如下内容：y‘= RHS
浏览 26提问于2019-08-05得票数 2
我想用python数值求解ODE，比如y'=f(x) (有边界条件y(0)=0)。我不知道函数f(x)的解析表达式是什么，相反，对于我想要解微分方程的域，我有一组函数的点(数据)。from scipy.int
浏览 2提问于2019-04-08得票数 1
回答已采纳
我正在尝试最小化两个有界函数之间的均方误差，curve_fit做得很好，但当两个函数之间的均方误差小于0.1时，我想停止计算。beta * y[0] * y[1] / N inf = -(sus + rec) return integrate.odeint(sir_model, (sus0, inf0, rec0), x, args=(b
浏览 65提问于2021-07-06得票数 1
我目前正在python下解决一个微分方程系统，使用来模拟带电粒子在一个场中(源代码来自这个)： returnnp.array([x[1], eta*Ez0(x[0])])
它工作得很好，但我想在x< 0时立即停止计算。: else:
return np.array([x[1], eta*Ez0(x[0
浏览 0提问于2010-11-20得票数 0
回答已采纳
4回答
由于odeint使用函数f作为参数，所以我将此函数编写为静态类型的Cython，并希望odeint的运行时间将显著减少。函数f包含在名为ode.pyx的文件中，如下所示：cimport numpy as np derivs[1] = -omega/Q + sin(theta) + d*cos(Omega*t)然后创建一个文件s
浏览 16提问于2017-03-16得票数 7
回答已采纳
1回答
本文给出的0波函数及其导数的值分别为0,1。这允许我们使用scipy.integrate.odeint函数来解决给定能量值的问题。现在的任务是找到能量本征值，给出进一步的边界条件，即L上的波函数是0，在python上使用根查找函数。我已经做了一些研究，只能找到一种叫“射击方法”的东西，我不知道如何实现。另外，我遇到了解BVP的枕函数，但是我似乎不明白这个函数的第二个输入(边界条件残差)到底是什么。*m_el*energy*y[0]&#
浏览 5提问于2019-12-17得票数 3
典型的纯python (C，Fortran)，如果它由简单的仿制品组成的紧密循环组成。例如，在这个中使用所描述的时，只需编写函数就可以在纯python中集成，如下所示： Si = y[0] Ri =它是否使用某种方法将函数f(y,t) 从python转换为本机代码？(比如f2py、scipy.weave、cython .)据我所知，odeint不需要任何C/C++或Fortran编译器，也
浏览 3提问于2014-03-24得票数 1
回答已采纳
1回答
= -(24/37)*np.cos(3*t)- (4/37)*np.sin(3*t) + (y0+24/37)*np.exp(0.5*t)
# plt.plot(t,y_ODEint)
当前的问题是，像<em
浏览 2提问于2022-03-08得票数 1
回答已采纳
当我运行代码时from scipy.integrate import odeint return y_out
X0 = [0,0,0.33206][
浏览 0提问于2014-09-20得票数 3
回答已采纳
1回答
scipy odeint)。然而，odeint迭代为负的时间值-超出了我正在寻找的实际解决方案的范围-这会导致错误，因为我的ODE依赖于时间的平方根，并且我的函数返回复数值而不是实数。有没有办法防止odeint变为负值？为什么它要尝试输入数组t之外的值 对于我的实际代码，我找到了一些方法来避免nan结果，比如在函数中添加if t<0: t=0 (在这里不起作用)，或者通过施加非常小的最大时间步长(hmax<tmin)，请注意，Scipy odeint Non
浏览 58提问于2019-06-18得票数 1
谁能提供一个在SciPy中为integrate.odeint函数提供雅可比矩阵的示例？我尝试在SciPy tutorial 中运行这段代码，但似乎从未调用过Dfun() (雅可比函数)。from numpy import * # addedfrom scipy.special import gamma, airyjacobian' # added x =
浏览 0提问于2013-07-09得票数 10
回答已采纳
1回答
给出指标的界外误差
、、、
我试着用Scipy的odeint函数来求解python中的微分方程。方程的形式是dy/dt = w(t)，其中w(t) = w1*(1+A*sin(w2*t))对于某些参数w1、w2和A。(f,0,t)import numpy as np return w[t0]
IndexError: index 1001 is out of bou
浏览 5提问于2015-09-18得票数 2
回答已采纳
我正在尝试创建一个名为calc(f,a,b)的函数，其中x是一个变量为f的方程，我想将此代码放入该函数中。def calc(f, a, b): integral = odeint(lambda y, x : f, 0, limits)return integral[1]
此函数使用内置的odeint函数获取积分。我的问题是，这个函数(x**2)需要传递给紧跟在y, x: f之后的odeint</em
浏览 1提问于2011-03-21得票数 1
回答已采纳
我正在将Python代码从2.7转换为Python3.8。2.代码y_trajectory = scipy.integrate.odeint(growth_a_derivs, y_start, t_array, atol=eps, args=myargs)
使用gro
浏览 4提问于2021-06-04得票数 0
回答已采纳
我的ODE为Mx''+Lx'+ f(x) =0，其中f(X)是多项式函数。请看我的完整代码，其中我定义了一个函数，即'diff‘的微分方程。然后，我使用'odeint‘，它调用'diff’和必要的参数，以解决微分方程。但是当f(x)是多项式时，参数列表就会变得太长。因此，我
浏览 0提问于2017-12-15得票数 1
回答已采纳
我对python和scipy比较陌生，因为我是MATLAB的转换者。我在对scipy.integrate中的odeint函数进行快速测试时，遇到了这个潜在的bug。考虑以下代码片段：from scipy.integrate import odeint
# ODE integration with u(t) = acel, a la
浏览 0提问于2011-04-13得票数 4
回答已采纳
1回答
最优误差
、、、
我想用Python解决ODEs的一个相当简单的参数估计问题。我一直在使用odeint函数来求解ODE，使用scipy.optimize库来查找参数。当我单独使用odeint函数时，没有问题，但是通过scipy.optimize，它给了我一个错误RuntimeError:函数返回的数组必须是一维的，但是得到了ndim=2。我不认为这是一个复制与，因为我可以独立运行odeint求解器.代码如下：from scipy.integrate
浏览 0提问于2018-07-29得票数 2
回答已采纳
云服务器
ICP备案
对象存储
云直播
实时音视频
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号