Relation between two Sinc-collocation methods for Volterra integral equations of the second kind and further improvement

Technische Zusammenfassung: Beziehung zwischen zwei Sinc-Kollokationsmethoden für Volterra-Integralgleichungen und weitere Verbesserungen

1. Problemstellung
Das Paper befasst sich mit der numerischen Lösung von Volterra-Integralgleichungen zweiter Art der Form:
$$u(t) - \int_a^t k(t, s)u(s) \, ds = g(t), \quad a \le t \le b.$$
Hierbei sind $k(t, s)$ (der Kern) und $g(t)$ gegebene stetige Funktionen, und $u(t)$ ist die zu bestimmende Lösung.
Es existieren bereits mehrere auf Sinc-Numerik basierende Methoden für diese Gleichungen:

• Die Sinc-Nyström-Methode (Muhammad et al.), die eine root-exponentielle Konvergenz erreicht.
• Die Sinc-Kollokationsmethode von Rashidinia und Zarebnia.
• Die Methode von Stenger, die ursprünglich für Anfangswertprobleme (Kern als Funktion einer Variablen) entwickelt wurde, aber auch auf den allgemeinen Fall (Kern als Funktion zweier Variablen) anwendbar scheint.

Bisher war die theoretische Beziehung zwischen Stengers Methode und der Methode von Rashidinia–Zarebnia unklar, und die Konvergenz der letzteren war nicht streng bewiesen. Zudem gibt es einen Trade-off zwischen Konvergenzgeschwindigkeit und Rechenaufwand bei den bestehenden Methoden.

2. Methodik
Die Studie nutzt Sinc-Numerik, die auf der Approximation von Funktionen durch Sinc-Funktionen und der Sinc-unbestimmten Integration basiert. Ein zentrales Element ist die Transformation des Integrationsintervalls $(a, b)$ auf die reelle Achse $\mathbb{R}$ mittels Variablentransformationen:

• SE-Transformation (tanh): $t = \psi_{SE}(x) = \frac{b-a}{2} \tanh(\frac{x}{2}) + \frac{b+a}{2}$.
• DE-Transformation (double-exponential): $t = \psi_{DE}(x) = \frac{b-a}{2} \tanh(\frac{\pi}{2} \sinh x) + \frac{b+a}{2}$.

Die Arbeit vergleicht drei Hauptansätze:

1. Sinc-Nyström: Approximation des Integraloperators direkt.
2. Sinc-Kollokation (Rashidinia–Zarebnia): Standard-Kollokation mit speziellen Randbedingungen.
3. Stengers Methode: Eine hybride Vorgehensweise, die zuerst eine Nyström-Lösung berechnet und diese dann durch eine verallgemeinerte Sinc-Approximation (Interpolation) darstellt.

3. Hauptbeiträge und theoretische Ergebnisse

• Beziehung zwischen Stenger und Rashidinia–Zarebnia:

  • Die Studie beweist theoretisch, dass Stengers Methode als eine weitere Sinc-Kollokationsmethode betrachtet werden kann.
  • Es wird gezeigt, dass die Lösungen beider Methoden ($v_N^{Stenger}$ und $w_N^{RZ}$) an den Kollokationspunkten übereinstimmen ($v_N^{Stenger}(t_i) = w_N^{RZ}(t_i)$), aber im Allgemeinen nicht äquivalent sind ($v_N^{Stenger} \neq w_N^{RZ}$), da sie unterschiedliche Approximationsoperatoren verwenden.
  • Stengers Methode wird theoretisch für den allgemeinen Fall (Kern $k(t,s)$ von zwei Variablen) gerechtfertigt, obwohl sie ursprünglich für einen Kern einer Variablen entwickelt wurde.

• Konvergenzanalyse:

  • Es wird streng bewiesen, dass beide Kollokationsmethoden (Stenger und Rashidinia–Zarebnia) die gleiche Konvergenzrate erreichen:
    $$O(\sqrt{N} \exp(-\sqrt{\pi d \alpha N}))$$
    wobei $d$ die Größe des analytischen Bereichs und $\alpha$ der Hölder-Exponent ist.
  • Obwohl die Konvergenzrate identisch ist, ist Stengers Methode in der Implementierung einfacher und robuster, insbesondere bei Systemen von Integralgleichungen.

• Verbesserung durch DE-Transformation:

  • Die Autoren entwickeln eine neue Sinc-Kollokationsmethode, die Stengers Ansatz mit der DE-Transformation (double-exponential) kombiniert.
  • Diese neue Methode erreicht eine signifikant höhere Konvergenzrate (fast exponentiell):
    $$O(\exp(-\pi d N / \log(2dN/\alpha)))$$
  • Diese Rate ist zwar theoretisch geringfügig niedriger als die der DE-Sinc-Nyström-Methode, aber in der Praxis überlegen.

4. Numerische Ergebnisse
Numerische Experimente mit zwei Testgleichungen bestätigen die theoretischen Vorhersagen:

• Genauigkeit vs. $N$: Die DE-Sinc-Nyström-Methode zeigt die beste Konvergenz in Bezug auf die Anzahl der Punkte $N$, gefolgt von der DE-Sinc-Kollokation. Die SE-basierten Methoden sind etwas langsamer.
• Genauigkeit vs. Rechenzeit: Dies ist der entscheidende Punkt. Die DE-Sinc-Kollokationsmethode (neu vorgeschlagen) ist der DE-Sinc-Nyström-Methode in Bezug auf die Rechenzeit überlegen.
  • Grund: Die Nyström-Methode erfordert die Berechnung der Sinus-Integral-Funktion ($Si(x)$) als Basisfunktion, was einen hohen Rechenaufwand verursacht. Die Kollokationsmethode vermeidet dies und nutzt nur Standard-Sinc-Funktionen.
• Vergleich der Kollokationsmethoden: Stengers Methode und die von Rashidinia–Zarebnia liefern fast identische Ergebnisse, wobei Stengers Ansatz aufgrund der einfacheren Implementierung bevorzugt wird.

5. Bedeutung und Fazit

• Theoretische Klärung: Das Paper schließt die Lücke im Verständnis der Beziehung zwischen Stengers und Rashidinia–Zarebnias Methoden und liefert den ersten strengen Konvergenzbeweis für beide.
• Praktische Empfehlung: Die Studie empfiehlt Stengers Sinc-Kollokationsmethode als die beste Wahl unter den drei verglichenen klassischen Methoden (besser als Nyström wegen geringerer Rechenkosten trotz ähnlicher Konvergenzrate).
• Innovation: Die vorgeschlagene DE-Sinc-Kollokationsmethode stellt einen signifikanten Fortschritt dar. Sie kombiniert die Einfachheit und Effizienz der Kollokation mit der extrem hohen Konvergenzgeschwindigkeit der DE-Transformation. Sie bietet das beste Verhältnis zwischen Genauigkeit und Rechenaufwand für Volterra-Integralgleichungen zweiter Art.

Zusammenfassend liefert das Paper einen umfassenden theoretischen Rahmen für Sinc-basierte Kollokationsmethoden und identifiziert eine optimierte Variante (DE-Kollokation nach Stenger), die sowohl theoretisch als auch numerisch überlegen ist.