有没有可能强制SymPy计算一个积分，或者告诉我为什么它不能计算积分？
SymPy是一个开源的Python库，用于符号数学计算。它提供了强大的符号计算功能，包括解方程、求导、求积分等。在SymPy中，可以使用sympy.integrate()函数来计算积分。
然而，对于某些积分，SymPy可能无法给出精确的解析解或者计算过程非常复杂，因此无法直接计算积分。对于这种情况，SymPy可能会返回一个未经计算的积分对象。
若要强制SymPy计算一个积分，可以使用sympy.integrate()函数的第二个参数，即积分变量。例如，要计算函数f(x)在区间[a, b]上的积分，可以使用以下代码：
import sympy
x = sympy.Symbol('x')
a = 0
b = 1
f = x**2
integral = sympy.integrate(f, (x, a, b))上述代码中，x是定义的符号变量，a和b是积分区间的上下界，f是要积分的函数。sympy.integrate()函数的第一个参数是要积分的函数，第二个参数是一个元组，其中第一个元素是积分变量，第二个和第三个元素是积分区间的上下界。
如果SymPy无法计算某个积分，可能是因为该积分没有解析解，或者计算过程过于复杂。在这种情况下，可以考虑使用数值积分的方法来近似计算积分结果。SymPy也提供了数值积分的函数sympy.nintegrate()，可以通过设置积分方法和积分精度来进行数值积分。
总结来说，SymPy可以计算许多常见的积分，但对于某些积分可能无法给出精确解析解或者计算过程非常复杂。在这种情况下，可以考虑使用数值积分来近似计算积分结果。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云云服务器（CVM）：弹性、稳定可靠的云服务器实例，可满足不同规模业务的需求。
- 腾讯云容器服务（TKE）：基于 Kubernetes 的高度可扩展的容器管理服务，可简化容器的部署、运维和扩缩容。
- 腾讯云函数计算（SCF）：基于事件驱动的无服务器计算服务，可弹性运行代码，免去服务器运维成本。
- 腾讯云对象存储（COS）：安全、稳定、低成本的云端存储服务，可用于存储和处理大规模非结构化数据。
- 腾讯云人工智能（AI）：提供了丰富的人工智能技术和服务，包括图像识别、语音识别、自然语言处理等，可帮助开发者构建智能化应用。
注意：以上提到的腾讯云产品仅为示例，并不代表其他品牌商的产品。
相关·内容
1回答
为什么SymPy不将标准日志-普通PDF集成到1？
我在Python3.x和SymPy 1.0.1中运行以下代码：
from sympy.stats import density, LogNormal
from sympy import Symbol, integrate, oo
mu, sigma = 0, 1
z = Symbol('z')
X = LogNormal('x', mu, sigma)
f = density(X)(z)
integrate(f, (z, 0, oo))
哪一个应该(？)返回1但产出：
sqrt(2)*Integral(ex
浏览 2提问于2016-11-06得票数 2
回答已采纳
为什么integral_0^1 log(X)/(x^2-1) dx在SymPy中不能工作？
AttributeError：“Not”对象没有属性“_eval_power”
(OK) Wolfram\Alpha示例：
∫%5B0,1%5D+log(x)%2F(x%5E2-1)+dx
integral_0^1 log(X)/(x^2-1) dx =π^2/8？
1.2337
(??)同情
from sympy import *
# var("x")
x = symbols('x', positive=True)
f=log(x)/(x^2-1)
print
浏览 0提问于2018-11-12得票数 0
回答已采纳
SymPy不能计算心脏的长度。如何编写代码来计算闭合线的长度？Python3.7.1 SymPy 1.3
from sympy import *
fi = symbols('fi')
cont = integrate(sqrt((1+sin(fi))**2+(cos(fi))**2), (fi, -pi/2, pi/2))
print(cont)
>>>
Integral(sqrt((sin(fi) + 1)**2 + cos(fi)**2), (fi, -pi/2, pi/2))
浏览 0提问于2019-02-04得票数 2
每当我尝试使用Sympy来集成我的函数时，它不会计算它，而只是返回乘以1.0的输入。
下面是我的代码：
import sympy as sp
a = sp.Symbol('a')
Ωm = 0.31
Ωr = 9*(10**(-5))
ΩΛ = 0.69
Ω0 = Ωm + Ωr + ΩΛ
I = sp.integrate(((Ωm*a**(-1)) + (Ωr*a**(-2)) + (ΩΛ*a**(2)) + (1-Ω0))**(-0.5),a)
print(I)
这是我得到的答案：
1.0*Integral((0.69*a**2 - 8.99999999999235
浏览 0提问于2017-10-21得票数 2
1回答
在SymPy中，是否可以对不定积分应用限制并对其求值？
import sympy
from sympy.abc import theta
y = sympy.sin(theta)
Y_indef = sympy.Integral(y)
Y_def = sympy.Integral(y, (theta, 0, sympy.pi / 2))
Y_def.evalf()生成一个数字。
我正在寻找像Y_indef.evalf((theta, 0, sympy.pi/2))这样的东西来获得相同的答案。
浏览 4提问于2013-02-08得票数 4
回答已采纳
2回答
我不能用Sympy来解这个不定积分。我查看了，它显然是收敛的。
import sympy as sp
a, b, C = sp.symbols("a, b, C", real=True)
E = sp.symbols("E", real=True, positive=True)
chi = C * sp.exp(-a * E) * sp.sinh(sp.sqrt(b * E))
sp.integrate(chi, E)
我试图用指数函数重写这个表达式，但没有结果：
sp.integrate(chi.rewrite(sp.exp).expand().powsimp
浏览 10提问于2020-03-02得票数 1
回答已采纳
我认为Sympy在计算trig函数的时出错了。例如：
from sympy import fourier_transform, sin
from sympy.abc import x, k
print fourier_transform(sin(x), x, k)
的预期答案是
但是Sympy返回0。有时，由于fourier_transform(Heaviside(t)*cos(t),t,omega)和fourier_transform(Heaviside(t)*sin(t),t,omega)返回正确的答案，函数工作得很好。我认为Sympy可能使用Laplace变换来计算傅里叶变换。
浏览 10提问于2014-01-18得票数 9
回答已采纳
2回答
我不明白为什么渐近不能计算从0到t的非常简单的积分。如何解决这个问题？ import sympy
from sympy import sin , cos , sqrt, asinh , log
t= sympy.symbols('t')
x = 't^2'
y = 't^3'
x_derivatives = sympy.diff(x , t)
y_derivatives = sympy.diff(y , t)
expression = x_derivatives**2 + y_derivatives**2
print(expression)
浏览 7提问于2020-10-22得票数 0
回答已采纳
1回答
如何提示逆高斯方程的积分
、、、
逆高斯PDF在Wolfram Alpha 上的符号积分
但是我不能用Sympy来复制它(这需要很长时间才能得到响应)：
from __future__ import division
from sympy import *
A, B, x = symbols('A B x')
ig_pdf = (sqrt(x)**(-3))*exp(-B*x-A/x)
pprint(ig_pdf, use_unicode=True)
ig_cdf = integrate(ig_pdf, x, conds='none')
pprint(ig_cdf, use_unicode=Tr
浏览 0提问于2017-08-23得票数 4
回答已采纳
1回答
假设我有
x=Symbol('x',commutative=False)
和一些表达式
expr = x**2+x
如果你试着把这个
integrate(expr,x)
出现以下错误：
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/usr/local/lib/python2.7/dist-packages/sympy/integrals/integrals.py", line 1295, in inte
浏览 13提问于2018-02-27得票数 2
回答已采纳
我用渐近法计算一些积分，如下所示。
#Calculate Calculus
import sympy
x = sympy.Symbol('x')
f = (6-x*x)**(1.5)
f.integrate()
这将失败并抛出例外情况，例如：
ValueError: gamma function pole
如果我只使用一个整数作为幂num，它就能正常工作。
f = (6-x*x)**(2)
#result: x**5/5 - 4*x**3 + 36*x
浏览 6提问于2016-05-04得票数 4
回答已采纳
我最近用Python代替了Matlab，因为我是真实的；Y被Sympy所兴奋。
但现在我遇到了以下问题：
我需要一种用Python可靠地计算连续傅里叶变换的方法。Sympy有问题的解决方案，包括Diracs (三角洲-函数)，因为他们发生，例如，三角函数等。
例如，如果我尝试
fourier_transform(cos(x),x,v)
输出是0，它应该基于狄拉克增量函数。
有谁知道，如果Sympy的这一部分要改进，还是有另一种方法可以用Python解析地找到fourier变换？
感谢您的回答或任何建议！
浏览 2提问于2015-10-19得票数 8
回答已采纳
1回答
我有以下微分方程，我想用SymPy来解决
这个微分方程有隐式解( h(0) = [0,1)和t= [0，inf])
但SymPy给出了
其他软件包，如Maxima，能够找到的。然而，使用SymPy我就做不到。有没有办法做到这一点？我的代码是
import sympy as sp
sp.init_printing(use_unicode=True)
h = sp.symbols('h', function=True)
t = sp.symbols('t')
eq = sp.Eq(sp.Derivative(h(t),t), (1 -
浏览 4提问于2017-07-26得票数 3
1回答
我是新使用Sympy的，人们告诉我，我可以用Sympy进行积分的数值逼近。我试过的是：
import numpy as np
from sympy.abc import x
import matplotlib.pyplot as plt
from sympy import *
init_printing(use_latex="mathjax")
x = symbols('x')
def Kpol(w,q,thresh,xf,rmf):
fw = 3.*(1.+w)/(5.+3.*w)
Krm = (thresh/fw/rmf/rmf)*((1+
浏览 11提问于2022-07-21得票数 -1
2回答
我使用符号积分来积分循环函数和幂函数的组合函数。
from sympy import *
import math
import numpy as np
t = Symbol('t')
integrate(0.000671813*(7/2*(1.22222222+sin(2*math.pi*t-math.pi/2))-6)**0.33516,t)
然而，当我完成输入时，它给出了一个奇怪的结果：
0.000671813*Integral((3.0*sin(6.28318530717959*t - 1.5707963267949) - 2.33333334)**0.33516, t
浏览 0提问于2016-12-09得票数 0
我有以下积分(受启发)，涉及两个实变量a和b (以及虚单位I)：
f = integrate(exp((-a*t**2+I*b*t)/(3*t**2+1)+I*t*x)*(x/(3*t**2+1)), (t, -oo, oo), (x, 0, oo))
我想迭代两个变量a和b，假设在-10到10之间，例如，for a in range(-10,11):后面跟着一个内环for b in range(-10,11):，然后计算积分，这允许我根据变量a和b绘制一个曲面。
这一做法如下：
from sympy.abc import a, b, x, t
from sympy import *
from
浏览 3提问于2021-12-30得票数 1
回答已采纳
1回答
我正在尝试运行下面的示例以表示同情
import sympy as sp
from sympy.parsing.sympy_parser import parse_expr
n, L, x = sp.symbols('n L x')
sp.integrate(sp.cos(sp.pi*x/L)**2, (x, -L/2, L/2))
由于某种原因，它实际上并不能解决这个积分。它只是打印出原来的积分。
>>> sp.integrate(sp.cos(sp.pi*x/L)**2, (x, -L/2, L/2))
L
-
2
/
|
|
浏览 2提问于2017-09-21得票数 2
回答已采纳
然而，我正在尝试使用sympy.intergrate()函数，但是在指定以下代码时，我一直得到一个TypeError：
import sympy as sp
a, b, z, x, c0 = sp.symbols('a,b,z,x,c0')
a = 0.1
b = 0.5
f = 0.147
c0 = 8.1
z = 1
l = (a * sp.exp(b*z) * c0*sp.exp(f*z))
sp.integrate(l (z, 1, 0))
TypeError Traceback (most
浏览 3提问于2014-03-31得票数 0
回答已采纳
1回答
Sympy中导数的积分函数
、、、、
我正在尝试用Python做一些变分演算，并且我遇到了一个问题，我不太明白它的起源。我正在使用Anaconda套件，木星笔记本和IPython版本5.3.0，与python3.6。
如果我写
from sympy import Funtion, symbols, integrate
from IPython.display import display
t, s = symbols('t s')
f = Function('f')
c = Function('c')
计算以下表达式(它们将是失败部分的元素)
display(f(t, s))
di
浏览 7提问于2017-04-03得票数 1
回答已采纳
考虑SymPy中两个高斯变量的求和：
from sympy import *
from sympy.stats import *
init_printing()
a = Normal('a', 0, Symbol('P', real=true))
b = Normal('b', 0, Symbol('Q', real=true))
ss = a + b
pprint(simplify(density(ss)))
我希望结果是一个具有更多方差的高斯的PDF。相反，结果如下：
⎧ ⎛
浏览 0提问于2018-07-03得票数 1
回答已采纳
我想我误解了什么。我认为这段代码的打印结果必须相等，但它们不是： import sympy
x = sympy.Symbol('x')
summation = sympy.summation(sympy.exp(-x), (x, 0, sympy.oo))
print(summation) # 1/(-exp(-1) + 1)
integration = sympy.integrate(sympy.exp(-x), (x, 0, sympy.oo))
print(integration) # 1 你能解释一下summation和integrate之间的区别吗
浏览 16提问于2019-03-22得票数 1
回答已采纳
我正在尝试使用Sympy来求一个带有复变元的双曲正切函数("Sech")的傅立叶变换。
import sympy as sy
C, t, T0, f, w = sy.symbols('C, T, T_0, f, omega', real=True)
Ut = sy.sech(t/T0) * sy.exp(-sy.sqrt(-1) * C / 2 * t * t / (T0 * T0))
Uf = sy.fourier_transform(Ut, t, f)
不幸的是，当我请求时，Sympy似乎只是挂起了。
这是一个bug吗?或者有没有更好的方法可以将请求呈现给渐
浏览 2提问于2017-09-11得票数 3
doit()函数在sympy中继续进行，并尽可能地计算表达式。例如：
from sympy import *
u = IndexedBase('u')
i = symbols('i')
test = Sum(u[i],(i,1,3))
test.doit()
会回来
另外：
from sympy import *
u,x = symbols('u, x')
test = Derivative(u,x)
test.doit()
创建符号微分，然后计算它。在这种情况下，评估结果为零：
但是如果我想要一个Derivativ
浏览 1提问于2017-09-03得票数 4
回答已采纳
1回答
在博士中，我们可以读到：
manualintegrate模块具有返回所使用的步骤的函数(有关更多信息，请参见模块docstring )。
但是打电话给help(sympy.integrals.manualintegrate)我们得到：
Help on function manualintegrate in module sympy.integrals.manualintegrate:
manualintegrate(f, var)
manualintegrate(f, var)
Compute indefinite integral of a single variable usi
浏览 0提问于2018-02-26得票数 6
回答已采纳
1回答
如果我做了
import sympy
k, V, Vprime = sympy.symbols('k, V, Vprime')
print sympy.diff(k + V(t), t)
我得到了Derivative(V(t), t)，正如我所期望的-导数分布，常数项有零导数。
但是，如果我通过替换构造一个等价表达式，simplify就不会分发这个导数。如何通过替换获得与直接计算表达式相同的结果？
sympy.diff(Vprime(t)).subs({Vprime(t): k + V(t)}).simplify()
返回Derivative(k + V(t), t)。
浏览 3提问于2015-08-19得票数 2
回答已采纳
1回答
我想用计算二进制判决的概率。相似函数用指数定义如下；
sA = b+exp(-|x-xA|)
sB = b+exp(-|x-xB|)
pA = 1/(1+(sB/sA))
为了得到损失，我们需要在各自的范围内集成pA和1-pA除以x。
loss = integrate(pA, (x,0,0.5)) + integrate(1-pA, (x,0.5,1))
当使用b=0 (0.5075)时，我得到了损失，但当b>0时出现了以下错误；
提高sympy.polys.polyerrors.PolynomialDivisionFailed：(f，g，K) PolynomialDivisionFa
浏览 3提问于2014-01-25得票数 1
回答已采纳
2回答
好的，我正在用SymPy编写一个程序。我使用sympy.Eq、sympy.solve和sympy.Symbol。当我编译我的程序并运行使用sympy.solve、sympy.Eq和sympy.Symbol的命令(这是一个用户界面)时。它在没有编译的情况下工作得很好，但是当我用py2exe (作为捆绑的可执行文件)编译它时，它给出了我的这个错误：
Traceback (most recent call last):
File "RateTime.py", line 46, in <module>
File "RateTime.py", line
浏览 0提问于2012-10-12得票数 0
回答已采纳
在不使用SymPy的情况下，我很难创建一个简单的for循环代码来集成给定的函数。我在考虑使用某种黎曼近似，但我不知道如何准确地做到这一点。到目前为止，我的代码是：
def xsquared(x):
n = 2
return x**n
def integral(fun, xmin, xmax):
total = 0
for a in range(xmin, xmax):
x = a
total += fun(x*1.235)
return total
print(integral(xsquared, 0, 4))
输出
浏览 4提问于2018-11-12得票数 1
回答已采纳
1回答
我认为SymPy在集成1/x时返回log(x)是好的，而不是log(abs(x))。这使它简单，让用户担心的标志。
然而，在某些情况下，它会导致错误的答案。例如,
from sympy import *
x = Symbol('x')
integrate(cos(x)**2/sin(x), x)
结果：
log(cos(x) - 1)/2 - log(cos(x) + 1)/2 + cos(x)
而它应该(很明显？)be
log(1 - cos(x))/2 - log(1 + cos(x))/2 + cos(x)
有什么办法能给出第二个答案吗？
浏览 1提问于2018-08-31得票数 1
回答已采纳
我开始学会同情了。我尝试了一个使用渐近的集成问题，我得到了这个异常：
File "/home/me/anaconda3/lib/python3.6/site-
packages/sympy/polys/domains/domain.py",
line 146,
in convert raise CoercionFailed("can't convert %s of type %s to %s" %
(element, type(element), self))
sympy.polys.polyerrors.
浏览 3提问于2017-03-29得票数 1
1回答
是渐近可以快速计算傅里叶变换，还是我需要另一个软件？因为我尝试了很多东西。不过，我对python和Sympy相当缺乏经验。有没有更快的方法用python进行这些Fourier变换？
(a +j \omega) / (4a^2 +(a+j \omega)^2)我尝试了这样的
from sympy import poly, pi, I, pow
from sympy.abc import a,f, n, k
n = poly(a**2 + t**2)
k = t / n
fourier_transform(k , t, 2*pi*f)
我收到这样的信息：
浏览 0提问于2020-10-07得票数 3
1回答
第二类广义积分
、、
我试图用渐近来计算第二类广义积分的精确值：
from sympy import integrate, log
from sympy.abc import x
print (integrate(log(x) * log(x) /(1+x*x), (x,0,1)))
这段代码返回了很多错误。也许我需要使用另一种方法？我已经尝试了Integral，但一无所获。
我想从德怀特表(863.61和863.10)中计算这些积分：
我可能会用数值方法来计算它们，但我更希望得到渐近的精确解。第二类广义积分有可能用渐近得到精确解吗？或者这些积分对于渐近来说太复杂了？
浏览 29提问于2017-12-22得票数 2
回答已采纳
2回答
在我的程序中，我必须重新创建galerkin-method，我必须用它做一些集成，现在我用sympy.integrate做了。 我的问题是，在方程包含多个三角函数或添加了一个e函数(它应该能够解决所有问题)的情况下，渐近将永远计算。我也在使用sympy.symbols，因为通过这些集成，我创建了一个必须求解的方程系统。我需要一个反导数，而不仅仅是一个值的解。 有没有一种数值积分方法或其他方法可以返回一个精确值？我尝试使用sympy.Integral(equation).evalf(1)，但是错误太高，或者返回的十进制数太长，这让我回到了运行时间太长的情况。 这些函数可能如下所示： (-4*x
浏览 20提问于2019-12-14得票数 1
回答已采纳
3回答
我想对表格的一个积分进行评估：
但是，使用python sympy无法做到这一点。这就是我所实施的：
import sympy as sym
r,a,b = sym.symbols('r,a,b')
I = r/sym.sqrt((a-r)*(r-b))
I_1 = sym.integrate (I,r)
print (I_1)
我得到的结果是：Integral(r/sqrt(-(-a + r)*(-b + r)), r)
因此，它不是求积分，而是简化方程。
有人能帮我吗？
浏览 7提问于2022-05-17得票数 2
我遇到了一些问题，下面是我的代码：
from sympy import integrate as inter
import sympy as sp
from sympy import symbols
def f(x):
f1 = (sp.exp(-x ** 2)) / (1 + x ** 2)
f2 = (2 * (sp.cos(x)) ** 2) / (1 + (x - 4) ** 2)
f = f1 + f2
return(f)
def integrate_f_from0(b):
x = symbols('x')
a =
浏览 4提问于2017-12-18得票数 3
回答已采纳
我目前正在使用Sympy来帮助我执行数学计算。现在，我正在尝试执行数值积分，但每次运行脚本时都会出现错误。下面是脚本：
from sympy import *
cst = { 'qe':1.60217646*10**-19, 'm0':N(1.25663706*10**-6) }
d = 3.6*10**-2
l = 20.3*10**-2
n = 217.0
I = 10.2
# Circum of loops
circ = l/n;
# Radius
r = N( circ/(2*pi) )
# Flux through a ring a distanc
浏览 0提问于2012-07-20得票数 3
回答已采纳
1回答
我想写一些代码来计算函数，关于x和y，使用二重积分。对于给定函数，内积分和外积分的下界均为-1，上界均为1。 下面是我尝试用来计算示例多项式的代码： import sympy
example_poly = sympy.Poly(x/2-y**2,x,y)
a = sympy.integrate(example_poly, (x, -1, 1))
b = sympy.integrate(a, (y,-1,1))
print(b) 然而，我只是得到了一个错误提示'GeneratorsNeeded: can't initialize from 'dict‘without
浏览 17提问于2021-05-05得票数 1
回答已采纳
1回答
我编写的代码看起来可能很愚蠢，因为它是一个派生函数的集成。因为它是我正在编写的关于声学分析的其他代码的基本基础。这个分析包含了乘法中不同的导数函数的积分。为此，我使用SciPy进行集成，使用渐近进行区分。但是它给出了一个显示TypeError的错误(“无法将表达式转换为浮点数”)。下面是我写的代码。希望能找到解决办法。
import sympy
from sympy import *
from scipy.integrate import quad
var('r')
def diff(r):
r=symbols('x')
Z = 64.25 *
浏览 3提问于2022-02-06得票数 0
1回答
我尝试了下面的代码来计算这个积分，但是输出是LaTeX格式的积分，而不是计算。
import sympy as smp
g, t,v = smp.symbols('g t v')
expr = smp.sqrt(1 - 1/(g*t +1)**2)
smp.Integral(expr,t)
如果SymPy不能计算这个积分，那么还有其他库可以计算它吗？
浏览 10提问于2022-10-25得票数 0
回答已采纳
我不明白为什么下面代码的答案是10而不是1。有人能帮我理解lambdify是怎么回事吗，或者是什么造成了错误的答案？
import sympy
from sympy.utilities.lambdify import lambdify
from sympy import Function
from sympy.abc import x, y
def kFct(xIndex,loc,k1,k2):
... if xIndex <= loc:
... return k1
... else:
... return k2
...
loc = 0.5
k1 = 1
k2
浏览 2提问于2013-01-07得票数 2
1回答
因此，我创建了这样一个分段函数：
x= sp.symbols('x')
f = sp.Piecewise(
(1, x==0),
(sp.sin(x)*(x+1)/x, True))
但是如果我用x代替0，我会得到nan
f.subs(x,0)
nan
那么第一个问题是为什么==不能很好地处理sympy.Piecewise？好吧，我把它改为：
f = sp.Piecewise(
(sp.sin(x)*(x+1)/x, sp.And(x < 0, x > 0)),
(1, True))
而且效果不错。
接
浏览 0提问于2018-11-07得票数 2
回答已采纳
1回答
我想使用SymPy正态分布，因为我用它做一些符号计算。
from sympy.stats import Normal, cdf
from scipy.stats import norm
norm.cdf(10, 10, 0.1)
def survival(x):
y = Normal("x", 10, 0.1)
return cdf(y)(x)
survival(10)
此代码将产生以下错误：
UnboundLocalError: local variable 'reprec' referenced before assignment
人们
浏览 1提问于2017-11-10得票数 2
回答已采纳
我正在使用Sympy包并试图进行集成，为此，我编写了一个小函数：
from sympy import *
from __future__ import division
init_printing()
x, Rn = symbols('x Rn')
def minimum(k):
expr = (x**(k/2-1)*exp(-x/2))/(2**(k/2)*gamma(k/2))
linkes_int = integrate(x*expr,(x,0,Rn))
rechtes_int = integrate(Rn*expr,(x,Rn,oo))
浏览 5提问于2020-06-17得票数 1
回答已采纳
3回答
其思想是计算以下向量场和曲线的线积分： ? ? 这是我尝试过的代码： import numpy as np
from sympy import *
from sympy import Curve, line_integrate
from sympy.abc import x, y, t
C = Curve([cos(t) + 1, sin(t) + 1, 1 - cos(t) - sin(t)], (t, 0, 2*np.pi))
line_integrate(y * exp(x) + x**2 + exp(x) + z**2 * exp(z), C, [x, y, z]) 但
浏览 68提问于2019-04-19得票数 1
回答已采纳
1回答
我试着找一个多项式函数的根，并使用下面的代码： import sympy
from sympy import Poly, roots
g=sympy.var("x")
p = Poly(x**25-96*x**12-4*x**3+2, gen=g)
print(roots(p)) 我不知道为什么它不能工作。如果我使用一个更简单的多项式函数，比如x**2-1，它就能工作
浏览 91提问于2021-03-22得票数 0
1回答
我无法让sympy将xδ(x)简化为0，其中δ为DiracDelta。
In : simplify(x*DiracDelta(x))
Out: x⋅δ(x)
有没有另一个简化的例程，可以计算狄拉克函数因子中的项，在其自变量的零点处，并检查它们是否消失？
这在构造分段函数时很有用，例如f(x)= x H(x) (其中H是Heaviside函数)，其导数应该是H(x)，但在sympy中是x⋅δ(x) + Heaviside(x)。Mathematica确实将xδ(x)简化为0。
浏览 4提问于2018-01-23得票数 4
回答已采纳
1回答
我是SymPy的新手，我希望有人能指出优化我的代码的方法。
我需要用非常高的十进制位(150-300)对一个涉及到的表达式进行数值计算，而且考虑到要计算的参数空间，每个参数集花费30秒或更长的时间，- which非常长。
我在集成处理中使用了lambdify和mpmath后端和meijerg=True，它显著降低了运行时间。还有其他可以使用的方法吗？理想情况下，将评估时间推到1秒以下是很好的。我的代码是：
import mpmath
from mpmath import mpf, mp
mp.dps = 150 # ideally would like to have this set to
浏览 15提问于2021-04-03得票数 2
2回答
我想知道为什么同情不能解决以下问题：
from sympy import *
ss = symbols('s', real = True)
a = symbols('a', real = True)
f = Function('f')
g = Function('g')
eq = Integral(a*g(ss) + f(ss),(ss,0,oo))
solve(eq, a)
返回是一个空的解决方案列表。我想告诉你足够多的东西，这样我就可以得到一个解决方案：
-1*Integral(f(ss),(ss,0,oo))/Integra
浏览 3提问于2022-11-22得票数 0
回答已采纳
1回答
我试图用Python来解决积分问题，我做了以下代码
from sympy import integrate
import sympy as sy
from scipy.integrate import quad
import math
k_2=Symbol('k_2')
k_1=Symbol('k_1')
t_d=Symbol('t_d')
C_0=Symbol('C_0')
y=Symbol('y')
f=lambda y:((k_2+C_0*sy.exp(-y/t_d))/(k_2+C_0))**(k_1*t_
浏览 10提问于2021-11-03得票数 0
1回答
如何使用python/sympy计算力对物体(在真空中)产生的动量和角动量？
我正在寻找一个库来做这件事，我并不坚持渐近，但一个带有渐近的答案就可以了(不要把问题弄得太宽泛)。
我都试了些什么？好吧，我试着阅读了关于渐近的文档，并意识到它不是针对程序员的，而是针对物理学生的。
浏览 1提问于2017-12-02得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
