Fresnel integral

Plots of Template:Math and Template:Math. The maximum of Template:Math is about Template:Math. If the integrands of Template:Mvar and Template:Mvar were defined using Template:Math instead of Template:Math, then the image would be scaled vertically and horizontally (see below).

The Fresnel integrals S(x) and C(x) are two transcendental functions named after Augustin-Jean Fresnel that are used in optics and are closely related to the error function (erf). They arise in the description of near-field Fresnel diffraction phenomena and are defined through the following integral representations:

$S(x)=\int _{0}^{x}\sin(t^{2})\,dt,\quad C(x)=\int _{0}^{x}\cos(t^{2})\,dt.$

The simultaneous parametric plot of S(x) and C(x) is the Euler spiral (also known as the Cornu spiral or clothoid). Recently, they have been used in the design of highways and other engineering projects.

Definition

Fresnel integrals with arguments Template:Pit²/2 instead of t² converge to 0.5.

The Fresnel integrals admit the following power series expansions that converge for all x:

${\begin{aligned}S(x)&=\int _{0}^{x}\sin(t^{2})\,dt=\sum _{n=0}^{\infty }(-1)^{n}{\frac {x^{4n+3}}{(2n+1)!(4n+3)}}.\\C(x)&=\int _{0}^{x}\cos(t^{2})\,dt=\sum _{n=0}^{\infty }(-1)^{n}{\frac {x^{4n+1}}{(2n)!(4n+1)}}.\end{aligned}}$

Some authors, including Abramowitz and Stegun, (eqs 7.3.1 – 7.3.2) use $\pi t^{2}/2$ for the argument of the integrals defining S(x) and C(x). This changes their limits at infinity from $\textstyle (1/2)\cdot {\sqrt {\pi /2}}$ to $1/2$ and the arc length for the first spiral turn from $(2\pi )$ to 2 (at $t=2$), all smaller by a factor $\textstyle {\sqrt {2/\pi }}$. The alternative functions are also called Normalized Fresnel integrals.

Euler spiral

Euler spiral $(x,y)=(C(t),S(t))$. The spiral converges to the centre of the holes in the image as $t$ tends to positive or negative infinity.

The Euler spiral, also known as Cornu spiral or clothoid, is the curve generated by a parametric plot of $S(t)$ against $C(t)$. The Cornu spiral was created by Marie Alfred Cornu as a nomogram for diffraction computations in science and engineering.

From the definitions of Fresnel integrals, the infinitesimals $dx$ and $dy$ are thus:

${\begin{aligned}dx&=C'(t)\,dt=\cos(t^{2})\,dt,\\dy&=S'(t)\,dt=\sin(t^{2})\,dt.\end{aligned}}$

Thus the length of the spiral measured from the origin can be expressed as

$L=\int _{0}^{t_{0}}{\sqrt {dx^{2}+dy^{2}}}=\int _{0}^{t_{0}}dt=t_{0}.$

That is, the parameter $t$ is the curve length measured from the origin $(0,0)$, and the Euler spiral has infinite length. The vector $(\cos(t^{2}),\sin(t^{2}))$ also expresses the unit tangent vector along the spiral, giving $\theta =t^{2}$. Since t is the curve length, the curvature $\kappa$ can be expressed as

$\kappa ={\frac {1}{R}}={\frac {d\theta }{dt}}=2t.$

And the rate of change of curvature with respect to the curve length is

${\frac {d\kappa }{dt}}={\frac {d^{2}\theta }{dt^{2}}}=2.$

An Euler spiral has the property that its curvature at any point is proportional to the distance along the spiral, measured from the origin. This property makes it useful as a transition curve in highway and railway engineering: If a vehicle follows the spiral at unit speed, the parameter $t$ in the above derivatives also represents the time. That is, a vehicle following the spiral at constant speed will have a constant rate of angular acceleration.

Sections from Euler spirals are commonly incorporated into the shape of roller-coaster loops to make what are known as clothoid loops.

Properties

• C(x) and S(x) are odd functions of x.
• Asymptotics of the Fresnel integrals as $x\to \infty$ are given by the formulas:

${\begin{aligned}S(x)&={\sqrt {\frac {\pi }{2}}}\left({\frac {{\mbox{sign}}{(x)}}{2}}-\left[1+O(x^{-4})\right]\left({\frac {\cos {(x^{2})}}{x{\sqrt {2\pi }}}}+{\frac {\sin {(x^{2})}}{x^{3}{\sqrt {8\pi }}}}\right)\right),\\C(x)&={\sqrt {\frac {\pi }{2}}}\left({\frac {{\mbox{sign}}{(x)}}{2}}+\left[1+O(x^{-4})\right]\left({\frac {\sin {(x^{2})}}{x{\sqrt {2\pi }}}}-{\frac {\cos {(x^{2})}}{x^{3}{\sqrt {8\pi }}}}\right)\right).\end{aligned}}$

