Fresnel integral

S(x) and C(x) The maximum of C(x) is about 0.977451424. If πt²/2 were used instead of t², then the image would be scaled vertically and horizontally (see below).

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

S(x)=\int_0^x \sin(t^2)\,\mathrm{d}t,\quad C(x)=\int_0^x \cos(t^2)\,\mathrm{d}t.

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.[1]

Definition

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

Normalised Fresnel integrals, S(x) and C(x). In these curves, the argument of the trigonometric function is πt2/2, as opposed to just t2 as above.
S(x)=\int_0^x \sin(t^2)\,\mathrm{d}t=\sum_{n=0}^{\infin}(-1)^n\frac{x^{4n+3}}{(2n+1)!(4n+3)}
C(x)=\int_0^x \cos(t^2)\,\mathrm{d}t=\sum_{n=0}^{\infin}(-1)^n\frac{x^{4n+1}}{(2n)!(4n+1)}

Some authors, including Abramowitz and Stegun, (eqs 7.3.1 7.3.2) use \frac{\pi}{2}t^2 for the argument of the integrals defining S(x) and C(x). To get these functions, multiply the above integrals by \sqrt{\frac{2}{\pi}} and multiply the argument x by \sqrt{\frac{\pi}{2}}.

Euler spiral

Main article: 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:

 \mathrm{d}x = C'(t)\,\mathrm{d}t = \cos(t^2)\,\mathrm{d}t \,
 \mathrm{d}y = S'(t)\,\mathrm{d}t = \sin(t^2)\,\mathrm{d}t \,

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

L = \int_0^{t_0} {\sqrt {\mathrm{d}x^2 + \mathrm{d}y^2}} = \int_0^{t_0}{\mathrm{d}t} = 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²), sin(t²)] also expresses the unit tangent vector along the spiral, giving θ = . Since t is the curve length, the curvature, \kappa can be expressed as:

 \kappa = \tfrac {1}{R} = \tfrac {\mathrm{d}\theta}{\mathrm{d}t} = 2t

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

\tfrac{\mathrm{d}\kappa}{\mathrm{d}t}=\tfrac {\mathrm{d}^2\theta}{\mathrm{d}t^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

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).
Complex Fresnel integral S(z)
Complex Fresnel integral C(z)
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].
or S(z) + i C(z) = \sqrt{\frac{\pi}{2}}\frac{1+i}{2} \operatorname{erf}\left(\frac{1+i}{\sqrt{2}}z\right)
\int_{0}^{\infty} \cos t^2\,\mathrm{d}t = \int_{0}^{\infty} \sin t^2\,\mathrm{d}t = \frac{\sqrt{2\pi}}{4} = \sqrt{\frac{\pi}{8}}.

Evaluation

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

The limits of C and S 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 x-axis, the bisector of the first quadrant y = x with x ≥ 0, and a circular arc of radius R centered at the origin.

As R goes to infinity, the integral along the circular arc tends to 0, the integral along the real axis tends to the half Gaussian integral

\int_0^\infty e^{-t^2} \mathrm{d}t = \frac{\sqrt{\pi}}{2},

and after routine transformations, the integral along the bisector of the first quadrant can be related to the limit of the Fresnel integrals.

Generalization

The integral \int x^m \exp(ix^n)\,\mathrm{d}x = \int\sum_{l=0}^\infty\frac{i^lx^{m+nl}}{l!}\,\mathrm{d}x
 = \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 [3]

\int x^m \exp(ix^n)\,\mathrm{d}x =\frac{x^{m+1}}{m+1}\,_1F_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),

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

\int x^m\sin(x^n)\,\mathrm{d}x = \frac{x^{m+n+1}}{m+n+1}
\,_1F_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

_1F_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)\,\mathrm{d}x=\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)\,\mathrm{d}x = \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)\,\mathrm{d}x = V_{n,m}(x)e^{ix^n}

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

Applications

The Fresnel integrals were originally used in the calculation of the field intensity in an environment related to the bending of light around opaque objects.[4] More recently, they have been used in the design of highways and railways, specifically their curvature transition zones[1] and roller coasters.[4] Another application is for calculating the transitions on a velodrome track to allow rapid entry to the bends and gradual exit.

See also

References

  1. 1.0 1.1 Stewart, James (2007). Essential Calculus. Belmont, Calif.: Thomson Brooks/Cole. p. 230. ISBN 0-495-01442-7.
  2. functions.wolfram.com, Fresnel integral S: Representations through equivalent functions and Fresnel integral C: Representations through equivalent functions. Note: Wolfram uses the Abramowitz & Stegun convention, which differs from the one in this article by factors of \sqrt{\pi/2}
  3. Mathar, R. J. (2012). "Series Expansion of Generalized Fresnel Integrals". arXiv:1211.3963.
  4. 4.0 4.1 Beatty, Thomas. "How to evaluate Fresnel Integrals" (PDF). FGCU MATH - SUMMER 2013. Retrieved 27 July 2013.

External links