Fehlerfunktion

Aus besserwiki.de
Fehlerfunktion
Plot of the error function
Darstellung der Fehlerfunktion
Allgemeine Informationen
Allgemeine Definition
AnwendungsgebieteWahrscheinlichkeitsrechnung, Thermodynamik
Bereich, Kodomäne und Bild
Bereich
Bild
Grundlegende Merkmale
ParitätUngerade
Spezifische Merkmale
Wurzel0
Abgeleitet
Antiderivative
Definition von Reihen
Taylor-Reihe

In der Mathematik ist die Fehlerfunktion (auch Gauß-Fehlerfunktion genannt), oft mit erf bezeichnet, eine komplexe Funktion einer komplexen Variablen, die wie folgt definiert ist:

Dieses Integral ist eine spezielle (nicht elementare) Sigmoidfunktion, die häufig in der Wahrscheinlichkeitsrechnung, der Statistik und bei partiellen Differentialgleichungen auftritt. In vielen dieser Anwendungen ist das Argument der Funktion eine reelle Zahl. Wenn das Funktionsargument reell ist, dann ist auch der Funktionswert reell.

In der Statistik hat die Fehlerfunktion für nicht-negative Werte von x die folgende Bedeutung: Für eine normalverteilte Zufallsvariable Y mit Mittelwert 0 und Standardabweichung 1/2 ist erf x die Wahrscheinlichkeit, dass Y in den Bereich [-x, x] fällt.

Zwei eng verwandte Funktionen sind die komplementäre Fehlerfunktion (erfc), definiert als

und die imaginäre Fehlerfunktion (erfi), definiert als

wobei i die imaginäre Einheit ist

Graph der Fehlerfunktion

Bezeichnung

Die Bezeichnung "Fehlerfunktion" und ihre Abkürzung erf wurden 1871 von J. W. L. Glaisher vorgeschlagen, weil sie mit der "Theorie der Wahrscheinlichkeit und insbesondere der Theorie der Fehler" zusammenhängen. Das Komplement der Fehlerfunktion wurde ebenfalls von Glaisher in einer separaten Veröffentlichung im selben Jahr diskutiert. Für das "Gesetz der Leichtigkeit" von Fehlern, deren Dichte gegeben ist durch

(die Normalverteilung) gegeben ist, berechnet Glaisher die Wahrscheinlichkeit eines Fehlers, der zwischen p und q liegt, als:

Anwendungen

Wenn die Ergebnisse einer Reihe von Messungen durch eine Normalverteilung mit der Standardabweichung σ und dem Erwartungswert 0 beschrieben werden, dann ist erf (a/σ 2) die Wahrscheinlichkeit, dass der Fehler einer einzelnen Messung zwischen -a und +a liegt, für ein positives a. Dies ist z. B. bei der Bestimmung der Bitfehlerrate eines digitalen Kommunikationssystems nützlich.

Die Fehler- und komplementären Fehlerfunktionen treten beispielsweise in Lösungen der Wärmegleichung auf, wenn die Randbedingungen durch die Heaviside-Stufenfunktion gegeben sind.

Die Fehlerfunktion und ihre Näherungen können verwendet werden, um Ergebnisse abzuschätzen, die mit hoher Wahrscheinlichkeit oder mit geringer Wahrscheinlichkeit gelten. Gegeben eine Zufallsvariable X ~ Norm[μ,σ] (eine Normalverteilung mit Mittelwert μ und Standardabweichung σ) und eine Konstante L < μ:

wobei A und B bestimmte numerische Konstanten sind. Wenn L ausreichend weit vom Mittelwert entfernt ist, nämlich μ - Lσln k, dann:

Die Wahrscheinlichkeit geht also gegen 0, wenn k → ∞.

Die Wahrscheinlichkeit, dass X im Intervall [La, Lb] liegt, lässt sich ableiten als

Die Fehlerfunktion hat eine gewisse Ähnlichkeit mit der Verteilungsfunktion der Normalverteilung. Sie hat jedoch eine Zielmenge von , während eine Verteilungsfunktion zwingend Werte aus dem Bereich annehmen muss.

