如何解决用sympy限制solveset的解的域
我试图对我想用 sympy 求解的变量施加约束。 (即设置一个需要尊重的域)。我想使用 solveset
(我已经看到其他使用 solve
的答案)。
假设我有以下 K_p
,我想为 K_p - x = 0
求解 t
:
from sympy import *
psi,mu_1,mu_2,tau,t,x = symbols("psi,x",positive = True,real = True)
K_p =(mu_1*mu_2*psi*t**2*tau**4 - 2*mu_1*mu_2*t*tau**2 - mu_1*psi*t*tau**2 + mu_1 - mu_2*psi*t*tau**2 + mu_2 + psi)/(mu_1*mu_2*t**2*tau**4 - mu_1*t*tau**2 - mu_2*t*tau**2 + 1)
我如何施加以下约束:
t < 1/(mu_1 * tau**2)
和 t < 1/(mu_2 * tau**2)
s_x = solveset(K_p - x,domain=S.Reals)
我已经试过了:
cset = ConditionSet(t,t < 1/(mu_1 * tau**2)).intersect(ConditionSet(t,t < 1/(mu_1 * tau**2)))
s_x = solveset(K_p - x,domain=cset)
int = Interval(-oo,1/(maximum(mu_1,mu_2)*tau**2))
s_x = solveset(K_p - x,domain=int)
但是 sympy 并没有真正解决它。
请注意,我设法通过这种愚蠢的蛮力方法获得了两种解决方案中的哪一种:
s_x = solveset(K_p - x,domain=S.Reals)
eq1 = s_x.args[0].args[1].args[0]
eq2 = s_x.args[0].args[1].args[1]
v_mu_1 = 2
v_mu_2 = 3
v_x = 50
v_psi = 10
v_tau =1
mint = min(1/(v_mu_1* v_tau**2),1/(v_mu_2* v_tau**2))
eq1.subs([(mu_1,v_mu_1),(mu_2,v_mu_2),(tau,v_tau),(x,v_x),(psi,v_psi)]) < mint
eq2.subs([(mu_1,v_psi)]) < mint
只有其中一个方程给了我应该的True
。我希望solveset的结果就是那个方程。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。