有没有办法让scipy.integrate.quad接受args中的数组？
scipy.integrate.quad是一个用于数值积分的函数，它的参数args用于传递给被积函数的额外参数。然而，args只能接受单个值，而不能接受数组。
如果想让scipy.integrate.quad接受args中的数组，可以使用lambda函数来包装被积函数，将args中的数组作为输入参数传递给lambda函数。下面是一个示例代码：
import numpy as np
from scipy import integrate
def integrand(x, a, b):
return a * x + b
a = np.array([1, 2, 3])
b = np.array([4, 5, 6])
result = integrate.quad(lambda x: integrand(x, a, b), 0, 1)
print(result)在这个示例中，我们定义了一个被积函数integrand，它接受两个参数a和b，并返回a * x + b的结果。然后，我们将a和b定义为数组，并使用lambda函数将它们作为参数传递给integrand。最后，我们调用scipy.integrate.quad函数进行数值积分。
需要注意的是，由于args参数无法接受数组，我们需要使用lambda函数来包装被积函数，并将数组作为参数传递给lambda函数。这样做可以实现让scipy.integrate.quad接受args中的数组。
相关·内容
2回答
集成Python中返回数组的函数
python、numpy、scipy
我有很多数据要集成，并希望找到一种方法来完成所有这一切只是矩阵，并愿意妥协的准确性，以提高性能。我的想法是这样的：
import numpy
import scipy
a = np.array([1,2,3])
def func(x):
return x**2 + x
def func2(x):
global a
return a*x
def integrand(x):
return func(x)*func2(x)
integrated = quad(integrand, 0, 1)
因此，我试图集成来自integrand的数组中的每个元素。
我知道
浏览 1提问于2015-12-18得票数 10
1回答
具有可变边界的超数组集成
python、numpy、scipy、integrate
我试图在一个点的列表上集成一个函数，并将整个数组传递给一个集成函数，以便将事情向量化。首先，打电话给scipy.integrate.quad太慢了，因为我有大约10000点要集成。使用scipy.integrate.romberg的速度要快得多，几乎是瞬间的，而四边形是慢的，因为你必须循环它或矢量化它。
我的函数很复杂，但是为了证明目的，假设我想从a到b集成x^2，但是x是一个计算x的标量数组。
import numpy as np
from scipy.integrate import quad, romberg
def integrand(x, y):
return x
浏览 7提问于2013-12-05得票数 1
1回答
在值数组上对两个变量的函数进行积分
python、arrays、integration
我目前正在尝试使用SciPy来解决这个积分问题： ? 我首先被建议使用插值，我尝试过，但由于某些原因无法解决，但可能是一个很好的方法。我找到了关于使用np.vectorize的this post，我认为它仍然可以工作，但我得到了一个错误。下面是我到目前为止编写的代码(还要注意，n和n，eq不是索引，它们只是变量名)： import numpy as np
from scipy import integrate
def K(x): #This is a function in the integral.
b = 0.252
return b*(((4/(x**3))
浏览 18提问于2021-10-17得票数 0
回答已采纳
1回答
如何将args()传递给integrate.quad？
python、lambda、scipy、parameter-passing
我在做一个数值积分，而不是用
scipy.integrate.dblquad
我尝试用以下方法直接建立双积分
scipy.integrate.quad
来自camz的回答：
a=1
b=2
def integrand(theta, t):
return sin(t)*sin(theta);
def theta_integral(t):
# Note scipy will pass args as the second argument
return integrate.quad(integrand, b-a, t , args=(t))[0]
integrate
浏览 5提问于2022-08-30得票数 0
回答已采纳
2回答
如何用Python的scipy.integrate.quad计算多元函数的单积分？
python、scipy
我正在尝试使用scipy.integrate.quad将一个函数集成到Python中。这个特殊的函数有两个参数。只有一个论点我想要积分。下面显示了一个示例。
from scipy import integrate as integrate
def f(x,a): #a is a parameter, x is the variable I want to integrate over
return a*x
result = integrate.quad(f,0,1)
这个示例不起作用(您可能很清楚)，因为Python在我尝试它时提醒我：
TypeError: f() takes e
浏览 2提问于2015-01-23得票数 14
回答已采纳
1回答
Wolfram和scipy.integrate.quad对同一个积分给出了不同的答案
python、math、numpy、scipy、wolframalpha
考虑以下职能：
import numpy as np
from scipy.special import erf
def my_func(x):
return np.exp(x ** 2) * (1 + erf(x))
当我使用-14的quad函数计算这个函数从scipy到-4的积分时，得到以下结果：
In [3]: from scipy import integrate
In [4]: integrate.quad(my_func, -14, -4)
/usr/local/lib/python2.7/dist-packages/scipy/integrate/quadpack.p
浏览 2提问于2015-08-30得票数 4
回答已采纳
4回答
如何使用数组输入调用python scipy quad
python、numpy、scipy、integrate、quad
我正在使用嵌套的scipy.integrate.quad调用来集成一个二维被积函数。被积函数是由numpy函数组成的，因此传递给它一个输入数组要比循环输入并为每个输入调用一次要高效得多-由于numpy的数组，它快了大约2个数量级。
然而……如果我只想在一个维度上积分被积函数--但是在另一个维度上有一组输入，事情就会失败--看起来'scipy‘quadpack包不能做numpy所做的处理阵列输入的事情。有没有其他人看到过这个-或者找到了修复它的方法-或者我误解了它。我从quad得到的错误是：
Traceback (most recent call last):
File "C
浏览 3提问于2012-08-10得票数 3
回答已采纳
1回答
numpy中向量化函数的性能行为
python、numpy、scipy
我想通过以下方式在python中执行数学积分：
1借助scipy.optimize.fsolve求解隐式方程，求被积函数的最大值
2在scipy.integrate.quad的帮助下，将被积函数的最大值移到零，并从-10到10进行积分
由于被积函数有一个自由参数xi，因此我想在numpy的帮助下使用一系列xi值来执行此操作，因此我使用numpy.vectorize。
现在有两种方法来向量化这个算法：
A分别向量化1和2，并将vec_1的结果作为输入提供给vec_2
B向量化一次执行1和2的函数
我注意到A比B快得多。下面是代码：
from scipy import optimize, inte
浏览 4提问于2013-01-21得票数 2
1回答
集特殊功能与绘图于一体
python、matplotlib、plot、integrate、bessel-functions
我正试图设计一个特别的整体(如。函数和我的最小代码如下所示。
#!/usr/bin/env python
import matplotlib.pyplot as plt
import numpy as np
import scipy.integrate as integrate
import scipy.special as sp
from scipy.special import jn
#x = np.arange(0.0, 10.0, 0.1)
U = np.linspace(0,10,1000)
#Delta = U**2
#Delta = U-4+8*integrate.qu
浏览 5提问于2017-07-09得票数 0
2回答
通过集成创建一个新的lambda函数
python、python-3.x、numpy、lambda、scipy
我使用lambda函数在Python中定义了一个现有函数。该函数表示概率分布的PDF。我想构造另一个lambda函数-一行，它代表CDF。
我不希望使用带有def关键字的单独函数定义。
下面是我一直使用的代码部分：
import numpy as np
import scipy.integrate as integrate
#define range
dx=0.01
X = np.arange(0,13,dx)
#define a piecewise function for the spline
ul = 1.0
f_pdf = lambda x: np.piecewise(x, [x
浏览 0提问于2017-11-08得票数 1
回答已采纳
1回答
Sympy:将符号表达式用作数值积分
python、sympy
我需要象征性地操作一个函数，然后对函数进行数值积分。如何在integrand函数中正确使用表达式f。如果这样做是明智的，我该如何正确地使用lambdify？非常感谢。
from sympy import *
import scipy.integrate as integrate
r = symbols('r') #define symbol
f = diff(r*r) #carry out symbolic manipulation
def integrand(x):
浏览 1提问于2016-10-22得票数 2
回答已采纳
2回答
如何实现从Numba调用的C与nquad的有效集成？
python、numpy、scipy、numba
我需要在python的6D中做一个数值积分。由于scipy.integrate.nquad函数的速度很慢，我目前正试图通过使用Numba将integrand定义为一个scipy.LowLevelCallable来加快速度。
通过复制给定的scipy.integrate.quad示例，我能够在1D使用完成此操作。
import numpy as np
from numba import cfunc
from scipy import integrate
def integrand(t):
return np.exp(-t) / t**2
nb_integrand = cfunc(&#
浏览 3提问于2017-08-22得票数 8
回答已采纳
1回答
振荡函数的积分在python中不收敛
python、scipy
我想用scipy.integrate.quad计算函数的积分。以下是integrand的外观：
我们可以注意到，integrand的大部分贡献将来自.1到4或5 ish。对于x= 10或更多，虽然函数是振荡的(在图片上很难分辨)，但它非常小，而且一直在变小。
下面是从0到某些上限的集成结果。积分的上界在x轴上。
在这里，虽然前一百倍的结果似乎是稳定的，但之后就不再是这样了，即使我期待一条直线.
我对python非常陌生，我不知道什么是最好的行动方针。我现在最好的猜测是，取一些小于100的值作为积分的上界，并放弃其他值，因为它只是integrate.quad的一个糟糕的收敛。
编辑：绘制了第
浏览 2提问于2017-07-18得票数 2
回答已采纳
1回答
使用SciPy集成一个本身涉及积分的函数
python、scipy、integral
所以我要做的是计算y上的积分，然后计算exp(-t*的解)，并把它积分到x上。
应该是这样的：
X (exp(-t* B))从0到Pi的积分
y(3.0*(sin(x)*sin(y)*sin(TM)+cos(x)*cos(TM))**2.0-1.0)**2.0上的B=积分从0到2Pi
我试着用but来做它，但是它不会在没有x的情况下对y进行集成。
到目前为止，这是我的代码：
from numpy import cos, sin, exp
import math
import scipy.integrate as integrate
t=0.0
TM=(54.74/180)*math.pi
浏览 1提问于2016-05-01得票数 2
回答已采纳
2回答
python中的双重积分
python、scipy
我想用python做一个二重积分。在的帮助下，我做了以下工作
import numpy as np
import matplotlib.pyplot as plt
from scipy import integrate
def func(x,a):
return a*x*x
def integrl(a):
return integrate.quad(func,0,10,args=(a))
print integrate.quad(lambda x: integrl(x),0,1)
但是我得到了以下错误
Traceback (most recent call last):
浏览 0提问于2012-10-14得票数 0
回答已采纳
2回答
python的quad函数中的参数"points“
python、scipy、integral
我想用Python的scipy包中的scipy.integrate.quad函数做一个主积分。例如，我想要计算(很抱歉使用latex代码)：
\begin{equation}
\lim_{\epsilon\rightarrow 0}
\left(\int_{-1}^{-\epsilon}dx\frac{1}{x}+\int_\epsilon^1 dx\frac{1}{x}\right)
\end{equation}
或
我从手册中发现，有一个可选的参数quad，points，它是有界积分区间中的一系列断点，在那里可能会出现被积函数的局部困难。points可以用来计算上面的表达式吗？如
浏览 2提问于2017-02-15得票数 0
1回答
黑体辐射与巨蟒的整合
python、physics、astronomy
我在试图绘制黑体辐射的整合图。但出于某种原因，图上有个小飞艇。我也收到警告了。
IntegrationWarning:算法不收敛。检测到舍入错误。
我认为这可能是高估了整合。但我不知道如何解决这个问题。
import scipy.integrate as integrate
import numpy as np #math module
import matplotlib #Plotting module
import matplotlib.pyplot as plt
#constants
h = 6.626e-34
浏览 4提问于2021-02-08得票数 0
1回答
有没有办法在jit_integrand_function和scipy.integrate.quad中使用scipy.optimize.fsolve？
scipy、jit、numba、quad
基于这里提供的解释，1，我正在尝试使用相同的想法来加速以下积分： import scipy.integrate as si
from scipy.optimize import root, fsolve
import numba
from numba import cfunc
from numba.types import intc, CPointer, float64
from scipy import LowLevelCallable
def integrand(t, *args):
a = args[0]
c = fsolve(lambda x: a * x**2 -
浏览 20提问于2020-07-16得票数 3
2回答
大数上的scipy.integrate.quad精度
python、numpy、scipy、calculus
我试图通过scipy.integrate.quad()计算这样一个积分(实际上是指数分布的cdf及其pdf)。
import numpy as np
from scipy.integrate import quad
def g(x):
return .5 * np.exp(-.5 * x)
print quad(g, a=0., b=np.inf)
print quad(g, a=0., b=10**6)
print quad(g, a=0., b=10**5)
print quad(g, a=0., b=10**4)
结果如下：
(1.0, 3.5807346295637055e
浏览 1提问于2016-09-15得票数 4
回答已采纳
2回答
用scipy积分方程
python、scipy
from scipy.integrate import quad
def integrand(a, b):
return a * x ** 2 + b
a = 2
b = 1
I = quad(integrand, 0, 1, args=(a,b))
I
这是我的程序。当我尝试运行它时，它显示错误：
integrand () takes 2 positional arguments but 3 were given ....
我不明白为什么当只有两个变量，即a和b时，它会要求3个参数。
有谁可以帮我？有人能澄清我的疑虑吗？
浏览 0提问于2018-08-22得票数 5
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
