当我使用poly()或.as_poly()时，Sympy的sqf()和sqf_list()给出了不同的结果
当使用poly()或.as_poly()函数时，Sympy的sqf()和sqf_list()函数会给出不同的结果。
首先，poly()函数用于创建多项式对象，而.as_poly()函数用于将表达式转换为多项式对象。这两个函数在处理多项式时非常有用。
接下来，sqf()函数是Sympy中的因式分解函数，它将多项式分解为平方因子的乘积。它返回一个表示分解结果的多项式对象。
而sqf_list()函数也是一个因式分解函数，但它返回一个包含分解结果的列表。列表的每个元素都是一个元组，其中包含一个因子和其对应的指数。
这两个函数的不同之处在于返回结果的形式。sqf()函数返回一个多项式对象，而sqf_list()函数返回一个列表。
在实际应用中，根据具体的需求选择使用哪个函数。如果只需要得到分解后的多项式对象，可以使用sqf()函数。如果需要得到分解结果的详细信息，可以使用sqf_list()函数。
以下是一些示例应用场景和推荐的腾讯云相关产品：
- 应用场景：
- 多项式因式分解
- 多项式求解
- 推荐的腾讯云相关产品：
- 云函数 SCF（https://cloud.tencent.com/product/scf）：用于实现无服务器计算，可以将多项式因式分解的计算逻辑部署为云函数，实现按需计算。
- 云数据库 CDB（https://cloud.tencent.com/product/cdb）：用于存储多项式数据，方便在多个应用程序之间共享和访问。
- 云服务器 CVM（https://cloud.tencent.com/product/cvm）：用于部署和运行多项式计算的应用程序，提供高性能的计算资源。
- 人工智能平台 AI Lab（https://cloud.tencent.com/product/ailab）：提供强大的人工智能算法和工具，可以应用于多项式计算中的模式识别和优化问题。
请注意，以上推荐的腾讯云产品仅供参考，具体选择应根据实际需求和项目要求进行评估。
相关·内容
我在做一些矩阵运算，有时涉及到其条目具有恒定值的矩阵。
但由于某些原因，即使结果很简单，我也无法得到矩阵运算将结果合并成一个多项式。例如，请考虑以下几点：
from sympy.matrices import *
import sympy
x= sympy.symbol.symbols('x')
Poly_matrix = sympy.Matrix([[sympy.Poly(x, x)],[sympy.Poly(x, x)]])
constant_matrix = sympy.Matrix([[0,1]])
constant_matrix_poly = constant_m
浏览 8提问于2019-09-23得票数 8
回答已采纳
我对Sympy Poly.all_coeffs()返回的结果的数据类型有疑问。我最近才开始使用Sympy。我的Sympy传递函数如下： ? 然后我运行以下代码： n,d = fraction(Gs)
num = Poly(n,s)
den = Poly(d,s)
num_c = num.all_coeffs()
den_c = den.all_coeffs() 我得到了： ? 然后我运行以下代码： from scipy import signal
#nu = [5000000.0]
#de = [4.99, 509000.0]
nu = num_c
de = den_c
sy
浏览 21提问于2020-01-10得票数 0
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回答
我想使用Sympy多项式，但我也想使用更高精度的系数。
这么做似乎给了我sympy.core.numbers.float系数的多项式。
import sympy
from sympy import Poly
from sympy.abc import x
from gmpy2 import mpfr, get_context
get_context().precision = 150
#float64 can't tell this from 1.0
one_and_change = mpfr('1.0000000000000000000000000000000000001
浏览 3提问于2015-02-20得票数 1
回答已采纳
如果我对我创建的一个很好的示例使用sqf()函数或sqf_list()函数，它会给我一个很好的结果。例如：
v = (x1 + 2) ** 2 * (x2 + 4) ** 5
sqf(v) = (1, [(x1 + 2, 2), (x2 + 4, 5)])
但是如果我第一次使用v = poly(v)，它只找到= (1, [(Poly(x1 + 2, x1, x2, domain='ZZ'), 2)])。
这是有意为之吗？我假设它是由
v = (x1 + 2) ^ 2 * (x2 + 4) ^ 5
转到
v = Poly(x1^2*x2^5 + 20*x1^2*x2^4 + 1
浏览 16提问于2019-08-17得票数 0
我试图找出为什么sympy.diff不像预期的那样区分sympy多项式。通常，如果定义了一个符号变量，并且没有使用sympy.diff定义多项式，那么sympy.Poly就可以很好地工作。但是，如果函数是使用sympy.Poly定义的，那么sympy.diff似乎并不实际计算导数。下面是一个代码示例，它展示了我的意思：
import sympy as sy
# define symbolic variables
x = sy.Symbol('x')
y = sy.Symbol('y')
# define function WITHOUT using sy.P
浏览 1提问于2015-06-29得票数 4
回答已采纳
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
我试图将dup_zz_mignotte_bound(f，K)的结果与其他边界进行比较，但当我添加dup_zz_mignotte_bound时，没有在所有边界中使用变量x。例如：
from sympy.polys import ring, ZZ
from sympy.abc import x
from sympy import factor
R, x = ring('x', ZZ)
poly = x**8 +8*x**7 +47*x**6 +136*x**5 +285*x**4 +171*x**3 - 20*x**2 - 21*x+2
print(factor(poly,
浏览 3提问于2020-04-29得票数 2
回答已采纳
我有以下几点
x=Symbol('x',commutative=False)
y=Symbol('y',commutative=False)
expr = 2*x + 87*x*y + 7*y
现在，这起作用了
integrate(expr,y,manual=True)
因为它给了
2*x*y + 87*x*y**2/2 + 7*y**2/2
但与x完全相同的是失败的：
integrate(expr,x,manual=True)
Traceback (most recent call last):
File "<stdin>",
浏览 0提问于2018-02-28得票数 1
回答已采纳
1回答
我试图实现一种算法，该算法涉及两个变量"x“和"y”中的多项式，但有时我需要将它们解释为单变量多项式(即，将x作为常数)，例如，为了使用函数 (扩展的欧几里得算法)。是否有一种简单的方法可以将"x“解释为常量而不是变量？
我试过以下几种方法：
import sympy
x = sympy.Symbol('x', constant=True)
y = sympy.Symbol('y')
f = sympy.Poly(x*y + y**2)
g = sympy.Poly(x+y)
(s, t, gcd) = sympy.gcdex(f,g)
浏览 3提问于2016-10-16得票数 3
回答已采纳
我想找到以下两个常量cons1和cons2的间隔，我编写了如下代码
from sympy import Poly
from sympy import Abs
from sympy.solvers.inequalities import solve_rational_inequalities
from sympy.abc import x
cons1=2*((x+2)**2)-Abs(x)-1
cons2=exp(2*x+1)-2.5
solve_rational_inequalities([[((Poly(cons1), Poly(1, x)), '<=')
浏览 0提问于2018-03-11得票数 2
回答已采纳
1回答
如何在SymPy中表示Galois字段？我在网上找不到这方面的任何文档，但是SymPy包含一个名为"galoistools“的模块，所以我想我应该试一试。我尝试了以下实验：
from sympy import *
x = symbols("x")
A = [LC(Poly(i*x, modulus=8) * Poly(j*x, modulus=8)) for i in range(1, 8) for j in range(1, i+1)]
B = [LC(Poly(i*x, domain=GF(8)) * Poly(j*x, domain=GF(8))) for i i
浏览 2提问于2018-05-05得票数 2
回答已采纳
我希望在Symbols中扩展SymPy类，以便可以添加一个布尔属性。我可以用一个符号来完成这个任务(请参阅我的问题，以及其他人的问题)。实现这一目标的代码如下所示：
from sympy.core.symbol import Symbol
class State(Symbol):
def __init__(self, name, boolean_attr):
self.boolean_attr = boolean_attr
super(State, self).__init__(name)
但是，这个解决方案的问题是，当我定义一个多项式或涉及多个Stat
浏览 1提问于2019-08-02得票数 0
回答已采纳
3回答
我有一个Sympy矩阵A和一个多项式表达式P，我想计算P(A)。
下面是一个示例：
x = Symbol('x')
P = x**2 - 3*x + 5
A = Matrix([ [1,3], [-1,2] ])
P.subs(x, A)
我期望Sympy计算A**2 - 3*A + 5*eye(2) (在示例中，结果是零矩阵)。
但是这个错误消息失败了：
AttributeError: ImmutableMatrix has no attribute as_coeff_Mul.
有什么办法能得到我想要的吗？
编辑：--我曾尝试将P转换为Sympy的多项式类，然后进行替换，但是
浏览 2提问于2015-09-07得票数 3
回答已采纳
我一直在尝试使用渐近的solvers.solve_poly_system，但当我使用sqrt函数求解一个方程时，我得到了以下错误：
sympy.polys.polyerrors.PolynomialError: sqrt(x1) contains an element of the set of generators.
下面是有问题的代码：
# An array of variables is declared as symbols in sympy. This has already worked to solve equations/derivatives.
for x in range
浏览 3提问于2021-04-15得票数 0
1回答
这学期我在研究生院上了一门密码学课程，有一次我们讨论的题目是NTRU。我正试图用纯Python编写这段代码，纯粹是作为一种爱好。当我试图找到多项式的逆模p(在本例中为p= 3)时，当我想要严格的正系数时，SymPy总是返回负系数。这是我的密码。我会解释我的意思。
import sympy as sym
from sympy import GF
def make_poly(N,coeffs):
"""Create a polynomial in x."""
x = sym.Symbol('x')
coeff
浏览 0提问于2018-12-13得票数 2
回答已采纳
我已经创建了一个函数，它在给定系数范围的情况下，用这些系数构造多项式，并输出所有根的列表。然而，Numba并不喜欢它。它是这样的：
import math
import numpy as np
import itertools
from numba import jit
from sympy.solvers import solve
from sympy import Symbol
from sympy import Poly
@jit
def polyn(ranges=[[-20,20],[-20,20],[-20,20],[-20,20]],step=4):
l = []
浏览 9提问于2017-12-16得票数 1
回答已采纳
1回答
我正在试图求解大型的任意方程组(对于流程流板应用程序来说，物理组件模型是用简单的数学编写的)。
在处理非整数时，我对渐近的solve_poly_system函数的行为感到困惑，并且无法通过文档找到它。这个简单的代码片段：
sympy.solve_poly_system([x-5,y-x])
给我(5,5)如我所期望的。但是，当我将数字更改为浮点数时：
sympy.solve_poly_system([x-5.1,y-x])
我得到了错误:引发DomainError(“不能在%s之上计算Groebner基”% opt.domain) sympy.polys.polyerrors.DomainEr
浏览 2提问于2013-12-10得票数 3
回答已采纳
1回答
我想知道在sympy中是否可以使用向量符号来求多项式和表达式的导数。例如，如果我有一个表达式作为两个坐标x1和x2的函数，我是否可以只调用一次diff(x)，其中x是x1和x2的向量，或者我是否需要对x1和x2进行两个单独的diff调用，并将它们堆叠在一个矩阵中？ 这说明了什么是有效的，而不是我想要的理想工作： import sympy
from sympy.matrices import Matrix
# I understand that this is possible:
x1 = sympy.symbol.symbols('x1')
x2 = sympy.sy
浏览 47提问于2019-08-30得票数 3
2回答
我不懂同情，但我想解决以下问题：
有多个不等式约束的形式
A < f(x) < B
C < g(x) < D
...
其中A，B，C，D只是数字。F和g是有理函数。
例如，我可以让以下内容发挥作用：
solve_rational_inequalities([[
((Poly(x-10000.00), Poly(1, x)), '>'),
((Poly(x-100000.00), Poly(1, x)), '<'),
((Poly((x/130000.00)-0.00), Poly(1, x)), '
浏览 0提问于2018-01-15得票数 0
回答已采纳
在使用sympy (当前版本)求解多项式方程时(在这种情况下，polynom为d 2)：
from sympy import solve_poly_system
solve_poly_system(4*d**2*sin(a)**2*sin(b)/cos(b)**2 - d*cos(a) + 4, d**2*sin(a)**2*sin(b)/cos(b)**2 - d*cos(a) + 8, 3*d**2*sin(a)**2*sin(b)/cos(b)**2 - d*cos(a) + 3 ,d ,a, b)
我收到以下错误：
PolynomialError: cos(a)包含发电机集的一个元素
浏览 1提问于2016-02-07得票数 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回答
在下面的代码中，我定义了实数符号来获取渐近极限，然后使用subs来获得数值结果。 #To take the limit with unknown values, I have to define real symbols
from sympy import *
x = symbols("x")
H = symbols("H")
p = symbols("p")
#-----
# Position_1
#-----
dict = {H: 5,p:3,x:1}
### Some operations
H = Symbol("H"
浏览 37提问于2020-04-03得票数 0
回答已采纳
我尝试为使用以下代码
# x, y -- x and y coordinates of point
# poly -- a list of tuples [(x, y), (x, y), ...]
def isPointInPath(x, y, poly):
num = len(poly)
i = 0
j = num - 1
c = False
for i in range(num):
if ((poly[i][1] > y) != (poly[j][1] > y))
浏览 2提问于2017-01-06得票数 0
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
回答已采纳
我有一个4x4的多项式渐近矩阵：
S = Matrix([[Poly(alpha**4 + alpha**3, alpha, modulus=2), Poly(alpha**4 + alpha**3 + alpha + 1, alpha, modulus=2),Poly(alpha**3 + alpha**2 + alpha, alpha, modulus=2), Poly(alpha**4 + alpha**3 + alpha**2 + alpha, alpha, modulus=2)],
[Poly(alpha**4 + alpha**3 + alpha + 1, alpha
浏览 22提问于2020-06-01得票数 0
1回答
我有一个CSV文件，包含80个特性(列)和4000个实例(行)。我想把这个文件转换成更高的维度(内核3)。新文件应该有88560列和4000行。我试着做一个只有3种价值的飞行员。我首先准备了符号a、b、c的标题，然后尝试将值放入，但无法用值替换字符串。代码如下所示。
import csv
fw = open('polyexp.csv', "w")
a,b,c= sympy.symbols("a b c")
formula = (a+b+c) ** 3
poly = str(formula.expand())
ls = poly.split(&#
浏览 15提问于2016-03-13得票数 1
回答已采纳
我有sympy.poly，我想共轭所有系数。这是很容易做到的，但是我如何从新的单调式和系数中重新组合多项式呢？
MWE：
from sympy import symbols, poly, I
x = symbols("x")
p = poly(x ** 2 + I, [x])
print(p.coeffs())
conj_coeffs = [c.conjugate() for c in p.coeffs()]
print(conj_coeffs)
# poly(p.coeffs(), p.monoms()) ??
[1, I]
[1, -I]
浏览 1提问于2020-06-11得票数 2
回答已采纳
我写了一个Python程序，根据给定的时变电场计算磁场。我使用的是sympy.vector的CoordSys3d
我今天下午把它修好了。它给出了正确的答案。然后，我升级到最新版本的Anaconda (更新了集成库)，现在当我调用SymPy ()时，它抛出了一个错误。
下面是错误回溯：
>>> integrate( jimmy, t )
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/opt/anaconda3/li
浏览 31提问于2021-09-26得票数 0
首先，
solve_poly_system( seq, *gens, **args),
有人确切地知道solve_poly_system的参数是什么意思吗？
我有系统，
rd = λk
ua = λk
k = ad
其中r，u，k是已知的，a，d，λ是变量。
当我使用solve( system，a，d，λ)时，我得到错误，
NotImplementedError: only zero-dimensional systems supported (finite number of solutions)
此错误在文件中引发
/usr/lib/python2.7/dist-packages/symp
浏览 1提问于2012-09-10得票数 1
回答已采纳
2回答
我需要显示一个周长。为了做到这一点，我认为我可以计算很多x的y的两个值，所以我这样做了：
import sympy as sy
from sympy.abc import x,y
f = x**2 + y**2 - 1
a = x - 0.5
sy.solve([f,a],[x,y])
这就是我得到的：
Traceback (most recent call last):
File "<input>", line 1, in <module>
File "/usr/lib/python2.7/dist-packages/sympy/sol
浏览 0提问于2012-10-19得票数 2
我正在实现一个算法，在它中，我需要得到因子，多项式f(x) = p(x)q(x)，其中p和q是相对素数。当然，我可以使用sympy.factor，但我希望构造一个函数，该函数可能返回sympy.factor结果中的因素列表。
例如,
x = sympy.Symbol('x')
y = sympy.Symbol('y')
g = sympy.Poly(y**4 + x*y**3 + y**2 + x*y, y)
sympy.factor(g)
抛出
y*(x + y)*(y**2 + 1)
因此，也许我可以构建一个解析器来解释字符串，并将其划分为圆括号，显示列表y
浏览 1提问于2016-10-25得票数 3
我考虑以下矩阵： M1 = Matrix([[1/7,2/7],[3/7,4/7]])
M2 = Matrix([[1,2],[3,4]])/7 它们显然是相同的，但当我确定它们的行列式时，我得到了不同的结果： print(M1.det())
print(M2.det()) 给出了以下结果： -0.0408163265306122
-2/49 我希望第一个结果表示为有理数，而不是浮点数。
浏览 13提问于2021-01-20得票数 0
回答已采纳
1回答
不知道如何乘两个不同变量的多项式。下面是我的代码，并在IPython上运行。
from sympy import *
from numpy import *
m1 = poly1d([1,0,0,1], variable = 'x')
m2 = poly1d([1,0], variable = 'y')
p=m1*m2
print(p)
期望结果是一个变量x和y的多项式，但下面是我的结果。
4
1 x + 1 x
浏览 2提问于2017-05-30得票数 0
回答已采纳
我试图为变量x (依赖于其他符号b, c, d)求解以下方程：
from sympy import symbols, sin, cos, solve, solveset, S
b, c, d, x = symbols('b c d x', real=True)
expr = sin(x)*sin(c - 2*d - x) - sin(2*b - 2*d - x)*sin(c - x) # find `x` such that `expr == 0`
还有适用于各种符号的其他限制因素，即：
0 < d < b < d+x < c
因此，我尝试使用
浏览 11提问于2022-04-27得票数 0
回答已采纳
1回答
我想要完全分解一个多边形，。
SymPy提供了factor来完成这一任务，但我非常惊讶的是，分解只能在整数根上完成，例如：
>>> from sympy import *
>>> z = symbols('z')
>>> factor(z**2 - 1, z)
(z - 1)*(z + 1)
>>> factor(z**2 + 1, z)
z**2 + 1
或
>>> factor(2*z - 1, z)
2*z - 1
>>> factor(2*z - 1, z, ex
浏览 1提问于2020-06-13得票数 0
回答已采纳
1回答
例子:让
M = Matrix([[1,2],[3,4]]) # and
p = Poly(x**3 + x + 1) # then
p.subs(x,M).expand()
给出错误：
<class'sympy.matrices.immutable.ImmutableDenseMatrix'>：不能添加TypeError和<class‘>--core.number.class>
这是很有道理的，因为第一项变成了矩阵，但最后一项(常量项)不是矩阵，而是标量。为了纠正这种情况，我将多项式改为
p = Poly(x**3 + x + x**0)
浏览 2提问于2021-01-11得票数 1
回答已采纳
1回答
我正在尝试线性化一个符号表达式，就像在中为collect函数给出的那样：a*x**2 + b*x**2 + a*x - b*x + c。我想用x代替x+dx，然后收集dx中线性的所有项。但是，当我用coeff函数尝试collect函数时，它没有给出所需的答案：
>>> from sympy import S, collect, expand, factor, Wild
>>> from sympy.abc import a, b, c, x, y, z
>>> collect(a*x**2 + b*x**2 + a*x - b*x + c,
浏览 0提问于2019-07-24得票数 0
我使用的是渐近性，并试图计算多项式的根数。
from sympy.abc import x
from sympy import Poly
p = Poly(x**4+0.1,x)
此时，p是具有区域'RR'：Poly(1.0*x**4 + 0.1, x, domain='RR')的多项式
如果我试图计算区间中的根数，就会得到：
p.count_roots(0,2)
TypeError: unsupported operand type(s) for *=: 'RealElement' and 'PythonRational'
但是
浏览 2提问于2016-10-19得票数 1
回答已采纳
1回答
为什么num_den_to_sympy(b, a)和sy.simplify(sy.expand(tf))不一样？扩展版本包含一个附加的s ** 4，我不知道它是从哪里来的。
import numpy as np
import sympy as sy
from scipy.signal import *
from IPython.display import display
sy.init_printing()
def num_den_to_sympy(num, den, symplify=True):
s = sy.Symbol('s')
G = sy.Po
浏览 5提问于2017-05-28得票数 2
回答已采纳
我遇到了这个问题，那就是我不能将整数加到BigDecimal中。在代码中，我查看了错误，它显示“二元运算符‘+’的操作数类型不正确”。如何将整数添加到BigDecimal？(在代码中表示为1+ sqf)
BigDecimal sqf = new BigDecimal(Math.sqrt(5));
sqf.setScale(100);
BigDecimal bd = new BigDecimal((1 + sqf) / 2); //Error here (1 + sqf)
bd.setScale(100);
System.out.println(isBuzzNumber(77707));
Sys
浏览 1提问于2014-09-11得票数 2
我想使用渐近来求解满足一组二次和线性方程的向量，例如x.T @ A @ x = x.T @ B @ x = 1和x.T @ C = 0。这是我的尝试： import numpy as np
from sympy.polys.polymatrix import PolyMatrix
from sympy import symbols
from sympy.solvers import solve_poly_system
x = PolyMatrix([symbols(f'x{k}') for k in range(2)])
A = np.eye(2)
B = np.array
浏览 12提问于2020-06-25得票数 1
回答已采纳
有没有办法把Sympy的解决方案和Astropy结合起来？例如，如果我想解开普勒关于所有其他变量的给定值的p的第三定律：
from sympy import solveset, Poly, Eq, Function, exp
from sympy import Symbol, pi
Ms = Symbol('Ms') # mass of sun
Mp = Symbol('Mp') # mass of planet
G = Symbol('G') # Gravitational consatant
a = Symbol('a') #
浏览 3提问于2020-12-29得票数 2
如何使用sympy.jl或另一个用于符号计算的Julia软件包获得多变量多项式的所有组合的系数？
下面是MATLAB的一个例子，
syms a b y
[cxy, txy] = coeffs(ax^2 + by, [y x], ‘All’)
cxy =
[ 0, 0, b]
[ a, 0, 0]
txy =
[ x^2y, xy, y]
[ x^2, x, 1]
我的目标是
[ x^2y, xy, y]
[ x^2, x, 1]
而不是x^2，y
我在和上也问了同样的问题，但是我想我应该问一下是否可以使用Sympy.py来完成这个任务。
利用朱莉
浏览 10提问于2022-11-05得票数 1
1回答
我有这个多项式，我需要找到它的根。我遇到的问题是，我的多项式有一个外部选择的参数，为了解释，我需要找到根。我的方法是在real_roots()方法中使用SymPy：
from sympy.abc import x, c
from sympy import Poly, real_roots
Poly(4-2*c*c*x+c*c*x**3,x,domain='RR[c]').real_roots()
但是我得到了错误NotImplementedError: sorted roots not supported over RR[c]，我很难解释这个错误。这是逼近这个多项式的根的正确方
浏览 11提问于2022-09-16得票数 1
回答已采纳
1回答
import sympy as S
F = S.FiniteField(101)
当我调用f = S.poly(y ** 2 - x ** 3 - x - 1,F)时，会得到以下错误：
'FiniteField‘对象没有属性'is_commutative’
但是有限域在定义上是可交换的！所以我不太清楚这个错误是什么意思！
以前有人遇到过这种情况吗？如何在有限域上声明多项式？
浏览 1提问于2015-05-01得票数 11
1回答
我正在尝试用符号来解决一个复数的多项式，以及它们与SymPy的共轭。我想我已经走了很长的路，但是solve没有给我任何解决方案，尽管多项式是可解的。
from sympy import *
# set up symbols
a, b = symbols("a b", real=True)
t = a+I*b
T = functions.conjugate(t)
# set up polynomial
a1=0.005+I*0.0009
a2=0.9+I*-0.9
a3=0.4+I*0.5
a4=8+I*-80
a5=284+I*-1.5
a6=27100+I*-11500
浏览 1提问于2013-04-05得票数 2
回答已采纳
3回答
我正在编写一个非常简单的脚本，它读取一个相当大的文件(3M行，1.1G文件)，其中包含多项式的litteral (str)表达式。然后，我使用Sympy进行一些符号计算，并将结果写入16个单独的文件。
我的脚本在运行时占用了更多的内存空间(> 20 Gb)，我不明白为什么。你能找到任何方法来改善这个脚本的内存使用吗？
from sympy import sympify
from sympy.abc import x,y
from sympy import degree
fin = open("out_poly","r")
A = fin.readline
浏览 1提问于2013-11-05得票数 4
回答已采纳
我希望在SymPy中生成一个多项式，但理想情况下，我希望尽可能通用地编写代码：
from sympy import *
sx, sa0, sa1, sa2 = symbols('x sa0 sa1 sa2')
sa = MatrixSymbol('a',3,1)
polynomial = sum([sa[i,0]*sx**i/factorial(i) for i in range(sa.shape[0])])
下面的表达式显式地给出了多项式，但代码并不像我所希望的那样一般：
poly2 = sx**2*sa2/2 + sx*sa1 + sa0
在这个表达式
浏览 2提问于2014-07-23得票数 0
回答已采纳
新的SymPy，并正在试验除法多项式。目前，我正在运行Python3.6.4和SymPy 1.1.1 (安装了pip)。在浏览了他们的文档之后，我在他们的上找到了以下有用的例子：
另一种选择是同时除以多个多项式。通常，输出不是唯一的，取决于除数的顺序和给定的单项顺序(如果指定的话)。
F= x*y + y*z + z*x >>> g1 =x+1 >>> g2 = 2*y +1 >>> q，r= div(f，g1，g2，x，y，z) >>> q y+ z，-1/2 +-2 >>> r3*z 1/2 -
浏览 0提问于2018-03-12得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
