python 二重积分 辛普森法_函数二重积分与采样数据Python

我正在寻找一种使用numpytrapz或scipy堆栈中的类似函数对采样数据执行双重积分的方法。在

我特别想计算一下函数：

其中f(x',y')是采样数组，F(x, y)是相同维度的数组。在

这是我的尝试，结果不正确：def integrate_2D(f, x, y):

def integral(f, x, y, x0, y0):

F_i = np.trapz(np.trapz(np.arcsinh(1/np.sqrt((x-x0+0.01)**2+(y-y0+0.01)**2)) * f, x), y)

return F_i

sigma = 1.0

F = [[integral(f, x, y, x0, y0) for x0 in x] for y0 in y]

return F

xlin = np.linspace(0, 10, 100)

ylin = np.linspace(0, 10, 100)

X,Y = np.meshgrid(xlin, ylin)

f = 1.0 * np.exp(-((X - 8.)**2 + (Y - 8)**2)/(8.0))

f += 0.5 * np.exp(-((X - 1)**2 + (Y - 9)**2)/(10.0))

F = integrate_2D(f, xlin, ylin)

输出数组似乎朝向生成网格的对角线，而它应该返回一个看起来像模糊输入数组的数组。在