MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: numerical inversion of laplace transform

  • To: mathgroup at smc.vnet.net
  • Subject: [mg74861] Re: numerical inversion of laplace transform
  • From: "Roman" <rschmied at gmail.com>
  • Date: Mon, 9 Apr 2007 06:09:11 -0400 (EDT)
  • References: <ev7j3g$l2j$1@smc.vnet.net>

Dan,
Since your p[t] is a step function, your problem can be solved
exactly.
   G[t] = Integrate[p[x]*F[t-x], {x,-Infinity,Infinity}]
Given your choice of p this becomes
   G[t] = Integrate[F[t-x], {x,-Infinity,a}]
where a=30. Substitute z=t-x, you get
   G[t] = Integrate[F[z], {z,t-a,Infinity}]
The first derivative of this is
   G'[t] = -F[t-a]
So you can find F simply from
   F[t] = -G'[t+a] = -G'[t+30]

Roman.

On Apr 7, 10:00 am, "dantimatter" <dantimat... at gmail.com> wrote:
> hello all,
>
> i have a function G(t) which is a convolution of F(t) and p(t):
> G=F*p .
> now p(t) is a step function (p(t)=UnitStep[30-t]) and i know what G(t)
> is, but i'd like to get at F.  what i've been trying to do is take the
> laplace transform of G and dividing by the laplace transform of p,
> then inverting that to get F. i'm having a lot of trouble.   the built
> in inverseLaplaceTransform function is unable to do it, and i haven't
> found any numerical routines that are reliable.  is there some reason
> why i shouldn't be able to do this?  perhaps there are better ways to
> get at F(t)?  your advice is much appreciated.
>
> thanks,
> dan




  • Prev by Date: Re: Default value
  • Next by Date: Re: numerical inversion of laplace transform
  • Previous by thread: Re: numerical inversion of laplace transform
  • Next by thread: Re: numerical inversion of laplace transform