Gaussian function

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 characteristic symmetric "bell shape curve" that quickly falls off towards plus/minus infinity. The parameter "a" is the height of the curve's peak, "b" is the position of the center of the peak, and "c" controls the width of the "bell".

Gaussian functions are widely used in statistics where they describe the normal distributions, in signal processing where they serve to define Gaussian filters, in image processing where two-dimensional Gaussians are used for Gaussian blurs, and in mathematics where they are used to solve heat equations and diffusion equations and to define the Weierstrass transform.

Properties

Gaussian functions arise by applying the exponential function to a general quadratic function. The Gaussian functions are thus those functions whose logarithm is a quadratic function.

The parameter "c" is related to the full width at half maximum (FWHM) of the peak according to

: mathrm{FWHM} = 2 sqrt{2 ln 2} c = 2.35482... cdot c.

Alternatively, the parameter "c" can be interpreted by saying that the two inflection points of the function occur at "x"="b"−"c" and "x"="b"+"c".

Gaussian functions are analytic, and their limit for "x"→±∞ is 0. Gaussian functions are among those functions that are elementary but lack elementary antiderivatives; the integral of the Gaussian function is the error function. Nonetheless their improper integrals over the whole real line can be evaluated exactly, using the Gaussian integral

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

and one obtains

:int_{-infty}^infty a e^{- { (x-b)^2 over 2 c^2 } },dx=accdotsqrt{2pi}

This integral is 1 if and only if "a" = 1/("c"√(2π)), and in this case the Gaussian is the probability density function of a normally distributed random variable with expected value μ="b" and variance σ2="c"2. These Gaussians are graphed in the accompanying figure.

Taking the Fourier transform of a Gaussian function with parameters "a", "b"=0 and "c" yields another Gaussian function, with parameters "ac", "b"=0 and 1/"c". So in particular the Gaussian functions with "b"=0 and "c"=1 are kept fixed by the Fourier transform (they are eigenfunctions of the Fourier transform with eigenvalue 1).

The product of two Gaussian functions is again a Gaussian, and the convolution of two Gaussian functions is again a Gaussian.

Two-dimensional Gaussian function

A particular example of a two-dimensional Gaussian function is :f(x,y) = A e^{- left(frac{(x-x_o)^2}{2sigma_x^2} + frac{(y-y_o)^2}{2sigma_y^2} ight)}.

Here the coefficient "A" is the amplitude, "x"o,yo is the center and σ"x", σ"y" are the "x" and "y" spreads of the blob. The figure on the left was created using "A" = 1, "x"o = 0, "y"o = 0, σ"x" = σ"y" = 1.

In general, a two-dimensional Gaussian function is expressed as

:f(x,y) = A e^{- left(a(x - x_o)^2 + 2b(x-x_o)(y-y_o) + c(y-y_o)^2 ight)}

where the matrix

:left [egin{matrix} a & b \ b & c end{matrix} ight]

is positive-definite.

Using this formulation, the figure on the left can be created using "A" = 1, ("x"o, "y"o) = (0, 0), "a" = "c" = 1, "b" = 0.

Meaning of parameters for the general equation

For the general form of the equation the coefficient "A" is the height of the peak and ("x"o, "y"o) is the center of the blob.

If we set

:a = frac{cos^2 heta}{2sigma_x^2} + frac{sin^2 heta}{2sigma_y^2}

:b = -frac{sin2 heta}{4sigma_x^2} + frac{sin2 heta}{4sigma_y^2}

:c = frac{sin^2 heta}{2sigma_x^2} + frac{cos^2 heta}{2sigma_y^2}

then we rotate the blob by an angle heta. This can be seen in the following examples:

Using the following MATLAB code one can see the effect of changing the parameters easily

A = 1;x0 = 0; y0 = 0;

sigma_x = 1;sigma_y = 2;

for theta = 0:pi/100:pia = cos(theta)^2/2/sigma_x^2 + sin(theta)^2/2/sigma_y^2;b = -sin(2*theta)/4/sigma_x^2 + sin(2*theta)/4/sigma_y^2 ;c = sin(theta)^2/2/sigma_x^2 + cos(theta)^2/2/sigma_y^2;

[X, Y] = meshgrid(-5:.1:5, -5:.1:5);Z = A*exp( - (a*(X-x0).^2 + 2*b*(X-x0).*(Y-y0) + c*(Y-y0).^2)) ;surf(X,Y,Z);shading interp;view(-36,36);axis equal;drawnowend

