Gaussian integral

Gaussian integral

The Gaussian integral, or probability integral, is the improper integral of the Gaussian function e^-x}^2} over the entire real line. It is named after the German mathematician and physicist Carl Friedrich Gauss, and the equation is:

:int_{-infty}^infty e^{-x^2},dx = sqrt{pi}.

This integral has wide applications including normalization in probability theory and continuous Fourier transform. It also appears in the definition of the error function.

Although no elementary function exists for the error function, as can be proven by the Risch algorithm, the Gaussian integral can be solved analytically through the tools of calculus. That is, there is no elementary "indefinite integral" for int e^{-x^2},dx, but the "definite integral" int_{-infty}^infty e^{-x^2},dx can be evaluated.

Computation

A standard way to compute the Gaussian integral is
* consider the function e^{-(x^2+y^2)}=e^{-r^2} on the plane mathbf{R}^2, and compute its integral two ways:
* on the one hand, by double integration in the Cartesian coordinate system, its integral is a square: left(int e^{-x^2},dx ight)^2;
* on the other hand, by shell integration (a case of double integration in polar coordinates), its integral is computed to be pi.

Comparing these two computations yields the integral, though one should take care about the improper integrals involved.

Brief proof

Briefly, one computes that on the one hand,:egin{align}int_{mathbf{R}^2} e^{-(x^2+y^2)},dx,dy &= int_{-infty}^infty int_{-infty}^infty e^{-(x^2+y^2)},dx,dy\ &= left ( int_{-infty}^infty e^{-x^2},dx ight ) cdot left ( int_{-infty}^infty e^{-y^2},dy ight )\ &= left ( int_{-infty}^infty e^{-x^2},dx ight )^2.end{align}

On the other hand,:egin{align}int_{mathbf{R}^2} e^{-(x^2+y^2)},dx,dy&= int_0^{2pi} int_0^{infin} re^{-r^2},dr,d heta\&= 2pi int_0^infty re^{-r^2},dr\&= 2pi int_{-infty}^0 frac{1}{2} e^s,ds= pi int_{-infty}^0 e^s,ds= pi (e^0 - e^{-infty})= pi (1 - 0)= pi,end{align}where the factor of r comes from the transform to polar coordinates (r,dr,d heta is the standard measure on the plane, expressed in polar coordinates), and the substitution involves taking s=-r^2, so ds=-2r,dr.

Combining these yields left ( int_{-infty}^infty e^{-x^2},dx ight )^2=pi, so int_{-infty}^infty e^{-x^2},dx=sqrt{pi}.

Careful proof

To justify the improper double integrals and equating the two expressions, we begin with an approximating function:

:I(a)=int_{-a}^a e^{-x^2}dx.

so that the integral may be found by

:lim_{a oinfty} I(a) = int_{-infty}^{infty} e^{-x^2}, dx.

Taking the square of "I" yields

:I^2(a)= left ( int_{-a}^a e^{-x^2}, dx ight )cdot left ( int_{-a}^a e^{-y^2}, dy ight )= int_{-a}^a left ( int_{-a}^a e^{-y^2}, dy ight ),e^{-x^2}, dx = int_{-a}^a int_{-a}^a e^{-(x^2+y^2)},dx,dy.

Using Fubini's theorem, the above double integral can be seen as an area integral int e^{-(x^2+y^2)},d(x,y), taken over a square with vertices {(−"a", "a"), ("a", "a"), ("a", −"a"), (−"a", −"a")} on the "xy"-plane.

Since the exponential function is greater than 0 for all real numbers, it then follows that the integral taken over the square's incircle must be less than I(a)^2, and similarly the integral taken over the square's circumcircle must be greater than I(a)^2. The integrals over the two disks can easily be computed by switching from cartesian coordinates to
polar coordinates: x=r,cos heta, y= r,sin heta, d(x,y) = r, d(r, heta):

