MathGroup Archive 2007

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

Search the Archive

Re: Hankel transform question

  • To: mathgroup at
  • Subject: [mg82496] Re: Hankel transform question
  • From: Jim Rockford <jim.rockford1 at>
  • Date: Mon, 22 Oct 2007 05:40:44 -0400 (EDT)
  • References: <ff9t5n$6cu$><ffcjm9$scg$>

On Oct 20, 6:01 am, d... at wrote:
> What I get in the development Mathematica kernel is simply an
> unevaluated Integrate.
> In[12]:= Integrate[s*BesselJ[v,w*s], {s,0,Infinity},
>   Assumptions->{v>0,w>0}]
> Out[12]= Integrate[s BesselJ[v, s w], {s, 0, Infinity},
>   Assumptions -> {v > 0, w > 0}]
> But the integral is in fact divergent. To see this, check the series
> expansion at infinity.
> In[13]:= InputForm[Series[s*BesselJ[v,w*s], {s,Infinity,2},
>   Assumptions->{v>0,w>0}]]
> Out[13]//InputForm=
> Cos[Pi/4 + (Pi*v)/2 - s*w]*
>   SeriesData[s, Infinity, {Sqrt[2/Pi]/Sqrt[w]}, -1, 3, 2] +
>   SeriesData[s, Infinity,
>     {(-1 + 4*v^2)/(4*Sqrt[2*Pi]*w^(3/2))}, 1, 3, 2]*
>     Sin[Pi/4 + (Pi*v)/2 - s*w]
> Notice that the "main" term is, up to scaling and phase change, of the
> form Cos[s]*Sqrt[s]. This diverges in classical integration.
> To get a regularized result, you can do as below.
> In[14]:= InputForm[Integrate[s*BesselJ[v,w*s], {s,0,Infinity},
>   Assumptions->{v>0,w>0}, GenerateConditions->False]]
> Out[14]//InputForm= v/w^2
> Offhand I do not see any reason to expect a delta function result.
> Certainly DiracDelta[s] would be incorrect, because s is the integral
> variable of (definite) integration, hence cannot appear in the result.

Thanks Daniel.  Let me add a few comments.  First, I was not expecting
an integral of f(s) to converge either, with f(s) of the form

f(s) = s^a BesselJ(v,w s)       (v is the Bessel order, w is a
transform variable)

if a > 1/2 , since we know that the Bessel function only decays like 1/
sqrt(s)  as s goes to infinity.

However, it was Mathematica itself (version 5.2) that suggested this
integral might make sense, given some sort of regularization (what is
this regularization, by the way?  Principal Value?)

In fact, from the mathworld page

one can see that the delta function   Delta[w]   can be written as the
Hankel Transform of unity  (I'm ignoring 2*pi conventions):

Delta[w] = Integrate[s BesselJ[0,w s],{s,0,Infinity}]

I'm certain that this is simply the analog of the usual representation
of the Delta function as a Fourier Transform (the Hankel version is
simply in polar coordinates).

In any case, it doesn't seem that Mathematica recognizes this variant
of the Delta function.


  • Prev by Date: Problem with Integral in mathematica 5.1
  • Next by Date: Re: Can Integrate[expr,{x,a,b}] give an incorrect result?
  • Previous by thread: Re: Hankel transform question
  • Next by thread: Controlling StepLeftButton and StepRightButton in Animator