Inverse Fourier Transform in 2-d

Click For Summary
SUMMARY

The discussion revolves around solving the inverse Fourier transform of a specific integral involving constants \(d\), \(k\), and \(\ell\). The integral is transformed into polar coordinates, leading to a more manageable form that incorporates Bessel functions. The final expression for the integral is derived as \(I = -i\frac{2\pi}{3} h k K_o(h\sqrt{k^2+\ell^2})\), where \(K_o\) is the modified Bessel function of the second kind. This transformation simplifies the computation significantly compared to the original form.

PREREQUISITES
  • Understanding of inverse Fourier transforms
  • Familiarity with polar coordinates in integration
  • Knowledge of Bessel functions, specifically \(J_0\) and \(K_0\)
  • Proficiency in complex analysis and Gaussian integrals
NEXT STEPS
  • Study the properties and applications of Bessel functions, particularly \(K_0\)
  • Learn techniques for transforming integrals into polar coordinates
  • Explore advanced topics in Fourier analysis, focusing on inverse transforms
  • Investigate Gaussian integrals and their applications in physics and engineering
USEFUL FOR

Mathematicians, physicists, and engineers working on Fourier analysis, particularly those dealing with complex integrals and Bessel functions in applied contexts.

nickthequick
Messages
39
Reaction score
0
Hi all,

I've been trying to solve the following

I = \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \frac{x}{(x^2+y^2+d^2)^{\frac{5}{2}}} e^{-i(kx+\ell y)} \ dx \ dy

where d,k,\ell are constants. I haven't been able to put this into a tractable analytic form and I figured I'd consult all of you experts for advice before I resorted to approximation methods. So does anyone see any obvious way of solving this?

Thanks!

Nick
 
Physics news on Phys.org
EDIT:
The angular integral is wrong!Go to polar coordinates. The integral over the angular coordinate is:
<br /> \int_{0}^{2 \pi}{\cos \phi \, e^{-i q \, \rho \, \cos (\phi - \phi_0)} \, d\phi}<br />
where q = \sqrt{k^2 + l^2}, and \tan \phi_0 = l/k. Due to the periodicity of the two functions, we can perform the translation \phi \rightarrow \phi + \phi_0, but keep the same range of integration. Then, we get:
<br /> \int_{0}^{2 \pi}{\cos (\phi + \phi_0) \, e^{-i q \, \rho \, \cos \phi} \, d\phi}<br />
Applying the addition theorem for the cosine, we get:
<br /> \cos \phi_0 \, \int_{0}^{2 \pi}{\cos \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi} - \sin \phi_0 \, \int_{0}^{2 \pi}{\sin \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi}<br />
The two integrals are:
<br /> \int_{0}^{2 \pi}{\cos \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi} = -2 \pi \, i \, J_1(q \, \rho)<br />
and
<br /> \int_{0}^{2 \pi}{\sin \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi} = 0<br />
(the last one may be seen by shifting \phi \rightarrow \phi + \pi. Then, we have an integral of an odd function over a symmetric interval.)

Then, the integral over the radial coordinate is:
<br /> I = -2 \pi \, i \, \cos \phi_0 \, \int_{0}^{\infty}{\frac{\rho^2}{(\rho^2 + d^2)^{5/2}} \, J_1(q \, \rho) \, d\rho}<br />
 
Last edited:
Perhaps you can make use of the identity

\frac{1}{z^\nu} = \frac{1}{\Gamma(\nu)} \int_0^\infty dt~t^{\nu-1} e^{-zt},
which holds for z, \nu &gt; 0.

You would then have a triple integral

I_{k\ell} = \frac{1}{\Gamma(5/2)}\int_{-\infty}^\infty dx \int_{-\infty}^\infty dy \int_0^\infty dt~x t^{3/2} \exp\left(-\left[ tx^2 + ikx + ty^2 + i\ell y + d^2t\right]\right).

This is Gaussian in each of x and y, so it looks like you should be able to at least reduce the integral to just one integral over t. Whether or not that will be easier than the integral Dickfore derived or doable at all, I don't know.
 
Hm. Thanks for the suggestions. Here's what I've come up with.

Rewrite I as (I've changed d to h here for clarity of presentation)

I = \frac{i}{3}\partial_k \partial_h\iint \frac{1}{(x^2+y^2+h^2)^{\frac{3}{2}}}e^{-i(kx+\ell y)} \ dx \ dy

=\frac{i}{3}\partial_k \partial_h\int_0^{2\pi}\int_0^{\infty} \frac{1}{(r^2+h^2)^{\frac{3}{2}}} e^{-i\vec{r}\cdot \vec{k}}r\ dr \ d\theta

WLOG, orient \vec{k}=(k,\ell) with the x-axis so that \vec{r}\cdot \vec{k} = |\vec{r}||\vec{k}| \cos \theta. Therefore, we have

I=-i\frac{2\pi}{3} \partial_k\partial_h\int_0^{\infty} \frac{J_o(|\vec{k}| r)}{(r^2+h^2)^{\frac{3}{2}}} r \ dr

where J_o is the Bessel function of the first kind. The integral on the RHS of the above equation can be solved in closed form, so that we find

I = -i\frac{2\pi}{3}\partial_k\partial_h K_o(h|\vec{k}|)= -i\frac{2\pi}{3} h k\ K_o(h\sqrt{k^2+\ell^2})

where K_o is a modified bessel function of the second kind. This might be as good as I can hope for. Also, I might have made some algebra mistakes along the way.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
11K
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
9K