Fehlerfunktion-Rechner

Was ist die Fehlerfunktion?

Die Fehlerfunktion, bezeichnet als erf(x), ist eine der wichtigsten speziellen Funktionen in der Mathematik und den angewandten Wissenschaften. Sie ist definiert als das Integral erf(x) = (2/√π) ∫₀ˣ e⁻t² dt. Der Faktor 2/√π normiert die Funktion so, dass erf(∞) = 1 gilt. Die Funktion liefert Werte zwischen −1 und +1, ist ungerade (erf(−x) = −erf(x)) und streng monoton steigend. Sie tritt natürlicherweise überall dort auf, wo Gauß-Integrale mit einer endlichen oberen Grenze vorkommen — etwa in der Wahrscheinlichkeitstheorie, der Wärmeleitung und der Diffusion.

Der Name „Fehlerfunktion“ stammt aus der Theorie der Beobachtungsfehler in Astronomie und Vermessungswesen, wo Gauß und Laplace sie verwendeten, um die Wahrscheinlichkeit zu beschreiben, dass ein Messfehler innerhalb eines bestimmten Bereichs liegt. Heute ist sie unverzichtbar in der Statistik (sie liefert die Verteilungsfunktion der Normalverteilung), der Physik (Diffusion, Wärmeleitungsgleichung) und der Ingenieurwissenschaft (Kommunikationstechnik, Signalerkennung).

Komplementäre Fehlerfunktion (erfc)

Die komplementäre Fehlerfunktion ist schlicht erfc(x) = 1 − erf(x). Sie stellt die Wahrscheinlichkeitsmasse im Randbereich (Tail) der Gauß-Verteilung jenseits von x dar. Während sich erf(x) auf die akkumulierte Wahrscheinlichkeit von 0 bis x konzentriert, betrachtet erfc(x) den verbleibenden Tail von x bis ∞. Besondere Werte sind erfc(0) = 1 und erfc(∞) = 0. In der Zuverlässigkeitstechnik taucht erfc oft direkt auf, weil man sich für Ausfallwahrscheinlichkeiten interessiert — die Wahrscheinlichkeit, dass eine Größe einen Schwellenwert überschreitet — statt für die akkumulierte Masse darunter.

Für große positive x fällt erfc(x) sehr schnell ab: erfc(1) ≈ 0,1573, erfc(2) ≈ 0,00468, erfc(3) ≈ 0,0000221. Dieser schnelle Abfall macht die direkte Berechnung für große x numerisch instabil, was die skalierte Version erfcx motiviert.

Wie man erf(x) berechnet — Abramowitz-&-Stegun-Näherung

Obwohl sich erf(x) nicht in geschlossener Form durch elementare Funktionen ausdrücken lässt, kann sie mit rationalen Polynomnäherungen sehr genau berechnet werden. Die am weitesten verbreitete Näherung stammt von Abramowitz und Stegun (Formel 7.1.26) und erreicht über den gesamten reellen Zahlenbereich einen maximalen absoluten Fehler von weniger als 1,5 × 10⁻⁷:

Man setzt t = 1/(1 + 0,3275911 |x|), berechnet dann das Polynom p(t) = t(0,254829592 + t(−0,284496736 + t(1,421413741 + t(−1,453152027 + t · 1,061405429)))) und schließlich erf(x) = sign(x) · (1 − p(t) · e⁻x²). Dieses fünfkoeffizientige Polynom in Horner-Form läuft in konstanter Zeit und ist die Standardwahl in den meisten wissenschaftlichen Rechenbibliotheken.

Häufige Werte von erf(x)

x	erf(x)	erfc(x)	erfcx(x)
0	0.000000	1.000000	1.000000
0.5	0.520500	0.479500	0.615080
1	0.842701	0.157299	0.427584
1.5	0.966105	0.033895	0.253994
2	0.995322	0.004678	0.134641
3	0.999978	0.000022	0.057289

Anwendungen der Fehlerfunktion

Normalverteilung und Wahrscheinlichkeit

Für eine standardnormalverteilte Variable Z ~ N(0,1) gilt P(−x ≤ Z ≤ x) = erf(x/√2). Äquivalent dazu ist die Verteilungsfunktion (CDF) der Normalverteilung Φ(x) = (1 + erf(x/√2))/2. Das bedeutet, dass sich jede Wahrscheinlichkeitsberechnung im Zusammenhang mit der Normalverteilung durch erf ausdrücken lässt. In Pythons SciPy ist scipy.special.erf direkt verfügbar; in C99 und C++11 ist std::erf eine Standardbibliotheksfunktion.

Wärmeleitung und Diffusion

Die Lösung der eindimensionalen Wärmeleitungsgleichung ∂T/∂t = α ∂²T/∂x² in einem halbunendlichen Medium mit sprunghafter Randbedingung lautet T(x,t) = T₀ erfc(x / (2√(αt))). Dabei ist α die Temperaturleitfähigkeit. Der Ausdruck x/(2√(αt)) ist die dimensionslose Ähnlichkeitsvariable für Diffusionsprobleme. Mit dieser Formel lässt sich berechnen, wie weit Wärme (oder ein diffundierender Stoff) zur Zeit t in ein Material eingedrungen ist — wesentlich in der Materialwissenschaft, Geophysik und Halbleiterfertigung.

Signalverarbeitung und Kommunikationstechnik

In der digitalen Kommunikationstechnik ist die Bitfehlerrate (BER) für binäre Phasenumtastung (BPSK), übertragen über einen additiven weißen gaußschen Rauschkanal (AWGN), gegeben durch BER = erfc(√(Eₛ/N₀))/2, wobei Eₛ/N₀ das Signal-Rausch-Verhältnis pro Bit ist. Die in der Kommunikationstechnik verwendete Q-Funktion lautet Q(x) = erfc(x/√2)/2. Systemdesigner nutzen erfc-Tabellen und -Rechner, um die minimale Sendeleistung zu bestimmen, die für eine gewünschte Fehlerrate benötigt wird — damit ist erfc grundlegend für das Design drahtloser Kommunikationssysteme.

erfc vs. erfcx — Wann welche Funktion verwenden

Für kleine bis moderate Werte von x (etwa |x| ≤ 3) ist die direkte Berechnung von erfc(x) = 1 − erf(x) numerisch stabil. Für große positive x wird erfc(x) jedoch extrem klein (z. B. erfc(5) ≈ 1,54 × 10⁻ῲῲ), sodass Auslöschungsfehler (Floating-Point-Kürzung) in 1 − erf(x) dominieren. Die skalierte Funktion erfcx(x) = e^(x²) · erfc(x) vermeidet dies: für x = 5 ist erfcx(5) ≈ 0,1107, eine perfekt darstellbare Zahl. Nachdem erfcx(x) ermittelt wurde, lässt sich erfc(x) = erfcx(x) · e^(−x²) bei Bedarf mit hochpräziser Arithmetik zurückgewinnen. Verwenden Sie erfcx, wenn Sie erfc für Argumente x > 3 in einem numerisch sensiblen Kontext berechnen.