12

我正在尝试在2D 区域上集成SciPy中的多变量函数。下面的Mathematica代码相当于什么?

In[1]:= F[x_, y_] := Cos[x] + Cos[y] 

In[2]:= Integrate[F[x, y], {x, -\[Pi], \[Pi]}, {y, -\[Pi], \[Pi]}]

Out[2]= 0

查看SciPy 文档,我只能找到对一维求积的支持。有没有办法在 SciPy 中进行多维积分?

4

2 回答 2

13

我认为它会像这样工作:

def func(x,y):
    return cos(x) + cos(y)

def func2(y, a, b):
    return integrate.quad(func, a, b, args=(y,))[0]

print integrate.quad(func2, -pi/2, pi/2, args=(-pi/2, pi/2))[0]

Wolfram|Alpha 同意

编辑:我刚刚发现dblquad似乎完全符合您的要求:

print integrate.dblquad(func, -pi/2, pi/2, lambda x:-pi/2, lambda x:pi/2)[0]
于 2010-03-03T02:43:16.560 回答
9

如果您想进行符号集成,请查看 sympy ( code.google.com/p/sympy ):

import sympy as s
x, y = s.symbols('x, y')
expr = s.cos(x) + s.sin(y)
expr.integrate((x, -s.pi, s.pi), (y, -s.pi, s.pi))
于 2010-03-03T20:27:55.813 回答