Die Fehlerfunktion kann verwendet werden, um mit Hilfe der Inversionsmethode normalverteilte Pseudozufallszahlen zu generieren.

Eigenschaften

Darstellungen in der komplexen Ebene
Integrand exp(-z2)
erf z

Die Eigenschaft erf (-z) = -erf z bedeutet, dass die Fehlerfunktion eine ungerade Funktion ist. Dies ergibt sich unmittelbar aus der Tatsache, dass der Integrand e-t2 eine gerade Funktion ist (die Antiderivative einer geraden Funktion, die im Ursprung Null ist, ist eine ungerade Funktion und umgekehrt).

Da die Fehlerfunktion eine ganze Funktion ist, die reelle Zahlen in reelle Zahlen überführt, gilt für jede komplexe Zahl z:

wobei z die komplex Konjugierte von z ist.

Die Integranden f = exp(-z2) und f = erf z sind in der komplexen z-Ebene in den Abbildungen rechts mit Bereichsfärbung dargestellt.

Die Fehlerfunktion bei +∞ ist genau 1 (siehe Gaußsches Integral). Auf der reellen Achse nähert sich erf z bei z → +∞ der Eins und bei z → -∞ der Eins. Auf der imaginären Achse tendiert sie zu ±i∞.

Taylor-Reihe

Die Fehlerfunktion ist eine ganze Funktion; sie hat keine Singularitäten (mit Ausnahme derjenigen im Unendlichen) und ihre Taylorentwicklung konvergiert immer, ist aber bekannt "[...] für ihre schlechte Konvergenz, wenn x > 1 ist."

Das definierende Integral kann nicht in geschlossener Form in Form von Elementarfunktionen ausgewertet werden, aber wenn man den Integranden e-z2 in seine Maclaurin-Reihe erweitert und Term für Term integriert, erhält man die Maclaurin-Reihe der Fehlerfunktion als:

die für jede komplexe Zahl z gilt. Die Nennerterme sind die Folge A007680 im OEIS.

Für die iterative Berechnung der obigen Reihe kann die folgende alternative Formulierung nützlich sein:

weil -(2k - 1)z2/k(2k + 1) den Multiplikator ausdrückt, mit dem der k-te Term in den (k + 1)-ten Term umgewandelt wird (wobei z als erster Term betrachtet wird).

Die imaginäre Fehlerfunktion hat eine sehr ähnliche Maclaurin-Reihe, die lautet:

die für jede komplexe Zahl z gilt.

Ableitung und Integral

Die Ableitung der Fehlerfunktion ergibt sich unmittelbar aus ihrer Definition:

Daraus ergibt sich auch unmittelbar die Ableitung der imaginären Fehlerfunktion:

Eine Gegenableitung der Fehlerfunktion, die man durch Integration durch Teile erhält, ist

Eine Antiderivative der imaginären Fehlerfunktion, die ebenfalls durch Integration durch Teile erhalten werden kann, ist

Ableitungen höherer Ordnung sind gegeben durch

wobei H die Hermite-Polynome der Physiker sind.

Bürmann-Reihe

Mit Hilfe des Satzes von Hans Heinrich Bürmann erhält man eine Entwicklung, die für alle reellen Werte von x schneller konvergiert als eine Taylor-Entwicklung:

wobei sgn die Vorzeichenfunktion ist. Behält man nur die ersten beiden Koeffizienten bei und wählt c1 = 31/200 und c2 = -341/8000, so zeigt die resultierende Näherung ihren größten relativen Fehler bei x = ±1,3796, wo er weniger als 0,0036127 beträgt:

Inverse Funktionen

Inverse Fehlerfunktion

Für eine komplexe Zahl z gibt es keine eindeutige komplexe Zahl w, die erf w = z erfüllt, so dass eine echte Umkehrfunktion mehrwertig wäre. Für -1 < x < 1 gibt es jedoch eine eindeutige reelle Zahl mit der Bezeichnung erf-1 x, die folgende Bedingung erfüllt

Die inverse Fehlerfunktion wird gewöhnlich mit dem Bereich (-1,1) definiert und ist in vielen Computeralgebrasystemen auf diesen Bereich beschränkt. Sie kann jedoch auf die Scheibe |z| < 1 der komplexen Ebene erweitert werden, indem man die Maclaurin-Reihe

mit c0 = 1 und

Wir haben also die Reihenentwicklung (gemeinsame Faktoren wurden aus Zähler und Nenner gestrichen):

(Nach der Streichung sind die Zähler-/Nennerbrüche die Einträge OEIS: A092676/OEIS: A092677 im OEIS; ohne Streichung sind die Zählerterme im Eintrag OEIS: A002067 angegeben). Der Wert der Fehlerfunktion bei ±∞ ist gleich ±1.

Für |z| < 1 gilt erf(erf-1 z) = z.

Die inverse komplementäre Fehlerfunktion ist definiert als

Für reelle x gibt es eine einzige reelle Zahl erfi-1 x, die erfi(erfi-1 x) = x erfüllt. Die inverse imaginäre Fehlerfunktion ist definiert als erfi-1 x.

Für jedes reelle x kann die Newton-Methode zur Berechnung von erfi-1 x verwendet werden, und für -1 ≤ x ≤ 1 konvergiert die folgende Maclaurin-Reihe:

wobei ck wie oben definiert ist.

Asymptotische Entwicklung

Eine nützliche asymptotische Entwicklung der komplementären Fehlerfunktion (und damit auch der Fehlerfunktion) für große reelle x ist

wobei (2n - 1)!! die doppelte Fakultät von (2n - 1) ist, die das Produkt aller ungeraden Zahlen bis zu (2n - 1) ist. Diese Reihe divergiert für jedes endliche x, und ihre Bedeutung als asymptotische Expansion ist, dass für jede ganze Zahl N ≥ 1 gilt

wobei der Rest in Landau-Schreibweise wie folgt lautet

wenn x → ∞.

In der Tat ist der genaue Wert des Rests

was sich leicht durch Induktion, Schreiben

und Integrieren nach Teilen.

Für ausreichend große Werte von x sind nur die ersten paar Terme dieser asymptotischen Entwicklung erforderlich, um eine gute Annäherung an erfc x zu erhalten (während für nicht allzu große Werte von x die obige Taylorentwicklung bei 0 eine sehr schnelle Konvergenz liefert).

Fortgesetzte Bruchentwicklung

Eine kontinuierliche Bruchentwicklung der komplementären Fehlerfunktion ist:

Integral der Fehlerfunktion mit Gaußscher Dichtefunktion

die mit Ng und Geller, Formel 13 in Abschnitt 4.3 verwandt erscheint, wobei die Variablen geändert werden.

Faktorielle Reihen

Die inverse faktorielle Reihe:

konvergiert für Re(z2) > 0. Hier bezeichnet

zn die steigende Fakultät und s(n,k) eine vorzeichenbehaftete Stirlingzahl erster Art. Es gibt auch eine Darstellung durch eine unendliche Summe, die die doppelte Fakultät enthält:

Numerische Approximationen

