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)

将方程重写为多项式

1 个答案: