Thêta algorithme

L'ε-algorithme et le ρ-algorithme, bien que issus de théories différentes et ayant des applications complémentaires, ont une structure algorithmique très proche. La seule différence entre les algorithmes est la valeur du numérateur dans la formule de récurrence (égale à 1 pour l'ε-algorithme et k+1 pour le ρ-algorithme). L’idée du θ-algorithme est de proposer un numérateur qui optimise l'accélération de la convergence, se calculant automatiquement avec l'évolution des termes de la suite à accélérer.

L'algorithme part de la structure générale:

${\displaystyle \theta _{-1}^{(n)}=0{\text{ }}\theta _{0}^{(n)}=S_{n}{\text{ }}\forall n}$

${\displaystyle \theta _{k+1}^{(n)}=\theta _{k-1}^{(n+1)}+{\frac {\omega _{k}^{(n)}}{\theta _{k}^{(n+1)}-\theta _{k}^{(n)}}}{\text{ }}\forall n,k}$

Ou S_n est la suite numérique dont on cherche à connaitre la limite S. Le numérateur ωⁿ_k est à déterminer pour optimiser au mieux la convergence. Brezinski a démontré que le numérateur optimal pour accélérer les colonnes paires était:

${\displaystyle \omega _{2k+1}^{(n)}=1{\text{ et }}\omega _{2k}^{(n)}=-\lim _{n\to \infty }{\frac {\theta _{2k}^{(n+2)}-\theta _{2k}^{(n+1)}}{D_{2k+1}^{(n+1)}-D_{2k+1}^{(n)}}}{\text{ }}\forall k}$

${\displaystyle {\text{ avec }}D_{k}^{(n)}={\frac {1}{\theta _{k}^{(n+1)}-\theta _{k}^{(n)}}}{\text{ }}\forall n,k}$

Ce choix nécessitant l'évaluation d'une limite dépendant de S_n, il est en pratique remplacé dans le θ-algorithme par sa meilleure approximation disponible, compte tenu des termes de S_n déjà calculés.

${\displaystyle \omega _{2k+1}^{(n)}=1{\text{ et }}\omega _{2k}^{(n)}=-{\frac {\theta _{2k}^{(n+2)}-\theta _{2k}^{(n+1)}}{D_{2k+1}^{(n+1)}-D_{2k+1}^{(n)}}}{\text{ }}\forall k,n}$

${\displaystyle {\text{ avec }}D_{k}^{(n)}={\frac {1}{\theta _{k}^{(n+1)}-\theta _{k}^{(n)}}}{\text{ }}\forall n,k}$

Finalement, le θ-algorithme consiste à calculer un tableau en initialisant la première colonne par la suite S_n, et en calculant les autres cases à l'aide des relations suivantes :

${\displaystyle \theta _{-1}^{(n)}=0{\text{ }}\theta _{0}^{(n)}=S_{n}{\text{ }}\forall n}$

${\displaystyle \theta _{2k+1}^{(n)}=\theta _{2k-1}^{(n+1)}+D_{2k}^{(n)}{\text{ }}\forall n,k}$

${\displaystyle \theta _{2k+2}^{(n)}={\frac {D_{2k+1}^{(n+1)}\theta _{2k}^{(n+1)}-D_{2k+1}^{(n)}\theta _{2k}^{(n+2)}}{D_{2k+1}^{(n+1)}-D_{2k+1}^{(n)}}}{\text{ }}\forall n,k}$

${\displaystyle {\text{ avec }}D_{k}^{(n)}={\frac {1}{\theta _{k}^{(n+1)}-\theta _{k}^{(n)}}}{\text{ }}\forall n,k}$

La formule de calcul du θ-algorithme permet de calculer les termes du tableau de haut en bas et de gauche à droite.

Les différentes colonnes d'indice k paires du tableau fournissent d'autres suites convergeant en général plus vite que la suite S_n d'origine. Les colonnes d'indice impaires sont considérés comme des intermédiaires de calcul. Il faut connaître 3 termes supplémentaires de la suite S_n pour passer d'une colonne paire du tableau à une autre. Ceci diffère de l'ε-algorithme et du ρ-algorithme, ou seulement 2 termes étaient nécessaires (dans le cas présent, un terme est utilisé pour estimer le facteur de relaxation).

• Mise en œuvre classique

