问题的动机
我正在尝试在所有空间上积分函数 f(x,y,z)。
我尝试使用 scipy.integrate.tplquad 和 scipy.integrate.nquad 进行积分,但这两种方法都将积分返回为 0(当积分应该是有限的时)。这是因为,随着积分量的增加,被积函数非零的区域被采样得越来越少。积分“错过”了这个空间区域。然而, scipy.integrate.quad 似乎能够通过执行变量的更改来处理 [-无穷大,无穷大] 的积分......
问题
是否可以使用 scipy.integrate.quad 3 次来执行三重积分。我想到的代码如下所示:
x_integral = quad(f, -np.inf, np.inf)
y_integral = quad(x_integral, -np.inf, np.inf)
z_integral = quad(y_integral, -np.inf, np.inf)
其中f是函数f(x, y, z),x_integral应该从x = [-无穷大,无穷大]积分,y_integral应该从y = [-无穷大,无穷大]积分,z_integral应该从z =积分[-无穷大,无穷大]。我知道quad想要返回一个 float ,因此不喜欢在x上积分函数f(x, y, z)以返回y和z的函数(如上面代码中的x_integral = ...行正在尝试做)。有没有办法实现上面的代码?
谢谢
最佳答案
下面是一个对 quad
进行嵌套调用的示例,该调用执行积分,给出球体体积的 1/8:
import numpy as np
from scipy.integrate import quad
def fz(x, y):
return quad( lambda z:1, 0, np.sqrt(x**2+y**2) )[0]
def fy(x):
return quad( fz, 0, np.sqrt(1-x**2), args=(x, ) )[0]
def fx():
return quad( fy, 0, 1 )[0]
fx()
>>> 0.5235987755981053
4/3*np.pi/8
>>> 0.5235987755982988
关于python - 使用 scipy.integrate.quad 执行 3D 积分,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/51302638/