如何在网格网格上应用积分函数-解网

问：

为了简单起见，假设我有一个形式的函数：

Equation link

我想在 x-y 网格上评估这个函数：

x = np.linspace(0,1,50)
y = np.linspace(0,1,50)
xx, yy = np.meshgrid(x, y)

where 用作积分的上限（x_2），但我不确定最佳方法。我想我正在寻找合适的集成功能来允许我做到这一点。我知道这个例子可以很容易地通过解析求解，然后在整个网格上执行，但我的真实例子更复杂，我想学习适当的数值方法。我真的很感激任何建议。提前致谢！xx

我尝试使用 scipy.integrate 中的四边形函数，但这只允许我计算一个确定的积分，而不是在整个网格上积分。

python numpy 数学

# vectorize, to implement function over xy meshgrid
@np.vectorize

# define function
def z(x,y):
    
    # define function to be integrated
    def integrand(x):
         return x**3
        
    # integrate function
    ans,error = quad(integrand,0,x)

    # calculate full function, with integrated terms
    z = ans + y
    
    # return solution
    return(G)

# define meshgrid
x = np.linspace(0,1,50)
y = np.linspace(0,1,50)

XX, YY = np.meshgrid(T,P)

# implement function across meshgrid
z_out = z(XX,YY)

谢谢！

上一个：如何从加速度估计角度和从三维角度估计加速度？

下一个：在 scipy 或 numpy 中结合置信区间和自定义泊松？

如何在网格网格上应用积分函数

How to apply an integration function over a meshgrid

评论