Après chaque calcul d'un nouveau terme de la suite initiale, une diagonale montante du tableau du θ-algorithme est initiée. L'ajout de trois termes de la suite initiale permet de finaliser deux diagonales montantes du tableau. Pour faire des économies de mémoire, il est inutile de mémoriser l'intégralité du tableau dans ce cas d'utilisation. Seules les deux dernières diagonales ascendantes sont indispensables, ainsi que l'usage de trois variables auxiliaires[2].

• Itération du θ-2

Chaque colonne paire du θ-algorithme génère une sous-suite qui converge vers la limite recherchée. Une possibilité est de partir d'une de ces sous-suites, et de lui appliquer de nouveau le θ-algorithme, de la même manière que pour la suite initiale. Plusieurs procédés peuvent être envisagés en fonction de la sous-suite choisie et l'utilisation de l'algorithme qu'on en fait. Le plus utilisé est l'itération de la 2ème colonne du θ-algorithme, appelé algorithme θ-2 (par analogie au Delta-2 d'Aitken, étant l'itéré de la 2ème colonne de l'ε-algorithme), ou transformation W itéré de Lubkin[3] . La procédure consiste à arrêter le développement du θ-algorithme à 2ème colonne, et de repartir de la suite de cette colonne comme point de départ pour appliquer à nouveau le θ-2. L'algorithme θ-2 bénéficie de plusieurs résultats théorique de convergence[4] - [5], dont le θ-algorithme est dépourvu, d'où son intérêt.

• optimisation par colonne

Lorsque la suite à accélérer à notre disposition comprend un nombre restreint de termes, il est utile d'optimiser au mieux leur usage pour approcher leur limite. Le θ-algorithme de base utilise 3 termes pour progresser de 2 colonnes (un des termes est sacrifié pour fournir une estimation du paramètre de relaxation ωⁿ_k). Dans l'idéal, le paramètre de relaxation optimal pour la colonne k est la limite de ωⁿ_k quand n tend vers l'infini. Une possibilité de modification est de calculer colonne par colonne. La première colonne étant calculée jusqu'au dernier terme disponible, colonne des facteurs de relaxation ωⁿ₁ en est déduite. Ceci permet d'estimer la limite de la suite ωⁿ₁, en retenant son dernier terme (ou en utilisant un accélérateur de convergence à cette suite). Cette valeur limite estimée est utilisée pour calculer la deuxième colonne de l'algorithme. Le calcul de cette colonne nécessite un terme de moins que dans la mise en œuvre classique. Les colonnes suivantes sont calculées de la même manière.

• Accélération de la convergence d'une suite numérique

La série suivante, série dite du problème de Bâle,

${\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n^{2}}}={\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+{\frac {1}{4^{2}}}+\cdots ={\frac {\pi ^{2}}{6}}\approx 1,64493406684822643...}$

converge très lentement. Le θ-algorithme peut être utilisé pour accélérer la convergence des sommes partielles de cette série. Voici le résultat :

La suite à accélérer est la somme partielle de la série, en ajoutant un terme à la fois. On constate que les colonnes paires convergent vers la limite, et ceci plus rapidement que la suite d'origine. De même, chaque colonne paire converge plus rapidement que la colonne paire précédente : le mécanisme d'estimation du facteur de relaxation s'avère performant dans ce cas.

Les colonnes impaires ne représentant que des intermédiaires de calcul, elles ont été remplacées dans le tableau par les colonnes des coefficients de relaxation déduits. Ceci permet d'exposer les rouages de l'algorithme et d’observer vers quelles valeurs les coefficients de relaxation convergent. Les valeurs des coefficients de relaxation semblent converger vers 2/1, 4/3, 6/5... Ceci suggère que le ρ-algorithme avec la suite auxiliaire 1,2,3,4,5,6... (le ρ-algorithme "classique" ) serait l'algorithme de choix pour accélérer cette suite. L'utilisation du θ-algorithme permet d'éviter à rechercher le bon algorithme avec la bonne suite auxiliaire.

En comparant les résultats obtenus avec l'ε-algorithme et le ρ-algorithme, on obtient:

Les meilleurs résultats sont obtenus avec le ρ-algorithme, ce qui était prévisible d'après les facteurs de relaxation du θ-algorithme. Celui-ci paye sa polyvalence en étant légèrement moins performant que le ρ-algorithme. L'ε-algorithme est inefficace à accélérer cette suite.

