Residuum (numerische Mathematik)

Vorlage:Hinweisbaustein Als Residuum bezeichnet man in der numerischen Mathematik die Abweichung vom gewünschten Ergebnis, welche entsteht, wenn in eine Gleichung Näherungslösungen eingesetzt werden. Angenommen, es sei eine Funktion <math>f</math> gegeben und man möchte ein <math>x</math> finden, so dass

Mit einer Näherung <math>\tilde{x}</math> an <math>x</math> ist das Residuum <math>r</math>

<math>r = f(x) - f(\tilde{x})= b - f(\tilde{x}).</math>

Der Fehler <math>e</math> zur Lösung hingegen ist

<math>e = x - \tilde{x}.</math>

Der Fehler ist in der Regel unbekannt, da <math>x</math> unbekannt ist, weswegen dieser als Abbruchkriterium in einem numerischen Verfahren nicht benutzbar ist. Das Residuum hingegen hängt nur von <math>\tilde{x}</math> ab.

Wenn das Residuum klein ist, folgt in vielen Fällen, dass der Fehler auch klein ist, also die Näherung nahe bei der Lösung liegt, das heißt der relative Fehler ist

<math>\frac{\|x - \tilde{x}\|}{\|x\|} \ll 1.~</math>

In diesen Fällen wird die zu lösende Gleichung als gut gestellt angesehen und das Residuum kann als Maß der Abweichung der Näherung von der exakten Lösung betrachtet werden. Bei linearen Gleichungssystemen, also <math>f(x)=Ax</math>, können sich die Norm des relativen Fehlers und die Norm des relativen Residuums um den Faktor der Kondition <math>\kappa</math> unterscheiden:

<math>\frac{\|x - \tilde{x}\|}{\|x\|} \le \kappa(A) \frac{\|Ax - A\tilde{x}\|}{\|Ax\|} = \kappa(A) \frac{\|b - A\tilde{x}\|}{\|b\|}. ~</math>

Residuum einer Operatorgleichung

Analog wird der Begriff des Residuums für Differential-, Integral- und Funktionalgleichungen verwendet, bei denen anstelle einer Zahl eine Funktion <math>u</math> gesucht ist, die eine Operatorgleichung

für alle Werte <math> x \in \mathcal X</math> erfüllt. Für eine Approximation <math>\tilde{u}</math> an <math>u</math> ist das Residuum die Funktion

<math>r = T(u) - T(\tilde{u}) = g - T(\tilde{u}).</math>

Als Maß für die Güte der Approximation kann dann zum Beispiel das Maximum des Residuums

<math>\max_{x\in \mathcal X} \|r(x)\| = \max_{x\in \mathcal X} \|g(x)-T(\tilde{u})(x)\|</math>

oder auch das gemittelte Residuum

<math>\int_{\mathcal X} \|g(x)-T(\tilde{u})(x)\|^2~{\rm d} x</math>

gewählt werden.

Literatur

C. T. Kelley: Iterative Methods for Linear and Nonlinear Equations. SIAM, ISBN 0-89871-352-8.
R. Schaback, H. Wendland: Numerische Mathematik. 5. Auflage, Springer, 2005.