Question

我使用sympy对一组幂律进行非线性求解，最后得到了一堆像

这样的解决方案

20394*sqrt(6)*x**(2/3)*y**(1/4)

我希望在暴露指数的合理性质的同时折叠数值系数，即我希望我的代码输出

49954.8938143201*x**(2/3)*y**(1/4)

然而，当我使用evalf（）时，我得到了

49954.8938143201*x**0.666666666666667*y**0.25

有没有办法让同事不评估指数，甚至只是提取答案的非符号部分，这样我就可以创建自己的打印功能，以我想要的方式格式化输出？

Answer 1

一种方法是选择所有Mul和Add子表达式，将它们的参数分成数字和符号（基于是否包含任何符号），并使用数字参数对该部分进行数值计算。

expr = 20394*sqrt(6)*x**(S(2)/3) + 3*log(3)*y**(S(3)/2) - 3*sqrt(7) + (3+sqrt(2))*x*y
eval_dict = {}
for a in expr.atoms(Mul, Add):
    numeric = [arg for arg in a.args if not arg.has(Symbol)]
    symbolic = [arg for arg in a.args if arg.has(Symbol)]
    eval_dict[a] = a.func(a.func(*numeric).evalf(), a.func(*symbolic)) 
print(expr.subs(eval_dict))

打印

49954.8938143201*x**(2/3) + 4.41421356237309*x*y + 3.29583686600433*y**(3/2) - 7.93725393319377

理性指数逃避数值评估，是Rational的实例。

在系统中评估系数但不是指数

1 个答案: