Función gaussiana

ImprimirCitar
Función matemática

En matemáticas, una función gaussiana, a menudo denominada simplemente gaussiana, es una función de la forma base

f()x)=exp⁡ ⁡ ()− − x2){displaystyle f(x)=exp(-x^{2}
f()x)=aexp⁡ ⁡ ()− − ()x− − b)22c2){displaystyle f(x)=aexp left(-{frac {(x-b)^{2}{2c^{2}}}right)}
abcabc

Las funciones gaussianas a menudo se usan para representar la función de densidad de probabilidad de una variable aleatoria normalmente distribuida con valor esperado μ = b y la varianza σ2 = c2. En este caso, la Gaussiana es de la forma

g()x)=1σ σ 2π π exp⁡ ⁡ ()− − 12()x− − μ μ )2σ σ 2).{fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {2fnMicrosoft Sans Serif} }}}exp left(-{2}{2}{frac {(x-mu)}{2}{sigma ^{2}}}derecha). }

Las funciones gaussianas se usan ampliamente en estadística para describir las distribuciones normales, en el procesamiento de señales para definir filtros gaussianos, en el procesamiento de imágenes donde se usan gaussianas bidimensionales para desenfoques gaussianos, y en matemáticas para resolver ecuaciones de calor y ecuaciones de difusión y para definir la transformada de Weierstrass.

Propiedades

Las funciones gaussianas surgen al componer la función exponencial con una función cuadrática cóncava:

f()x)=exp⁡ ⁡ ()α α x2+β β x+γ γ ),{displaystyle f(x)=exp(alpha x^{2}+beta x+gamma),}

  • α α =− − 1/2c2,{displaystyle alpha =-1/2c^{2},}
  • β β =b/c2,{displaystyle beta =b/c^{2},}
  • γ γ =In⁡ ⁡ a− − ()b2/2c2).{displaystyle gamma =ln a-(b^{2}/2c^{2}). }

(Nota: en In⁡ ⁡ a,a=1/()σ σ 2π π ){displaystyle ln a,a=1/(sigma {sqrt {2pi}}}}, no estar confundido con α α =− − 1/2c2,{displaystyle alpha =-1/2c^{2},})

Las funciones gaussianas son, pues, aquellas funciones cuyo logaritmo es una función cuadrática cóncava.

El parámetro c está relacionado con el ancho total a la mitad del máximo (FWHM) del pico según

FWHM=22In⁡ ⁡ 2c.. 2.35482c.{displaystyle {text{FWHM}=2{sqrt {2ln 2},capprox 2.35482,c}

La función se puede expresar en términos de FWHM, representada por w:

f()x)=ae− − 4()In⁡ ⁡ 2)()x− − b)2/w2.{displaystyle f(x)=ae^{-4(ln2)(x-b)^{2}/w^{2}}

Alternativamente, el parámetro c puede interpretarse diciendo que los dos puntos de inflexión de la función ocurren en x = b ± c.

El ancho total a la décima parte del máximo (FWTM) para un gaussiano podría ser de interés y es

FWTM=22In⁡ ⁡ 10c.. 4.29193c.{displaystyle {text{FWTM}=2{sqrt {2ln 10},capprox 4.29193,c}

Las funciones gaussianas son analíticas y su límite como x → ∞ es 0 (para el caso anterior de b = 0).

Las funciones gaussianas se encuentran entre aquellas funciones que son elementales pero carecen de antiderivadas elementales; la integral de la función gaussiana es la función de error:

∫ ∫ e− − x2dx=π π 2er⁡ ⁡ x+C.{displaystyle int e^{-x^{2},dx={frac {sqrt {pi} - ¿Qué?

Sin embargo, sus integrales impropias sobre toda la línea real se pueden evaluar exactamente, utilizando la integral de Gauss

∫ ∫ − − JUEGO JUEGO JUEGO JUEGO e− − x2dx=π π ,{displaystyle int _{-infty ¿Qué?
∫ ∫ − − JUEGO JUEGO JUEGO JUEGO ae− − ()x− − b)2/()2c2)dx=ac⋅ ⋅ 2π π .{displaystyle int _{-infty }{infty }ae^{-(x-b)^{2}/(2c^{2})},dx=accdot {sqrt {2pi}}}

Curvas Gaussianas normalizadas con valor esperado μ y diferencia σ2. Los parámetros correspondientes son a=1σ σ 2π π {textstyle a={tfrac {1}{sigma {sqrt {2pi} }, b = μ y c = σ.

Esta integral es 1 si y sólo si a=1c2π π {textstyle a={tfrac {1}{c{sqrt {2pi}}}} {c} {c}} {c}} {c}}}} {c}}}}} {c}}}} {c}}}}} {c}} {c}} {c} {c}}} {c}}}}}}}}}}}}} {c}}}}}}}}}}}}}}}} {c}}}}}}}}}}} {c}}}}} {c}} {c}}}}}}}}} {c}}}}}}}}}}} {c}}}}}}}}}}}}}}}} {c}}}}}}}}}} {c}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}} (la constante normalizadora), y en este caso el Gaussian es la función de densidad de probabilidad de una variable aleatoria normalmente distribuida con valor esperado μ = b y diferencia σ2 = c2:

g()x)=1σ σ 2π π exp⁡ ⁡ ()− − ()x− − μ μ )22σ σ 2).{fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {2fnMicrosoft Sans Serif} }}}exp left({frac {-x-mu)}{2sigma ^{2}}}right). }

Estas gaussianas se trazan en la figura adjunta.

Las funciones gaussianas centradas en cero minimizan el principio de incertidumbre de Fourier.

El producto de dos funciones gaussianas es un Gausiano, y la convolución de dos funciones gaussianas es también un Gausiano, con varianza siendo la suma de las diferencias originales: c2=c12+c22{displaystyle ¿Qué?. El producto de dos funciones de densidad de probabilidad gaussiana (PDF), sin embargo, no es en general un PDF Gausiano.

Tomando la transformación de Fourier (convención universal de frecuencia angular) de una función Gausiana con parámetros a = 1, b = 0 y c produce otra función Gaussiana, con parámetros c{displaystyle c}, b = 0 y 1/c{displaystyle 1/c}. Así que en particular las funciones gaisianas con b = 0 y c=1{displaystyle c=1} se mantienen fijos por la transformación Fourier (son eigenfunctions de la transformación Fourier con eigenvalue 1). Una realización física es la del patrón de difracción: por ejemplo, una diapositiva fotográfica cuya transmisión tiene una variación gausiana es también una función gausiana.

El hecho de que la función gaussiana sea una función propia de la transformada continua de Fourier nos permite derivar la siguiente identidad interesante de la fórmula de suma de Poisson:

.. k▪ ▪ Zexp⁡ ⁡ ()− − π π ⋅ ⋅ ()kc)2)=c⋅ ⋅ .. k▪ ▪ Zexp⁡ ⁡ ()− − π π ⋅ ⋅ ()kc)2).{displaystyle sum _{kin mathbb {Z} }exp left(-pi cdot left({frac {k}{c}right)^{2}right)=ccdot sum _{kin mathbb {Z}exp left(-pi cdot (kc)^{2}right). }

Integral de una función gaussiana

La integral de una función gaussiana arbitraria es

∫ ∫ − − JUEGO JUEGO JUEGO JUEGO ae− − ()x− − b)2/2c2dx=2aSilenciocSilencioπ π .{displaystyle int _{-infty }{infty }a,e^{-(x-b)^{2}/2c^{2}},dx={sqrt {2}a, sometida a tortura,{sqrt {f} }}

Una forma alternativa es

∫ ∫ − − JUEGO JUEGO JUEGO JUEGO ke− − fx2+gx+hdx=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO ke− − f()x− − g/()2f))2+g2/()4f)+hdx=kπ π fexp⁡ ⁡ ()g24f+h),{displaystyle int _{-infty }{infty }k,e^{-fx^{2}+gx+h},dx=int {f}f} {f}f}f}f}f}f}f}f}f}f}}f}}f}}f}f}, dx=k,{f} {f} {f}} {f}}}f}f}f}f} {f}f}f}f}f}f}f}f}f}f}f}f}f}f}f}f} {f}
f

Relación con la integral gaussiana estándar

La integral