:int_0^{2pi}int_0^a re^{-r^2},dr,d heta < I^2(a) < int_0^{2pi}int_0^{asqrt{2 re^{-r^2},dr,d heta.

(See to polar coordinates from cartesian coordinates for help with polar transformation.)

Integrating,

: pi (1-e^{-a^2}) < I^2(a) < pi (1 - e^{-2a^2}).

By the squeeze theorem, this gives the Gaussian integral

:int_{-infty}^infty e^{-x^2}, dx = sqrt{pi}.

Relation to the gamma function

Since the integrand is an even function,

:int_{-infty}^{infty} e^{-x^2} dx = 2 int_0^infty e^{-x^2} dx

which, after a change of variable, turns into the Euler integral

:int_0^infty e^{-t} t^{-1/2} dt , = , Gammaleft(frac{1}{2} ight)

where Γ is the gamma function. This shows why the factorial of a half-integer is a rational multiple of sqrt pi. More generally,

:int_0^infty e^{-ax^b} dx = a^{-1/b} , Gammaleft(1+frac{1}{b} ight).

Generalizations

The integral of a Gaussian function

The integral of any Gaussian function is reducible in terms of the Gaussian integral

:int_{-infty}^{infty} ae^{-(x+b)^2/c^2},dx.

The constant "a" can be factored out of the integral. Replacing "x" with "y" - "b" yields

:aint_{-infty}^infty e^{-y^2/c^2},dy.

Substituting "y" with "cz" gives

:a |c| int_{-infty}^infty e^{-z^2},dz:=a |c| sqrt{pi}.

"n"-dimensional and functional generalization

Suppose "A" is a symmetric positive-definite invertible covariant tensor of rank two. Then,

:int expleft( - frac 1 2 A_{ij} x^i x^j ight) d^nx=sqrt{frac{(2pi)^n}{det A

where the integral is understood to be over R"n". This fact is applied in the study of the multivariate normal distribution.

Also,

:int x^{k_1}cdots x^{k_{2N , expleft( - frac 1 2 A_{ij} x^i x^j ight) d^nx=sqrt{frac{(2pi)^n}{det A , frac{1}{2^N N!} , sum_{sigma in S_{2N(A^{-1})^{k_{sigma(1)}k_{sigma(2) cdots (A^{-1})^{k_{sigma(2N-1)}k_{sigma(2N)

where σ is a permutation of {1, ..., 2"N"} and the extra factor on the right-hand side is the sum over all combinatorial pairings of {1, ..., 2"N"} of "N" copies of "A"−1.

Alternatively,

:int f(vec x) , expleft( - frac 1 2 A_{ij} x^i x^j ight) d^nx=sqrt{(2pi)^nover det{A , left. expleft({1over 2}(A^{-1})^{ij}{partial over partial x^i}{partial over partial x^j} ight)f(vec{x}) ight|_{vec{x}=0}

for some analytic function "f", provided it satisfies some appropriate bounds on its growth and some other technical criteria. (It works for some functions and fails for others. Polynomials are fine.) The exponential over a differential operator is understood as a power series.

While functional integrals have no rigorous definition (or even a nonrigorous computational one in most cases), we can "define" a Gaussian functional integral in analogy to the finite-dimensional case. There is still the problem, though, that (2pi)^infty is infinite and also, the functional determinant would also be infinite in general. This can be taken care of if we only consider ratios:

:frac{int f(x_1)cdots f(x_{2N}) e^{-iint frac{A(x_{2N+1},x_{2N+2}) f(x_{2N+1}) f(x_{2N+2})}{2} d^dx_{2N+1} d^dx_{2N+2 mathcal{D}f}{int e^{-iint frac{A(x_{2N+1},x_{2N+2}) f(x_{2N+1}) f(x_{2N+2})}{2} d^dx_{2N+1} d^dx_{2N+2 mathcal{D}f},

:=frac{1}{2^N N!}sum_{sigma in S_{2NA^{-1}(x_{sigma(1)},x_{sigma(2)})cdots A^{-1}(x_{sigma(2N-1)},x_{sigma(2N)}).

In the deWitt notation, the equation looks identical to the finite-dimensional case.

"n"-dimensional with linear term

Once again A is a symmetric matrix, then:int e^{-A_{ij} x^i x^j+B_i x_i} d^nx=sqrt{ frac{pi^n}{det{A }e^{frac{1}{4}B^TA^{-1}B}.

Integrals of similar form

:egin{align}&int_0^infty x^{2n} e^{-x^2/a^2},dx &&= sqrt{pi} frac{(2n)!}{n!} left(frac{a}{2} ight)^{2n+1} \&int_0^infty x^{2n+1}e^{-x^2/a^2},dx &&= frac{n!}{2} a^{2n+2}end{align}

Beyond Gaussian Integrals

Exponentials of other even polynomials can easily be solved using series. For example the solution to the integral of the exponential of a quartic polynomial is:

int_{-infty}^{infty}{exp({a x^4+b x^3+c x^2+d x+f})dx} = e^fsum_{n,m,p=0}^{infty}{ frac{ b^{4n{(4n)!}frac{c^{2m{(2m)!}frac{d^{4p{(4p)!} frac{ Gamma(3n+m+p+frac14) }{a^{3n+m+p+frac14} } }.

These integrals turn up in subjects such as quantum field theory.

References

*MathWorld|title=Gaussian Integral|urlname=GaussianIntegral
*David Griffiths. Introduction to Quantum Mechanics. 2nd Edition back cover.


Wikimedia Foundation. 2010.

Игры ⚽ Поможем решить контрольную работу

Look at other dictionaries:

  • Integral — This article is about the concept of integrals in calculus. For the set of numbers, see integer. For other uses, see Integral (disambiguation). A definite integral of a function can be represented as the signed area of the region bounded by its… …   Wikipedia

  • Gaussian function — In mathematics, a Gaussian function (named after Carl Friedrich Gauss) is a function of the form::f(x) = a e^{ { (x b)^2 over 2 c^2 } }for some real constants a > 0, b , c > 0, and e ≈ 2.718281828 (Euler s number).The graph of a Gaussian is a… …   Wikipedia

  • Gaussian orbital — In computational chemistry and molecular physics, Gaussian orbitals (also known as Gaussian type orbitals, GTOs or Gaussians) are functions used as atomic orbitals in the LCAO method for the computation of electron orbitals in molecules.… …   Wikipedia

  • Integral de Gauss — Función gaussiana . El área encerrada bajo esa curva con el eje x es . En matemáticas la integral de Gauss, integral gaussi …   Wikipedia Español

  • Gaussian integer — In number theory, a Gaussian integer is a complex number whose real and imaginary part are both integers. The Gaussian integers, with ordinary addition and multiplication of complex numbers, form an integral domain, usually written as Z[i]. The… …   Wikipedia

  • Gaussian quadrature — In numerical analysis, a quadrature rule is an approximation of the definite integral of a function, usually stated as a weighted sum of function values at specified points within the domain of integration.(See numerical integration for more on… …   Wikipedia

  • Gaussian surface — A Gaussian surface is a closed two dimensional surface through which a flux or electric field is to be calculated. The surface is used in conjunction with Gauss s law (a consequence of the divergence theorem), allowing one to calculate the total… …   Wikipedia

  • Gaussian curvature — In differential geometry, the Gaussian curvature or Gauss curvature of a point on a surface is the product of the principal curvatures, κ 1 and κ 2, of the given point. It is an intrinsic measure of curvature, i.e., its value depends only on how… …   Wikipedia

  • Integral domain — In abstract algebra, an integral domain is a commutative ring that has no zero divisors,[1] and which is not the trivial ring {0}. It is usually assumed that commutative rings and integral domains have a multiplicative identity even though this… …   Wikipedia

  • Integral element — In commutative algebra, an element b of a commutative ring B is said to be integral over its subring A if there are such that That is to say, b is a root of a monic polynomial over A.[1] If B consists of elements that are integral over A, then B… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”