Numerische Behandlung Volterrascher Integralgleichungen 2. Art

Gesucht ist die Lösung der Gleichung

(11.78)

für x aus dem Intervall . Numerische Lösungsansätze bestehen darin, das Integral durch eine Quadraturformel zu approximieren:

(11.79a)

Das Integrationsintervall und somit die Quadraturformel sind von x abhängig. Das wird durch den Index von zum Ausdruck gebracht. Man erhält als Näherungsausdruck für (11.78)

(11.79b)

Die Funktion ist eine Näherung für die Lösung von (11.78). Die Anzahl und Lage der Stützstellen der Quadraturformel ist von x abhängig, wodurch deren Wahl stark eingeschränkt ist. Ist eine Stützstelle von , so müssen und insbesondere bekannt sein. Dies erfordert aber zuvor eine Auswertung der rechten Seite von (11.79b) für , was einer Quadratur über dem Intervall entspricht. Aus diesem Grund ist die Verwendung der häufig bevorzugten GAUSSschen Quadraturformeln nicht möglich. Man löst das Problem durch die Wahl von Stützstellen und verwendet Quadraturformeln Q[a,xn] mit Stützstellen . Die Funktionswerte in den Stützstellen werden abkürzend bezeichnet durch Für erhält man (vgl. Integralgleichungen mit ausgearteten Kernen)

(11.79c)

und damit:

(11.79d)

Dabei hat Q[a,x1] die Stützstellen x0 und x1 und folglich die Gestalt

(11.79e)

mit geeigneten Koeffizienten w0 und . Setzt man dieses Verfahren fort, kann man die nacheinander aus der allgemeinen Beziehung

(11.79f)

bestimmen. Die Quadraturformeln Q[a,xk] haben folgende Form:

(11.79g)

Damit lautet (11.79f):

(11.79h)

Die einfachste Quadraturformel ist die linksseitige Rechteckformel. Dabei ist

(11.79i)

Man erhält damit das System

=  
= (11.80a)
=  


und allgemein
(11.80b)

Eine etwas genauere Approximation des Integrals gewährleistet die Trapezformel. Die Stützstellen seien zur Vereinfachung äquidistant,

(11.80c)

Angewandt auf (11.79f) ergibt das

(11.80d)
(11.80e)

Die jeweils zu berechnende Größe kommt dabei auch auf der rechten Seite vor. Die Gleichungen sind aber leicht nach den gesuchten Funktionswerten umzustellen.

Hinweis: Mit der angeführten Methode können auch nichtlineare Integralgleichungen näherungsweise gelöst werden. In diesem Fall wird bei Anwendung der Trapezformel zur Bestimmung der jedesmal die Lösung einer nichtlinearen Gleichung erforderlich sein. Dies kann man umgehen, wenn man die Trapezformel nur auf das Intervall [a,xk-1] anwendet und das Intervall [xk-1,xk] mit der linksseitigen Rechteckformel behandelt. Ist h genügend klein, wird dieser Quadraturfehler die Lösung nicht sehr beeinflussen.

Beispiel


Die Integralgleichung soll nach der Vorschrift (11.79f) mit der linksseitigen Rechteckformel näherungsweise gelöst werden. Als Stützstellen werden die äquidistanten Werte zugrunde gelegt, d.h. .






In der folgenden Tabelle sind zum Vergleich die Werte der exakten Lösung sowie der Näherungslösungen, die mittels linksseitiger Rechteckformel und Trapezformel ermittelt wurden, aufgeführt. Die Berechnung erfolgte mit der Schrittweite h=0,1.
    x          exakt          Rechteckformel     Trapezformel    
0,2 2,0401 2,0602 2,0401
0,4 2,1621 2,2030 2,1620
0,6 2,3709 2,4342 2,3706
0,8 2,6749 2,7629 2,6743
1,0 3,0862 3,2025 3,0852