Sympy - solve不能完全解决方程系统(尽管方程系统很容易求解)
Sympy是一个Python库，用于符号计算和数学建模。它提供了一组功能强大的工具，可以处理代数方程、微积分、离散数学等数学问题。
在解决方程系统时，Sympy的solve函数可以用于求解代数方程组。然而，有时候即使方程系统很容易求解，Sympy的solve函数可能无法完全解决方程系统。这可能是由于以下原因之一：
- 复杂度：某些方程系统可能非常复杂，包含大量的未知数和方程。在这种情况下，Sympy的solve函数可能无法找到解析解或者需要很长的计算时间。
- 数值解：有些方程系统可能只能通过数值方法求解，而不是解析方法。Sympy的solve函数主要用于求解解析解，对于数值解可能不太适用。
- 特殊情况：某些方程系统可能具有特殊的性质或约束条件，这可能导致Sympy的solve函数无法找到所有解或者找到错误的解。
在这种情况下，可以考虑使用其他数学建模工具或数值计算库来解决方程系统。例如，可以使用NumPy、SciPy或者其他数值计算库来进行数值求解。另外，还可以考虑使用专门的数学建模软件，如MATLAB或Mathematica。
总结起来，尽管Sympy的solve函数在大多数情况下可以有效地解决方程系统，但在某些复杂或特殊情况下可能无法完全解决。在这种情况下，可以考虑使用其他数学建模工具或数值计算库来解决方程系统。
相关·内容
我正在使用Sympy (版本0.7.3，Python 2.7.5，Mac OS X)求解器来求解一些矩阵方程，我遇到了不一致的行为，这会破坏我的整个结果。准确地说，结果似乎取决于应该求解的方程的顺序。
以下代码生成了一个最小的工作示例：
from sympy import *
axx, bxx, byy = symbols('axx bxx byy')
此命令
solve([axx - bxx, byy])
产生{axx: bxx, byy: 0}作为结果，而切换方程的顺序
solve([byy, axx - bxx])
给出了{byy: 0, bxx: axx}，这在数学上当
浏览 1提问于2013-08-29得票数 1
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
回答已采纳
我有一个问题，就是渐近模没有解下面的方程。
import sympy as sy
from sympy import Symbol
from sympy import *
deld = Symbol('deld')
sinh = sy.sinh
a = (-0.00026478*deld - 1.75012e-9*sinh(0.62831 - 62831.85307*deld) + 1.750123e-9*sinh(62831.85307*deld - 0.62831) + 2.64786e-9)
b = (0.00026478*deld + 1.75012e-9*sinh(
浏览 10提问于2022-07-20得票数 1
回答已采纳
5回答
我正在尝试使用SymPy在python中求解一个方程。我有一个生成的方程式(类似于function = y(8.0-(y**3.0))，我在SymPy中使用它来创建一个新的方程式，如下所示:输出y(8.0-(y**3.0)) == 2的eq = sympy.Eq(function, 2)。但是sympy.solve(eq)似乎不起作用。
>>> from sympy import Eq, Symbol as sym, solve
>>> y = sym('y')
>>> eqa = Eq(y(8.0-(y**3.0)), 8
浏览 1提问于2010-12-15得票数 7
1回答
我使用下面的代码来求解一个三次方程。
from sympy.solvers import solve
from sympy import Symbol
x = Symbol('x')
print(solve(-0.0643820896731566*x**3 + 0.334816369385245*x**2 + 1.08104426781115*x - 2.05750838005246,x))
由于它是一个具有实数系数的三次方程，所以不可能有三个不同的复根。但它给出了以下结果。
[-3.19296319480108 - 0.e-22*I, 1.43925417946882 + 0
浏览 2提问于2019-04-12得票数 4
回答已采纳
给出了这个例子：
from sympy.core.symbol import symbols
from sympy.solvers.solveset import nonlinsolve
x, y, z = symbols('x, y, z', real=True)
nonlinsolve([x*y - 1, 4*x**2 + y**2 - 5], [x, y])
{(-1, -1), (-1/2, -2), (1/2, 2), (1, 1)}
但即使是在他们网站上的实时外壳中，这也会导致一个错误：
>>> from sympy.solvers.solvese
浏览 2提问于2016-11-24得票数 8
回答已采纳
1回答
我试着用渐近来解决下面的三角方程，但是它给了我这个错误：
NotImplementedError: could not solve b*tan(3*atan(6/b)/2) - 9
有没有办法让我用渐近解出下面的方程？
import sympy
from sympy import *
a = sympy.symbols("a")
b = sympy.symbols("b")
eq1 = sympy.Eq(b*tan(2*a), 6)
eq2 = sympy.Eq(b*tan(3*a), 9)
result = sympy.solve([eq1, eq2]
浏览 47提问于2021-02-19得票数 1
回答已采纳
2回答
我正试图解出下列方程式：
cos(a) - cos(b) -1=0
和
sin(a) - sin(b) =0
使用SymPy，我得不到(a,b)的解决方案。
from sympy import symbols, sin, cos, solve, Eq
a, b = symbols("a b")
eq1 = Eq(cos(a) - cos(b) - 1, 0)
eq2 = Eq(sin(a) - sin(b), 0)
solve([eq1, eq2], [a, b])
然而，至少有一个解决办法：
A= 60°
B= 120°
为什么SymPy不能解决这组三角方程？
浏览 19提问于2022-08-16得票数 1
回答已采纳
我正在使用python2.7和同情。我想要解一个方程，但考虑到变量的值是受约束的(并不是整数变量的所有值)，例如-25 <= a <= 25。
我该怎么做？
我的代码
from sympy import *
from sympy import sin
from sympy.solvers import solve
a, b, c, d = symbols('a b c d')
print solve(cos(a)*sin(b) - sin(a)*cos(b)*cos(d - c), a, b, c, d)
浏览 4提问于2015-07-03得票数 0
我有下面的方程系统
from sympy.abc import x,y,z
from sympy import solve, Eq
eqs = [Eq(x,y), Eq(y,z)]
result = solve(eqs, x)
现在我希望能够为z替换一个值，并为x求值。我该怎么做呢？
浏览 0提问于2014-11-07得票数 1
我想用数值方法解决这个方程
beta.ppf(x，a，b)-c=0
其中a，b，c是已知的常数。当我试着
from sympy.solvers import solve
from sympy import Symbol
from scipy.stats import beta
x=Symbol('x')
solve(beta.ppf(x,a,b)-c,x)
它还给了我
TypeError:无法确定关系的真值
我怎么才能修好它？
浏览 2提问于2016-08-28得票数 0
我有个方程要解。方程可描述为上述公式。N和S是常数，例如N= 201和S= 0.5。我在蟒蛇中用同情来解决这个问题。python脚本如下所示：
from sympy import *
x=Symbol('x')
print solve( (((1-x)/200) **(1-x))* x**x - 2**(-0.5), x)
然而，有一个RuntimeError:在__instancecheck__中超过了最大递归深度。
我也尝试过使用数学，它可以输出0.963的结果。
任何建议都是欢迎的。谢谢。
浏览 4提问于2017-10-24得票数 0
回答已采纳
我正在试着用SymPy的dsolve用ics解下面的微分方程 ? from sympy import Function, Derivative, dsolve, symbols
t, k = symbols('t,k', real=True)
M = Function('M')(t)
M_ = Derivative(M,t)
Eqn = M_ - k*M
sol = dsolve(Eqn, ics={M.subs(t,0): 100, M.subs(t, 6): 97}) 我得到以下错误： File "..\Python\Python39\l
浏览 116提问于2020-12-26得票数 0
1回答
解决()执行中止
、、
我很难用Sympy来解方程。例如，当我运行代码时：
print(correction(10))
我希望它打印一个数字，它是f。相反，它给了我错误:执行中止。
def correction(r):
from sympy import cosh, log, exp, symbols, solve
f = symbols('f')
def equation():
return cosh(((r - 1.0)/(r + 1.0))*(log(2.0)/f)) - 0.5*exp(log(2.0)/f)
correction = so
浏览 2提问于2016-12-14得票数 2
回答已采纳
1回答
我想通过使用r重新排列变量的下面的等式。
P = C * ((1-(1+r)**(-n)))/r + fv*(1+r)**(-n)
to
r = blabla...
我明白，同情与这样的重新安排任务有关。所以我写了下面的代码。
# Solve the equation for r
import sympy
from sympy import symbols
P, C, r, n, fv = sympy.symbols('P C r n fv')
eq = sympy.Eq(P, C * ((1-1/(1+r)**n))/r + fv/(1+r)**n)
sympy.solve
浏览 2提问于2020-04-13得票数 1
回答已采纳
1回答
当我用交感解数学问题时遇到了一些问题，我的代码是这样的，
eq1=Eq(0, -sympy.sqrt(x**2) + sympy.sqrt((x + 4)**2))
sympy.solve(eq1,x)
这个方程的答案应该是-2，但我得到了一个空的列表。
浏览 0提问于2018-09-13得票数 2
回答已采纳
2回答
我需要解一个方程，但是蜘蛛返回了这个错误：
import sympy as sym
import sympy as sympy
Re=100
Epsilon=0.00000075
D=0.01
f=symbols('f')
eq=(((1/(-2* sympy.log((Epsilon/D/ 3.7)+( 2.51 / (Re*(f**0.5))), 10)))**2 -f),0)
sym.solve((eq), (f))
print(f)
NotImplementedError: could not solve -4*f*log(20270270270270
浏览 1提问于2022-03-17得票数 0
我正在试着解x和y的联立方程，我没有得到任何结果(代码一直在运行)。我觉得这个错误与在方程式中使用sqrt有关，但不确定。有人能帮我解决这个问题吗？
from __future__ import division
from sympy import Symbol,sqrt,solve
x = Symbol('x')
y = Symbol('y')
z = Symbol('z')
a = Symbol('a')
b = Symbol('b')
c = Symbol('c')
d = Symbol(&#
浏览 0提问于2018-07-19得票数 4
我想用符号系数来求解非常简单的方程：
from sympy import *
a, b = symbols('a b', commutative=False)
x = IndexedBase('x')
basis = [a, b, a * b - b * a]
el = b * a - a * b
coefs = [x[k] for k in range(len(basis))]
eq = el - sum([c * bel for c, bel in zip(coefs, basis)])
solve(eq.expand(), coefs)
方程是-x[0]
浏览 1提问于2018-06-20得票数 4
我试图通过修改找到的代码来解决问题。
import sympy as sp
t = sp.symbols('t', real=True)
x = sp.Function('x', real=True)
diffeq = sp.Eq(x(t).diff(t), sp.sqrt(1 - x(t)**2))
res = sp.dsolve(diffeq, ics={x(0): 0})
math.se文章建议，分段函数满足方程，但也许我的代码没有添加所有必要的假设。我得到以下回溯：
Traceback (most recent call last):
File
浏览 11提问于2022-10-22得票数 0
回答已采纳
2回答
我有两个耦合的非线性微分方程组。
为此，我编写了以下代码，以使用Sympy获取两个函数的值：
from sympy import *
t = symbols('t')
c, b, B, alp, mu = symbols('c b B alp mu', integer=True)
f, g = symbols('f g', cls=Function)
print(solve([Eq(f(t).diff(t)+c*f(t)+0.5*B*f(t)**2-b*g(t), 0), Eq(g(t).diff(t)+b*g(t)-c*f(t), (1-
浏览 1提问于2017-01-14得票数 0
1回答
当我试图用RSymPy库求解一个方程时，我得到了NotImplementedError:多元方程误差。
我不知道我错了什么地方，或者我应该使用其他的多元方程库，或者有一个解决办法来使它工作？
此外，如果不对其进行硬编码，我也无法将变量形式的值传递给方程参数(ymax、ymin和ec50)。
library(rSymPy)
ymax <- Var("ymax")
ymin <- Var("ymin")
ec50 <- Var("ec50")
sympy("var('x y')")
sympy(
浏览 4提问于2017-02-24得票数 0
回答已采纳
2回答
问题是:我试图解一个二阶的形式的矩阵方程：
其中X(查找)和C(已知)的大小为nxn。(n为1000左右)。C是已知的对称协方差矩阵。(X也应该是对称的)
这是我的代码：
from sympy import solve
from sympy import Indexed, IndexedBase, Tuple
import numpy as np
X = IndexedBase('X',shape=(n,n))
eqs = Tuple(np.dot(X,X)-np.dot(C,X)-np.eye(n))
solve(eqs, X)
这样做对吗？我的密码需要很长时间。我正
浏览 6提问于2017-12-02得票数 1
回答已采纳
3回答
当使用solve计算二次方程的根时，SymPy返回可以简化的表达式，但我无法使其简化。一个最小的例子如下所示：
from sympy import *
sqrt(-24-70*I)
在这里，SymPy只返回sqrt(-24-70*I)，而Mathematica或Maple将回答与5-7*I等价的问题。
我知道有两个平方根，但是这种行为意味着SymPy将返回非常复杂的解决方案
z = symbols("z")
solve(z ** 2 + (1 + I) * z + (6 + 18 * I), z)
同时，Maple和Mathematica都会很高兴地给出两个高斯整数来解决这个方
浏览 8提问于2016-12-25得票数 6
回答已采纳
我想用交感来进行电子设计计算。我有一个方程，我对某些电阻值进行求解。我只能使用标准电阻值，所以我再次在方程中设置了选择值，并得到了最终结果。
'''Example calculations for LMZ22010 switching regulator'''
from IPython import get_ipython
get_ipython().magic('reset -sf')
from sympy.solvers import solve
from sympy import Symbol, Eq, symbols,
浏览 2提问于2018-04-13得票数 1
回答已采纳
1回答
SymPy为我提供了复杂的解决方案，即使我将变量约束为实数，就像在this answer中一样。 这里有一个寻找单位球面和平面X=1的交集的例子。应该恰好有一个实数解，(1,0,0)。 >>> import sympy
>>> x,y,z = sympy.symbols('x,y,z', real=True)
>>> list(sympy.solve([x**2 + y**2 + z**2 - 1, x - 1], [x,y,z]))
[(1, -I*z, z), (1, I*z, z)] 这些都是正确的，但不是真的。 我
浏览 22提问于2020-08-26得票数 0
我要用“n-解数”函数来解这个方程。不幸的是，出现了此错误：
ValueError: Could not find root within given tolerance. (435239733.760000060718 > 2.16840434497100886801e-19)
Try another starting point or tweak arguments.
守则是：
import sympy
d=[0.3, 32.6, 33.4, 241.7, 396.2, 444.4, 480.8, 588.9, 1043.9, 1136.1, 1288.1, 1408.1, 143
浏览 6提问于2022-06-07得票数 0
2回答
SymPy可以方便地求解简单系数较短的二次方程。例如：
from pprint import pprint
from sympy import *
x,b,f,Lb,z = symbols('x b f Lb z')
eq31 = Eq((x*b + f)**2, 4*Lb**2*z**2*(1 - x**2))
pprint(eq31)
sol = solve(eq31, x)
pprint(sol)
但系数稍大一点，就不能：
from pprint import pprint
from sympy import *
c3,b,f,Lb,z = symbols('c3
浏览 1提问于2014-07-22得票数 1
正如推荐的那样，我尝试使用sympy.solveset而不是sympy.solve来解决一组两个方程。 from sympy import *
a1, a1dot, pa1, pa1dot, a2, a2dot, pa2, pa2dot, m, l = symbols('a1, a1dot, pa1, pa1dot, a2, a2dot, pa2, pa2dot, m, l')
eq1 = 2*m*l**2*a1dot + m*l**2*a2dot*cos(a1 - a2) - pa1 #equation one
eq2 = m*l**2*a2dot + m*l**2*a1
浏览 27提问于2021-08-06得票数 1
回答已采纳
在使用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
回答已采纳
2回答
我试图让SymPy解决一个方程系统，但它给我一个错误，如下所示：
NotImplementedError: could not solve 3*sin(3*t0/2)*tan(t0) + 2*cos(3*t0/2) - 4
有没有别的方法可以让我解方程组：
sin(x)+(y-x)cos(x) = 0
-1.5(y-x)sin(1.5x)+cos(1.5x) = 2
我使用：
from sympy import *
solve([sin(x)+(y-x)cos(x), -1.5(y-x)sin(1.5x)+cos(1.5x)-2], x, y)
浏览 19提问于2018-01-12得票数 3
回答已采纳
我试图得到弗洛伊登斯坦方程(Psi)的一个角度：
k1 * cos(psi) - k2 * cos(fi) + k3 - cos(psi - fi) =0
我有k1，k2，k3和fi的值。我尝试了以下几种方法：
from sympy import Symbol, solve, cos
x = Symbol('x')
realPsi = solve(k1 * cos(x) - k2 * cos(fi) + k3 - cos(x - fi), x)
我收到这个错误：
File "/usr/lib/python2.7/dist-packages/sympy/solvers/
浏览 1提问于2013-05-28得票数 1
回答已采纳
我想用sympy.solve来解一个变量的方程。方程中有‘np.deg 2rad’，'np.tan‘和'np.pi’。例如，我希望解a的方程'a-np.tan(np.deg2rad(b))*np.pi‘，这是我的代码：
import numpy as np
import math
import sympy as sy
a,b = sy.symbols('a b')
sy.solve(a-np.tan(np.deg2rad(b))*np.pi, a)
当我尝试时，我会得到以下错误：
AttributeError: 'Symbol' obj
浏览 10提问于2022-11-07得票数 0
所以我和Sympy一起努力构建一个通用的物理问题求解器/生成器。其中一个部分是，我要找一个函数，它将取kwargs，根据它得到的，重新排列方程，并替换它中的值。幸亏如此，我终于找到了我需要的东西。
但是.我尝试过将sympy.solve放在一个for循环中来生成所有这些表达式，然后我遇到了.某物。
import sympy
R, U, I, eq = sympy.symbols('R U I eq')
eq = R - U/I
for x in 'RUI':
print(x)
print(sympy.solve(eq, x))
输出？
R
[U/
浏览 2提问于2016-03-28得票数 3
回答已采纳
我有一个联立方程计算器，它包含两个方程的字符串，它应该输出到X和Y的解。这是它的代码。
x, y = symbols('x,y')
transformations=(standard_transformations + (implicit_multiplication_application,))
eqs_sympy = [Eq(parse_expr(e.split('=')[0], transformations=transformations),
parse_expr(e.split('=')[1], transfo
浏览 2提问于2021-08-08得票数 1
回答已采纳
2回答
我有一个方程，我需要解一个分数。我有更复杂的问题要解决，但下面是一个很小的例子:以下面的简单函数Y = X*a为例。我想为Y/X解决问题，所以我期望Y/X =a。这是代码，它生成一组空的答案
from sympy import *
X,Y,a = symbols('X Y a')
testEq = Eq(Y,X*a)
solve(testEq,Y/X)
我想我是误会了，任何帮助都很感激！
浏览 3提问于2017-12-17得票数 1
回答已采纳
1回答
使用Python和SymPy，我试图为b求解这个方程：
我的代码：
From sympy import *
b = Symbol('b')
x = ((b/(2*math.pi*math.e))*((math.pi*b)**(1/b)))**(1/(2*(b-1)))
solve(x-1.0034,b)
我得到了一个错误: NotImplementedError:实现了多个generators...No算法来求解方程。
你知道哪里有错误吗？或者说，这个等式是如此的困难以至于Python无法解决它？谢谢
浏览 0提问于2020-03-24得票数 0
回答已采纳
1回答
求解变量子集的多元方程
、、、
我正在使用渐近来解决一些方程，我遇到了一个问题。我对许多方程都有这个问题，但我会用一个例子来说明。我有一个有多个变量的方程，我想用所有变量来求解这个方程，但其中一个变量被排除了。例如方程0 = 2^n*(2-a) - b + 1。这里有三个变量a、b和n。我想要获取a和b的值，而不是n的值，这样a和b就不能包含n。 2^n*(2-a) - b + 1 = 0
# Since we don't want to solve in terms of n we know that (2 - a)
# has to be zero and -b + 1 has to be zero.
2
浏览 22提问于2019-01-02得票数 3
1回答
我正在尝试求解一组方程，然后将结果存储在变量中，以便绘制或进一步计算。
%pylab inline
from sympy import init_printing;init_printing()
from sympy import *
d,vf,a,vi,t,x,h,g,theta,ans=symbols('d vf a vi t x h g theta ans')
equations=[Eq(sin(theta),(0.5*g*t**2+h)/(vi*t)),Eq(cos(theta),x/(vi*t))]
ans=solve(equations,[h,t],dict=Tru
浏览 0提问于2013-07-23得票数 0
我试图解一个带有对数的简单方程，但是当我试图执行代码时，我得到了一个RuntimeError: maximum recursion depth exceeded。以下是我正在做的事情：
import sympy as sp
import numpy as np
pH = sp.Symbol("pH")
pCO2 = sp.Symbol("pCO2")
HCO3 = sp.Symbol("HCO3")
myhco3 = 10.0**(7.0 - 6.1 + np.log10(0.03 * 44.0))
hh = sp.Eq(pH, 6.1
浏览 1提问于2016-07-04得票数 4
回答已采纳
3回答
我试图使用下面的python代码(当然是使用SymPy )来求解一组方程：
def Solve(kp1, kp2):
a, b, d, e, f = S('a b d e f'.split())
equations = [
Eq(a+b, 2.6),
Eq(2*a + b + d + 2*f, 7),
Eq(d + e, 2),
Eq(a*e,kp2*b*d),
Eq( ((b * (f**0.5))/a)*((a+b+d+e+f+13.16)**-0.5), kp1)
]
return
浏览 4提问于2015-04-13得票数 5
回答已采纳
1回答
我是个蟒蛇初学者，我正试着解决一个方程式。下面是方程，我正试图为v求解：
**fSw - ((1)/(1+(Kro*Meu_w)/(Krw*Meu_nw))) = 0**
我想做的是设置一个Meu_w和Meu_nw值，并为Kro和Krw值进行求解。我尝试过制作Krw和Kro数组，但是当我运行我的代码时，我得到的唯一返回是"[]“。我认为这是一个问题，因为我试图将Kro和Krw的值集提供给方程，即，我想要的是在将Krw和Krw的数组输入方程之后得到多个fSw值，但我不知道如何做到这一点。
我的代码如下，
from sympy.solvers import solve
from sympy
浏览 6提问于2022-11-08得票数 0
回答已采纳
在二维空间中，我有一个泊松方程，如下所示：
下面是我解决这个问题的尝试：
import sympy as sp
x, y = sp.symbols('x, y')
f = sp.Function('f')
u = f(x, y)
eq = sp.Eq(u.diff(x, 2) + u.diff(y, 2), u)
print(sp.pdsolve(eq))
它给出了一个错误：
psolve: Cannot solve -f(x, y) + Derivative(f(x, y), x, x) + Derivative(f(x, y), y, y)
这样
浏览 0提问于2018-10-31得票数 0
回答已采纳
1回答
from sympy import solve, symbols
from math import e
a = symbols('a')
r = solve(1604+14*a**(-1)-14x1604*(1-e**(-1604*a))**-1, a)
print(r)
当我试图解决这个方程时，程序似乎一直在运行(持续了4个小时)，没有解。有人能帮忙吗？
浏览 2提问于2022-06-03得票数 -1
1回答
我试图在python中求解一个方程组(带有渐近性)，但我得到了"TypeError:无效NaN比较“。这些方程是“几乎线性的”(线性方程组，但也包含一个abs部分)。
对于函数f，我试图找到两点x0，x1，使得f( x0 ) = x1和f(x1) =x0(因此这两个方程形成一个系统)。
我有：
from sympy import *
c = 2
x = Symbol('x', real=True)
f1 = c * (1/2) * (Abs(2*x-1)) # define a function with variable 'x', also Abs
浏览 6提问于2022-10-27得票数 0
回答已采纳
动机加泰罗尼亚数的生成函数满足二次方程是众所周知的。我想要一个函数的前几个系数，由一个代数方程(不一定是二次方程！)隐式定义。
示例.
import sympy as sp
sp.init_printing() # math as latex
from IPython.display import display
z = sp.Symbol('z')
F = sp.Function('F')(z)
equation = 1 + z * F**2 - F
display(equation)
solution = sp.solve(equation
浏览 0提问于2017-09-26得票数 2
回答已采纳
1回答
我正在试着解决关于x的如下方程：
这里给出了alpha's和K，N将大于1000。有没有一种方法可以指定LHS，给出一个使用alpha的np.array？我的希望是定义：
eqn = Eq(LHR - K)
solve(eqn,x)
通过告诉渐近那个LHS=和( a_i + x)。
任何关于求解器的技巧都会以最快的速度做到这一点，也将不胜感激。谢谢!
我希望是这样的：
from sympy import Symbol, symbols, solve, summation, log
import numpy as np
N=10
K=1
alpha=np.random.rand
浏览 2提问于2013-02-05得票数 2
2回答
我有一个方程式如下
y = (x^3-1)((exp(0.019/x))*1.11 -1)
现在，对于y的已知值，如何得到x值？
我试着喜欢这个
from sympy import *
x=Symbol('x')
y=((x**3)-1)*((2.718**(0.019/x))*1.1 -1)
z=solve(y - 0.11,x)
错误：
NotImplementedError: multiple generators [x, 1359**(1/(1000*x)), 500**(1/(1000*x))]
No algorithms are implemented t
浏览 2提问于2020-10-01得票数 0
回答已采纳
新来的同情，我真的看不出我做错了什么。我正在解三角形的面积，已知的高度，但要求是定义两个符号。
def solve_math_problem():
from sympy import symbols, Eq, solve
from fractions import Fraction
base, height = symbols('base, height')
eq1 = Eq(Fraction(1, 2) * base * height, 24)
eq2 = eq1.subs(height, 8)
result = solve([eq1, eq2], base)
print
浏览 3提问于2021-06-16得票数 0
回答已采纳
1回答
我玩过的游戏有一个谜语，它涉及到解下列方程式：
x*411 + y*295 + z*161 = 3200
我不想认为我刚刚把它塞进了sympy，我还没有真正用过它：
>>> from sympy import *
>>> x, y, z = symbols('x y z', integer=True, positive=True)
>>> solve(x*411 + y*295 + z*161 - 3200, [x, y, z])
[{x: -295*y/411 - 161*z/411 + 3200/411}]
嗯，这只给了
浏览 1提问于2017-09-02得票数 6
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
