Question

当我使用函数f2时，积分求解（相当快），但包含一些浮点常量。函数f1不使用浮点指数，但无法计算积分。它反而重新显示要解决的积分（经过很长一段时间）。

所以作为SymPy的新用户，我想知道1）我是否在f1中使用了一些错误的命令？ 2）是否可以更快地执行SymPy（因为它实际上并不与Maple的速度相比）。

from sympy import *
from IPython.display import display

init_printing(use_unicode=False, wrap_line=False, no_global=True)

def f1():
    x, y = symbols('x y')
    w, h = symbols('w h', real=True, nonzero=True, positive=True)
    result = Integral((1/sqrt(((y-x)**2 + h**2)**3)), (x,0,w), (y,0,w))
    display(result)
    result = integrate((1/sqrt(((y-x)**2 + h**2)**3)), (x,0,w), (y,0,w))
    display(result)

def f2():
    x, y = symbols('x y')
    w, h = symbols('w h', real=True, nonzero=True, positive=True)
    result = Integral(1/(((y-x)**2 + h**2)**1.5), (x,0,w), (y,0,w))
    display(result)
    result = integrate(1/(((y-x)**2 + h**2)**1.5), (x,0,w), (y,0,w))
    display(result)

Sympy版

>>> sympy.__version__
>>> '0.7.6.1'

Answer 1

你写的两个积分并不完全相同。你可以做到

In [22]: integrate((1/sqrt(((y-x)**2 + h**2))**3), (x,0,w), (y,0,w))
Out[22]:
          ________
         ╱      2
        ╱      w
  2⋅   ╱   1 + ──
      ╱         2
    ╲╱         h     2
- ──────────────── + ─
         h           h

请注意**3的展示位置的差异。原因是

In [25]: sqrt(x**3)
Out[25]:
   ____
  ╱  3
╲╱  x

In [26]: sqrt(x)**3
Out[26]:
 3/2
x

对于一般x，两者并不相等。对于你的情况，它们实际上是相同的，因为根中的表达式是正的，但SymPy没有注意到这一点。

SymPy无法解决重写积分问题

1 个答案: