A note on improvement by iteration for the approximate solutions of second kind Fredholm integral equations with Green's kernels

Titel: Eine Anmerkung zur Verbesserung durch Iteration für die Näherungslösungen von Fredholm-Integralgleichungen zweiter Art mit Greenschen Kernen

1. Problemstellung

Der Artikel befasst sich mit der numerischen Lösung linearer Operatorgleichungen der Form $x - Kx = f$, wobei $K$ ein Fredholm-Integraloperator mit einem Kern $\kappa(s, t)$ vom Typ der Greenschen Funktion ist. Solche Kerne sind typischerweise nicht glatt (nicht $C^\infty$), da sie eine Stetigkeit der Funktion selbst aufweisen, aber ihre Ableitungen an der Diagonalen $s=t$ einen Sprung aufweisen können.

Das Ziel ist die Approximation der eindeutigen Lösung $\phi$ durch diskrete Methoden. Der Fokus liegt auf der Analyse der Konvergenzordnung (Convergence Order) verschiedener numerischer Verfahren:

• Kollokationsmethode ($\phi^C_n$)
• Iterierte Kollokationsmethode ($\phi^S_n$)
• Modifizierte Kollokationsmethode ($\phi^M_n$)
• Iterierte modifizierte Kollokationsmethode ($\tilde{\phi}^M_n$)

Ein zentrales Problem in der Literatur ist, dass bei Verwendung von nicht-Gaußschen Kollokationspunkten (z. B. äquidistante Punkte) und stückweisen Polynomen geraden Grades die Iteration bei glatten Kernen oft keine Verbesserung der Konvergenzordnung bringt. Die Autoren untersuchen, ob dies auch für Greensche Kerne gilt und ob durch Modifikationen der Methode eine Verbesserung erzielt werden kann.

2. Methodik

Die Autoren verwenden einen funktionalanalytischen Ansatz basierend auf Interpolationsprojektionen:

• Raum und Partition: Der Definitionsbereich $[0, 1]$ wird in $n$ gleich große Intervalle unterteilt ($h = 1/n$). Der Approximationsraum $X_n$ besteht aus stückweisen Polynomen vom Grad $\le 2r$ bezüglich dieser Partition.
• Kollokationspunkte: In jedem Teilintervall werden $2r+1$ Kollokationspunkte gewählt. Diese sind nicht die Gauß-Punkte (Nullstellen der Legendre-Polynome), sondern äquidistante Punkte (oder allgemein $2r+1$ beliebige Punkte pro Intervall).
• Projektionsoperator: $P_n$ ist der interpolierende Projektionsoperator auf $X_n$ an den gewählten Kollokationspunkten.
• Kern-Eigenschaften: Der Kern $\kappa(s, t)$ ist vom Greenschen Typ, d.h., er ist auf den Dreiecken $\Omega_1 = \{0 \le t \le s \le 1\}$ und $\Omega_2 = \{0 \le s \le t \le 1\}$ durch glatte Funktionen $\kappa_1$ bzw. $\kappa_2$ definiert, die auf der Diagonalen stetig sind ($\kappa_1(s,s) = \kappa_2(s,s)$).
• Analyse-Werkzeuge:
  • Verwendung von dividierten Differenzen (divided differences) zur Abschätzung der Fehlerterme.
  • Herleitung von Lemmas, die die Beschränktheit von dividierten Differenzen des Operators $Kx$ an den Kollokationspunkten zeigen (Lemma 2.1 und 2.2).
  • Anwendung von Taylor-Entwicklungen und Mittelwertsätzen, um die Fehlerterme in Abhängigkeit von der Schrittweite $h$ zu schätzen.

3. Wichtige Beiträge und Ergebnisse

Die Hauptergebnisse werden in den Sätzen 4 und 5 zusammengefasst und stellen eine signifikante Verbesserung gegenüber früheren Ergebnissen für nicht-Gaußsche Punkte dar:

A. Kollokations- und Iterierte Kollokationsmethode

Für die Standard-Kollokationslösung $\phi^C_n$ und die iterierte Lösung $\phi^S_n = K\phi^C_n + f$ gilt unter der Annahme, dass die exakte Lösung $\phi \in C^{2r+2}([0,1])$ ist:

• Fehler der Kollokationslösung: $\|\phi - \phi^C_n\|_\infty = O(h^{2r+1})$
• Fehler der iterierten Lösung: $\|\phi - \phi^S_n\|_\infty = O(h^{2r+2})$

Bedeutung: Im Gegensatz zu Fällen mit glatten Kernen, wo die Iteration bei nicht-Gaußschen Punkten oft keine Verbesserung bringt, zeigt sich hier eine Verbesserung der Konverdenzordnung um einen Faktor $h$ (von $2r+1$ auf $2r+2$).

B. Modifizierte und Iterierte Modifizierte Kollokationsmethode

Für die modifizierte Kollokationsmethode (eingeführt von Kulkarni), definiert durch den Operator $K_n^M = P_n K + K P_n - P_n K P_n$, und die iterierte modifizierte Lösung $\tilde{\phi}^M_n$:

• Fehler der modifizierten Lösung:
  • Für $r \ge 1$: $\|\phi - \phi^M_n\|_\infty = O(h^{2r+2})$
  • Für $r = 0$: $\|\phi - \phi^M_n\|_\infty = O(h^3)$
• Fehler der iterierten modifizierten Lösung:
  • Für $r \ge 1$: $\|\phi - \tilde{\phi}^M_n\|_\infty = O(h^{2r+3})$
  • Für $r = 0$: $\|\phi - \tilde{\phi}^M_n\|_\infty = O(h^4)$

Bedeutung: Die modifizierte Methode liefert bereits eine höhere Konvergenzordnung als die Standard-Kollokation. Die Iteration verbessert die Ordnung weiter (um jeweils einen weiteren Grad $h$), was bei Greenschen Kernen und nicht-Gaußschen Punkten bisher nicht in dieser Form nachgewiesen war.

4. Signifikanz und Schlussfolgerung

• Überwindung einer Limitierung: Bisherige Arbeiten zeigten, dass bei nicht-Gaußschen Punkten die Iteration oft keine Superkonvergenz (Superconvergence) liefert. Dieser Artikel widerlegt dies für den wichtigen Fall von Greenschen Kernen.
• Theoretischer Durchbruch: Die Autoren beweisen, dass die spezifische Struktur der Greenschen Kerne (Stetigkeit, aber Unstetigkeit der Ableitungen auf der Diagonalen) es erlaubt, durch geschickte Wahl der Projektion und Iteration die Konvergenzordnung zu steigern, selbst wenn die Kollokationspunkte nicht optimal (Gauß) gewählt sind.
• Praktische Relevanz: Da die Wahl von äquidistanten oder einfachen Kollokationspunkten in der Praxis oft einfacher zu implementieren ist als die Berechnung von Gauß-Punkten, bietet diese Arbeit eine theoretische Grundlage für effizientere Algorithmen, die dennoch hohe Genauigkeit erreichen.
• Gegenbeispiel: Im Abschnitt 4 wird ein Beispiel angeführt, das zeigt, dass bei sehr einfachen Fällen (stückweise konstante Funktionen, $r=0$, bestimmte nicht-Gaußsche Punkte) die Superkonvergenz tatsächlich versagen kann. Dies unterstreicht die Notwendigkeit der hier entwickelten Analyse für den Fall $r \ge 0$ mit stückweisen Polynomen höheren Grades.

Zusammenfassend demonstriert der Artikel, dass für Fredholm-Integralgleichungen mit Greenschen Kernen die Iteration der Kollokations- und modifizierten Kollokationsmethoden zu einer messbaren Steigerung der Konvergenzordnung führt, selbst wenn nicht-Gaußsche Kollokationspunkte verwendet werden. Dies ist ein wichtiger Schritt zur Optimierung numerischer Verfahren für eine breite Klasse von Integralgleichungen.