问题描述
我试图对我想用 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 (将#修改为@）