ZeroDivisionError в SymPy определенный интеграл; границы интегрирования переменных - PullRequest
Новая
       104

ZeroDivisionError в SymPy определенный интеграл; границы интегрирования переменных

1 голос
/ 17 апреля 2020

Я использую Julia и вызываю SymPy для интеграции двумерного нормального PDF в две переменные, a и n. Есть третья переменная p, которую я хочу включить в границы интеграции, но SymPy выдает ошибку ZeroDivisionError, когда я пытаюсь это сделать.

Выполнение следующего, 

using SymPy
σ = 0.2
ρ = -0.9
f(a,n) = (1/(2*π*(σ^2)*sqrt(1-ρ^2)))*exp(-(1/(2*(1-ρ^2)))*((a^2)/σ^2 + (n^2)/σ^2 -((2*ρ*a*n)/σ^2)))
@vars a n p
integrate(a*f(a,n), (n, -oo, p*a), (a, -oo, oo))

Выход 

ERROR: PyError ($(Expr(:escape, :(ccall(#= C:\Users\Anshu\.julia\packages\PyCall\zqDXB\src\pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'ZeroDivisionError'>

Странно, следующие попытки выполняются без проблем:

Просто интегрируем f, а не a * f: 

integrate(f(a,n), (n, -oo, p*a), (a, -oo, oo))

Замена p на постоянную 1: 

integrate(a*f(a,n), (n, -oo, a), (a, -oo, oo))

1 Ответ

1 голос
/ 17 апреля 2020

После определения сигмы и ро, я получаю: 

julia> integrate(a*f(a,n), (n, -oo, p*a), (a, -oo, oo)) |> string
"0.159154943091895*Integral(a*exp(-a^2/(-2*ρ^2*σ^2 + 2*σ^2))*Integral(exp(-n^2/(-2*ρ^2*σ^2 + 2*σ^2))*exp(2*a*n*ρ/(-2*ρ^2*σ^2 + 2*σ^2)), (n, -oo, a*p)), (a, -oo, oo))/(σ^2*sqrt(-(ρ - 1)*(ρ + 1)))"

Возможно, это ваша версия, у меня есть 

julia> sympy.__version__
"1.5.1"
...