Python Sympy符号

如何解决Python Sympy符号

当我使用“ x”和“ z”作为符号时，此代码没有问题：

from sympy import *
x,z = symbols('x z')
y = -6*x**2 + 2*x*z**0.5 + 50*x - z
solve((diff(y,x),diff(y,z)))
y.subs({x: 5,z: 25})

但是当我使用“ q”和“ a”时，solve不会给我任何解决方案。

q,a = symbols('q a')
y = -6*q**2 + 2*q*a**0.5 + 50*q - a
solve((diff(y,q),a)))
y.subs({q: 5,a: 25})

如您所见，我使用“ subs”检查目标函数中是否没有错字。

更新：我使用“ Symbol”分别设置每个变量，但是再次使用“ q”和“ a”不起作用。

# This works
x = Symbol('x')
z = Symbol('z')
y = -6*x**2 + 2*x*z**0.5 + 50*x - z
solve((diff(y,z)))

# This does not work
q = Symbol('q')
a = Symbol('a')
y = -6*q**2 + 2*q*a**0.5 + 50*q-a
solve((diff(y,a)))

谢谢。

解决方法

知道了！

这完全取决于变量的字母顺序。

在您的第一个示例中，如果用x代替z，用z代替x，

内部求解将表达式发送到_solve中的函数sympy.solvers，该函数然后尝试求解您的方程式并且失败多次。

最后，最后一次尝试是尝试通过内部函数-sqrt(a) + q选取x - sqrt(z)或_ok_syms的符号，并使用按字母顺序对它们进行排序的参数来解决reversed或x - sqrt(z)即使没有此参数，它仍然可以，但是如果用x: sqrt(z)包裹起来，它会神奇地使您的示例以完全相反的方式工作）。

因此它确实将-sqrt(a) + q解析为a: q**2，并将50 - 10*sqrt(z)解析为-12*q + 2*sqrt(q**2) + 50。

在第一种情况下，它以易于解决的sqrt(q**2)结尾，在第二种情况下，它由于无法简化INSERT INTO而在executemany上丢失了。

来源： 在以下方面进行了大量测试： https://github.com/sympy/sympy/blob/master/sympy/solvers/solvers.py