如何解决如何使sympi排列公式使分子中的最高阶系数为1
我正在使用 sympy 和 python 来解决一些方程。我试图安排一个方程,使最高程度的系数为+1。例如下面的代码产生输出
Km,Kb,L,s,R,J,b = sym.symbols("Km,b")
G1 = Km / (L * s + R)
G2 = 1 / (J * s + b)
Msys = (G1 * G2) / (1 + G1 * G2 * Kb)
Msys = sym.expand(Msys)
Msys = sym.simplify(Msys)
Msys = sym.collect(Msys,s)
print(Msys)
#outputs Km/(J*L*s**2 + Kb*Km + R*b + s*(J*R + L*b))
我需要去掉 s**2 前面的系数(即设置为 1)。这通常是通过将顶部和底部除以 J*L 来完成的。我发现了一个名为 monic 的 sympi 函数,它应该可以做到这一点,但它只适用于分子。
这个问题的最佳解决方案是什么?
提前谢谢
解决方法
这不是一个理想的解决方案,因为它以我需要的形式获取公式。基本上,我将公式分为顶部和底部,然后乘以 1 / (J * L) 最高阶变量前面的系数。
Km,Kb,L,R,J,b,Gc,Ain,Amp,Hoi,s = sym.symbols("Km,s")
G1 = 1/ (L * s + R)
print(G1)
G2 = 1 / (J * s + b)
print(G2)
Msys = (Km * G1 * G2) / (1 + Km * G1 * G2 * Kb)
Msys = sym.expand(Msys)
Msys = sym.simplify(Msys)
Msys = sym.collect(Msys,s)
g3 = Msys * ((1/(J*L))/(1/(J*L)))
print(g3)
top,bot = Msys.as_numer_denom()
top = top * (1/(J*L))
print(top)
bot = bot * (1/(J*L))
bot = sym.expand(bot)
bot = sym.collect(bot,s)
print(bot)
Msys = top * 1/bot
虽然此解决方案有效,但由于其复杂性,它并不理想。我可能可以进一步开发它并使其成为我可以重用的功能,但必须有更简单的方法
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。