Approximation mit elementaren Funktionen

  • Abramowitz und Stegun geben mehrere Näherungen mit unterschiedlicher Genauigkeit an (Gleichungen 7.1.25-28). Dies ermöglicht es, die schnellste Annäherung zu wählen, die für eine bestimmte Anwendung geeignet ist. In der Reihenfolge der zunehmenden Genauigkeit sind dies:
    (maximaler Fehler: 5×10-4) wobei a1 = 0,278393, a2 = 0,230389, a3 = 0,000972, a4 = 0,078108
    (maximaler Fehler: 2,5×10-5) wobei p = 0,47047, a1 = 0,3480242, a2 = -0,0958798, a3 = 0,7478556
    (maximaler Fehler: 3×10-7) wobei a1 = 0,0705230784, a2 = 0,0422820123, a3 = 0,0092705272, a4 = 0,0001520143, a5 = 0,0002765672, a6 = 0,0000430638
    (maximaler Fehler: 1,5×10-7) mit p = 0,3275911, a1 = 0,254829592, a2 = -0,284496736, a3 = 1,421413741, a4 = -1,453152027, a5 = 1,061405429 Alle diese Näherungen gelten für x ≥ 0. Um diese Näherungen für negative x zu verwenden, nutzt man die Tatsache, dass erf x eine ungerade Funktion ist, also erf x = -erf(-x).
  • Exponentialschranken und eine reine Exponentialapproximation für die komplementäre Fehlerfunktion sind gegeben durch
  • Die obigen Ausführungen wurden auf Summen von N Exponentialen mit zunehmender Genauigkeit in Bezug auf N verallgemeinert, so dass erfc x durch 2(2x) genau angenähert oder begrenzt werden kann, wobei
    Insbesondere gibt es eine systematische Methodik zur Lösung der numerischen Koeffizienten {(an,bn)}N
    n = 1, die eine Minimax-Approximation oder -Schranke für die eng verwandte Q-Funktion ergeben: Q(x) ≈ (x), Q(x) ≤ (x), oder Q(x) ≥ (x) für x ≥ 0. Die Koeffizienten {(an,bn)}N
    n = 1 für viele Variationen der exponentiellen Näherungen und Grenzen bis zu N = 25 wurden als umfassender Datensatz für den freien Zugang freigegeben.
  • Eine enge Annäherung der komplementären Fehlerfunktion für x ∈ [0,∞) wird von Karagiannidis & Lioumpas (2007) gegeben, die für die geeignete Wahl der Parameter {A,B} zeigten, dass
    Sie ermittelten {A,B} = {1.98,1.135}, was eine gute Näherung für alle x ≥ 0 ergibt. Es stehen auch alternative Koeffizienten zur Verfügung, um die Genauigkeit für eine bestimmte Anwendung anzupassen oder den Ausdruck in eine enge Schranke zu verwandeln.
  • Eine untere Schranke für einen einzelnen Term ist
    wobei der Parameter β so gewählt werden kann, dass der Fehler im gewünschten Approximationsintervall minimiert wird.
  • Eine weitere Annäherung wird von Sergei Winitzki mit Hilfe seiner "globalen Padé-Annäherungen" gegeben:
    wobei
    Diese Näherung ist so konzipiert, dass sie in der Nähe von 0 und in der Nähe von unendlich sehr genau ist, und der relative Fehler ist kleiner als 0,00035 für alle reellen x. Bei Verwendung des alternativen Wertes a ≈ 0,147 verringert sich der maximale relative Fehler auf etwa 0,00013. Diese Näherung kann invertiert werden, um eine Näherung für die inverse Fehlerfunktion zu erhalten:
  • Eine Näherung mit einem maximalen Fehler von 1,2×10-7 für jedes reelle Argument ist:
    mit
    und

Tabelle der Werte

x erf x 1 - erf x
0 0 1
0.02 0.022564575 0.977435425
0.04 0.045111106 0.954888894
0.06 0.067621594 0.932378406
0.08 0.090078126 0.909921874
0.1 0.112462916 0.887537084
0.2 0.222702589 0.777297411
0.3 0.328626759 0.671373241
0.4 0.428392355 0.571607645
0.5 0.520499878 0.479500122
0.6 0.603856091 0.396143909
0.7 0.677801194 0.322198806
0.8 0.742100965 0.257899035
0.9 0.796908212 0.203091788
1 0.842700793 0.157299207
1.1 0.880205070 0.119794930
1.2 0.910313978 0.089686022
1.3 0.934007945 0.065992055
1.4 0.952285120 0.047714880
1.5 0.966105146 0.033894854
1.6 0.976348383 0.023651617
1.7 0.983790459 0.016209541
1.8 0.989090502 0.010909498
1.9 0.992790429 0.007209571
2 0.995322265 0.004677735
2.1 0.997020533 0.002979467
2.2 0.998137154 0.001862846
2.3 0.998856823 0.001143177
2.4 0.999311486 0.000688514
2.5 0.999593048 0.000406952
3 0.999977910 0.000022090
3.5 0.999999257 0.000000743

Verwandte Funktionen

Komplementäre Fehlerfunktion

Die komplementäre Fehlerfunktion, bezeichnet als erfc, ist definiert als

Damit ist auch erfcx definiert, die skalierte komplementäre Fehlerfunktion (die anstelle von erfc verwendet werden kann, um einen arithmetischen Unterlauf zu vermeiden). Eine andere Form von erfc x für x ≥ 0 ist nach ihrem Entdecker als Craigsche Formel bekannt:

Dieser Ausdruck gilt nur für positive Werte von x, aber er kann in Verbindung mit erfc x = 2 - erfc(-x) verwendet werden, um erfc(x) für negative Werte zu erhalten. Diese Form hat den Vorteil, dass der Bereich der Integration fest und endlich ist. Eine Erweiterung dieses Ausdrucks für erfc der Summe zweier nichtnegativer Variablen lautet wie folgt:

Imaginäre Fehlerfunktion

Die imaginäre Fehlerfunktion, bezeichnet als erfi, ist definiert als

wobei D(x) die Dawson-Funktion ist (die anstelle von erfi verwendet werden kann, um einen arithmetischen Überlauf zu vermeiden).

Trotz des Namens "imaginäre Fehlerfunktion" ist erfi x reell, wenn x reell ist.

Wenn die Fehlerfunktion für beliebige komplexe Argumente z ausgewertet wird, wird die resultierende komplexe Fehlerfunktion gewöhnlich in skalierter Form als Faddeeva-Funktion bezeichnet:

Kumulative Verteilungsfunktion

Die Fehlerfunktion ist im Wesentlichen identisch mit der kumulativen Standardnormalverteilungsfunktion Φ, die in einigen Softwaresprachen auch als norm(x) bezeichnet wird, da sie sich nur durch Skalierung und Übersetzung unterscheiden. In der Tat,

oder umgeordnet für erf und erfc:

Folglich ist die Fehlerfunktion auch eng mit der Q-Funktion verwandt, die die Schwanzwahrscheinlichkeit der Standardnormalverteilung darstellt. Die Q-Funktion lässt sich in Bezug auf die Fehlerfunktion wie folgt ausdrücken

Die Umkehrung von Φ ist als Normalquantilfunktion oder Probitfunktion bekannt und kann in Form der umgekehrten Fehlerfunktion wie folgt ausgedrückt werden

Die Standard-Normalquantilfunktion wird häufiger in der Wahrscheinlichkeitsrechnung und Statistik verwendet, während die Fehlerfunktion häufiger in anderen Bereichen der Mathematik verwendet wird.

Die Fehlerfunktion ist ein Spezialfall der Mittag-Leffler-Funktion und kann auch als konfluente hypergeometrische Funktion (Kummersche Funktion) ausgedrückt werden:

Sie hat einen einfachen Ausdruck in Form des Fresnel-Integrals.

In Bezug auf die regularisierte Gammafunktion P und die unvollständige Gammafunktion,

sgn x ist die Vorzeichenfunktion.

Verallgemeinerte Fehlerfunktionen

Grafik der verallgemeinerten Fehlerfunktionen En(x):
Graue Kurve: E1(x) = 1 - e-x/π
rote Kurve: E2(x) = erf(x)
grüne Kurve: E3(x)
blaue Kurve: E4(x)
goldene Kurve: E5(x).

Einige Autoren erörtern die allgemeineren Funktionen:

Bemerkenswerte Fälle sind:

  • E0(x) ist eine gerade Linie durch den Ursprung: E0(x) = x/eπ
  • E2(x) ist die Fehlerfunktion, erf x.

Nach der Division durch n! sehen alle En für ungerade n ähnlich (aber nicht identisch) aus. In ähnlicher Weise sehen die En für gerade n nach einer einfachen Division durch n! ähnlich (aber nicht identisch) aus. Alle verallgemeinerten Fehlerfunktionen für n > 0 sehen auf der positiven x-Seite des Diagramms ähnlich aus.

Diese verallgemeinerten Funktionen können für x > 0 äquivalent durch die Gammafunktion und die unvollständige Gammafunktion ausgedrückt werden:

Daher können wir die Fehlerfunktion mit Hilfe der unvollständigen Gammafunktion definieren:

Iterierte Integrale der komplementären Fehlerfunktion

Die iterierten Integrale der komplementären Fehlerfunktion sind definiert durch

Die allgemeine Rekursionsformel lautet

Sie haben die Potenzreihe

aus denen die Symmetrieeigenschaften folgen

und

Implementierungen

Als reelle Funktion eines reellen Arguments

  • In Posix-kompatiblen Betriebssystemen muss der Header math.h die Funktionen erf und erfc (doppelte Genauigkeit) sowie ihre Gegenstücke mit einfacher und erweiterter Genauigkeit erff, erfl und erfcf, erfcl deklarieren und die mathematische Bibliothek libm muss sie bereitstellen.
  • Die GNU Scientific Library bietet erf, erfc, log(erf) und skalierte Fehlerfunktionen.

Als komplexe Funktion eines komplexen Arguments

  • libcerf, numerische C-Bibliothek für komplexe Fehlerfunktionen, bietet die komplexen Funktionen cerf, cerfc, cerfcx und die reellen Funktionen erfi, erfcx mit etwa 13-14 Stellen Genauigkeit, basierend auf der Faddeeva-Funktion, wie sie im MIT Faddeeva-Paket implementiert ist

Verwendung

Wärmeleitungsgleichung

Die Fehlerfunktion und die komplementäre Fehlerfunktion kommen beispielsweise in Lösungen der Wärmeleitungsgleichung vor, wenn Randwertbedingungen durch die Heaviside-Funktion vorgegebenen sind.

Wertetabelle

0,00 0,0000000 1,0000000 1,30 0,9340079 0,0659921
0,05 0,0563720 0,9436280 1,40 0,9522851 0,0477149
0,10 0,1124629 0,8875371 1,50 0,9661051 0,0338949
0,15 0,1679960 0,8320040 1,60 0,9763484 0,0236516
0,20 0,2227026 0,7772974 1,70 0,9837905 0,0162095
0,25 0,2763264 0,7236736 1,80 0,9890905 0,0109095
0,30 0,3286268 0,6713732 1,90 0,9927904 0,0072096
0,35 0,3793821 0,6206179 2,00 0,9953223 0,0046777
0,40 0,4283924 0,5716076 2,10 0,9970205 0,0029795
0,45 0,4754817 0,5245183 2,20 0,9981372 0,0018628
0,50 0,5204999 0,4795001 2,30 0,9988568 0,0011432
0,55 0,5633234 0,4366766 2,40 0,9993115 0,0006885
0,60 0,6038561 0,3961439 2,50 0,9995930 0,0004070
0,65 0,6420293 0,3579707 2,60 0,9997640 0,0002360
0,70 0,6778012 0,3221988 2,70 0,9998657 0,0001343
0,75 0,7111556 0,2888444 2,80 0,9999250 0,0000750
0,80 0,7421010 0,2578990 2,90 0,9999589 0,0000411
0,85 0,7706681 0,2293319 3,00 0,9999779 0,0000221
0,90 0,7969082 0,2030918 3,10 0,9999884 0,0000116
0,95 0,8208908 0,1791092 3,20 0,9999940 0,0000060
1,00 0,8427008 0,1572992 3,30 0,9999969 0,0000031
1,10 0,8802051 0,1197949 3,40 0,9999985 0,0000015
1,20 0,9103140 0,0896860 3,50 0,9999993 0,0000007