Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Sympy computing the inverse laplace transform

I am having some trouble computing the inverse laplace transform of a symbolic expression using sympy. In matlab and in the book I am working from the expression s/(s^2 + w^2) transforms to cos(wt).

When I attempt to do this using sympy like so:

expression = s/(s**2+w**2)
Answer = sympy.inverse_laplace_transform(expression, s, t)

I get that

Answer = (-I*exp(2*t*im(w))*sin(t*re(w)) + exp(2*t*im(w))*cos(t*re(w)) + I*sin(t*re(w)) + cos(t*re(w)))*exp(-t*im(w))*Heaviside(t)/2

What am I doing wrong?

like image 733
SomeRandomPhysicist Avatar asked Dec 08 '16 15:12

SomeRandomPhysicist


1 Answers

Sympy assumes that w is complex-valued. The simpler approach is to provide the option real=True in the definition of the symbol.

s, t = sp.symbols('s, t')
w = sp.symbols('w', real = True)
expression = s/(s**2+w**2)

sympy.inverse_laplace_transform(expression, s, t)

cos(t*w)*Heaviside(t)

like image 190
Stelios Avatar answered Nov 09 '22 10:11

Stelios