How to Integrate the Sinc Function

Reviewed by Joseph Meyer

Last Updated: August 7, 2020

The cardinal sine function, also known as the sinc function, is the function

$\operatorname {sinc}x={\begin{cases}{\dfrac {\sin x}{x}}&{\text{if }}x\neq 0,\\\ \ \ 1&{\text{if }}x=0.\end{cases}}$

This function frequently pops up first as an example of evaluation of limits, and it is well-known that $\lim _{{x\to 0}}{\frac {\sin x}{x}}=1;$ hence, why the function at 0 is defined to be that limiting value. However, this function primarily finds wider applicability in signal analysis and related fields. For example, the Fourier transform of a rectangular pulse is the sinc function.

Evaluating the integral of this function is rather difficult because the antiderivative of the sinc function cannot be expressed in terms of elementary functions. This means that we cannot directly apply the fundamental theorem of calculus. We will instead employ Richard Feynman's trick of differentiating under the integral. We will also show a more general solution using residue theory.

Method 1

Method 1 of 2:

Differentiation Under the Integral

Download Article

{"smallUrl":"https:\/\/www.wikihow.com\/images\/thumb\/8\/8b\/1000px-Si_sinc.png\/460px-1000px-Si_sinc.png","bigUrl":"\/images\/thumb\/8\/8b\/1000px-Si_sinc.png\/703px-1000px-Si_sinc.png","smallWidth":460,"smallHeight":366,"bigWidth":704,"bigHeight":560,"licensing":"<div class=\"mw-parser-output\"><p>License: <a target=\"_blank\" rel=\"nofollow noreferrer noopener\" class=\"external text\" href=\"https:\/\/creativecommons.org\/licenses\/by-sa\/3.0\/\">Creative Commons<\/a><br>Details: Georg-Johann<br>\n<\/p><\/div>"}

1
Begin with the integral to be evaluated. We are evaluating over the entire real line, so the limits will be positive and negative infinity. Above is a visualization of the function with both definitions - unnormalized (in red) and normalized (in blue). We will be evaluating the unnormalized sinc function.
- $\int _{{-\infty }}^{{\infty }}{\frac {\sin x}{x}}{\mathrm {d}}x$
- We see from the graph that ${\frac {\sin x}{x}}$ is an even function, which can be confirmed by looking at the function above. Then, we can factor out a 2.
  - $2\int _{{0}}^{{\infty }}{\frac {\sin x}{x}}{\mathrm {d}}x$
- The integral above with bounds of 0 to infinity is also known as the Dirichlet integral.
2
Define a function $G(t)$ . The purpose of defining such a function with an argument $t$ is so that we can work with an integral that is easier to evaluate, whilst meeting the conditions of the sinc integral for appropriate values of $t.$ In other words, putting the $e^{{-tx}}$ term inside the integral is valid, since the integral converges for all $t\geq 0,$ while setting $t=0$ recovers the original integral. This reformulation means that we are ultimately evaluating $G(0).$
- $G(t)=\int _{{0}}^{{\infty }}{\frac {\sin x}{x}}e^{{-tx}}{\mathrm {d}}x$
Advertisement
3
Differentiate under the integral. We can move the derivative under the integration sign because the integral is being taken with respect to a different variable. While we do not justify this operation here, it is widely applicable for a great many functions. Keep in mind that $t$ is to be treated as a variable throughout the evaluation, not a constant.
- ${\begin{aligned}{\frac {{\mathrm {d}}G}{{\mathrm {d}}t}}&={\frac {{\mathrm {d}}}{{\mathrm {d}}t}}\int _{{0}}^{{\infty }}{\frac {\sin x}{x}}e^{{-tx}}{\mathrm {d}}x\\&=\int _{{0}}^{{\infty }}{\frac {\partial }{\partial t}}e^{{-tx}}{\frac {\sin x}{x}}{\mathrm {d}}x\\&=-\int _{{0}}^{{\infty }}e^{{-tx}}\sin x{\mathrm {d}}x\end{aligned}}$
4
Evaluate ${\frac {{\mathrm {d}}G}{{\mathrm {d}}t}}$ . This is, in fact, the evaluation for the Laplace transform of $\sin x.$ The most basic way to evaluate this integral is by using integration by parts, which we work out below. See the tips for a more powerful way to integrate this. Pay attention to the signs.
- ${\begin{aligned}{\frac {{\mathrm {d}}G}{{\mathrm {d}}t}}&=e^{{-tx}}\cos x{\Big |}_{{0}}^{{\infty }}+\int _{{0}}^{{\infty }}te^{{-tx}}\cos x{\mathrm {d}}x\\&=e^{{-tx}}\cos x{\Big |}_{{0}}^{{\infty }}+\left[te^{{-tx}}\sin x{\Big |}_{{0}}^{{\infty }}+\int _{{0}}^{{\infty }}t^{{2}}e^{{-tx}}\sin x{\mathrm {d}}x\right]\\&=e^{{-tx}}\cos x{\Big |}_{{0}}^{{\infty }}+te^{{-tx}}\sin x{\Big |}_{{0}}^{{\infty }}-t^{{2}}{\frac {{\mathrm {d}}G}{{\mathrm {d}}t}}\end{aligned}}$
- ${\frac {{\mathrm {d}}G}{{\mathrm {d}}t}}(1+t^{{2}})=(0-1)+(0-0)$
- ${\frac {{\mathrm {d}}G}{{\mathrm {d}}t}}=-{\frac {1}{1+t^{{2}}}}$
5
Integrate both sides with respect to $t$ . This recovers $G(t)$ under a different variable. Since the integrand is the differential of a well-known function, this evaluation is trivial.
- $-\int {\frac {1}{1+t^{{2}}}}{\mathrm {d}}t=-\tan ^{{-1}}t+C$
- Here, we recognize that $G(t)=0$ as $t\to \infty$ for both this integral and the one defined in step 2. However, $\lim _{{t\to \infty }}\tan ^{{-1}}t={\frac {\pi }{2}},$ so $C={\frac {\pi }{2}}$ as well.
- Therefore, $G(t)=-\tan ^{{-1}}t+{\frac {\pi }{2}}.$
6
Evaluate the sinc integral. Now that we have $G(t),$ where $t\geq 0,$ we can substitute 0 for $t$ and find that $G(0)={\frac {\pi }{2}}.$
- $G(0)=\int _{{0}}^{{\infty }}{\frac {\sin x}{x}}{\mathrm {d}}x={\frac {\pi }{2}}$
- Finally, we recall that to integrate over all the reals, we simply multiply by 2, as $\operatorname {sinc}x$ is an even function.
  - $\int _{{-\infty }}^{{\infty }}{\frac {\sin x}{x}}{\mathrm {d}}x=\pi$
