Message Boards Message Boards


Reproducing Hankel transform pairs?

Posted 1 month ago
1 Reply
1 Total Likes

There is a classic Hankel transform pair from optics that I am trying to reproduce with Mathematica, without luck, so I can confidently generate an unknown case. The function in the spatial frequency domain, H, is:


with rho < 1. There are several versions of the normalized inverse Hankel transform of the frequency domain function available in the literature with hNorm[r] = h[r]/h[0] where h is the inverse hankel transform of H. I'd be happy to find any of them. For discussion sake one is published to be:


When I try the inversion, however, I only get a result for the rho*Sqrt[1-rho^2] term. The inverse transform for the ArcCos term (hStar) is returned unchanged:




I've tried to enforce various appropriate assumptions and have tried other Bessel function orders in the transform without luck. Any ideas?


InverseHankelTransform[]'s support is still somewhat limited,so don't be surprised if some things don't work yet.

I have only workaround: InverseHankelTransform[ArcCos[rho], rho, r] is:

    InverseHankelTransform[FourierCosTransform[ArcCos[a rho], a, s], 
      rho, r, Assumptions -> {s > 0}][[1]], s, a] /. a -> 1 // Expand

(*(\[Pi] BesselJ[0, r/2] BesselJ[1, r/2])/(2 r) - (
I Sqrt[\[Pi]] MeijerG[{{}, {1/2}}, {{0, 1}, {0}}, r^2/4])/r^2*)

but solution is expressed by MeijerG function.

  FullSimplify[HankelTransform[%, r, rho] // FunctionExpand, 
   Assumptions -> 0 < rho < 1]
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract