In der Mathematik ist eine Gauß-Funktion , oft einfach als Gauß-Funktion bezeichnet , eine Funktion der Grundform mit parametrischer Erweiterung für beliebige reelle Konstanten a , b und c ≠ 0. Sie ist nach dem Mathematiker Carl Friedrich Gauß benannt . Der Graph einer Gauß-Funktion hat die charakteristische symmetrische Form einer Glockenkurve . Der Parameter a ist die Höhe des Scheitelpunkts der Kurve, b die horizontale Position des Scheitelpunkts und c (die Standardabweichung , auch als Gaußsche Effektivwertbreite bezeichnet ) die Breite der Glockenkurve.

Gauß-Funktionen werden häufig verwendet, um die Wahrscheinlichkeitsdichtefunktion einer normalverteilten Zufallsvariablen mit Erwartungswert μ = b und Varianz σ² = c² darzustellen . In diesem Fall hat die Gauß-Funktion die Form [ 1 ] .

Gauß-Funktionen finden breite Anwendung in der Statistik zur Beschreibung von Normalverteilungen , in der Signalverarbeitung zur Definition von Gauß-Filtern , in der Bildverarbeitung, wo zweidimensionale Gauß-Funktionen für Gaußsche Unschärfen verwendet werden , und in der Mathematik zur Lösung von Wärmeleitungs- und Diffusionsgleichungen sowie zur Definition der Weierstrass-Transformation . Sie werden außerdem häufig in der Quantenchemie zur Bildung von Basissätzen eingesetzt .

Eigenschaften

Gaußsche Funktionen entstehen durch die Komposition der Exponentialfunktion mit einer konkaven quadratischen Funktion : wobei (Anmerkung: in , nicht zu verwechseln mit )

Die Gaußschen Funktionen sind also diejenigen Funktionen, deren Logarithmus eine konkave quadratische Funktion ist.

Der Parameter c steht in Zusammenhang mit der Halbwertsbreite (FWHM) des Peaks gemäß folgender Formel:

Die Funktion kann dann mithilfe der FWHM, dargestellt durch w , ausgedrückt werden :

Alternativ kann der Parameter c auch so interpretiert werden, dass die beiden Wendepunkte der Funktion bei x = b ± c auftreten .

Die volle Breite bei einem Zehntel des Maximums (FWTM) einer Gauß-Funktion könnte von Interesse sein und ist

Gaußsche Funktionen sind analytisch , und ihr Grenzwert für x → ∞ ist 0 (für den obigen Fall b = 0 ).

Gaußsche Funktionen gehören zu den Funktionen, die elementar sind , aber keine elementaren Stammfunktionen besitzen ; das Integral der Gaußschen Funktion ist die Fehlerfunktion :

Dennoch lassen sich ihre uneigentlichen Integrale über die gesamte reelle Zahlengerade mithilfe des Gaußschen Integrals exakt auswerten, und man erhält:

Normalisierte Gaußsche Kurven mit Erwartungswert μ und Varianz σ² . Die entsprechenden Parameter sind b = μ und c = σ .

Dieses Integral ist genau dann 1, wenn (die Normierungskonstante ), und in diesem Fall ist die Gauß- Funktion die Wahrscheinlichkeitsdichtefunktion einer normalverteilten Zufallsvariablen mit Erwartungswert μ = b und Varianz σ 2 = c 2 :

Diese Gauß-Funktionen sind in der beigefügten Abbildung dargestellt.

Das Produkt zweier Gauß-Funktionen ist wieder eine Gauß-Funktion, und auch die Faltung zweier Gauß-Funktionen ist wieder eine Gauß-Funktion, deren Varianz der Summe der ursprünglichen Varianzen entspricht. Das Produkt zweier Gaußscher Wahrscheinlichkeitsdichtefunktionen (PDFs) ist jedoch im Allgemeinen keine Gaußsche PDF.

Das Fourier- Unschärfeprinzip wird genau dann zu einer Gleichung, wenn (modulierte) Gauß-Funktionen betrachtet werden. [ 2 ]

