python - 极坐标中数据点的Python数值积分

Question

我想使用 scipy 对一组 3D 数据点进行数值积分。其中一个积分 d(phi)f(x)*cos(phi)来自0 to 2.pi，其中我有 f(x) 的值，x 可以表示为 2 个向量K和L的大小以及 phi（这两个向量之间的角度）x = sqrt(K^2 + L^2 +2KLcos(phi))，即 K 将用于另一个积分，而 L 将保持不变

到目前为止我尝试的是 scipy.simps：

my_Phifunc = np.cos(xnew) * interpol_func
area2=simps(my_Phifunc, xnew)

我将我的数据点作为 interpol_func 插值在数组 xnew 上。

但是，这给我的结果与预期的答案相去甚远，我有点不确定我的想法是否正确

标签： pythonnumerical-integration