∫ ∫ − − JUEGO JUEGO JUEGO JUEGO ae− − ()x− − b)2/2c2dx{displaystyle int _{-infty }{infty }ae^{-(x-b)^{2}/2c^{2}},dx}
abcaxSí. = xb
a∫ ∫ − − JUEGO JUEGO JUEGO JUEGO e− − Sí.2/2c2dSí.,{displaystyle aint _{-infty - Sí.
z=Sí./2c2{displaystyle z=y/{sqrt {2c^{2}}}
a2c2∫ ∫ − − JUEGO JUEGO JUEGO JUEGO e− − z2dz.{displaystyle a{sqrt {2c^{2}int} ¿Qué?

Entonces, usando la identidad integral de Gauss

∫ ∫ − − JUEGO JUEGO JUEGO JUEGO e− − z2dz=π π ,{displaystyle int _{-infty } {infty }e^{-z^{2},dz={sqrt {pi}}}}

tenemos

∫ ∫ − − JUEGO JUEGO JUEGO JUEGO ae− − ()x− − b)2/2c2dx=a2π π c2.{displaystyle int _{-infty }{infty }ae^{-(x-b)^{2}/2c^{2}},dx=a{sqrt {2pi c^{2}}}}}

Función gaussiana bidimensional

Parcela 3d de una función Gausiana con un dominio bidimensional

Forma básica:

f()x,Sí.)=exp⁡ ⁡ ()− − x2− − Sí.2){displaystyle f(x,y)=exp(-x^{2}-y^{2}

En dos dimensiones, la potencia a la que se eleva e en la función gaussiana es cualquier forma cuadrática definida negativa. En consecuencia, los conjuntos de niveles de la Gaussiana siempre serán elipses.

Un ejemplo particular de una función gaussiana bidimensional es

f()x,Sí.)=Aexp⁡ ⁡ ()− − ()()x− − x0)22σ σ X2+()Sí.− − Sí.0)22σ σ Y2)).{displaystyle f(x,y)=Aexp left(-left({frac {(x-x_{0}^{2}{2sigma {fnK} {fnK} {fnK}} {2}}}{2sigma ¿Sí?

Aquí el coeficiente A es la amplitud, x0, y0 es el centro, y σx, σy son los diferenciales x e y del blob. La figura de la derecha se creó usando A = 1, x0 = 0, y0 = 0, σx = σy = 1.

El volumen bajo la función gaussiana viene dado por

V=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO ∫ ∫ − − JUEGO JUEGO JUEGO JUEGO f()x,Sí.)dxdSí.=2π π Aσ σ Xσ σ Y.{displaystyle V=int _{-infty }{infty }int _{-infty }{infty }f(x,y),dx,dy=2pi Asigma _{X}sigma _{Y}

En general, una función gaussiana elíptica bidimensional se expresa como

f()x,Sí.)=Aexp⁡ ⁡ ()− − ()a()x− − x0)2+2b()x− − x0)()Sí.− − Sí.0)+c()Sí.− − Sí.0)2)),{displaystyle f(x,y)=Aexp {Big (}-{big (}a(x-x_{0})^{2}+2b(x-x_{0})(y-y-y_{0})+c(y-y-y_{0})}{big)}{Big)}}}}} {
[abbc]{displaystyle {begin{bmatrix}a Pulbb

Usando esta formulación, la figura de la derecha se puede crear usando A = 1, (x0, y0) = (0, 0), a = c = 1/2, b = 0.

Significado de parámetros para la ecuación general

Para la forma general de la ecuación, el coeficiente A es la altura del pico y (x0, y0) es el centro de la mancha.

Si configuramos

a=#2⁡ ⁡ Silencio Silencio 2σ σ X2+pecado2⁡ ⁡ Silencio Silencio 2σ σ Y2,b=− − pecado⁡ ⁡ 2Silencio Silencio 4σ σ X2+pecado⁡ ⁡ 2Silencio Silencio 4σ σ Y2,c=pecado2⁡ ⁡ Silencio Silencio 2σ σ X2+#2⁡ ⁡ Silencio Silencio 2σ σ Y2,{displaystyle {begin{aligned}a ventaja={frac {cos ^{2}theta }{2sigma ¿Qué? }{2sigma ¿Por qué? }{4sigma ¿Qué? }{4sigma ¿Por qué? }{2sigma ¿Qué? }{2sigma ¿Qué?
Silencio Silencio {displaystyle theta }b


Para recuperar los coeficientes Silencio Silencio {displaystyle theta }, σ σ X{displaystyle sigma _{X} y σ σ Y{displaystyle sigma _{Y} desde a{displaystyle a}, b{displaystyle b} y c{displaystyle c} uso

Silencio Silencio =12arctan⁡ ⁡ ()2ba− − c),Silencio Silencio ▪ ▪ [− − 45,45],σ σ X2=12()a⋅ ⋅ #2⁡ ⁡ Silencio Silencio +2b⋅ ⋅ #⁡ ⁡ Silencio Silencio pecado⁡ ⁡ Silencio Silencio +c⋅ ⋅ pecado2⁡ ⁡ Silencio Silencio ),σ σ Y2=12()a⋅ ⋅ pecado2⁡ ⁡ Silencio Silencio − − 2b⋅ ⋅ #⁡ ⁡ Silencio Silencio pecado⁡ ⁡ Silencio Silencio +c⋅ ⋅ #2⁡ ⁡ Silencio Silencio ).{displaystyle {begin{aligned}theta {frac {1}{2}}arctan left({frac {2b}{a-c}right),quad theta in [-45,45],\sigma {2cdot cos ^{2}theta +2bcdot cos theta cos theta theta theta theta theta thetatheta ¿Qué? {fnMicrosoft Sans Serif} {2cdot sin ^{2}theta -2bcdot cos theta sin theta {}}


Ejemplos de rotaciones de blobs gaussianos se pueden ver en los siguientes ejemplos:

Silencio Silencio =0{displaystyle theta =0}
Silencio Silencio =− − π π /6{displaystyle theta =-pi /6}
Silencio Silencio =− − π π /3{displaystyle theta =-pi /3}

Usando el siguiente código de Octave, uno puede ver fácilmente el efecto de cambiar los parámetros:

A = 1;x0 = 0; Y0 = 0;sigma_X = 1;sigma_ Y = 2;[X, Y] = meshgrid()-5:.1:5, -5:.1:5);para theta = 0:pi/100:pi a = #()theta)^2 / ()2 * sigma_X^2) + pecado()theta)^2 / ()2 * sigma_ Y^2); b = pecado()2 * theta) / ()4 * sigma_X^2) - pecado()2 * theta) / ()4 * sigma_ Y^2); c = pecado()theta)^2 / ()2 * sigma_X^2) + #()theta)^2 / ()2 * sigma_ Y^2); Z = A * exp()-()a * ()X - x0).^2 + 2 * b * ()X - x0) . * ()Y - Y0) + c * ()Y - Y0).^2)); surf()X, Y, Z); desintegración interp; vista()-36, 36) Espera a Buttonpressfinal

Estas funciones se usan a menudo en el procesamiento de imágenes y en modelos computacionales de la función del sistema visual; consulte los artículos sobre el espacio de escala y la adaptación de formas afines.

Véase también distribución normal multivariante.

Función gaussiana o supergaussiana de orden superior

Se puede tomar una formulación más general de una función gausiana con una caída plana y gaisiana al elevar el contenido del exponente a una potencia P{displaystyle P}:

f()x)=Aexp⁡ ⁡ ()− − ()()x− − x0)22σ σ X2)P).{displaystyle f(x)=Aexp left(-left({frac {(x-x_{0}^{2}}{2sigma Está bien. }

Esta función se conoce como función supergaussiana y se usa a menudo para la formulación de haces gaussianos. Esta función también se puede expresar en términos del ancho total a la mitad del máximo (FWHM), representado por w:

f()x)=Aexp⁡ ⁡ ()− − In⁡ ⁡ 2()4()x− − x0)2w2)P).{displaystyle f(x)=Aexp left(-ln 2left(4{frac {(x-x_{0}}{2}}{2}}right)^{P}right). }

En una formulación bidimensional, una función gaissa a lo largo x{displaystyle x} y Sí.{displaystyle y} se puede combinar con potencialmente diferente PX{displaystyle P_{X} y PY{displaystyle P_{Y} para formar una distribución triangular rectangular:

f()x,Sí.)=Aexp⁡ ⁡ ()− − ()()x− − x0)22σ σ X2)PX− − ()()Sí.− − Sí.0)22σ σ Y2)PY).{displaystyle f(x,y)=Aexp left(left({frac {(x-x_{0}^{2}{2sigma ¿Por qué? ¿Sí?
f()x,Sí.)=Aexp⁡ ⁡ ()− − ()()x− − x0)22σ σ X2+()Sí.− − Sí.0)22σ σ Y2)P){displaystyle f(x,y)=Aexp left(-left({frac {(x-x_{0}^{2}{2sigma {fnK} {fnK} {fnK}} {2}}}{2sigma ¿Qué?

Función gaussiana multidimensional

En un n{displaystyle n}- espacio dimensional una función gaissa se puede definir como

f()x)=exp⁡ ⁡ ()− − xTCx),{displaystyle f(x)=exp(-x^{mathsf {T}Cx),}
x=[x1⋯ ⋯ xn]{displaystyle x={begin{bmatrix}x_{1} ¿Qué?n{displaystyle n}C{displaystyle C}n× × n{displaystyle ntimes n}T{displaystyle {} {fnMithsf}}

La parte integral de esta función gaisiana sobre todo n{displaystyle n}- espacio dimensional se da como

∫ ∫ Rnexp⁡ ⁡ ()− − xTCx)dx=π π nDetC.{displaystyle int _{mathbb [R] ^{n}exp(-x^{mathsf {T}Cx},dx={sqrt {frac {pi ^{n}{det - Sí.

Se puede calcular fácilmente diagonalizando la matriz C{displaystyle C} y cambiar las variables de integración a los eigenvectores C{displaystyle C}.

Más generalmente, una función gaussiana desplazada se define como

f()x)=exp⁡ ⁡ ()− − xTCx+sTx),{displaystyle f(x)=exp(-x^{mathsf {T}Cx+s^{mathsf {T}x),}
s=[s1⋯ ⋯ sn]{displaystyle s={begin{bmatrix}s_{1} {fn}}C{displaystyle C}CT=C{displaystyle C^{mathsf {T}=C}
∫ ∫ Rne− − xTCx+vTxdx=π π nDetCexp⁡ ⁡ ()14vTC− − 1v)↑ ↑ M.{displaystyle int _{mathbb {fnh}e^{-x^{mthsf {T}}Cx+v^{mathsf {T}x},dx={sqrt {frac {pi ^{n}{det {C}}exp left({frac} {1}{4}v^{mathsf {T}C^{-1}vright)equiv {mathcal {M}}
∫ ∫ Rne− − xTCx+vTx()aTx)dx=()aTu)⋅ ⋅ M,Dondeu=12C− − 1v.{displaystyle int _{mathbb {fnh}e^{-x^{mthsf {T}}Cx+v^{mathsf {T}x}(a^{mathsf {T}x),dx=(a^{T}u)cdot {mathcal {M},{text{ Where }u={frac C^{-1}v.}
∫ ∫ Rne− − xTCx+vTx()xTDx)dx=()uTDu+12tr⁡ ⁡ ()DC− − 1))⋅ ⋅ M.{displaystyle int _{mathbb {fnh}e^{-x^{mthsf {T}}Cx+v^{mathsf {T}x}(x^{mathsf {T}Dx),dx=left(u^{mathsf {f}Du+frac {1}{2}operatorname {tr} (DC^{-1})derecha)cdot {mathcal] {M}}
∫ ∫ Rne− − xTC.x+s.Tx()− − ∂ ∂ ∂ ∂ x▪ ▪ ∂ ∂ ∂ ∂ x)e− − xTCx+sTxdx=()2tr⁡ ⁡ ()C.▪ ▪ CB− − 1)+4uTC.▪ ▪ Cu− − 2uT()C.▪ ▪ s+C▪ ▪ s.)+s.T▪ ▪ s)⋅ ⋅ M,{displaystyle {begin{aligned} ¿Qué? {fnh}e^{-x^{mthsf {T}}C'x+s'^{mathsf {T}x}left(-{frac {partial }{partial # Lambda {frac {partial }{partial x}right)e^{-x^{mathsf Cx+s. {T}x},dx\\\fnMicrosoft} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {f}}}}Lambda Cu-2u} {Mathsf} {cHFF}(C'Lambda s+CLambda s')+s'^{mathigned}Lambda send}
u=12B− − 1v,v=s+s.,B=C+C..{textstyle u={2}B^{-1}v, v=s+s', B=C+C'}

Estimación de parámetros

Varios campos como la fotometría estelar, la caracterización de haz Gaussian y el trabajo de espectroscopia de línea de emisión/absorción con funciones Gaussianas demostradas y necesitan calcular con precisión los parámetros de altura, posición y anchura de la función. Hay tres parámetros desconocidos para una función gaisiana 1D (a, b, c) y cinco para una función de Gaussian 2D ()A;x0,Sí.0;σ σ X,σ σ Y){displaystyle (A;x_{0},y_{0};sigma _{X},sigma _{Y})}.

El método más común para estimar los parámetros gaussianos es tomar el logaritmo de los datos y ajustar una parábola al conjunto de datos resultante. Si bien esto proporciona un procedimiento de ajuste de curva simple, el algoritmo resultante puede estar sesgado al ponderar excesivamente los valores de datos pequeños, lo que puede producir grandes errores en la estimación del perfil. Uno puede compensar parcialmente este problema a través de la estimación de mínimos cuadrados ponderados, reduciendo el peso de los valores de datos pequeños, pero esto también puede estar sesgado al permitir que la cola de la Gaussiana domine el ajuste. Para eliminar el sesgo, se puede usar un procedimiento de mínimos cuadrados iterativamente reponderados, en el que los pesos se actualizan en cada iteración. También es posible realizar una regresión no lineal directamente sobre los datos, sin involucrar la transformación logarítmica de datos; para más opciones, vea ajuste de distribución de probabilidad.

Precisión de parámetros

Una vez que se tiene un algoritmo para estimar los parámetros de la función gaussiana, también es importante saber cuán precisas son esas estimaciones. Cualquier algoritmo de estimación de mínimos cuadrados puede proporcionar estimaciones numéricas para la varianza de cada parámetro (es decir, la varianza de la altura, la posición y el ancho estimados de la función). También se puede utilizar la teoría del límite de Cramér-Rao para obtener una expresión analítica para el límite inferior de las varianzas de los parámetros, dadas ciertas suposiciones sobre los datos.

  1. El ruido del perfil medido es i.i.d. Gaussian, o el ruido es distribuido por Poisson.
  2. El espaciado entre cada muestreo (es decir, la distancia entre píxeles que miden los datos) es uniforme.
  3. El pico es "bien muestreado", por lo que menos del 10% de la zona o volumen bajo el pico (rea si un Gaussian 1D, volumen si un Gaussian 2D) se encuentra fuera de la región de medición.
  4. El ancho del pico es mucho más grande que la distancia entre las ubicaciones de muestra (es decir, los píxeles detectores deben ser al menos 5 veces más pequeños que el Gaussian FWHM).

Cuando estas suposiciones estén satisfechas, la siguiente matriz de covariancia K aplica para los parámetros de perfil 1D a{displaystyle a}, b{displaystyle b}, y c{displaystyle c} bajo i.i.d. Ruido gausiano y bajo ruido Poisson:

KGauss=σ σ 2π π δ δ XQ2()32c0− − 1a02ca20− − 1a02ca2),KPoiss=12π π ()3a2c0− − 120ca0− − 120c2a),{displaystyle mathbf {K} _{text{Gauss}={frac {sigma ^{2}{sqrt {pi}delta {X}Q^{2} {begin{pmatrix}{frac} {3}{2c} {2c} {2c}} {0} {2c} {2c} {2}}}} {0\\\fn} {f}} {f} {f}}} {f}}}}} {\f} {f}}}} {f}} {f}}}}}}}}}}}} {\\\f}}}}}}}}}}}}}\\\\f} {\f}}\\\\\\\\\f}}}}}\\\\\\\\\\\\fn}}\\\\\\\\fn}}}}}}}}}}}\\\\\\\\\\fn}}}}}}}}} {-1}{a}} {2c}}end{pmatrix}qquad {fnK} {f} {fnMicroc {1}{sqrt {2pi} # {begin{pmatrix}{frac} {3a}{2c} {0} {2}} {c}{a} {0} {0} {f} {0f} {f}} {f}f}} {f} {f}f} {f} {f} {f}} {f}f}}f} {f}f}}}}}}}f}f} {f} {f}}f}f}f}f}f}f}f}f}f} {f}f}f}f}f}f}f}f}f}f}f} {f}f}f}f}f}f}f}f}f}f}f}f}f}f}f}fn}f}f}f}f}f}f}f}f}f}f {1}{2} {0} {c}}end{2a} {c}}f}}}} {f}}}} {f}}}}}} {f}}}}} {f}}}}}}}}} {}}} {}}}}} {}}}}} {}}}}}}}} {}}}}}} {}}}}}}}}}}}}}}}}}} {}}}}}}}}}}}}}}}}}}}}}}}}}}}}} {}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}} {}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}} {}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}
δ δ X{displaystyle delta ¿Qué?Q{displaystyle Q}σ σ {displaystyle sigma }
Var⁡ ⁡ ()a)=3σ σ 22π π δ δ XQ2cVar⁡ ⁡ ()b)=2σ σ 2cδ δ Xπ π Q2a2Var⁡ ⁡ ()c)=2σ σ 2cδ δ Xπ π Q2a2{displaystyle {begin{aligned}operatorname {var} (a) {3sigma ^{2}{2{sqrt {pi },delta ¿Por qué? {2sigma ^{2}c}{delta ¿Qué? ###operatorname {var} (c) âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa {var âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa âTMa â {2sigma ^{2}c}{delta ¿Qué? },Q^{2}}end{aligned}}

y en el caso del ruido de Poisson,

Var⁡ ⁡ ()a)=3a22π π cVar⁡ ⁡ ()b)=c2π π aVar⁡ ⁡ ()c)=c22π π a.{displaystyle {begin{aligned}operatorname {var} (a) {3a}{2{sqrt {2pi} - ¿Qué? {C}{sqrt {2pi} ####,a}\operatorname {var} (c) Sentirme={frac {c}{2{sqrt {2pi} },a}end{aligned}}

Para los parámetros de perfil 2D que dan la amplitud A{displaystyle A}, posición ()x0,Sí.0){displaystyle (x_{0},y_{0}}, y ancho ()σ σ X,σ σ Y){displaystyle (sigma _{X},sigma - Sí. del perfil, se aplican las siguientes matrices de covariancia:

KGauss=σ σ 2π π δ δ Xδ δ YQ2()2σ σ Xσ σ Y00− − 1Aσ σ Y− − 1Aσ σ X02σ σ XA2σ σ Y000002σ σ YA2σ σ X00− − 1Aσ σ Sí.002σ σ XA2σ σ Sí.0− − 1Aσ σ X0002σ σ YA2σ σ X)KPoisson=12π π ()3Aσ σ Xσ σ Y00− − 1σ σ Y− − 1σ σ X0σ σ XAσ σ Y00000σ σ YAσ σ X00− − 1σ σ Y002σ σ X3Aσ σ Y13A− − 1σ σ X0013A2σ σ Y3Aσ σ X).{displaystyle {begin{aligned}mathbf {K} _{text{Gauss}={frac {sigma ^{2}{pi delta ¿Por qué? ################################################################################################################################################################################################################################################################ {2}{sigma _{X}sigma - ¿Qué? {-1}{Asigma - ¿Qué? {-1}{Asigma - ¿Por qué? {fnK}{2}sigma ################################################################################################################################################################################################################################################################ - ¿Qué? - ¿Qué? {-1}{Asigma {2sigma} {fnK}{2}sigma ¿Qué? {-1}{Asigma - ¿Qué? {2sigma - ¿Qué? ¿Por qué? {Poisson}={frac {1}{2pi} {fnMicroc {3A}{sigma} _{X}sigma ################################################################################################################################################################################################################################################################ - ¿Qué? {-1}{sigma - ¿Por qué? {sigma {X}{Asigma ¿Qué? {sigma {Y}{Asigma - ¿Qué? {-1}{sigma {2sigma} {X}{3Asigma - ¿Qué? {1}{3A}\\\fnMicroc {-1}{sigma ################################################################################################################################################################################################################################################################ {2sigma {Y}{3Asigma {fnK}end{pmatrix}}}end{aligned}}

Gaussiana discreta

El núcleo discreto Gaussiano (sólido), en comparación con el núcleo Gaussiano muestrado (revisado) para escalas t=0.5,1,2,4.{displaystyle t=0.5,1,2,4}

Se puede pedir un análogo discreto al Gaussiano; esto es necesario en aplicaciones discretas, particularmente en el procesamiento de señales digitales. Una respuesta simple es muestrear la Gaussiana continua, produciendo el kernel gaussiano muestreado. Sin embargo, esta función discreta no tiene los análogos discretos de las propiedades de la función continua y puede provocar efectos no deseados, como se describe en el artículo Implementación del espacio de escala.

Un enfoque alternativo es utilizar el núcleo gaussiano discreto:

T()n,t)=e− − tIn()t){displaystyle T(n,t)=e^{-t}I_{n}(t)}
In()t){displaystyle I_{n}(t)}

Este es el análogo discreto de la Gaussiana continua en que es la solución a la ecuación de difusión discreta (espacio discreto, tiempo continuo), así como la Gaussiana continua es la solución a la ecuación de difusión continua.

Aplicaciones

Las funciones gaussianas aparecen en muchos contextos en las ciencias naturales, las ciencias sociales, las matemáticas y la ingeniería. Algunos ejemplos incluyen:

  • En las estadísticas y la teoría de probabilidad, las funciones gausianas aparecen como la función de densidad de la distribución normal, que es una distribución de probabilidad límite de sumas complicadas, según el teorema límite central.
  • Las funciones gausianas son la función de Green para la ecuación de difusión (homogénea e isotrópica) (y para la ecuación de calor, que es lo mismo), una ecuación diferencial parcial que describe la evolución del tiempo de una densidad de masa bajo la difusión. Específicamente, si la densidad de masa a tiempo t=0 es dado por un delta Dirac, lo que significa esencialmente que la masa se concentra inicialmente en un solo punto, luego la distribución masiva a la vez t será dada por una función Gaussiana, con el parámetro a estar relacionado linealmente con 1/t y c estar relacionado linealmente con t; este Gaussian que varia tiempo es descrito por el núcleo de calor. Más generalmente, si la densidad de masa inicial es φ(x), entonces la densidad de masa en tiempos posteriores se obtiene tomando la convolución de φ con una función Gausian. La convolución de una función con un Gaussian también se conoce como una transformación Weierstrass.
  • Una función Gaussiana es la función de onda del estado del suelo del oscilador armónico cuántico.
  • Los orbitales moleculares utilizados en la química computacional pueden ser combinaciones lineales de funciones gausianas llamadas orbitales gaussianas (ver también el conjunto de bases (química)).
  • Matemáticamente, los derivados de la función Gaussiana pueden ser representados usando funciones Hermite. Para la varianza unitaria, n- el derivado del Gaussian es la función gaussiana en sí multiplicada por el n- polinomio Hermite, hasta escala.
  • En consecuencia, las funciones gauchas también están asociadas con el estado de vacío en la teoría del campo cuántico.
  • Las vigas gausianas se utilizan en sistemas ópticos, sistemas de microondas y láser.
  • En la representación espacial de escala, las funciones gaussianas se utilizan como núcleos de suavizado para generar representaciones multiescala en la visión informática y el procesamiento de imágenes. Específicamente, los derivados de Gaussianos (funciones hermitas) se utilizan como base para definir un gran número de tipos de operaciones visuales.
  • Las funciones gausianas se utilizan para definir algunos tipos de redes neuronales artificiales.
  • En la microscopía de fluorescencia se utiliza una función Gaussiana 2D para aproximar el disco Airy, describiendo la distribución de intensidad producida por una fuente de puntos.
  • En el procesamiento de señales sirven para definir filtros gaussianos, como en el procesamiento de imágenes donde se utilizan Gaussians 2D para bordas gausianas. En el procesamiento digital de señales, se utiliza un núcleo Gaussiano discreto, que puede definirse mediante muestreo de un Gausian o de otra manera.
  • En geoestadística se han utilizado para comprender la variabilidad entre los patrones de una imagen compleja de entrenamiento. Se utilizan con métodos del núcleo para agrupar los patrones en el espacio de características.

Contenido relacionado

Grupo nilpotente

En matemáticas, específicamente en teoría de grupos, un grupo nilpotente G es un grupo que tiene una serie central superior que termina en G. De manera...

Esponja Menger

En matemáticas, la esponja de Menger es una curva fractal. Es una generalización tridimensional del conjunto de Cantor unidimensional y la alfombra de...

Tomas Fincke

Thomas Fincke fue un matemático y físico danés, y profesor en la Universidad de Copenhague durante más de 60...
Más resultados...
Tamaño del texto:
Copiar