Die Fouriertransformation (unitär, Kreisfrequenzkonvention) einer Gauß-Funktion mit den Parametern a = 1 , b = 0 und c ergibt eine andere Gauß-Funktion mit den Parametern , b = 0 und ⁠1/C[ 3 ] Insbesondere bleiben die Gauß-Funktionen mit b = 0 und c = a durch die Fouriertransformation fixiert (sie sind Eigenfunktionen der Fouriertransformation mit Eigenwert 1). Ein physikalisches Beispiel hierfür ist das Beugungsmuster : Beispielsweise ist ein fotografischer Dia,dessen Transmission eine Gaußsche Variation aufweist, ebenfalls eine Gauß-Funktion.

Die Tatsache, dass die Gauß-Funktion eine Eigenfunktion der kontinuierlichen Fourier-Transformation ist, erlaubt es uns, aus der Poisson-Summationsformel die folgende interessante Identität abzuleiten :

Integral einer Gauß-Funktion

Das Integral einer beliebigen Gauß-Funktion ist

Eine alternative Form ist diejenige, bei der f strikt positiv sein muss, damit das Integral konvergiert.

Beziehung zum Standard-Gauß-Integral

Das Integral für reelle Konstanten a , b , c > 0 lässt sich berechnen, indem man es in die Form eines Gauß-Integrals bringt . Zunächst kann die Konstante a einfach aus dem Integral ausgeklammert werden. Anschließend wird die Integrationsvariable von x auf y = x − b und dann auf y = ...

Anschließend verwendet man die Identität des Gaußschen Integrals.

Wir haben

Zweidimensionale Gauß-Funktion

3D-Darstellung einer Gauß-Funktion mit einem zweidimensionalen Definitionsbereich

Grundform:

In zwei Dimensionen ist das Argument der Exponentialfunktion in der Gauß-Funktion eine beliebige negativ definite quadratische Form . Folglich sind die Niveaumengen der Gauß-Funktion stets Ellipsen.

Ein konkretes Beispiel für eine zweidimensionale Gauß-Funktion ist

Hierbei ist der Koeffizient A die Amplitude, x₀ , y₀ der Mittelpunkt und σₓ , σᵧ die Ausdehnung des Bereichs in x- bzw. y- Richtung  .  Die Abbildung rechts wurde mit A = 1, x₀ = 0 , y₀ = 0 und σₓ = σᵧ = 1 erstellt .

Das Volumen unter der Gauß-Funktion ist gegeben durch

Im Allgemeinen wird eine zweidimensionale elliptische Gauß-Funktion wie folgt ausgedrückt: wobei die Matrix positiv definit ist .

Mit dieser Formulierung kann die Figur rechts mit A = 1 , ( x 0 , y 0 ) = (0, 0) , a = c = ⁠ erstellt werden.1/2, b = 0.

Bedeutung der Parameter für die allgemeine Gleichung

Bei der allgemeinen Form der Gleichung ist der Koeffizient A die Höhe des Gipfels und ( x 0 , y 0 ) der Mittelpunkt des Flecks.

Setzen wir dann , dann drehen wir den Blob um einen positiven Winkel gegen den Uhrzeigersinn (für eine negative Drehung im Uhrzeigersinn kehren wir die Vorzeichen des b- Koeffizienten um). [ 4 ]

Um die Koeffizienten und aus und zurückzuerhalten , verwenden Sie

Beispiele für Rotationen von Gaußschen Flächen können in den folgenden Beispielen gesehen werden:

Anhand des folgenden Octave -Codes lässt sich der Effekt einer Parameteränderung leicht erkennen:

A = 1 ; x0 = 0 ; y0 = 0 ;sigma_X = 1 ; sigma_Y = 2 ;[ X , Y ] = meshgrid ( - 5 :. 1 : 5 , - 5 :. 1 : 5 );für theta = 0 : pi / 100 : pi a = 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 );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 ); waitforbuttonpress ; end

Solche Funktionen werden häufig in der Bildverarbeitung und in Computermodellen der Funktion des visuellen Systems verwendet – siehe die Artikel über Skalenraum und affine Formadaptation .

Siehe auch multivariate Normalverteilung .

Gaußsche Funktion höherer Ordnung oder Super-Gaußsche Funktion oder verallgemeinerte Gaußsche Funktion

Eine allgemeinere Formulierung einer Gauß-Funktion mit flachem Maximum und Gaußschem Abfall erhält man durch Potenzieren des Exponenten :

Diese Funktion ist als Super-Gauß-Funktion bekannt und wird häufig zur Formulierung von Gaußschen Strahlen verwendet. [ 5 ] Diese Funktion kann auch durch die Halbwertsbreite (FWHM), dargestellt durch w , ausgedrückt werden :