Complex Fresnel integral S(z)

• Using the power series expansions above, the Fresnel integrals can be extended to the domain of complex numbers, and they become analytic functions of a complex variable.
• The Fresnel integrals can be expressed using the error function as follows:

Complex Fresnel integral C(z)

${\begin{aligned}S(z)&={\sqrt {\frac {\pi }{2}}}{\frac {1+i}{4}}\left[\operatorname {erf} \left({\frac {1+i}{\sqrt {2}}}z\right)-i\operatorname {erf} \left({\frac {1-i}{\sqrt {2}}}z\right)\right],\\C(z)&={\sqrt {\frac {\pi }{2}}}{\frac {1-i}{4}}\left[\operatorname {erf} \left({\frac {1+i}{\sqrt {2}}}z\right)+i\operatorname {erf} \left({\frac {1-i}{\sqrt {2}}}z\right)\right].\end{aligned}}$

or

${\begin{aligned}C(z)+iS(z)&={\sqrt {\frac {\pi }{2}}}{\frac {1+i}{2}}\operatorname {erf} \left({\frac {1-i}{\sqrt {2}}}z\right),\\S(z)+iC(z)&={\sqrt {\frac {\pi }{2}}}{\frac {1+i}{2}}\operatorname {erf} \left({\frac {1+i}{\sqrt {2}}}z\right).\end{aligned}}$

• C and S are entire functions.

Limits as x approaches infinity

The integrals defining C(x) and S(x) cannot be evaluated in the closed form in terms of elementary functions, except in special cases. The limits of these functions as x goes to infinity are known:

$\int _{0}^{\infty }\cos t^{2}\,dt=\int _{0}^{\infty }\sin t^{2}\,dt={\frac {\sqrt {2\pi }}{4}}={\sqrt {\frac {\pi }{8}}}\approx 0.6267.$

The sector contour used to calculate the limits of the Fresnel integrals

The limits of Template:Math and Template:Math as the argument tends to infinity can be found by the methods of complex analysis. This uses the contour integral of the function

$e^{-t^{2}}$

around the boundary of the sector-shaped region in the complex plane formed by the positive Template:Math-axis, the bisector of the first quadrant Template:Math with Template:Math, and a circular arc of radius Template:Math centered at the origin.

As Template:Math goes to infinity, the integral along the circular arc $\gamma _{2}$ tends to Template:Math

$\left|\int _{\gamma _{2}}e^{-t^{2}}dt\right|\leq \int _{\gamma _{2}}|e^{-t^{2}}|dt=R\int _{0}^{\pi /4}e^{-R^{2}\cos 2t}dt\leq R\int _{0}^{\pi /4}e^{-R^{2}(1-{\frac {4}{\pi }}t)}dt={\frac {\pi }{4R}}\left(1-e^{-R^{2}}\right),$

where polar coordinates $z=Re^{it}$ were used and Jordan's inequality was utilised for the second inequality. The integral along the real axis $\gamma _{1}$ tends to the half Gaussian integral

$\int _{0}^{\infty }e^{-t^{2}}\,dt={\frac {\sqrt {\pi }}{2}}.$

Note too that because the integrand is an entire function on the complex plane, its integral along the whole contour is zero. Overall, we must have

$\int _{0}^{\infty }e^{-t^{2}}\,dt=\int _{\gamma _{3}}e^{-t^{2}}\,dt,$

where $\gamma _{3}$ denotes the bisector of the first quadrant, as in the diagram. To evaluate the right hand side, parametrize the bisector as

$t=re^{\pi i/4}={\frac {\sqrt {2}}{2}}(1+i)r$

where r ranges from 0 to $+\infty$. Note that the square of this expression is just $+ir^{2}$. Therefore, substitution gives the right hand side as

$\int _{0}^{\infty }e^{-ir^{2}}{\frac {\sqrt {2}}{2}}(1+i)\,dr.$

Using Euler's formula to take real and imaginary parts of $e^{-ir^{2}}$ gives this as

${\begin{aligned}&\int _{0}^{\infty }(\cos(r^{2})-i\sin(r^{2})){\frac {\sqrt {2}}{2}}(1+i)\,dr\\={}&{\frac {\sqrt {2}}{2}}\int _{0}^{\infty }\left[\cos(r^{2})+\sin(r^{2})+i\left(\cos(r^{2})-\sin(r^{2})\right)\right]\,dr={\frac {\sqrt {\pi }}{2}}+0i,\end{aligned}}$

where we have written $0i$ to emphasize that the original Gaussian integral's value is completely real with zero imaginary part. Letting $I_{C}=\int _{0}^{\infty }\cos(r^{2})\,dr,I_{S}=\int _{0}^{\infty }\sin(r^{2})\,dr$ and then equating real and imaginary parts produces the following system of two equations in the two unknowns $I_{C},I_{S}$:

${\begin{aligned}I_{C}+I_{S}&={\sqrt {\frac {\pi }{2}}},\\I_{C}-I_{S}&=0.\end{aligned}}$

Solving this for $I_{C}$ and $I_{S}$ gives the desired result.

Generalization

The integral

$\int x^{m}\exp(ix^{n})\,dx=\int \sum _{l=0}^{\infty }{\frac {i^{l}x^{m+nl}}{l!}}\,dx=\sum _{l=0}^{\infty }{\frac {i^{l}}{(m+nl+1)}}{\frac {x^{m+nl+1}}{l!}}$

is a confluent hypergeometric function and also an incomplete gamma function

${\begin{aligned}\int x^{m}\exp(ix^{n})\,dx&={\frac {x^{m+1}}{m+1}}\,_{1}F_{1}\left({\begin{array}{c}{\frac {m+1}{n}}\\1+{\frac {m+1}{n}}\end{array}}\mid ix^{n}\right)\\&={\frac {1}{n}}i^{(m+1)/n}\gamma \left({\frac {m+1}{n}},-ix^{n}\right),\end{aligned}}$

which reduces to Fresnel integrals if real or imaginary parts are taken:

$\int x^{m}\sin(x^{n})\,dx={\frac {x^{m+n+1}}{m+n+1}}\,_{1}F_{2}\left({\begin{array}{c}{\frac {1}{2}}+{\frac {m+1}{2n}}\\{\frac {3}{2}}+{\frac {m+1}{2n}},{\frac {3}{2}}\end{array}}\mid -{\frac {x^{2n}}{4}}\right)$.

The leading term in the asymptotic expansion is

$_{1}F_{1}\left({\begin{array}{c}{\frac {m+1}{n}}\\1+{\frac {m+1}{n}}\end{array}}\mid ix^{n}\right)\sim {\frac {m+1}{n}}\,\Gamma \left({\frac {m+1}{n}}\right)e^{i\pi (m+1)/(2n)}x^{-m-1},$

and therefore

$\int _{0}^{\infty }x^{m}\exp(ix^{n})\,dx={\frac {1}{n}}\,\Gamma \left({\frac {m+1}{n}}\right)e^{i\pi (m+1)/(2n)}.$

For m = 0, the imaginary part of this equation in particular is

$\int _{0}^{\infty }\sin(x^{a})\,dx=\Gamma \left(1+{\frac {1}{a}}\right)\sin \left({\frac {\pi }{2a}}\right),$

with the left-hand side converging for a > 1 and the right-hand side being its analytical extension to the whole plane less where lie the poles of $\Gamma (a^{-1})$.

The Kummer transformation of the confluent hypergeometric function is

$\int x^{m}\exp(ix^{n})\,dx=V_{n,m}(x)e^{ix^{n}},$

with

$V_{n,m}:={\frac {x^{m+1}}{m+1}}\,_{1}F_{1}\left({\begin{array}{c}1\\1+{\frac {m+1}{n}}\end{array}}\mid -ix^{n}\right).$

Numerical approximation

For computation to arbitrary precision, the power series is suitable for small argument. For large argument, asymptotic expansions converge faster. Continued fraction methods may also be used.

For computation to particular target precision, other approximations have been developed. Cody developed a set of efficient approximations based on rational functions that give relative errors down to Template:Val. A FORTRAN implementation of the Cody approximation that includes the values of the coefficients needed for implementation in other languages was published by van Snyder. Boersma developed an approximation with error less than Template:Val.

Applications

The Fresnel integrals were originally used in the calculation of the electromagnetic field intensity in an environment where light bends around opaque objects. More recently, they have been used in the design of highways and railways, specifically their curvature transition zones, see track transition curve. Other applications are roller coasters or calculating the transitions on a velodrome track to allow rapid entry to the bends and gradual exit.Template:Citation needed

See also

• Böhmer integral
• Fresnel zone
• Track transition curve
• Euler spiral
• Zone plate
• Dirichlet integral

References

• Template:Cite book
• Template:AS ref
• Template:Cite journal
• Template:Cite journal
• Template:Cite web (Uses πt²/2 instead of t².)
• Template:Dlmf
• Template:Cite journal
• Template:Cite arxiv

External links

• Cephes, free/open-source C++/C code to compute Fresnel integrals among other special functions. Used in SciPy and ALGLIB.
• Faddeeva Package, free/open-source C++/C code to compute complex error functions (from which the Fresnel integrals can be obtained), with wrappers for Matlab, Python, and other languages.
• Template:Springer
• Template:Cite web
• Template:Mathworld
• Template:Mathworld