Sympy -power 变量方程失败债券收益率

如何解决Sympy -power 变量方程失败债券收益率

我是 sympy 的新手，正试图让它为我解方程。我一直在使用完整方程时遇到问题，因此已将其简化，我相信挑战在于使用 -ve 变量。例如

y=symbols('y')
eqn=Eq(4*2.71828**(y*5),200)
solve(eqn,y)

以上工作正常，并给出了 0.782405127372565 的答案

然而，我真正想做的是有一个 -y，如下所示。这没有给我答案。有什么想法吗？

y=symbols('y') 
eqn=Eq(4*2.71828**(-y*5),y)

解决方法

2.71828**(y*5) 给一个搜索精确解的符号库看起来很糟糕。请使用exp(-y*5)。一般来说，最好尽可能avoid floats。

为了避免复数作为解决方案，您可以将 y 限制为 real。 只要有可能，sympy 就会返回精确的符号解。您可以通过 N() 或 .evalf() 将它们转换为数字。

from sympy import Eq,symbols,solve,exp,N

y = symbols('y',real=True)
eqn = Eq(4 * exp(-y * 5),200)
sol = solve(eqn,y)
print(sol,[N(s) for s in sol])

输出：

[log(2000**(1/5)/10)] [-0.782404601085629]

没有 real=True，您将获得所有 5 个复杂的解决方案。