Loading [MathJax]/jax/output/HTML-CSS/jax.js

First time here? Check out the FAQ!

Ask Your Question
0

inverse laplace transforms of shifts

asked 12 years ago

calc314 gravatar image

I've been trying to get Sage to solve an ODE using Laplace transforms. Unfortunately, shifted functions appear to cause a problem for the inverse_laplace command. David Joyner's ODE book (dated 2008) says that Sage does not have the functionality to do this. Has this situation changed? Do I need to do something differently here to get this to work?

Here is a snippet of something I tried.

u=piecewise([[(0,3),0],[(3,infinity),exp(-2*(t-3))]],t)
u.laplace(t,s)
inverse_laplace(_,s,t)

The first two lines work properly, but the last line returns a formal result:

ilt(e^(-3*s)/(s + 2), s, t)
Preview: (hide)

2 Answers

Sort by » oldest newest most voted
1

answered 8 years ago

mforets gravatar image

updated 8 years ago

A simple implementation using the local heaviside function (see also generalized functions) is written in this question, and is such that:

var('t s')
u=piecewise([[(0,3),0],[(3,infinity),exp(-2*(t-3))]], var=t)
F(s) = u.laplace(t,s);
InverseLaplace(F, s, t)

produces t  e(2t+6)H(t3).

Preview: (hide)
link

Comments

InverseLaplace is not currently a built-in function.

calc314 gravatar imagecalc314 ( 8 years ago )
0

answered 12 years ago

kcrisman gravatar image

I guess Maxima can't do it.

sage: inverse_laplace?
def inverse_laplace(ex, t, s):
    r"""
    Attempts to compute the inverse Laplace transform of
    ``self`` with respect to the variable `t` and
    transform parameter `s`. If this function cannot find a
    solution, a formal function is returned.

(%i1) display2d:false;
(%o1) false
(%i2) f(s):=%e^(-3*s)/(s+2);
(%o2) f(s):=%e^-(3*s)/(2+s)
(%i4) ilt(f(s),t,s);    
(%o4) 'ilt(%e^-(3*s)/(s+2),t,s)

Does that answer things? If you want to file an enhancement request to Maxima, please do, maybe it's not so hard for them to implement; however, presumably there are arbitrarily hard ones like this that no system could do...

Preview: (hide)
link

Comments

Thanks for the help. It seems like the t-domain shifts and s-domain shifts that we teach in a standard ODE course should be reasonable to implement. Web searches indicate that there are some folks who have implemented their own code in maxima to handle this.

calc314 gravatar imagecalc314 ( 12 years ago )

Ooh, maybe that could be implemented upstream, then. Can you post some links?

kcrisman gravatar imagekcrisman ( 12 years ago )

The following appears to be a deeper patch for maxima: http://comments.gmane.org/gmane.comp.mathematics.maxima.general/33278

calc314 gravatar imagecalc314 ( 12 years ago )

And, this link appears to be for an approach that might be useful upstream: http://www.ma.utexas.edu/pipermail/maxima/2007/008424.html

calc314 gravatar imagecalc314 ( 12 years ago )

I'm not sure about the robustness of either of the approaches above, but they are interesting.

calc314 gravatar imagecalc314 ( 12 years ago )

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 12 years ago

Seen: 1,398 times

Last updated: Feb 19 '17