python-3.x - 如何使用 sympy 的 nonlinsolve 符号式地找到稳态

Question

我有一个标准的经济学问题，我想用 sympy 的 nonlinsolve 来解决。

我想象征性地找到模型的稳态并为此求解非线性方程组。

我很高兴看到 sympy 有 nonlinsolve 应该能够解决这些问题，但我到目前为止还不能这样做。

在下文中，您会找到无法自行解决的问题部分：

    from sympy.solvers.solveset import nonlinsolve
    from sympy import var

    var('c k mc n r_k eta w alpha', Rational = True, positive=True)

    eq1 = -r_k + alpha * mc * k**(-1 + alpha) * n**(1 - alpha)
    eq2 = -w + mc * (1 - alpha) * k**alpha * n**(-alpha)
    eq3 = c**-1 * w - n**eta

    eq_system = [eq1,eq2,eq3]

    nonlinsolve(eq_system , [n,k,w])

alpha 和 eta 是参数，而其他符号是变量。我知道它们是有理正数，但由于我还不明白的原因，代码抛出：

ValueError：alpha：指数必须是正整数

我也不明白为什么 alpha 必须是整数，如果我这样声明 alpha 也没有帮助。

关于如何使用 sympy 获得有意义的解决方案的任何提示？

非常感谢，

托尔

标签： python-3.xsympysymbolic-mathsolvereconomics