Hi,

could someone help me find a closed form for the inverse fourier transform of

Namely

sech is the hyperbolic secant :

and m>0

Thanks for you help

AlexisM

Printable View

- Jul 19th 2010, 04:15 AMAlexisMInverse fourier transform of sech(mt/2)^(1/m)
Hi,

could someone help me find a closed form for the inverse fourier transform of

Namely

sech is the hyperbolic secant :

and m>0

Thanks for you help

AlexisM - Jul 19th 2010, 06:46 AMAckbeet
Could you please be more explicit about what is in the argument of the sech function? Do you mean

or

You have to be very careful with parentheses! - Jul 19th 2010, 08:43 AMAlexisM
- Jul 19th 2010, 09:05 AMAckbeet
WolframAlpha balks at that integral (as does my local copy of Mathematica, as you'd expect). Are you sure the inverse FT has a closed-form solution? It sure doesn't look like it.

Also, I should point out that the two integrals you claim are equal are only equal if the first one is real. If you take the real part of the first integral, you get the second one.

[EDIT]: I think you also inadvertently changed the limits of integration in the second integral without multiplying by two to take advantage of the symmetry. - Jul 19th 2010, 09:22 AMAlexisM
yes i forgot the factor 2 (i'll edit the previous post to correct that). This, I think possible, because the function sech is even. So if i split the integral on both side of zero, replaces x with -x, dx with -dx and invert the sens of integration I can factorize and get a factor exp(itx)+exp(-itx) = 2*cos(tx).

I'm not sure it has a closed form (may be it's not the good expression but I am french so please forgive me if I misuse some expressions) , but i'm about sure you can give an answer in the form of a serie or something like that. - Jul 19th 2010, 09:38 AMAckbeet
[EDIT]: You can try a term-by-term series expansion integration. However, perhaps a more fundamental question is this: does the inverse FT integral you've given us even converge?

- Jul 19th 2010, 11:24 AMAlexisM
it converges (it is computable with closed form) for m->0, m=1, m->infty.

Here is a demonstration for the FT of using residues.

I am not familiar with residues but I wonder if the same method can't be used to compute the inverse FT of using the similarity between

and

- Jul 19th 2010, 01:58 PMAckbeet
[EDIT]: This post is mostly bosh.

Computable with closed form is not the same thing as convergent. In this case, by convergent, I mean that the integral is finite and "settles down" to one specific value.

Let's look at the asymptotics here. Suppose we look at

in the integral

Then

This implies that

If we take half the integral, namely,

we see that this integral does not converge.

These are not rigorous arguments. However, the problem I see with your inverse FT is that you've got oscillatory behavior (which is good for convergence, usually) multiplied by an exponential function with a positive exponent (not good). So it's not merely that you've got a very large function: you've got huge alternating areas "canceling" each other out. Basically, it's a divergent integral with no defined value, so far as I can see.

Are you sure that you can even take the inverse FT of this function? - Jul 19th 2010, 02:52 PMAlexisM
I think we have a confusion about the définition of

This function "behaves" well, since it decreses exponentially.

And when , then

and

So we have an oscillatory function damped with an exponentially decreasing function. I can't see why it would not converge. - Jul 19th 2010, 04:23 PMAckbeet
Oops. You're right. I was thinking cosh. Sech is a very nicely behaved function, I agree. And I agree that the integral should converge. However, it does not seem to have a nice antiderivative. I would probably go with the series solution.

- Jul 20th 2010, 02:01 AMAlexisM
Someone pointed to me this result from "A course on modern analysis, of Whittaker and Watson" :

for where is the beta function.

If I do the subsitutions :

this gives :

for

Thus :

This holds for a limited range of values of . Is there any reason that this integral could not exists for other values of ? Or does it mean that we can have different forms of the solution for different domains ? - Jul 20th 2010, 06:34 AMAckbeet
Wow. That's some pretty heavy analysis you have there. Good work. On the face of it, it seems like that approach could work just fine. If you look here, you'll see that the real parts of both arguments to the beta function must be positive. You have that here, since in the OP you mentioned that m>0. There doesn't seem to me any a priori reason why you couldn't have any values for .

- Jul 22nd 2010, 03:28 AMAlexisM
Yes, I think that the is just there to ensure that the real part of arguments is positive. So here I don't have to care about it. Would you know where I can find the demonstration for the formula, or some hints to do it ?

- Jul 22nd 2010, 03:36 AMAckbeet
Other than the Whittaker and Watson book you mentioned, I wouldn't know for certain where you could find that result and its proof. Some options you could try, though, would be books on mathematical methods for physicists. Those kinds of books always have significant sections on special functions. You could try Boas, or Arfken and Weber, or Courant and Hilbert.

Good luck! - Jul 22nd 2010, 09:04 AMAlexisM
Ok, I've finally found out myself :

I use mainly these two equivalences :

The second equality can be symmetrized using :

With a notation change (t->x) and applying this formula to our problem we get :

Now if i make the variable change :

we have

And we use the fact that the function we integrate is symmetric in u :

And the condition is slightly different (less strict) from the one in Whittaker & Watson, namely which gives no constraint on the imaginary part.