7

我有以下表达式:

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)

目前我手动替换bjcj在简化u0表达式中得到:

u0 = a2*(v*sin(bj) + cj*vt*cos(bj) + 2*cj*w*cos(bj))

是否可以使用 SymPy 来实现这一点,避免手动替换?

4

2 回答 2

11

我猜你缺少的是subs它将替换任意表达式,而不仅仅是符号

>>> print simplify(u0).subs({pi*j*(phi + t)/(2*phi): bj, j*pi/(2*phi): cj})
a2*(pi*j*vt*cos(bj) + 2*pi*j*w*cos(bj) + 2*phi*v*sin(bj))/(2*phi)

(我之所以使用simplify它,是因为这导致了pi*j*(phi + t)/(2*phi)而不是pi*j/2 + pi*j*t/(2*phi),但不是其他要求)

阅读http://docs.sympy.org/0.7.3/tutorial/basic_operations.html#substitution有关替换和替换的更多信息。如果你想做更高级的替换,看看replace方法。

于 2013-08-02T19:23:28.967 回答
1

cse您可以在例程中找到常见的子表达式。

于 2013-08-02T08:51:46.553 回答