scipy 相关问题

SciPy是一个用于Python编程语言的算法和数学工具的开源库。

使用 scipy.integrate.quad 积分 gamma 函数

当 a = 3/2 使用 scipy.integrate.quad 时,如何求 Gamma(x) = x^(a-1) * e^-x 从 0 到无穷大的积分? 到目前为止,我已经尝试了一些代码,但我被告知“quad:第一个参数......

回答 1 投票 0

如何使用 numba 使用数组中的值重新定义积分来加速积分

我尝试使用 python numba 更快地计算积分。尽管使用 numba 的计时速度对于单次计算来说几乎快了 10 倍,但当我循环 redefini 的过程时......

回答 1 投票 0

scipy dblquad 在简单二重积分中提供错误的结果

我正在尝试在Python中计算一个简单的双定积分:平方[0,1] x [0,1]中的函数Max(0, (4-12x) + (6-12y))。 我们可以用 Mathematica 来做并得到准确的结果...

回答 1 投票 0

复函数的二重积分

按照这里的方法, 我们可以这样做吗: ∫∫ (x+y)+i(x+y) dxdy = ∫ ∫(x+y)dx + i∫(x+y)dx dy = ∫∫(x+y)dxdy + i∫∫(x+ y)dx dy 问题是,我正在尝试计算

回答 1 投票 0

如何正确使用integrate.cumtrapz?

我尝试用 f(x)=x^2 进行梯形积分,我知道反导数的样子,所以 F(x) = (1/3)x^3 这是我的代码,就像我尝试过的那样: 像这样导入 scipy.integrate p=x**2 l=它。

回答 1 投票 0

以数组为上限对函数进行数值积分

我正在尝试对 1/H(z)dz 进行数值积分,其中函数 H(z) 由下式给出 def Hz_th(theta,z): aa,H0,omegam0=theta 返回 np.sqrt(H0**2*(1 + z)**3*omegam0 + 3*H0**2*(1 + z)**3*aa*o...

回答 1 投票 0

scipy的solve_ivp有一个关键字events。我的问题是关于这个关键字的

据我了解solve_ivp文档,events=event所需的函数event(t, y, arg)必须是这样的,即 $event(t_0, y_0, args) = 0. 如果事件已发生。 solve_ivp 然后尝试

回答 1 投票 0

是否有可能以某种方式获取这个积分?

我有以下代码,但显示错误: IntegrationWarning:已达到最大细分数 (50)。 如果增加限制没有改善,建议分析...

回答 1 投票 0

如何计算新点位于 Voronoi 图的哪个位置?

我编写了一个小脚本,用于显示本教程中 M 点的 voronoi 图。我使用 scipy.spatial。 我想给出一个新的平面点,并说出这个点在 voronoi 图的哪个位置。 ...

回答 2 投票 0

为什么我的 scipy.optimize.minimize(method="newton-cg") 函数停留在局部最大值?

我想找到依赖于2个变量的函数的局部最小值。为此,我的计划是使用 scipy.optimize.minimize 函数和“newton-cg”方法,因为我可以计算...

回答 1 投票 0

在 scipy 中将重叠数组拼接在一起

给定两个 numpy 数组(矩阵) A = np.linspace(1,9,9).reshape(3,3) B = np.linspace(10,18,9).reshape(3,3) 我们可以通过以下方式将它们组合成块对角矩阵: 从 scipy.linalg 导入块...

回答 0 投票 0

通过python生成具有一定概率密度函数的随机变量

我想生成具有麦克斯韦能量(例如f(ene))或任意分布的随机变量。 概率密度函数如下: def f(烯): 勒=3 返回 2(ene/pi/le**3)*np.e...

回答 3 投票 0

确保找到矩阵方程的平凡解

我正在尝试用 numpy 求解矩阵方程 $Ax = b$ (实际上我正在使用 scipy.sparse 但我相信问题仍然相同)。在我的设置中,$b$ 是某个函数的导数。 在我的...

回答 1 投票 0

在Python中使用输出数组计算数组的积分

我想计算以下形式的积分 我希望将结果作为数组(最终将它们绘制为欧米茄的函数)。我有 将 numpy 导入为 np 将 pylab 导入为 plt 来自...

回答 1 投票 0

如何使用 scipy.integrate 积分 $\int_0^T K_{\lambda}(t,s) \sigma(s)ds$

我需要整合表达式 使用 scipy 的集成模块。 我努力了 定义 K(t,s): 返回 np.exp(lbda*(T-s))*np.exp(lbda*t)/(1-np.exp(lbda*T)) 定义西格玛: 返回 np.sin(2*np.pi...

回答 1 投票 0

最快的双积分方法

我正在使用 scipy 的双积分 dblquad,我正在尝试提高速度。我已经检查了网上提出的解决方案,但无法使它们发挥作用。 为了缓解这个问题,我准备了......

回答 1 投票 0

在 Python 中将数据拟合到具有多个变量的互补误差函数

我无法将实验数据拟合到 Python 3.7.4 中的互补误差函数。 将 matplotlib.pyplot 导入为 plt 导入数学 将 numpy 导入为 np 从 scipy import 优化 来自 scipy

回答 2 投票 0

计算 2D 插值积分时出现错误。比较 numpy 数组

我的优化任务涉及计算以下积分并找到 xl 和 xu 的最佳值: 迭代花费的时间太长,因此我决定通过计算...的积分来加快迭代速度

回答 2 投票 0

scipy正交积分精度警告和nan差值

我正在尝试使用 scipy.integration.quadrature 集成复杂的函数组合,它会抛出准确性警告,并且“最新差异”值是(吞下)nan: 'C:\程序

回答 1 投票 0

SciPy:`scipy.special.fresnel(x[, out1, out2])`中的参数是什么?

我正在尝试使用 SciPy 绘制回旋函数。这是菲涅尔积分的语法,但我无法理解 scipy.special.fresnel(x[, out1, out...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.