Question

在我的代码中的某个时刻，我在一个sympy符号上添加了一个上限。 sympy.refine(x, sympy.Q.positive(upper_bound - x))

我会从符号“x”中检索programmaticaly上位数（和其他约束）的值。任何想法如何实现这一目标？

欢呼声

Answer 1

细化函数不起作用 - 它只是简化给定约束的表达式。例如：

In [54]: import sympy

In [55]: x = sympy.S('x')

In [56]: expr = sympy.S('sqrt(x**2)')

In [57]: sympy.refine(expr, sympy.Q.positive(x))
Out[57]: x

In [58]: sympy.refine(expr, sympy.Q.negative(x))
Out[58]: -x

In [59]: sympy.refine(expr, sympy.Q.real(x))
Out[59]: Abs(x)

不幸的是，使用不等式似乎没有做任何有用的事情：

In [62]: sympy.refine(expr, sympy.Q.is_true(x>0))
Out[62]: sqrt(x**2)

你可以用solveset做一些有用的事情：

In [68]: expr = sympy.S('A * x**2 + B * x + C')

In [69]: sympy.solveset(expr, x, sympy.Interval(1,10))
Out[69]: ConditionSet(x, Eq(A*x**2 + B*x + C, 0), [1, 10])

或许是一个更有用的例子：

In [19]: a = sympy.S('(x**2)*(sin(x)+x)')

In [20]: x = sympy.S('x')

In [25]: b = sympy.solveset(a,x,sympy.Interval(-2,2))

In [26]: b
Out[26]: ConditionSet(x, Eq(x + sin(x), 0), [-2, 2])

In [34]: b.base_set
Out[34]: [-2, 2]

如何查询sympy约束？

1 个答案: