Fourier Transform of a Sinc function

A few days ago, I was trying to do the convolution between a Sinc function and a Gaussian function. But I got stuck from the first step, when I tried to solve that by using the convolution theorem, namely the Fourier transform of the Sinc(x), although I knew it is very easy to find the right answer by Googling or Mathematica. But it worth a try to be done by hand. So, here we go:
We know that the Fourier transform of Sinc(z) is,

    \[\int_{-\infty}^{+\infty} {{\sin(z)}\over{z}} e^{-i \omega z} dz\]

and

    \[\int_{-\infty}^{+\infty} {{\sin(z)}\over{z}} e^{-i \omega z} dz = \int_{-\infty}^{+\infty} {{e^{iz}-e^{-iz}}\over{2iz}} e^{-i \omega z} dz.\]

So,

(1)   \begin{equation*} \int_{-\infty}^{+\infty} {{\sin(z)}\over{z}} e^{-i \omega z} dz={1 \over {2i}} \int_{-\infty}^{+\infty} {{e^{i(1-\omega)z}}\over{z}} dz + {1 \over {2i}} \int_{+\infty}^{-\infty} {{e^{-i(1+\omega)z}}\over{z}} dz. \end{equation*}

Let us consider the first item, when \omega < 1, namely (1-\omega)>0, we can choose the path below to do the contour integration.

Using the method of complex residues, we take the contour with no singular point, separating the path into four parts, namely A, B, C and D shown as the red letters in the figure. Thus,

(2)   \begin{equation*} \begin{split} \oint {{e^{i(1-\omega)z}}\over{z}} dz=\lim_{\substack{{\rho \to 0}\\{R\to\infty}}} \int_{-R}^{-\rho}{{e^{i(1-\omega)z}}\over{z}} dz+ \lim_{\rho \to 0} \int_{B}{{e^{i(1-\omega)z}}\over{z}} dz +\lim_{\substack{{\rho \to 0}\\{R\to\infty}}} \int^{R}_{\rho}{{e^{i(1-\omega)z}}\over{z}} dz\\ +\lim_{R \to \infty} \int_{D}{{e^{i(1-\omega)z}}\over{z}} dz = 2\pi i {\mathrm{Res}} f(z)=0. \end{split} \end{equation*}

Also if (1-\omega)>0, namely ~\omega < 1, and when ~{R \to \infty} and ~{\rho \to 0}, it is easy to proof using Jordan’s lemma, that the fourth term is equal to zero. Thus we can the rewrite the rest items as,

    \[\int_{-\infty}^{0_{-}}{{e^{i(1-\omega)z}}\over{z}} dz+\int^{+\infty}_{0_{+}}{{e^{i(1-\omega)z}}\over{z}} dz =\int_{-\infty}^{+\infty}{{e^{i(1-\omega)z}}\over{z}} dz =-\lim_{\rho \to 0} \int_{B}{{e^{i(1-\omega)z}}\over{z}} dz.\]

In order to calculate the value of the item on the right side, I use z=\rho e^{i\theta}, so dz=i \rho e^{i\theta}d\theta, then we have,

    \[\lim_{\rho \to 0} \int_{B}{{e^{i(1-\omega)z}}\over{z}} dz=\lim_{\rho \to 0}\int_{\pi}^0 \exp[i (1-\omega) \rho e^{i \theta}] i d\theta=\int_{\pi}^0id\theta=i(0-\pi)=-i \pi.\]

So, in the end, we get the value of the first integration item shown in the Eq.(1),

(3)   \begin{equation*} \int_{-\infty}^{+\infty}{{e^{i(1-\omega)z}}\over{z}} dz = i\pi,\quad \omega<1. \end{equation*}

When \omega>1, according to the Jordan’s lemma, we can not get the fourth item in the Eq.(2) to be zero. So we have to choose another path for finalizing the integration.
So I chose the path below to deal with such a situation,

It is not hard to calculate the same way as the previous one,

(4)   \begin{equation*} \begin{split} \oint {{e^{i(1-\omega)z}}\over{z}} dz=\lim_{\substack{{\rho \to 0}\\{R\to\infty}}} \int_{R}^{\rho}{{e^{i(1-\omega)z}}\over{z}} dz+ \lim_{\rho \to 0} \int_{B}{{e^{i(1-\omega)z}}\over{z}} dz +\lim_{\substack{{\rho \to 0}\\{R\to\infty}}} \int^{-R}_{-\rho}{{e^{i(1-\omega)z}}\over{z}} dz\\ +\lim_{R \to \infty}\int_{D}{{e^{i(1-\omega)z}}\over{z}} dz = 2\pi i {\mathrm{Res}} f(z)=0, \end{split} \end{equation*}

and here again I use the same trick, z=R e^{i\theta}, so dz=i R e^{i\theta}d\theta. Because \omega>1,

    \[\lim_{R \to \infty} \int_{B}{{e^{i(1-\omega)z}}\over{z}} dz=\lim_{R \to \infty}\int_{\pi}^0 \exp[i (1-\omega) R e^{i \theta}] i d\theta=\int_{\pi}^0 0 \times id\theta=0,\]

    \[\lim_{\rho \to 0} \int_{B}{{e^{i(1-\omega)z}}\over{z}} dz=\lim_{\rho \to 0}\int^{\pi}_0 \exp[i (1-\omega) \rho e^{i \theta}] i d\theta=\int^{\pi}_0 i d\theta=i(\pi-0)=i \pi.\]

So, in the end,

(5)   \begin{equation*} \int_{-\infty}^{+\infty}{{e^{i(1-\omega)z}}\over{z}} dz = i\pi,\quad \omega>1. \end{equation*}

Now, let’s move on to the second part in the Eq.(1). Following the same steps, using contour above real axis I get,

(6)   \begin{equation*} \int_{+\infty}^{-\infty}{{e^{-i(1+\omega)z}}\over{z}} dz = \lim_{\rho \to 0}\int_{\pi}^0 \exp[-i (1+\omega) \rho e^{i \theta}] i d\theta = -i\pi,\quad \omega< -1; \end{equation*}

and using contour blow real axis I get,

(7)   \begin{equation*} \int_{+\infty}^{-\infty}{{e^{-i(1+\omega)z}}\over{z}} dz = \lim_{\rho \to 0}\int^{\pi}_0 \exp[-i(1+\omega) \rho e^{i \theta}] i d\theta = i\pi,\quad \omega>-1. \end{equation*}

Finally we can write the answer as, if -1< \omega<1

    \[\int_{-\infty}^{+\infty} {{\sin(z)}\over{z}} e^{-i \omega z} dz={1 \over {2i}} \int_{-\infty}^{+\infty} {{e^{i(1-\omega)z}}\over{z}} dz + {1 \over {2i}} \int_{+\infty}^{-\infty} {{e^{-i(1+\omega)z}}\over{z}} dz={{\pi}\over2}+{{\pi}\over2}=\pi;\]

and if \omega>1 and \omega<-1

    \[\int_{-\infty}^{+\infty} {{\sin(z)}\over{z}} e^{-i \omega z} dz={1 \over {2i}} \int_{-\infty}^{+\infty} {{e^{i(1-\omega)z}}\over{z}} dz + {1 \over {2i}} \int_{+\infty}^{-\infty} {{e^{-i(1+\omega)z}}\over{z}} dz={{\pi}\over2}-{{\pi}\over2}=0.\]

(Note: if we define the Fourier Transform as \int_{-\infty}^{+\infty} {{\sin(z)}\over{z}} e^{-2\pi i \omega z} dz, Then the result could be a little different.)

    Share it ! (点击分享)

2 Thoughts on “Fourier Transform of a Sinc function

  1. Marvin Chrome 19.0.1049.3 Windows 7 on 2012/05/15 at 20:03 said:

    我发现我有很多的流量是从你这里来的~除了豆瓣就是这里了~

发表评论

电子邮件地址不会被公开。 必填项已用*标注

您可以使用这些HTML标签和属性: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>

WP-SpamFree by Pole Position Marketing

Post Navigation