Such functions are often used in image processing and in computational models of visual system function -- see the articles on scale space and affine shape adaptation.

Also see multivariate normal distribution.

Applications

Gaussian functions appear in many contexts in the natural sciences, the social sciences, mathematics, and engineering. Some examples include:
* In statistics and probability theory, Gaussian functions appear as the density function of the normal distribution, which is a limiting probability distribution of complicated sums, according to the central limit theorem.
* Gaussian functions are closely related to the (homogeneous and isotropic) diffusion equation (and, which is the same thing, to the heat equation), a partial differential equation that describes the time evolution of a mass-density under diffusion. Specifically, if the mass-density at time "t"=0 is given by a Dirac delta, which essentially means that the mass is initially concentrated in a single point, then the mass-distribution at time "t" will be given by a Gaussian function, with the parameter "a" being linearly related to 1/√"t" and "c" being linearly related to √"t". More generally, if the initial mass-density is φ("x"), then the mass-density at later times is obtained by taking the convolution of φ with a Gaussian function. The convolution of a function with a Gaussian is also known as a Weierstrass transform.
* A Gaussian function is the wave function of the ground state of the quantum harmonic oscillator.
* The molecular orbitals used in computational chemistry can be linear combinations of Gaussian functions called Gaussian orbitals (see also basis set (chemistry)).
* Mathematically, the derivatives of the Gaussian function are used to define Hermite polynomials.
* Consequently, Gaussian functions are also associated with the vacuum state in quantum field theory.
* Gaussian beams are used in optical and microwave systems,
* Gaussian functions are used as smoothing kernels for generating multi-scale representations in computer vision and image processing -- see the article on scale space representation. Specifically, derivatives of Gaussians are used as a basis for defining a large number of types of visual operations.
* Gaussian functions are used to define some types of artificial neural networks.
* In fluorescence microscopy a 2D Gaussian function is used to approximate the Airy disk, describing the intensity distribution produced by a point source.

ee also

*Lorentzian function
*Multivariate normal distribution

External links

* [http://mathworld.wolfram.com/GaussianFunction.html Mathworld, includes a proof for the relations between c and FWHM]


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • Gaussian blur — Gaussian blur, named for Carl Friedrich Gauss, describes blurring an image by a Gaussian function. It is a widely used effect in graphics software, typically to reduce image noise and reduce detail. The visual effect of this blurring technique is …   Wikipedia

  • Gaussian broadening — refers to broadening effects in spectral lines, these can be produced by Doppler broadening and natural broadening.This effect is similar to Gaussian blur effect in image processing produced by convolution with the Gaussian function.The term is… …   Wikipedia

  • Gaussian beam — In optics, a Gaussian beam is a beam of electromagnetic radiation whose transverse electric field and intensity (irradiance) distributions are well approximated by Gaussian functions. Many lasers emit beams that approximate a Gaussian profile, in …   Wikipedia

  • 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 {… …   Wikipedia

  • Gaussian filter — In electronics and signal processing, a Gaussian filter is a filter whose filter window is the Gaussian function:x mapsto sqrt{frac{a}{picdot e^{ a cdot x^2}or with the standard deviation as parameter:x mapsto frac{1}{sqrt{2cdotpi}cdotsigma}cdot… …   Wikipedia

  • Gaussian noise — is noise that has a probability density function (abbreviated pdf) of the normal distribution (also known as Gaussian distribution). In other words, the values that the noise can take on are Gaussian distributed. It is most commonly used as… …   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 network model — The Gaussian network model (GNM), one of many things named after Carl Gauss, is a representation of a biological macromolecule as an elastic mass and spring network to study, understand, and characterize mechanical aspects of its long scale… …   Wikipedia

  • Gaussian adaptation — Articleissues citations missing = July 2008 COI = y expert = Mathematics notability = July 2008 jargon = July 2008 OR = September 2007 primarysources = July 2008 technical = July 2008Gaussian adaptation (GA) is an evolutionary algorithm designed… …   Wikipedia

  • Gaussian period — In mathematics, in the area of number theory, a Gaussian period is a certain kind of sum of roots of unity. They permit explicit calculations in cyclotomic fields, in relation both with Galois theory and with harmonic analysis (discrete Fourier… …   Wikipedia

Share the article and excerpts

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