- It is worth memorizing this answer, as it can pop up in multiple contexts.

Method 2

Method 2 of 2:

Residue Theory

Download Article

1
Consider the integral below. Recall that $\sin x$ is simply the imaginary part of the exponential function $e^{{ix}}.$ This integral is continuous except for the singularity at $x=0.$
- $\int _{{-\infty }}^{{\infty }}{\frac {e^{{ix}}}{x}}{\mathrm {d}}x$
2
Consider the contour integral with an indented contour. The easiest improper integrals evaluated using residue theory use a semicircular arc that traces the real line from some boundary $-a$ to $a$ and arcs counterclockwise back to $-a,$ while $a\to \infty .$ However, we cannot use this because of the pole at the origin. The solution is to use an indented contour that goes around the pole.
- $\int _{{\Gamma }}{\frac {e^{{iz}}}{z}}{\mathrm {d}}z$
- The contour $\Gamma$ is split into four parts. We begin from $-a$ and traverse the real line to some small number $-\epsilon .$ Then a semicircular arc $\gamma _{{\epsilon }}$ with radius $\epsilon$ goes clockwise to $\epsilon$ on the real axis. This contour then goes to $a,$ from which a semicircular arc $\gamma _{{a}}$ with radius $a$ goes counterclockwise and back to $-a.$ The important thing to note here is that this integral does not have any singularities within the contour, and is therefore 0. We can therefore write the following.
- $\int _{{-a}}^{{-\epsilon }}{\frac {e^{{ix}}}{x}}{\mathrm {d}}x+\int _{{\epsilon }}^{{a}}{\frac {e^{{ix}}}{x}}{\mathrm {d}}x=-\int _{{\gamma _{{\epsilon }}}}{\frac {e^{{iz}}}{z}}{\mathrm {d}}z-\int _{{\gamma _{{a}}}}{\frac {e^{{iz}}}{z}}{\mathrm {d}}z$
3
Use Jordan's lemma to evaluate the $\gamma _{{a}}$ integral. Typically, for this integral to vanish, the degree of the denominator must be at least two greater than the degree of the numerator. Jordan's lemma implies that if such a rational function is multiplied by an $e^{{iz}}$ term, then the degree of the denominator need only be at least one greater. Therefore, this integral vanishes.
- $\int _{{\gamma _{{a}}}}{\frac {e^{{iz}}}{z}}{\mathrm {d}}z=0$
4
Evaluate the $\gamma _{{\epsilon }}$ integral.
- If you are familiar with contour integrals of ${\frac {1}{z}}$ involving circular arc contours, the example involves the fact that the integral depends on the angle that the arc traverses. In our example, the arc is being integrated from the angle $\pi$ to $0$ in a clockwise fashion. Such an integral will therefore equal $-\pi i.$
- We can generalize this result to arcs of any angle, but more importantly, for residues. See the tips for the theorem that this step uses. The residue at the origin is easily found to be $\operatorname {Res}(0)=1.$
  - $\lim _{{\epsilon \to 0^{{+}}}}\int _{{\gamma _{{\epsilon }}}}{\frac {e^{{iz}}}{z}}{\mathrm {d}}z=-\pi i\operatorname {Res}(0)=-\pi i$
5
Arrive at the answer to our integral. Because $\epsilon \to 0$ and $a\to \infty ,$ negate our result (see step 2) to arrive at our answer.
- $\int _{{-\infty }}^{{\infty }}{\frac {e^{{ix}}}{x}}{\mathrm {d}}x=\int _{{-a}}^{{-\epsilon }}{\frac {e^{{ix}}}{x}}{\mathrm {d}}x+\int _{{\epsilon }}^{{a}}{\frac {e^{{ix}}}{x}}{\mathrm {d}}x=\pi i$
6
Consider the imaginary part of the integral above. The above result really gives us two real results. First of all, the integral of the sinc function immediately follows.
- $\operatorname {Im}\int _{{-\infty }}^{{\infty }}{\frac {e^{{ix}}}{x}}{\mathrm {d}}x=\int _{{-\infty }}^{{\infty }}{\frac {\sin x}{x}}{\mathrm {d}}x=\pi$
- Second of all, the principal-valued integral of a related function ${\frac {\cos x}{x}}$ follows as well if we take the real part of our result, which is 0.

Community Q&A

Search

Add New Question

Question

Hi guys. The last sentence in Method 2.6 has to be clarified, I believe. cos(x)/x diverges for x->0 in a non-integrable fashion. To say something vanishs assumes it is well-defined

Community Answer

Yes, you are correct. The integral vanishes only in a principal-valued sense.

Thanks! We're glad this was helpful.
Thank you for your feedback.
If wikiHow has helped you, please consider a small contribution to support us in helping more readers like you. We’re committed to providing the world with free how-to resources, and even $1 helps us in our mission. Support wikiHow
Yes No

Not Helpful 0 Helpful 0

Ask a Question

200 characters left

Include your email address to get a message when this question is answered.

Submit

Tips

Some conventions define $\operatorname {sinc}x$ for $x\neq 0$ to be $\operatorname {sinc}x={\frac {\sin \pi x}{\pi x}}$ instead. This definition has a pretty normalization such that $\int _{{-\infty }}^{{\infty }}{\frac {\sin \pi x}{\pi x}}{\mathrm {d}}x=1.$ Verify this answer by redoing the steps above with this normalized function.

Thanks
Helpful 0 Not Helpful 0
In step 4 of method 1, we integrated via integration by parts. However, this process can be a bit slow. Another way is to use complex numbers by "complexifying" the integral.
- Recall Euler's formula relating exponents to rotation $e^{{ix}}=\cos x+i\sin x.$ We can rewrite the integral as the imaginary part of the complex exponential $\operatorname {Im},$ where for any complex number $z=a+bi,$ $\operatorname {Im}(z)=b.$
  - ${\frac {{\mathrm {d}}G}{{\mathrm {d}}t}}=-\operatorname {Im}\int _{{0}}^{{\infty }}e^{{(i-t)x}}{\mathrm {d}}x$
- We now have an exponential integral instead of a product of two functions. Now we simply do the integral and pick out the imaginary part.
  - ${\begin{aligned}-\operatorname {Im}\int _{{0}}^{{\infty }}e^{{(i-t)x}}{\mathrm {d}}x&=-\operatorname {Im}{\frac {e^{{(i-t)x}}}{i-t}}{\Bigg |}_{{0}}^{{\infty }}\\&=-\operatorname {Im}\left({\frac {1}{t-i}}\right)\\&=-{\frac {1}{t^{{2}}+1}}\end{aligned}}$
- This is a much more applicable and powerful method of integration. An advantage of this method is that we obtain the real part for free.
Thanks
Helpful 0 Not Helpful 0
In step 4 of method 2, we use a theorem that relates the contour integral around simple poles to the angle of the arc and the residue of the pole. Specifically, the theorem states this: If $f(z)$ has a simple pole at $z=z_{{0}}$ and $\gamma _{{\epsilon }}$ is a circular arc defined by $\gamma _{{\epsilon }}:z(t)=z_{{0}}+\epsilon e^{{it}},\ t_{{1}}\leq t\leq t_{{2}},$ then the following integral holds.
- $\lim _{{\epsilon \to 0^{{+}}}}\int _{{\gamma _{{\epsilon }}}}f(z){\mathrm {d}}z=i(t_{{2}}-t_{{1}})\operatorname {Res}(f(z);z_{{0}})$
- In our example, our contour was clockwise and integrated over $\pi ,$ so the result would be $-\pi i\operatorname {Res}(0),$ as used.
Thanks
Helpful 0 Not Helpful 0