scipy 相关问题

如何使用数组输入调用 python scipyquad

我正在使用嵌套的 scipy.integrate.quad 调用来积分二维被积函数。被积函数由 numpy 函数组成 - 因此向它传递一个输入数组比......更有效

Python+Scipy+Integration：处理带有尖峰的函数中的精度误差

我正在尝试使用 scipy.integrate.quad 在非常大的范围（0..10,000）内集成函数。该函数在其大部分范围内为零，但在非常小的范围内有一个尖峰（例如 1,602.....

scipy.linprog 可行性的错误？ （A_ub @ x0 <= b_ub).all() is True ---but--- linprog(np.zeros_like(x0), A_ub=A_ub, b_ub=b_ub) infeasible

使用 numpy、scipy 版本 numpy 1.25.0 scipy 1.11.0 以下 scipy.optimize.linprog 调用， 将 numpy 导入为 np 从 scipy.optimize 导入 linprog A_ub = np.array( [[-0.15729144, 0.

处理odeint中极小/极大的值

我正在尝试使用 scipy.integrate.odeint 求解这组耦合微分方程。 我写了这段Python代码： def odes(x, r): #常量 阿尔法 = 1.473 贝塔 = 52.46 伽玛...

scipy 四边形似乎提供不准确的积分值

我对 scipy 中的integrate.quad 即使对于简单的积分也返回不准确的结果感到困惑。下面我将举一个简单的例子，很容易就能跨越...

存储数组中的积分值，然后使用该新数组

我是Python新手，所以我真的很挣扎。我想定义一个函数，对不同值的数组进行一定的计算，将这些新计算的值存储在一个新的

Python 中的演化函数

更新的问题 继我原来的帖子之后，通过使用 @Attack68 的代码，我创建了一个程序，该程序通过选择乘法函数基成功地演化了该函数......

以常量变量为输入的积分和绘图函数

我正在做一个python项目，我需要计算摩托艇的速度及其位置。 速度的公式是 v = v0*e^-bt/m ，位置是 p = (mv0/b)(1-e^-bt/m) 其中 b 是...

使用 nquad 进行二重积分

这里有问题。到目前为止，这是我的代码： 从 scipy import 集成 导入数学 将 numpy 导入为 np a = 0.250 s02 = 214.0 a_s = 0.0163 def 被积函数(r, R, s02, a_s, a): 返回 2.0 *...

python 四元集成似乎不准确

我是Python的新手，正在尝试对函数进行数值积分。一切似乎都有效，但我得到的结果与我在 Mathematica 中得到的结果有很大不同（我......

类型错误：使用 trapz 方法时只能将 size-1 数组转换为 Python 标量

这是我的Python代码 将 numpy 导入为 np 从 scipy.integrate 导入quad、trapz、simps 导入数学 定义 f(x)： 返回 math.exp(-(x**2/2))/math.sqrt(2*math.pi) 结果，错误 =quad(f, 0, np.inf)

使用Integrate.quad进行积分时出现无效浮动错误

我有一个清单： q1 = [0,1.9488e-06,1.5473e-05,5.1829e-05,0.00012192,0.00023633,0.00040526,0.00063862,0.00094596,0.0013365,0.0018192,0.0024025, 0.0030949,0.0039041,0.0048379,0.0059036] 我是...

output[0] = y0 ValueError：使用序列设置数组元素

我已经为此苦苦挣扎了几天。我正在尝试估计分段高斯函数的密度。谁能告诉我为什么我现在收到错误 类型错误：输出[0] = y0

Python：计算数组的积分

导入numpy 将 matplotlib.pyplot 导入为 plt 从 scipy import 集成 定义 f(x,y)： 返回 x*y + x**2 def 积分(x,y): I = 整合.quad(f, 0, x, args=(y,))[0] 返回我 绝对 ga...

有没有办法用 scipy.fft 在傅立叶空间中进行数值积分？

我有兴趣在使用 scipy 对一些数据进行 fft 处理后集成到傅立叶空间中。我一直在关注傅立叶空间中的堆栈交换后数值积分......

如何使用 scipy.integrate.nquad 并在使用 LowLevelCallable cfunc 时向其传递参数

我非常绝望地尝试实现一个集成，我想将其结果用作偏微分方程的初始数据。积分本身是标量值，需要在整个空间 $\mathbb...

scipysolve_ivp events：导数积分中的事件

将 numpy 导入为 np 从 scipy.optimize 导入 fsolve 从 scipy.integrate 导入四元组 从 scipy.integrate 导入solve_ivp def dYdW(W, Y): X, y = Y dXdW = (kprime/Fa0)*(1-X)*y*(1/(1+epsilon*X...

用于 SciPy 集成和插值的 Numba

我使用 Numba 来加速我的代码。它效果很好，可提供 2-3 倍的改进。然而，我的代码中花费的主要时间（大约 90%）是在 scipy 四边形积分和插值上......

整合混沌系统

我有3个方程如下 u'=-a*(u-v) v'= c*u-v-u*w w'=-b*w+u*v a=5.0 b=0.9 且 c=8.2 我尝试使用 scipy.integrate.odeint 解决从 t=0 到 t=10 我的初始条件是 u(0)=0, v(0)...

视线上的密度分布积分

我的问题是这样的： 我知道密度是球体半径的函数。假设密度 rho(1000) 和半径(1000) 已经通过数值计算。我想找到积分...