Dirichlet integral

Peter Gustav Lejeune Dirichlet

In mathematics, there are several integrals known as the Dirichlet integral, after the German mathematician Peter Gustav Lejeune Dirichlet, one of which is the improper integral of the sinc function over the positive real number line.

$${\int }_0^{\mathrm{\infty }}\frac{\sin x}{x}dx=\frac{\pi }{2}.$$

This integral is not absolutely convergent, meaning $\left|\frac{\sin x}{x}\right|$ has an infinite Lebesgue or Riemann improper integral over the positive real line, so the sinc function is not Lebesgue integrable over the positive real line. The sinc function is, however, integrable in the sense of the improper Riemann integral or the generalized Riemann or Henstock–Kurzweil integral. This can be seen by using Dirichlet's test for improper integrals.

It is a good illustration of special techniques for evaluating definite integrals, particularly when it is not useful to directly apply the fundamental theorem of calculus due to the lack of an elementary antiderivative for the integrand, as the sine integral, an antiderivative of the sinc function, is not an elementary function. In this case, the improper definite integral can be determined in several ways: the Laplace transform, double integration, differentiating under the integral sign, contour integration, and the Dirichlet kernel. But since the integrand is an even function, the domain of integration can be extended to the negative real number line as well.

Evaluation

Laplace transform

Let $f(t)$ be a function defined whenever Template:Tmath. Then its Laplace transform is given by $$ℒ\{f(t)\}=F(s)={\int }_0^{\mathrm{\infty }}{e}^{-st}f(t)dt,$$ if the integral exists.

A property of the Laplace transform useful for evaluating improper integrals is $$ℒ\left[\frac{f(t)}{t}\right]={\int }_s^{\mathrm{\infty }}F(u)du,$$ provided $\underset{t\to 0}{lim}\frac{f(t)}{t}$ exists.

In what follows, one needs the result Template:Tmath, which is the Laplace transform of the function $\sin t$ Template:See below as well as a version of Abel's theorem (a consequence of the final value theorem for the Laplace transform).

Therefore, $$\begin{array}{cc}{\int }_0^{\mathrm{\infty }}\frac{\sin t}{t}dt& =\underset{s\to 0}{lim}{\int }_0^{\mathrm{\infty }}{e}^{-st}\frac{\sin t}{t}dt=\underset{s\to 0}{lim}ℒ\left[\frac{\sin t}{t}\right]\\ & =\underset{s\to 0}{lim}{\int }_s^{\mathrm{\infty }}\frac{du}{u^2+1}=\underset{s\to 0}{lim}\arctan u{|}_s^{\mathrm{\infty }}\\ & =\underset{s\to 0}{lim}[\frac{\pi }{2}-\arctan (s)]=\frac{\pi }{2}.\end{array}$$

Double integration

Evaluating the Dirichlet integral using the Laplace transform is equivalent to calculating the same double definite integral by changing the order of integration, namely, $$(I_1={\int }_0^{\mathrm{\infty }}{\int }_0^{\mathrm{\infty }}{e}^{-st}\sin tdtds)=(I_2={\int }_0^{\mathrm{\infty }}{\int }_0^{\mathrm{\infty }}{e}^{-st}\sin tdsdt),$$ $$(I_1={\int }_0^{\mathrm{\infty }}\frac{1}{s^2+1}ds=\frac{\pi }{2})=(I_2={\int }_0^{\mathrm{\infty }}\frac{\sin t}{t}dt),\text{ provided }s>0.$$ The change of order is justified by the fact that for all Template:Tmath, the integral is absolutely convergent.

