0

我有以下三重积分:

积分

数据为:t 一维数组(大小 124)、pr 一维数组(大小 10)、lat 一维数组(241)和 lon 一维数组(大小 480) V 四维数组(124、10、241、480)

我想在每个 lat 上整合 v over t、pr 和 lon 所以,我使用了以下代码:

def M(T, lam, P, V, phi):

    return integrate.tplquad(V*R*np.cos(phi), 0, T,
                         lambda T: 0, lambda T: lam,
                         lambda T,lam: 0, lambda T,lam: P)[0]

for i in range(lat.shape[0]):
    a = M(t, lon, pr, v[:, :, i, :], lat[i])

但我得到了错误:

ValueError:具有多个元素的数组的真值不明确。使用 a.any() 或 a.all()

tplquad 的工作原理仍然让我感到困惑。对使用 tplquad 或其他功能进行积分有任何帮助吗?

4

0 回答 0