2D интегралы в SciPy

Question

2D интегралы в SciPy

Я пытаюсь интегрировать многовариантную функцию в SciPy по 2D области. Каков был бы эквивалент следующего кода 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?

12

python integration wolfram-mathematica scipy multidimensional-array

задан Peter Mortensen 10 December 2010 в 16:23

2 ответа

Если вы хотите сделать символьную интеграцию, посмотрите на 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))

9

ответ дан 2 December 2019 в 07:02

Другие вопросы по тегам:

python integration wolfram-mathematica scipy multidimensional-array

Похожие вопросы:

score 13 · Accepted Answer

Думаю, это будет работать примерно так:

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 согласен

edit: I just discovered dblquad that seems to do exactly what you want:

print integrate.dblquad(func, -pi/2, pi/2, lambda x:-pi/2, lambda x:pi/2)[0]