Differentiation under the integral sign (Feynman's trick)

First rewrite the integral as a function of the additional variable Template:Tmath, namely, the Laplace transform of Template:Tmath. So let $$f(s)={\int }_0^{\mathrm{\infty }}{e}^{-st}\frac{\sin t}{t}dt.$$

In order to evaluate the Dirichlet integral, we need to determine Template:Tmath. The continuity of $f$ can be justified by applying the dominated convergence theorem after integration by parts. Differentiate with respect to $s>0$ and apply the Leibniz rule for differentiating under the integral sign to obtain $$\begin{array}{cc}\frac{df}{ds}& =\frac{d}{ds}{\int }_0^{\mathrm{\infty }}{e}^{-st}\frac{\sin t}{t}dt={\int }_0^{\mathrm{\infty }}\frac{\partial }{\partial s}{e}^{-st}\frac{\sin t}{t}dt\\ & =-{\int }_0^{\mathrm{\infty }}{e}^{-st}\sin tdt.\end{array}$$

Now, using Euler's formula Template:Tmath, one can express the sine function in terms of complex exponentials: $$\sin t=\frac{1}{2i}(e^{it}-e^{-it}).$$

Therefore, $$\begin{array}{cc}\frac{df}{ds}& =-{\int }_0^{\mathrm{\infty }}{e}^{-st}\sin tdt=-{\int }_0^{\mathrm{\infty }}{e}^{-st}\frac{e^{it}-e^{-it}}{2i}dt\\ & =-\frac{1}{2i}{\int }_0^{\mathrm{\infty }}[e^{-t(s-i)}-e^{-t(s+i)}]dt\\ & =-\frac{1}{2i}{[\frac{-1}{s-i}{e}^{-t(s-i)}-\frac{-1}{s+i}{e}^{-t(s+i)}]}_0^{\mathrm{\infty }}\\ & =-\frac{1}{2i}[0-(\frac{-1}{s-i}+\frac{1}{s+i})]=-\frac{1}{2i}(\frac{1}{s-i}-\frac{1}{s+i})\\ & =-\frac{1}{2i}\left(\frac{s+i-(s-i)}{s^2+1}\right)=-\frac{1}{s^2+1}.\end{array}$$

Integrating with respect to $s$ gives $$f(s)=\int \frac{-ds}{s^2+1}=A-\arctan s,$$

where $A$ is a constant of integration to be determined. Since Template:Tmath, Template:Tmath, using the principal value. This means that for $s>0$ $$f(s)=\frac{\pi }{2}-\arctan s.$$

Finally, by continuity at Template:Tmath, we have Template:Tmath, as before.

Complex contour integration

Consider $$f(z)=\frac{e^{iz}}{z}.$$

As a function of the complex variable Template:Tmath, it has a simple pole at the origin, which prevents the application of Jordan's lemma, whose other hypotheses are satisfied.

Define then a new function $$g(z)=\frac{e^{iz}}{z+i\epsilon }.$$

The pole has been moved to the negative imaginary axis, so $g(z)$ can be integrated along the semicircle $\gamma$ of radius $R$ centered at $z=0$ extending in the positive imaginary direction, and closed along the real axis. One then takes the limit Template:Tmath.

The complex integral is zero by the residue theorem, as there are no poles inside the integration path Template:Tmath: $$0={\int }_{\gamma }g(z)dz={\int }_{-R}^R\frac{e^{ix}}{x+i\epsilon }dx+{\int }_0^{\pi }\frac{e^{i(R{e}^{i\theta }+\theta )}}{R{e}^{i\theta }+i\epsilon }iRd\theta .$$

The second term vanishes as $R$ goes to infinity. As for the first integral, one can use one version of the Sokhotski–Plemelj theorem for integrals over the real line: for a complex-valued function Template:Mvar defined and continuously differentiable on the real line and real constants $a$ and $b$ with $a<0<b$ one finds $$\underset{\epsilon \to 0^{+}}{lim}{\int }_a^b\frac{f(x)}{x\pm i\epsilon }dx=\mp i\pi f(0)+𝒫{\int }_a^b\frac{f(x)}{x}dx,$$

where $𝒫$ denotes the Cauchy principal value. Back to the above original calculation, one can write $$0=𝒫\int \frac{e^{ix}}{x}dx-\pi i.$$

By taking the imaginary part on both sides and noting that the function $\frac{\sin x}{x}$ is even, we get $${\int }_{-\mathrm{\infty }}^{+\mathrm{\infty }}\frac{\sin (x)}{x}dx=2{\int }_0^{+\mathrm{\infty }}\frac{\sin (x)}{x}dx.$$

Finally, $$\underset{\epsilon \to 0}{lim}{\int }_{\epsilon }^{\mathrm{\infty }}\frac{\sin (x)}{x}dx={\int }_0^{\mathrm{\infty }}\frac{\sin (x)}{x}dx=\frac{\pi }{2}.$$

Alternatively, choose as the integration contour for $f$ the union of upper half-plane semicircles of radii $\epsilon$ and $R$ together with two segments of the real line that connect them. On one hand the contour integral is zero, independently of $\epsilon$ and Template:Tmath; on the other hand, as $\epsilon \to 0$ and $R\to \mathrm{\infty }$ the integral's imaginary part converges to $2I+\mathrm{\Im }(\ln 0-\ln (\pi i))=2I-\pi$ (here $\ln z$ is any branch of logarithm on upper half-plane), leading to Template:Tmath.

Dirichlet kernel

Consider the well-known formula for the Dirichlet kernel: $$D_n(x)=1+2\sum _{k=1}^n\cos (2kx)=\frac{\sin [(2n+1)x]}{\sin (x)}.$$

It immediately follows that: $${\int }_0^{\frac{\pi }{2}}{D}_n(x)dx=\frac{\pi }{2}.$$

Define $$f(x)=\{\begin{array}{cc}\frac{1}{x}-\frac{1}{\sin (x)}& x\ne 0\\ 0& x=0\end{array}$$

Clearly, $f$ is continuous when Template:Tmath; to see its continuity at Template:Tmath apply L'Hopital's Rule: $$\underset{x\to 0}{lim}\frac{\sin (x)-x}{x\sin (x)}=\underset{x\to 0}{lim}\frac{\cos (x)-1}{\sin (x)+x\cos (x)}=\underset{x\to 0}{lim}\frac{-\sin (x)}{2\cos (x)-x\sin (x)}=0.$$

Hence, $f$ fulfills the requirements of the Riemann-Lebesgue Lemma. This means: $$\underset{\lambda \to \mathrm{\infty }}{lim}{\int }_0^{\pi /2}f(x)\sin (\lambda x)dx=0⟹\underset{\lambda \to \mathrm{\infty }}{lim}{\int }_0^{\pi /2}\frac{\sin (\lambda x)}{x}dx=\underset{\lambda \to \mathrm{\infty }}{lim}{\int }_0^{\pi /2}\frac{\sin (\lambda x)}{\sin (x)}dx.$$

We would like to compute: $$\begin{array}{cc}{\int }_0^{\mathrm{\infty }}\frac{\sin (t)}{t}dt=& \underset{\lambda \to \mathrm{\infty }}{lim}{\int }_0^{\lambda \frac{\pi }{2}}\frac{\sin (t)}{t}dt\\ =& \underset{\lambda \to \mathrm{\infty }}{lim}{\int }_0^{\frac{\pi }{2}}\frac{\sin (\lambda x)}{x}dx\\ =& \underset{\lambda \to \mathrm{\infty }}{lim}{\int }_0^{\frac{\pi }{2}}\frac{\sin (\lambda x)}{\sin (x)}dx\\ =& \underset{n\to \mathrm{\infty }}{lim}{\int }_0^{\frac{\pi }{2}}\frac{\sin ((2n+1)x)}{\sin (x)}dx\\ =& \underset{n\to \mathrm{\infty }}{lim}{\int }_0^{\frac{\pi }{2}}{D}_n(x)dx=\frac{\pi }{2}\end{array}$$

However, we must justify switching the real limit in $\lambda$ to the integral limit in Template:Tmath, which will follow from showing that the limit does exist.

Using integration by parts, we have: $${\int }_a^b\frac{\sin (x)}{x}dx={\int }_a^b\frac{d(1-\cos (x))}{x}dx={\frac{1-\cos (x)}{x}|}_a^b+{\int }_a^b\frac{1-\cos (x)}{x^2}dx$$

Now, as $a\to 0$ and $b\to \mathrm{\infty }$ the term on the left converges with no problem. See the list of limits of trigonometric functions. We now show that ${\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\frac{1-\cos (x)}{x^2}dx$ is absolutely integrable, which implies that the limit exists.

First, we seek to bound the integral near the origin. Using the Taylor-series expansion of the cosine about zero, $$1-\cos (x)=1-\sum _{k\ge 0}\frac{(-1{)}^{(k+1)}{x}^{2k}}{(2k)!}=\sum _{k\ge 1}\frac{(-1{)}^{(k+1)}{x}^{2k}}{(2k)!}.$$

Therefore, $$\left|\frac{1-\cos (x)}{x^2}\right|=|-\sum _{k\ge 0}\frac{x^{2k}}{2(k+1)!}|\le \sum _{k\ge 0}\frac{|x{|}^k}{k!}=e^{|x|}.$$

Splitting the integral into pieces, we have $${\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\left|\frac{1-\cos (x)}{x^2}\right|dx\le {\int }_{-\mathrm{\infty }}^{-\epsilon }\frac{2}{x^2}dx+{\int }_{-\epsilon }^{\epsilon }{e}^{|x|}dx+{\int }_{\epsilon }^{\mathrm{\infty }}\frac{2}{x^2}dx\le K,$$

for some constant Template:Tmath. This shows that the integral is absolutely integrable, which implies the original integral exists, and switching from $\lambda$ to $n$ was in fact justified, and the proof is complete.

See also

• Dirichlet distribution
• Dirichlet principle
• Sinc function
• Fresnel integral

External links

• Template:MathWorld