[Date Index]
[Thread Index]
[Author Index]
Re: numerical inversion of laplace transform
*To*: mathgroup at smc.vnet.net
*Subject*: [mg75143] Re: numerical inversion of laplace transform
*From*: dantimatter <dantimatter at gmail.com>
*Date*: Wed, 18 Apr 2007 05:05:58 -0400 (EDT)
*References*: <evfk8c$5bk$1@smc.vnet.net><evksdj$ppp$1@smc.vnet.net>
Hi Roman et al.,
I don't seem to be getting the same answer that you indicated I should
get in the April 12th response. Is there something about Fourier
transforms I'm not quite 'getting'?? The Fourier transform of the step
function p[t]=UnitStep[33.6-t] is
FTp = -((I*Cos[33.6*=CF=89])/(Sqrt[2*Pi]*=CF=89)) + Sqrt[Pi/2]*DiracDelta=
[=CF=89] +
Sin[33.6*=CF=89]/(Sqrt[2*Pi]*=CF=89)
The function that I've fit to my data is g[t]:
g[t] = 6.06 - 4.17*Cos[(Pi*t)/84] + 1.19*Cos[(Pi*t)/42] -
2=2E95*Sin[(Pi*t)/84] + 0.71*Sin[(Pi*t)/42]
so that the Fourier transform of g is
FTg = (-439.34 - 310.621*I)*DiracDelta[Pi - 84*=CF=89] + (62.45 +
37.59*I)*DiracDelta[Pi - 42*=CF=89] +
15.20*DiracDelta[=CF=89] + (62.45 - 37.59*I)*DiracDelta[Pi + 42*=CF=89] -
(439.34 - 310.62*I)*DiracDelta[Pi + 84*=CF=89]
I want to invert FTg/FTp, but when I use the InverseFourierTransform
routine, I get
-0.29*Cos[(Pi*t)/84] + 0.24*Cos[(Pi*t)/42] - 0.38*Sin[(Pi*t)/84] -
0=2E10*Sin[(Pi*t)/42]
What is it that I'm doing wrong? I don't get an answer with a
constant term. Overall, your answer looks a lot better than mine. I
suspect the problem comes partly from me having a DiracDelta in my
expression for FTp, or - and this is more likely - I'm not using the
Fourier transform properly. Can you see what I'm doing wrong?
As always, your help is greatly appreciated.
Dan
On Apr 12, 3:59 am, "Roman" <rschm... at gmail.com> wrote:
> Dan,
>
> Given your G(t), what I wrote can be simplified drastically. First,
> transform to exponential writing using
> TrigToExp[G[t]]
> Then, you can manually do the transformation I wrote at 3) on the
> prefactors of the exponentials. For example, the term with Exp[-i*Pi*t/
> 84], which has w=Pi/84, acquires an extra factor of -i*w/(1-
> Exp[i*w*a]) = 0.0257383 - 0.0187*i. The constant term gets a prefactor
> of 1/a.
>
> Putting everything together, you get
>
> F[t] = 303/(50*a)
> - Cos[(Pi*t)/84]*((59*Pi)/3360 + (139*Pi*Cot[(a*Pi)/168])/5600)
> + Cos[(Pi*t)/42]*((71*Pi)/8400 + (17*Pi*Cot[(a*Pi)/84])/1200)
> + ((139*Pi)/5600 - (59*Pi*Cot[(a*Pi)/168])/3360)*Sin[(Pi*t)/84]
> - ((17*Pi)/1200 - (71*Pi*Cot[(a*Pi)/84])/8400)*Sin[(Pi*t)/42]
>
> With a = 33.6 as you specify, you get
>
> F[t] = 0.18035714285714283
> - 0.16249350579592095*Cos[(Pi*t)/84]
> + 0.04101478009014459*Cos[(Pi*t)/42]
> + 0.0020508865613427935*Sin[(Pi*t)/84]
> - 0.035877998487864125*Sin[(Pi*t)/42]
>
> Roman.
>
> On Apr 11, 8:12 am, "dantimatter" <dantimat... at gmail.com> wrote:
>
> > hi dimitris,
>
> > sorry about that posting. here we go again:
>
> > G(t)=6.06 - 4.17*Cos[(Pi*t)/84] + 1.19*Cos[(Pi*t)/42] -
> > 2.95*Sin[(Pi*t)/84] + 0.71*Sin[(Pi*t)/42]
> > p(t) = UnitStep[33.6 - t]
>
> > so that the Laplace-transformed function that I'd like to invert, G(s)/
> > p(s) is
>
> > G(s)/p(s) =
> > (3.83*^7*E^(33.6*s)*(0.0029 - 0.049*s + s^2)*(0.0054 + 0.031*s + s^2))/
> > ((-1. + E^(33.6*s))*(Pi^2 + 1764.*s^2)*(Pi^2 + 7056.*s^2))
>
> > thanks!
> > dan
Prev by Date:
**Re : Re: Interpolation**
Next by Date:
**Re: LegendreP Evaluation Mystery**
Previous by thread:
**Re: numerical inversion of laplace transform**
Next by thread:
**GenerateConditions setting**
| |