对于使用Sympy在python中的一个大学项目,我需要评估部分替换后的2var符号函数的输出,例如:
F= (x-1)/sqrt((x-10)^2+(y-10)^2) /x=1 => f=0
如果我使用的是1 var func,它可以很好地工作:
>>> from sympy import *
>>> x = symbols('x')
>>> F = 1-x
>>> F.evalf(subs={x: 1})
0.e-125但是在2var函数的情况下,subs似乎没有做正确的替换:
>>> from sympy import *
>>> x,y = symbols('x y')
>>> F = (x-1)/sqrt((x-10)**2+(y-10)**2)
>>> F.evalf(subs={x: 1})
(x - 1.0)*((x - 10.0)**2 + (y - 10.0)**2)**(-0.5)我期望F=0,因为x=1
发布于 2018-10-10 16:14:12
我不是同情专家,但当y的值使表达式为非零时,我发现了一个反例:
>>> F.evalf(subs={x:1,y:float("nan")})
nan当y为nan时,结果不是零。因此,不知道y值的表达式是不可能简化的(注意:交感愉快地计算F.evalf(subs={x:1,y:oo})为0,即使结果也应该是nan )。
发布于 2018-10-10 19:23:26
如前所述,如果y= nan创建一个数学异常,那么替换就不会发生。但我发现如果我用潜艇我就能得到我想要的:
>>> F.subs(x,1)
0现在,如果我明确地替换y= nan,那么我就得到了nan的答案。
>>> F.subs(y,nan)
nan我猜evalf只是在评估替换之前检查更多的异常。
https://stackoverflow.com/questions/52744190
复制相似问题