Re: Hankel transform question
- To: mathgroup at smc.vnet.net
- Subject: [mg82496] Re: Hankel transform question
- From: Jim Rockford <jim.rockford1 at gmail.com>
- Date: Mon, 22 Oct 2007 05:40:44 -0400 (EDT)
- References: <ff9t5n$6cu$1@smc.vnet.net><ffcjm9$scg$1@smc.vnet.net>
On Oct 20, 6:01 am, d... at wolfram.com 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
http://mathworld.wolfram.com/HankelTransform.html
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.
J