在Sympy中收集多元多项式函数的子表达式
在Sympy中，可以使用subs方法来收集多元多项式函数的子表达式。subs方法用于替换表达式中的变量或者子表达式。以下是一个完整的答案：
在Sympy中，可以使用subs方法来收集多元多项式函数的子表达式。subs方法用于替换表达式中的变量或者子表达式。具体用法如下：
- 首先，我们需要导入Sympy库：
import sympy as sp - 定义多元多项式函数：
f = sp.sin(x) + sp.cos(y) - 定义要收集的子表达式：
sub_expr = sp.sin(x) - 使用
subs方法替换子表达式：substituted_expr = f.subs(sub_expr, sp.Symbol('a'))- 这里我们将
sub_expr中的sp.sin(x)替换为一个新的符号sp.Symbol('a')
- 这里我们将
- 最终的结果将会是
substituted_expr：sp.cos(y) + a
这样，我们就成功收集了多元多项式函数f中的子表达式sp.sin(x)。
对于Sympy中的多元多项式函数的子表达式收集，我们可以举一个具体的应用场景：在数学建模中，有时候需要对多个变量进行组合运算，而其中某些组合可能会多次出现。如果我们能够将这些组合运算进行收集，并用一个新的符号代替，可以简化计算过程，提高效率。
推荐的腾讯云相关产品是腾讯云计算服务（Tencent Cloud Computing Services），该服务为企业和个人提供了丰富的云计算产品和解决方案。具体产品介绍和链接如下：
- 云服务器（CVM）：腾讯云提供的弹性云服务器，可满足不同规模和需求的业务场景。
- 云数据库 MySQL：基于MySQL的关系型云数据库，提供高性能、高可靠的数据库服务。
- 对象存储（COS）：腾讯云提供的高可用、高性能、弹性扩展的对象存储服务。
- 人工智能引擎（AI Engine）：腾讯云提供的全场景人工智能服务，包括图像识别、语音识别、自然语言处理等功能。
通过使用这些腾讯云的产品，您可以快速搭建和部署云计算应用，提高工作效率和灵活性。
希望以上信息对您有所帮助，如果还有其他问题，请随时提问。
相关·内容
我试图合成两个函数，并得到一个双值结果‘’
#!/usr/bin/python
from sympy import *
init_printing(use_unicode=True)
x= symbols('x')
f = x/(x+1);
g = x/(x+2);
print(compose(f,g))
这表明: x/((x + 1)*(x + 2))
应该是x/(2x+2)
我还是不明白。有人有主意吗？
谢谢
浏览 5提问于2022-09-19得票数 1
回答已采纳
我想找到以下两个常量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回答
我想写一些代码来计算函数，关于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
回答已采纳
2回答
收集相似的条件
、、、
我正在用解决SymPy：
solveset(Eq((e-m)/(e+2*m) = n*(a-m)/(a+2*m)), m)
是否有可能通过将类似的术语组合在一起来简化解决方案，就像Mathematica一样？
我发现最近的菜谱是，但它并不直接适用于这里。
我已经进步到：
但是我怎样才能得到更紧凑的形式呢？
浏览 19提问于2022-01-19得票数 2
回答已采纳
我试图将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
回答已采纳
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
回答已采纳
2回答
如何求取渐近表达式中的系数
b = sympy.Symbol('b')
a = sympy.Symbol('a')
x = sympy.Symbol('x')
y = 2*x**3 + a*x**2 +b*x
其结果应该是：
coeff for x**3 =2 coeff for x**2 = a coeff for x = b
浏览 6提问于2020-03-05得票数 0
回答已采纳
如果我们知道变量满足某些方程，是否有可能简化SymPy中的表达式？例如，在Mathematica中，我们可以这样写： Simplify[a+b-c, a+b==c] 当然，在这种情况下，可以求解a并进行替换。但是，对于长表达式，进行全局替换可能没有意义。如果目标是生成尽可能短的表达式，则可能需要对某些项应用替换，其余项保持不变，或者为b而不是a进行求解。 我认为sympy.assumptions模块不能对这几个变量相互限制。是否有可能在SymPy中以任何其他方式实现Mathematica的Simplify[expr, assum]的功能？ 或者，有没有其他开源项目可以做这样的事情？
浏览 23提问于2021-04-23得票数 2
回答已采纳
1回答
在尝试求解四阶多项式x**4+b*x+c=0时，我尝试了一种渐近性的方法，得到了一个意想不到的结果。当b=0,c=-1时，它的根是{−1,1,−,}。然而，这正是我从python3.8.1中得到的：
from sympy import *
init_printing()
from sympy import __version__ as sympy_version
sympy_version
>>'1.6.2'
x,b,c,r,r1,e,e1=symbols('x,b,c,r,r1,e,e1')
e=sympify('x**4+b*x+c
浏览 3提问于2020-11-11得票数 2
回答已采纳
我有以下代码：
from sympy import *
a = Symbol('a')
p = Poly(sqrt(a), domain=QQ)
p.eval(a,2)
在eval之后，我期望的是2的平方根，但是我得到的是：
ValueError: tuple.index(x): x not in tuple
我是不是误会了什么？
浏览 1提问于2017-04-27得票数 1
回答已采纳
我正在实现一个算法，在它中，我需要得到因子，多项式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
1回答
积分交感
、、、、
我试图使用sympy执行以下集成；
x = Symbol('x')
expr = (x+3)**5
integrate(expr)
我期待的答案是：
但被归还的是：
以下代码在MATLAB中工作：
syms x
y = (x+3)^5;
int(y)
我不知道我做错了什么，以便使用sympy来执行这个操作。
浏览 0提问于2015-11-26得票数 4
回答已采纳
1回答
在我正在做的一个项目中，我正在处理包含复杂指数的表达式，比如下面的一个，我的目标是尽可能地简化它： from sympy import Rational, exp, I, pi, pretty, cos, sin
E = Rational(1,20) + (Rational(1,4) + exp(2*I*pi/5)/4)*exp(-4*I*pi/5)/5 + exp(-2*I*pi/5)/20 + (exp(4*I*pi/5)/4 + exp(2*I*pi/5)/4)*exp(2*I*pi/5)/5 + (exp(-2*I*pi/5)/4 + exp(-4*I*pi/5)/4)*exp(4
浏览 15提问于2020-07-16得票数 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
回答已采纳
我试图在python中求解下面的多项式。其中一些常量(2.11和.125567481)是动态的，并将在程序中更改。我需要一种方法来计算这个方程在python中的解。我试过使用Sympy，但它一直在执行，但始终没有想出解决方案。Wolfram给出了4种解决方案，其中3种是复杂的，1种是真实的。我只需要真正的那个。在python中有什么方法可以有效地求解这个方程吗？在mathematica中很容易做到，但我无法找到一种在python3中高效地实现它的方法。
(要将它与symPy的函数结合使用，我必须操纵方程，使其等于零。它仍然没有找到孤子。)
浏览 4提问于2017-10-23得票数 2
回答已采纳
至少它看起来是这样的。下面的代码对于3x3和6x6矩阵的行为是正确的。
deter = mat.det('method':'berkowitz')
#self.resultFileHandler.writeLogStr(str(deter))
a = sy_polys.Poly(deter, k)
对于3x3，执行这段代码需要0.8s，对于6x6则需要288 s( det函数只有650 it，保利的其余部分)。对于10x10，要么是复杂性以惊人的速度上升，要么是其他原因阻止了它从保利电话中返回(我等了一个星期)。不会引发任何异常。
行列式
浏览 6提问于2016-06-30得票数 3
回答已采纳
1回答
我在找一个Haskell库来处理多元多项式。我找到了，但是没有例子，对我来说，文档太复杂了。我想举几个例子：
构造多元多项式
两个多元多项式的相加与相乘
用这个图书馆。
另外，我想知道这个库是否允许处理符号系数。
(1+a)*x1**2 + x1*x2 + a
其中a是一个文字变量(我可以用SymPy和朱莉娅中的SymPy来实现这一点)。也就是说，多项式在域(或环？) Q(a)或Z(a)中有系数。
浏览 5提问于2022-07-23得票数 1
回答已采纳
假设我在SymPy中有一个类似这样的表达式： s = sp.symbols('s')
Y = s**3 + 3*s**2 + 2*s
Y ? 现在，我想获得该表达式的“订单列表”。下面是我想要实现的一个示例:我想创建一个tf()对象。其中，tf()是一个类似Matlab语言的函数，它接受两个列表，一个分母和一个分子，并从中创建一个传递函数： tf([1], [1, 3, 2, 0]) 我想从我的SymPy表达式中提取这两个列表，以便像在tf()函数中一样使用它们。我该怎么做呢？
浏览 19提问于2020-12-31得票数 0
回答已采纳
1回答
from sympy import Symbol
from sympy.solvers import solveset
x = Symbol('x')
equation = x**2 - 2*x + 1
result = solveset(equation, x) 在这里，result的计算结果为FiniteSet(1)。 因为这是一个二次方程，所以必须有2个根。 而且，由于两个根是相同的(在本例中为1)，因此结果只返回一次。 为了解决线性递归，我需要知道重复的根，以及每个根重复了多少次。 如何做到这一点？
浏览 7提问于2021-10-01得票数 1
在Mathematica中，您可以使用函数FullSimplifyexpression，假设来简化使用假设的表达式。例如，如果我这样做：
FullSimplifyx^2-y^2，x^2-y^2==1，
那么结果将是1，因为这是‘最简单的表达式，等同于我给软件的函数。
现在我需要用Python做同样的事情，但是我不知道该怎么做。我已经在此页面上看到了有关命令sympy.refine的文档：
然而，我不能使用等式作为假设。似乎不可能用Q函数来假设这样的事情。
我试着这样做：
import sympy as sp
x,y=sp.symbols('x y')
sp.refine(x*
浏览 14提问于2021-06-25得票数 1
我有一张自动生成方程式的纸条。
方程是用渐近符号构造的。
我想知道这是否是一种方法，以检验这些方程在某些变量方面是否是线性的。
例如：
a, b, c, d = sympy.symbols('a, b, c, d')
eq1 = c*b*a + b*a + a + c*d
检查以下内容：eq1在a、d方面是线性的吗
True
浏览 1提问于2016-03-29得票数 7
回答已采纳
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
回答已采纳
我一直在尝试使用渐近的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
2回答
我先从sympy python lib开始。
如果，我有这样的表达式
from sympy.abc import a,b,c,p,q
e = p * ( a + b ) + q * ( a + c )
如何使用a,b,c作为因子？喜欢
a(p+q) + b*p + c*q
浏览 0提问于2013-06-02得票数 3
回答已采纳
我有以下身份： 1 = a + b + c 假设我有这样的表达式： expr = x*(a + b + c) 它可以简化为x。 有没有一种方法可以将它声明为SymPy，以便它可以简化它们？实际上我是男性化地做这项工作： >>> import sympy
>>> sympy.vars("x a b c")
>>> expr = x*(a + b + c)
>>> expr.subs(a + b + c, 1)
x
浏览 10提问于2021-10-06得票数 1
回答已采纳
让我们假设我们有一个表达式：
x*sin(y) + y*cos(x) + x*cos(z)
从上面的表达式中，我想得到一个新的多项式，它只有系数x的项，即：
x*sin(y) + x*cos(z)
或
x*(sin(y) + cos(z))
显然，设置y=0不是一个解决方案，因为x*sin(y)将变为0。如何滤除SymPy中多项式的项，使其只选择给定系数的项？
浏览 0提问于2018-09-04得票数 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
回答已采纳
我在做一些矩阵运算，有时涉及到其条目具有恒定值的矩阵。
但由于某些原因，即使结果很简单，我也无法得到矩阵运算将结果合并成一个多项式。例如，请考虑以下几点：
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
回答已采纳
2回答
符号方程求解器的设计
、、、
对于一个类，我被分配到一个项目中去做一个基本的方程求解器。它需要解线性方程组。例如，我应该能够解的一些方程：
a*x +b=c
a_x +b= c_x +d
16*x - 9*x =a
X/16-x/9=a
8*(x+3) - 5*(x + 4) = 12
x^2 + 2*x +1=0
(x + 1) *(x-1)=0
(x + 1) (x-1)= 2(x-1)
X- 2*sqrt(x) +1=0
日志3*x - log 6=1
日志3*x - log 6*x^2 =a
2^(x-1) *a= 2^(3*x+1)
我看了看，只找到了几个有
浏览 1提问于2011-11-30得票数 3
2回答
Sympy:在多项式中删除高阶项
、、、、
使用Sympy，假设我们有一个表达式f，它是符号"x“的多项式(以及其他可能的符号)。
我想知道，如果有一种有效的方法，把所有的项都降到比整数n还要大的f。
作为一个特例，我有一个非常复杂的函数，但我只想把项保持在x中的二阶以下，那么有效的方法是什么呢？
最明显的，也不是很有效的方法是，对于小于n的每m，取m导数，设x=0，得到x^m的系数，然后重建多项式。但接受衍生品并不是最有效率的事情。
浏览 7提问于2014-07-28得票数 11
回答已采纳
1回答
我正在尝试学习凯克卡克哈希函数的密码分析。关于零和判别器的一篇论文讨论了逆度为3的Keccak。
我无法弄清楚:你怎么得到3的逆度？例如，Keccak中的一个排列Chi在GF(2)中进行乘法。所以我猜这意味着，Keccak的图像变换有2级，但我不知道它们是如何推导出3级的，用于前像、逆或者向后。
如果有人能向我指出一个资源，这将是很大的帮助，因为我找不到一个。
浏览 0提问于2013-06-25得票数 6
回答已采纳
在使用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回答
当我运行下面的代码时，我得到
TypeError: can't multiply sequence by non-int of type "Add'
有没有人能解释一下为什么会出现这个错误？
from sympy.core.symbol import symbols
from sympy.solvers.solveset import nonlinsolve
x, y, z, r, R, a, m, n, b, k1, k2 = symbols('x,y,z,r,R,a,m,n,b,k1,k2', positive=True)
f1 = r * x
浏览 2提问于2018-02-27得票数 0
from sympy.parsing.sympy_parser import parse_expr
import sympy as sp
def differentiate(exp, n):
parse = parse_expr(exp)
diff = sp.diff(parse, 'x' , n)
answer = sp.expand(diff)
return answer
print(differentiate("x**5 + 4*x**4 + 3*x**2 + 5", 3))
print(differentiate(
浏览 2提问于2020-05-22得票数 0
回答已采纳
1回答
对于某个项目，我使用渐近来计算表达式模块化另一个函数。这些函数都有二进制系数(所以x^2 + 2x = x^2$)。他们的应用是在Galois领域。
我的问题是，当使用带逆的渐近rem函数(例如，x**-1)时，只是返回数字的逆(因此在本例中答案是1/x)，而不是返回模逆。
由于下面的评论，这里有一些进一步的澄清。我正在做的一个过于简化的版本是：
from sympy import *
x = symbols('x')
f = raw_input() #here f = '(x^3 + x)*(x + 1)^2 + (x^2 + x)/(x^3) + (x)^-1
浏览 3提问于2016-05-02得票数 2
回答已采纳
2回答
我想用SymPy获取多变量函数的导数，然后对a)要打印的符号结果，然后b)在要打印的点上的导数的结果。我使用以下代码
import math as m
import numpy
import scipy
#define constants
lambdasq = 0.09
Ca = 3
qOsq = 2
def f1(a,b,NN,ktsq,x):
return NN*x**(-a)*ktsq**b*m.exp(m.sqrt(16*Ca/
浏览 0提问于2018-01-20得票数 2
回答已采纳
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
sympy.line.equation()
同样的价值和类型，但有什么区别呢？“
我怎么才能修好它？
z1和z2有什么区别？
from sympy import *
var('x y z1 z2')
z1=8*x+6*y+48
print("#z1",z1,type(z1))
z2=Line(Point(-6,0),Point(0,-8)).equation()
print("#z2",z2,type(z2))
if type(z1)==type(z2):
print("#","type==")
els
浏览 6提问于2022-07-16得票数 0
回答已采纳
1回答
我有一个需要解的方程组
我知道以下变量：Sale Price、n、BR和每个t_i的不同值。此外，我知道p_i的所有值都是相等的。
我想为p_i解决这个问题。我怎么能这么做？
在下一部分中，我希望根据上面的结果来设置p_i，并更改Sale Price变量，并观察BR在这方面的变化。
这和Sympy有关系吗？
浏览 5提问于2015-05-01得票数 1
我正在尝试测试任何子表达式是否包含sin函数(或任何trig函数)。
from sympy import sin, symbols, Wild
A, B, x, y = symbols('A, B, x, y')
W1=Wild('W1')
我可以这么做：
>> (A*sin(x)+B*sin(y)).has(sin(x))
out: True
但这是行不通的：
>>: (A*sin(x)+B*sin(y)).has(sin(W1))
out: False
我如何测试一个或多个sin函数，而不考虑争论？
浏览 0提问于2014-09-10得票数 1
回答已采纳
1回答
Sympy简化逻辑
、、、
我最近开始使用Sympy，从一开始它的行为在我看来就很奇怪。
考虑一个变量的两个多项式，fx和tn_shifted。我从另一个减法得到多义p2。然后我对它进行简化，得到一个完全不同的表达式，如下所示。
import sympy as spy
from sympy import Symbol
x = Symbol('x')
a, b = 1, 3
n = 3
tn = lambda x: x**3 - 3*x
fx = spy.simplify(x**3 - 2*x**2 + x + 1)
tn_shifted = (b - a)**n / (2**(2*n - 1)) *
浏览 6提问于2016-04-06得票数 1
回答已采纳
我目前正在尝试使用SymPy来生成函数及其梯度，并对其进行数值计算。为简单起见，我将使用以下函数作为示例(请记住，实际函数要长得多)：
import sympy as sp
def g(x):
return sp.cos(x) + sp.cos(x)**2 + sp.cos(x)**3
对这个函数及其导数进行数值计算非常容易：
import numpy as np
g_expr = sp.lambdify(x,g(x),modules='numpy')
dg_expr = sp.lambdify(x,sp.diff(g(x)),modules='numpy
浏览 1提问于2015-06-10得票数 6
2回答
我刚刚开始使用，我对它的一些行为感到有点惊讶，例如，这不是我期望的结果：
>>> import sympy as s
>>> (-1)**s.I == s.E**(-1* s.pi)
False
>>> s.I**s.I == s.exp(-s.pi/2)
False
为什么这些返回值是False，有没有办法让它从一种写复数的方式转换成另一种方式？
浏览 0提问于2009-12-24得票数 4
回答已采纳
laguerre多项式似乎没有按某些顺序收敛，运行下面的代码就可以证明这一点。
import numpy as np
from sympy import mpmath as mp
from scipy.special import genlaguerre as genlag
from sympy.mpmath import laguerre as genlag2
from matplotlib import pyplot as plt
def laguerre(x, r_ord, phi_ord, useArbitraryPrecision=False):
if (r_ord <
浏览 4提问于2014-10-21得票数 0
回答已采纳
我们需要两个微分运算符[B]和[C]矩阵，如下所示：
B = sympy.Matrix([[ D(x), D(y) ],
[ D(y), D(x) ]])
C = sympy.Matrix([[ D(x), D(y) ]])
ans = B * sympy.Matrix([[x*y**2],
[x**2*y]])
print ans
[x**2 + y**2]
[ 4*x*y]
ans2 = ans * C
print ans2
[2*x, 2*y]
[4*y, 4*x]
这也可以用于计算向量场
浏览 4提问于2013-03-18得票数 4
回答已采纳
1回答
我想用Python找到任意阶Chebysev多项式的根。我见过勒让德多项式的。但是，我使用定义为的方法构造了多项式
import numpy as np
import sympy as sp
f0 = lambda x: chebyt(0,x)
f1 = lambda x: chebyt(1,x)
f2 = lambda x: chebyt(2,x)
f3 = lambda x: chebyt(3,x)
f4 = lambda x: chebyt(4,x)
plot([f0,f1,f2,f3,f4],[-1,1])
我尝试过使用np.roots(f4)，但是我收到了以下错误：TypeEr
浏览 1提问于2019-12-02得票数 2
回答已采纳
1回答
SymPy符号数学库()为求解奇异方程和方程组提供了一个求解模块。下面是一个用法的例子：
import sympy as sym
x,y,z = sym.symbols('x,y,z')
c1 = sym.Symbol('c1')
eq1 = sym.Eq(2*x**2+y+z,1)
eq2 = sym.Eq(x+2*y+z,c1)
eq3 = sym.Eq(-2*x+y,-z)
result = sym.solve([eq1,eq2,eq3],(x,y,z))
print(result)
在尝试了一系列自定义方法来解决以下非线性系统后，我使用了solve()，
浏览 0提问于2020-12-04得票数 2
回答已采纳
2回答
我使用Anaconda包中的SymPy 1.4来评估该工具。我有以下几首颂歌要解：
在使用SymPy时，我有以下代码和结果：
import sympy as sp
x=sp.Symbol('x')
y=sp.Function('y')(x)
diffeq=sp.Eq(sp.diff(y,x), (5*x**2-x*y+y**2)/x**2)
sp.dsolve(diffeq,y)
如果我使用WolframAlpha online，我会得到：
这是软件中的错误吗？如果是这样，我可以在哪里报告它？
浏览 3提问于2019-08-28得票数 1
1回答
我有以下Sympy表达式
expr=b0*d0*u0 - b0*d1*u1 - b1*d0*u1 - b1*d1*u0 + d0*b0*u0 - d0*b1*u1 - d1*b0*u1 - d1*b1*u0
我想知道，例如，产品是否
d0*u0
在这个表达式中。为此，我使用
print(expr.has(d0*u0))
但结果是
False
然而，如果我不询问它是否在表达式中就替换了这个子表达式，Sympy就会毫不费力地这样做，例如
print(expr.subs(d0*u0,x0))
b0*x0 - b0*d1*u1 - b1*d0*u1 - b1*d1*u0 + d0*b0*u0 - d0
浏览 0提问于2018-01-27得票数 2
回答已采纳
3回答
我有几个关于向量u的范数或范数平方的表达式。我想简化这些表达式，用一个已知的值来代替u的范数。然而，似乎显然涉及规范的简单倍数的表达式并没有得到简化。
例如，此代码执行我所期望的操作：
import sympy as sp
u1,u2,u3 = sp.symbols('u_1, u_2, u_3',real=True,positive=True)
utu = u1**2 + u2**2 + u3**2
print("Ex. 1")
print(utu.subs(utu,1))
这将产生预期的输出。
Ex. 1
1
然而，2*utu并没有像我所期望的那样简化
浏览 13提问于2022-05-05得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
