极坐标数据点的Python数值积分

极坐标数据点的Python数值积分,python,numerical-integration,Python,Numerical Integration,我想使用scipy在3D中对一组数据点进行数值积分。其中一个积分是 d(phi)f(x)*cos(phi)从0到2.pi,其中f(x)和x的值可以表示为两个向量的大小K和L和phi(这两个向量之间的角度),即x=sqrt(K^2+L^2+2KLcos(phi))和K将在另一个积分中使用,而我将保持不变 到目前为止,我尝试的是scipy.simps: my_Phifunc = np.cos(xnew) * interpol_func area2=simps(my_Phifunc, xnew) 我

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

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

my_Phifunc = np.cos(xnew) * interpol_func
area2=simps(my_Phifunc, xnew)
我将数据点插值为数组xnew上的interpol_func


然而,这给了我一个与预期答案相去甚远的结果,我有点不确定我的想法是否正确

这个问题有点令人困惑。三维数据点?x是两个向量的大小,或者,等等,一个角度?最好把你的问题归结起来,然后编辑这篇文章。这个问题有点让人困惑。三维数据点?x是两个向量的大小,或者,等等,一个角度?最好把你的问题归结起来,然后编辑这篇文章。