In einer zweidimensionalen Formulierung kann eine Gauß-Funktion entlang und mit potenziell unterschiedlichen und kombiniert werden , um eine rechteckige Gauß-Verteilung oder eine elliptische Gauß-Verteilung zu bilden.

Mehrdimensionale Gauß-Funktion

In einem -dimensionalen Raum kann eine Gaußsche Funktion definiert werden als , wobei eine Spalte von Koordinaten, eine positiv definite Matrix und die Matrix-Transposition bezeichnet .

Das Integral dieser Gauß-Funktion über den gesamten -dimensionalen Raum ist gegeben als

Sie lässt sich leicht berechnen, indem man die Matrix diagonalisiert und die Integrationsvariablen durch die Eigenvektoren von ersetzt .

Allgemeiner ist eine verschobene Gauß-Funktion definiert als , wobei der Verschiebungsvektor ist und die Matrix als symmetrisch, und positiv definit angenommen werden kann. Die folgenden Integrale mit dieser Funktion können mit derselben Methode berechnet werden:

Wo

Schätzung der Parameter

Zahlreiche Forschungsbereiche wie die Sternphotometrie , die Charakterisierung von Gaußschen Strahlen und die Emissions-/Absorptionslinienspektroskopie arbeiten mit abgetasteten Gauß-Funktionen und benötigen eine genaue Bestimmung der Parameter Höhe, Position und Breite dieser Funktion. Eine eindimensionale Gauß-Funktion weist drei unbekannte Parameter ( a , b , c ) auf, eine zweidimensionale hingegen fünf .

The most common method for estimating the Gaussian parameters is to take the logarithm of the data and fit a parabola to the resulting data set.[7][8] While this provides a simple curve fitting procedure, the resulting algorithm may be biased by excessively weighting small data values, which can produce large errors in the profile estimate. One can partially compensate for this problem through weighted least squares estimation, reducing the weight of small data values, but this too can be biased by allowing the tail of the Gaussian to dominate the fit. In order to remove the bias, one can instead use an iteratively reweighted least squares procedure, in which the weights are updated at each iteration.[8] It is also possible to perform non-linear regression directly on the data, without involving the logarithmic data transformation; for more options, see probability distribution fitting.

Parameter precision

Once one has an algorithm for estimating the Gaussian function parameters, it is also important to know how precise those estimates are. Any least squares estimation algorithm can provide numerical estimates for the variance of each parameter (i.e., the variance of the estimated height, position, and width of the function). One can also use Cramér–Rao bound theory to obtain an analytical expression for the lower bound on the parameter variances, given certain assumptions about the data.[9][10]

  1. The noise in the measured profile is either i.i.d. Gaussian, or the noise is Poisson-distributed.
  2. The spacing between each sampling (i.e. the distance between pixels measuring the data) is uniform.
  3. The peak is "well-sampled", so that less than 10% of the area or volume under the peak (area if a 1D Gaussian, volume if a 2D Gaussian) lies outside the measurement region.
  4. The width of the peak is much larger than the distance between sample locations (i.e. the detector pixels must be at least 5 times smaller than the Gaussian FWHM).

Sind diese Annahmen erfüllt, gilt die folgende Kovarianzmatrix K für die 1D-Profilparameter , , und unter iid-Gaußschem Rauschen und unter Poisson-Rauschen: [ 9 ] Dabei ist die Breite der zur Abtastung der Funktion verwendeten Pixel, die Quanteneffizienz des Detektors und die Standardabweichung des Messrauschens. Die einzelnen Varianzen der Parameter sind somit im Fall von Gaußschem Rauschen:

und im Fall des Poisson-Rauschens,

Für die 2D-Profilparameter, die Amplitude , Position und Breite des Profils angeben, gelten die folgenden Kovarianzmatrizen: [ 10 ]

wobei die Varianzen der einzelnen Parameter durch die Diagonalelemente der Kovarianzmatrix gegeben sind.

Diskrete Gaußsche

Der diskrete Gaußsche Kern (durchgezogene Linie) im Vergleich zum abgetasteten Gaußschen Kern (gestrichelte Linie) für verschiedene Skalen

