使用 sympy 求解表达式

如何解决使用 sympy 求解表达式

我想用 sympy 解决下面的表达式。 我已经定义了所有的符号和函数。 c,x,R 是函数，f 是我想为 p 求解的方程。我将有两个可能的 p ，我想绘制如下所示。 p 也在 0 和 1 之间，我没有放那个，也不知道如何定义。 x 是 c 函数 wrt k 的导数，我把它的倒数放在 f 中。

代码：

import sympy as sym
sym.init_printing()

p,q,a,c0,c1,z,k = sym.symbols('p,k')

c = sym.Eq(c0+0.5*c1*(1/k**2))

x = sym.diff(c,k)

R = sym.Eq(z*(a+x**(-1)))

f = sym.Eq((1-p+P*q)*(a+x**(-1))-c-(1-p)*R)

sym.solve([f],(p))

解决方法

也许看到一个具体的工作示例会帮助您理清问题陈述。您正在求解的表达式在 p 中是线性的，因此只有一个解。让 e1 和 e2 是 x（而不是 p）的两个函数。如果你想知道它们在哪里交叉，你需要知道什么时候 e1 == e2 或者等价地，什么时候 e1 - e2 = 0:

>>> from sympy import solve,solveset,Interval
>>> from sympy.abc import x
>>> e1 = 1/x
>>> e2 = (5*x + 9)/(25*x - 5)
>>> solve(e1 - e2)
[8/5 - sqrt(39)/5,sqrt(39)/5 + 8/5]  <--- two solutions but one is larger than 1
>>> solveset(e1-e2,x,Interval(0,1))
FiniteSet(8/5 - sqrt(39)/5)  <--- single solution in range of interest