将方程重写为多项式

Question

from sympy import *
K, T, s = symbols('K T s')
G = K/(1+s*T)
Eq1 =Eq(G+1,0)

我想用 sympy 重写等式 Eq1 作为多项式：1+K+T*s==0

我该怎么做？

我花了几个小时搜索和尝试简化方法，但找不到优雅、简短的解决方案。

SymPy 中的实际问题：

from IPython.display import display 
import sympy as sp
sp.init_printing(use_unicode=True,use_latex=True,euler=True)
Kf,Td0s,Ke,Te,Tv,Kv,s= sp.symbols("K_f,T_d0^',K_e,T_e,T_v,K_v,s")
Ga= Kf/(1+s*Tv)
Gb= Ke/(1+s*Te)
Gc= Kf/(1+s*Td0s)
G0=Ga*Gb*Gc
G1=sp.Eq(G0+1,0)
display(G1)

如何告诉 Sympy 将等式 G1 重写为形状为 s^3*(...)+s^2*(...)+s*(...)+(...) 的多项式=...？

课本上的实际问题：http://i.imgur.com/J1MYo9H.png

它应该是什么样子：http://i.imgur.com/RqEDo7H.png

这两个方程是等价的。

Answer 1

您可以执行以下操作。

import sympy as sp
Kf,Td0s,Ke,Te,Tv,Kv,s= sp.symbols("K_f,T_d0^',K_e,T_e,T_v,K_v,s")
Ga= Kf/(1+s*Tv)
Gb= Ke/(1+s*Te)
Gc= Kf/(1+s*Td0s)
G0=Ga*Gb*Gc

丢掉分母

eq = (G0 + 1).as_numer_denom()[0]

展开方程并收集 s 次方的项。

eq = eq.expand().collect(s)

最终方程

Eq(eq, 0)
Eq(K_e*K_f**2 + T_d0^'*T_e*T_v*s**3 + s**2*(T_d0^'*T_e + T_d0^'*T_v + T_e*T_v) + s*(T_d0^' + T_e + T_v) + 1, 0)

将方程重写为多项式

Rewrite equation as polynomial

python

sympy