我有以下表達式:SymPy,簡單化/置換使用已知的圖案或子表達式
from sympy import pi, sin, cos, var, simplify
var('j,u,v,w,vt,wt,a2,t,phi')
u0 = v*a2*sin(pi*j/2 + pi*j*t*phi**(-1)/2) + pi*vt*a2*cos(pi*j/2 + pi*j*t*phi**(-1)/2)*j*phi**(-1)/2 + pi*w*a2*cos(pi*j/2 + pi*j*t*phi**(-1)/2)*j*phi**(-1)
哪些可以簡化爲:
print simplify(u0)
#a2*(pi*j*vt*cos(pi*j*(phi + t)/(2*phi)) + 2*pi*j*w*cos(pi*j*(phi + t)/(2*phi)) + 2*phi*v*sin(pi*j*(phi + t)/(2*phi)))/(2*phi)
鑑於子表達式:
bj = pi*j*(phi + t)/(2*phi)
cj = j*pi/(2*phi)
目前我手動替換SIM卡中的bj
和cj
plified u0
表達獲得:
u0 = a2*(v*sin(bj) + cj*vt*cos(bj) + 2*cj*w*cos(bj))
是否有可能使用SymPy實現這一目標,避免了手動替換?
這太棒了!謝謝! –