sympy有一个内置的方法来设置多个多项式彼此相等吗



这里有一个简单的例子,说明我如何开始尝试在sympy中求解2个方程组。我的目标是首先设定彼此相等的方程式:

import sympy as sp
x,y,z = sp.symbols("X, Y, X")
print(sp.Poly((x-y),sp.Derivative).coeffs())
print(sp.Poly(sp.Eq(x,y),sp.Derivative).coeffs())

然而,我需要设置彼此相等,而不仅仅是用两个方程。以下是我尝试过但不起作用的:

import sympy as sp
x,y,z = sp.symbols("X, Y, X")
print(sp.Poly(sp.Eq(x,y,z),sp.Derivative).coeffs())

类型错误:new((从2到3个位置参数中获取,但有4个被赋予

对于类似x-y的更代数的方法,我不知道如何用3个或更多的方程来实现。

如何将我的示例扩展到三个或多个方程?这可能与症状有关吗?如果没有,是否有其他类似的库?如有任何帮助,我们将不胜感激。

编辑:因此,当我有两个多项式(或阶乘方程(方程,除了给定的β常数之外,它们的结构非常相似,理论上它们的一阶微分应该相等。这模拟了maclaurin级数,其中我的2个阶乘方程试图收敛到一个一般的未知函数。我会使用sp.idiff((来隐式地区分这两个方程。sp.idiff((将隐式求解dx0,其中我仍然需要求解dx1。我求解dx1如下:

equation1 = sp.idiff(expr1,array-like,respective_beta)
equation2 = sp.idiff(expr2,array-like,respective_beta)
coeffs = sp.Poly((equation1-equation2),sp.Derivative).coeffs()
dx1 = coeffs[-1]/coeffs[0]

在得到dx1之后,我可以反解dx0,最后用我通过实验设计得到的所有符号来模拟一个叫做多元delta定理的东西。

现在的问题是,我不仅有expr1,expr2,还有expr3。将来,我还会有expr4、expr5、expr6。我需要从这组方程中提取符号系数,这些方程彼此相等,才能同时反解dx5、dx4、dx3、dx2、dx1。

sp.Poly(sp.Eq(expr1,expr2,expr3,expr4,expr5,expr6),sp.Derivative).coeffs

上述方式显然会抛出一个错误:

TypeError:new((接受2到3个位置参数,但;N〃;被给予

所以问题又来了,会有一个数学解决方案或从sympy到"的方法吗;使相等;不求解的表达式集。sympy支持这一点吗?

与其声明它们相等,不如从另一个中减去一个,并将其用作solve()列表中的新方程!

>>> from sympy import *
>>> x, y, z = symbols("x y z")
>>> expr1 = x + y**2
>>> expr2 = z + y**3
>>> expr3 = expr1 - expr2
>>> expr3
x - y**3 + y**2 - z
>>> solve([expr3], x)
{x: y**3 - y**2 + z}

最新更新