Improvement of conformal maps combined with the Sinc approximation for derivatives over infinite intervals

Titre de l'étude

Amélioration des applications conformes combinées à l'approximation Sinc pour les dérivées sur des intervalles infinis.

1. Problématique

L'article s'attaque au problème de l'approximation numérique des dérivées de fonctions définies sur des intervalles infinis (notamment $(0, \infty)$ et $(-\infty, \infty)$).

• Contexte : F. Stenger a proposé des formules d'approximation efficaces combinant l'approximation Sinc et des applications conformes adaptées. Ces méthodes permettent d'atteindre une convergence de type "racine-exponentielle" ($O(\exp(-c\sqrt{n}))$).
• Limitation identifiée : Pour certains cas spécifiques (cas 2 et 4 de la classification de Stenger), les applications conformes originales ($\psi_2$ et $\psi_4$) limitent la largeur du domaine de régularité analytique ($d$) et le taux de décroissance ($\mu$) des fonctions transformées. Cela contraint la constante $c$ dans l'erreur d'approximation, limitant ainsi la vitesse de convergence.
• Objectif : Améliorer le taux de convergence en concevant de nouvelles applications conformes qui élargissent le domaine d'analyticité et optimisent la décroissance des fonctions transformées.

2. Méthodologie

Les auteurs proposent une modification des applications conformes utilisées dans la formule d'approximation de Stenger pour les dérivées d'ordre $m$.

• Approche de base : L'approximation de la dérivée $f^{(m)}(t)$ est obtenue en différenciant une série Sinc pondérée par une fonction de suppression de divergence $g_i(t)$ et une application conforme $\phi_i$.
  $$f^{(m)}(t) \approx \sum_{k=-M}^{N} \frac{f(\phi_i(kh))}{g_i(\phi_i(kh))} \left(\frac{d}{dt}\right)^m \left\{ g_i(t) S(k, h)(\phi_i^{-1}(t)) \right\}$$
• Nouvelles applications conformes : Au lieu d'utiliser les fonctions originales de Stenger, les auteurs introduisent :
  • Pour le cas $(0, \infty)$ avec décroissance exponentielle (remplacement de $\psi_2(x) = \text{arsinh}(e^x)$) :
    $$\phi_2(x) = \log(1 + e^x)$$
  • Pour le cas $(-\infty, \infty)$ avec décroissance algébrique à gauche et exponentielle à droite (remplacement de $\psi_4(x) = \sinh(\log(\text{arsinh}(e^x)))$) :
    $$\phi_4(x) = 2 \sinh(\log(\log(1 + e^x)))$$
• Analyse théorique : Les auteurs établissent de nouveaux théorèmes de convergence (Théorèmes 3.1 et 3.2). La preuve repose sur :
  1. La démonstration que les fonctions transformées appartiennent à l'espace de Hardy $H^1$ sur un domaine strip $D_d$ plus large.
  2. L'estimation de l'erreur de discrétisation (via le théorème de Okayama et Tanaka) et de l'erreur de troncature (via le lemme de Stenger).
  3. L'utilisation de la formule de Faà di Bruno pour gérer les dérivées d'ordre supérieur des compositions de fonctions.

3. Contributions Clés

• Amélioration des bornes de convergence : La contribution principale est l'augmentation des paramètres $d$ (largeur du domaine d'analyticité) et $\mu$ (taux de décroissance) dans l'expression de l'erreur.
  • Pour le cas 2, la borne supérieure de $d$ passe de $\pi/2$ (Stenger) à $\pi$ (Nouvelle méthode).
  • Pour le cas 4, non seulement $d$ augmente, mais $\mu$ peut également être plus grand car la condition de décroissance exponentielle est moins restrictive avec la nouvelle transformation.
• Nouveaux théorèmes de convergence : Les auteurs prouvent rigoureusement que la nouvelle formule conserve la convergence racine-exponentielle $O(n^{(m+1)/2} \exp(-\sqrt{\pi d \mu n}))$ mais avec une constante $c' = \sqrt{\pi d \mu}$ potentiellement supérieure à celle de Stenger.
• Preuves techniques : Développement de lemmes spécifiques pour borner les dérivées des nouvelles applications conformes inverses et leurs compositions, garantissant la validité des estimations d'erreur.

4. Résultats

Des expériences numériques ont été menées en double précision (C) et avec une précision arbitraire (Mathematica) pour valider la théorie.

• Cas test 1 ($t \in (0, \infty)$) : Fonction $f(t) = \sqrt{\frac{t}{1+t}} e^{-t}(1-e^{-t})^2$.
  • Comparaison des erreurs maximales pour $f(t)$, $f'(t)$ et $f''(t)$.
  • Résultat : La formule améliorée converge plus rapidement que celle de Stenger. Les graphiques montrent une pente plus raide pour l'erreur en fonction du nombre de points d'évaluation.
• Cas test 2 ($t \in (-\infty, \infty)$) : Fonction $f(t) = \frac{1}{(4+t^2)(1+e^{\pi t/2})}$.
  • Comparaison similaire pour les dérivées jusqu'à l'ordre 2.
  • Résultat : Là encore, la nouvelle formule surpasse celle de Stenger. De plus, le profil de convergence de la méthode de Stenger apparaît plus "irrégulier" (bumpy) que celui de la méthode améliorée, bien que la raison mécanique exacte de cette irrégularité ne soit pas élucidée.
• Observation générale : Dans tous les cas, l'erreur absolue diminue plus vite avec la nouvelle méthode, confirmant l'augmentation effective de la constante de convergence.

5. Signification et Perspectives

• Impact : Cette étude offre des outils numériques plus précis et efficaces pour la différentiation numérique sur des domaines infinis, ce qui est crucial pour la résolution numérique d'équations différentielles (EDP) et d'intégrales.
• Optimisation des paramètres : En augmentant le paramètre $d$ (la largeur du domaine de régularité), la méthode permet d'atteindre une précision donnée avec moins de points d'évaluation ($N$), réduisant ainsi le coût computationnel.
• Perspectives futures : Les auteurs notent que les transformations utilisées sont de type "exponentielle simple". Ils suggèrent que l'utilisation de transformations "double-exponentielles" (DE), déjà efficaces pour les intégrales et les approximations sur intervalles finis, pourrait permettre d'améliorer encore davantage ces formules pour les dérivées sur intervalles infinis.

En résumé, ce papier démontre théoriquement et numériquement qu'une modification astucieuse des applications conformes dans le cadre de l'approximation Sinc permet d'optimiser significativement la vitesse de convergence pour le calcul de dérivées sur des intervalles infinis.