我有以下三重积分:
数据为: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 或其他功能进行积分有任何帮助吗?