Le problème de Bâle est un cas particulier de la fonction zêta de Riemann, correspondant à ζ(2). Pour les abscisses fractionnaires de cette fonction, le ρ-algorithme classique n'est plus efficace: il nécessite une suite auxiliaire spécifique, ou une variante[6]. Par exemple, pour ζ(1,5),

${\displaystyle \zeta (1,5)=\sum _{n=1}^{\infty }{\frac {1}{n^{1,5}}}={\frac {1}{1^{1,5}}}+{\frac {1}{2^{1,5}}}+{\frac {1}{3^{1,5}}}+{\frac {1}{4^{1,5}}}+\cdots \approx 2,612375348685...}$

En comparant les résultats des trois algorithmes, on obtient:

Le seul algorithme accélérant significativement la série est le θ-algorithme.

Pour la série de Leibniz :

${\displaystyle {\frac {\pi }{4}}=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{2n+1}}\approx 0,7853981633974...}$

Cette fois ci, l'ε-algorithme donne de bons résultats, de même que le θ-algorithme.

Ces exemples exposent la polyvalence du θ-algorithme par rapport aux deux algorithmes dont il est issu.

• Accélération de la convergence des séries de Fourier et de Tchebychev

Les sommes partielles des séries de Fourier ou de Tchebychev peuvent fournir une suite qui peut être accélérée par le θ-algorithme. Cependant, une manière plus efficace de procéder est de transformer, par un changement de variable, ces séries en une série entière. Les sommes partielles de ces séries entières seront les donnés d'entrée pour le θ-algorithme.

La série de Fourier :

${\displaystyle f(\theta )={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }(a_{n}\cos {n\theta }+b_{n}\sin {n\theta })}$

ou celle de Tchebychev:

${\displaystyle f(x)={\frac {c_{0}}{2}}+\sum _{n=1}^{\infty }c_{n}T_{n}(x)}$

se transforment chacune en une série entière (d'une variable complexe)

${\displaystyle F(z)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }(a_{n}-\mathrm {i} b_{n})z^{n}{\text{ avec }}z=\mathrm {e} ^{\mathrm {i} \theta }}$

ou

${\displaystyle F(z)={\frac {c_{0}}{2}}+\sum _{n=1}^{\infty }c_{n}z^{n}{\text{ avec }}z=\mathrm {e} ^{\mathrm {i} \theta };\theta =Arcos(x)}$

dont la partie réelle est la série d'origine. Les T_n(x) représentent les polynômes de Tchebychev de 1ère espèce.

Par exemple, la fonction Arc sinus sur [–1, +1], dont la série de Tchebychev vaut :

${\displaystyle {\frac {\pi }{4}}Arcsin(x)=\sum _{n=0}^{\infty }{\frac {T_{2n+1}(x)}{(2n+1)^{2}}}}$ devient après changement de variable la série entière :

${\displaystyle F(z)=\sum _{n=1}^{\infty }{\frac {z^{2n+1}}{(2n+1)^{2}}}{\text{ avec }}z=\mathrm {e} ^{\mathrm {i} \theta };\theta =Arcos(x)}$

dont on pourra accélérer la convergence en utilisant le θ-algorithme (avec des nombres complexes) sur les sommes partielles de cette série. La partie réelle du résultat sera la valeur accélérée de la série de Tchebychev d'origine.

Le graphique ci dessous compare la fonction Arc sinus avec sa série de Tchebychev en retenant les 7 premiers termes (série de Tchebychev de degré 13), ainsi que la valeur accélérée de cette série partielle par l'ε-algorithme et le θ-algorithme. Avec 7 termes, l'ε-algorithme a pu être poussé jusqu'à la 6ème colonne, et le θ-algorithme, jusqu'à la 4ème colonne.

Accélération de la convergence de série de Tchebychev de l'Arc sinus(x) par le thêta algorithme

La série de Tchebychev converge assez lentement (coefficients décroissant en 1/n²), du fait que la fonction Arc sinus n'est pas dérivable aux extrémités de l'intervalle. L'ε-algorithme accélère significativement la convergence, sauf dans les zones proches des extrémités. Le θ-algorithme donne de bons résultats, sur tout l'intervalle, incluant les extrémités.