Man könnte nach einem diskreten Analogon zur Gauß-Funktion fragen; dies ist in diskreten Anwendungen, insbesondere in der digitalen Signalverarbeitung , notwendig . Eine einfache Lösung besteht darin, die kontinuierliche Gauß-Funktion abzutasten, wodurch sich der abgetastete Gauß-Kern ergibt . Diese diskrete Funktion besitzt jedoch nicht die diskreten Analoga der Eigenschaften der kontinuierlichen Funktion und kann, wie im Artikel „ Skalenraumimplementierung“ beschrieben, zu unerwünschten Effekten führen .

Ein alternativer Ansatz ist die Verwendung des diskreten Gaußschen Kernels : [ 11 ] wobei die modifizierten Bessel-Funktionen ganzzahliger Ordnung bezeichnen.

Dies ist das diskrete Analogon der kontinuierlichen Gauß-Verteilung, da sie die Lösung der diskreten Diffusionsgleichung (diskreter Raum, kontinuierliche Zeit) darstellt, genau wie die kontinuierliche Gauß-Verteilung die Lösung der kontinuierlichen Diffusionsgleichung ist. [ 11 ] [ 12 ]

Anwendungen

Gaußsche Funktionen treten in vielen Kontexten in den Naturwissenschaften , den Sozialwissenschaften , der Mathematik und den Ingenieurwissenschaften auf . Einige Beispiele hierfür sind:

See also

References

  1. ^Squires, G. L. (2001-08-30). Practical Physics (4 ed.). Cambridge University Press. doi:10.1017/cbo9781139164498. ISBN 978-0-521-77940-1.
  2. ^Folland, Gerald B.; Sitaram, Alladi (1997). "The uncertainty principle: A mathematical survey". The Journal of Fourier Analysis and Applications. 3 (3): 207–238. Bibcode:1997JFAA....3..207F. doi:10.1007/BF02649110. ISSN 1069-5869.
  3. ^Weisstein, Eric W. "Fourier Transform – Gaussian". MathWorld. Retrieved 19 December 2013.
  4. ^Nawri, Nikolai. "Berechnung von Kovarianzellipsen"(PDF). Archived from the original(PDF) on 2019-08-14. Retrieved 14 August 2019.
  5. ^Parent, A., M. Morin, and P. Lavigne. "Propagation of super-Gaussian field distributions". Optical and Quantum Electronics 24.9 (1992): S1071–S1079.
  6. ^"GLAD optical software commands manual, Entry on GAUSSIAN command"(PDF). Applied Optics Research. 2016-12-15.
  7. ^ Caruana, Richard A.; Searle, Roger B.; Heller, Thomas; Shupack, Saul I. (1986). „Schneller Algorithmus zur Spektrenauflösung“. Analytical Chemistry . 58 (6). American Chemical Society (ACS): 1162–1167 . doi : 10.1021/ac00297a041 . ISSN 0003-2700 . 
  8. ^ a b Hongwei Guo, „Ein einfacher Algorithmus zur Anpassung einer Gauß-Funktion“, IEEE Sign. Proc. Mag. 28(9): 134-137 (2011).
  9. ^ a b N. Hagen, M. Kupinski und EL Dereniak, „Gaussian profile estimation in one dimension“, Appl. Opt. 46:5374–5383 (2007)
  10. ^ a b N. Hagen und EL Dereniak, „Gaussian profile estimation in two dimensions“, Appl. Opt. 47:6842–6851 (2008)
  11. ^ a b Lindeberg, T., "Scale-space for discrete signals," PAMI(12), Nr. 3, März 1990, S. 234–254.
  12. ^ Campbell, J, 2007, Das SMM-Modell als Randwertproblem unter Verwendung der diskreten Diffusionsgleichung , Theor Popul Biol. 2007 Dez;72(4):539–46.
  13. ^ Haddad, RA und Akansu, AN, 1991, Eine Klasse schneller Gaußscher Binomialfilter für die Sprach- und Bildverarbeitung , IEEE Trans. on Signal Processing, 39-3: 723–727
  14. ^ Honarkhah, M und Caers, J, 2010, Stochastische Simulation von Mustern mittels distanzbasierter Mustermodellierung , Mathematical Geosciences, 42: 487–517

Weiterführende Literatur

  • Haberman, Richard (2013). „10.3.3 Inverse Fouriertransformation einer Gauß-Funktion“. Angewandte partielle Differentialgleichungen . Boston: PEARSON. ISBN 978-0-321-79705-6Die