Pulsars
0 %
Log inSign up

Contractions et suites itérées

Suites itérées et convergence géométrique

La suite des itérées

Étant donné ff et un point de départ u0u_0, la suite des itérées est définie par récurrence :

un+1=f(un)u_{n+1} = f(u_n)

Si cette suite converge vers un réel cc et si ff est continue, alors en passant à la limite dans la relation de récurrence on obtient f(c)=cf(c) = c : la limite est nécessairement un point fixe. Toute la question est donc de garantir la convergence.

Décroissance géométrique de l'erreur

Supposons ff kk-contractante et admettant un point fixe cc. En appliquant la définition au couple (un,c)\left(u_n, c\right), et puisque f(c)=cf(c) = c :

un+1c=f(un)f(c)kunc\lvert u_{n+1} - c \rvert = \lvert f(u_n) - f(c) \rvert \leq k\,\lvert u_n - c \rvert

Une récurrence immédiate donne alors la majoration fondamentale :

uncknu0c\lvert u_n - c \rvert \leq k^n\,\lvert u_0 - c \rvert

Comme 0k<10 \leq k < 1, on a kn0k^n \to 0, donc uncu_n \to c : la convergence est géométrique, et elle a lieu quel que soit le point de départ.

Un exemple à suivre pas à pas

Itérons f(x)=cosxf(x) = \cos x à partir de u0=0u_0 = 0. Le point fixe, appelé nombre de Dottie, vaut c0,739085c \approx 0{,}739085.

nn unu_n unc\lvert u_n - c \rvert
00 00 0,73910{,}7391
22 0,5403020{,}540302 0,19880{,}1988
44 0,6542900{,}654290 0,08480{,}0848
66 0,7013690{,}701369 0,03770{,}0377
88 0,7221020{,}722102 0,01700{,}0170
1010 0,7314040{,}731404 0,00770{,}0077

L'erreur est divisée par un peu plus de 22 toutes les deux étapes, ce qui correspond bien à k20,71k^2 \approx 0{,}71 par itération de rang pair :

012345678910erreur |u_n - c|rang n

Estimations d'erreur utilisables

En pratique on ne connaît pas cc, donc la majoration en u0c\lvert u_0 - c \rvert ne sert à rien telle quelle. Deux variantes se calculent, elles, à partir des seules valeurs observées :

unckn1ku1u0etunck1kunun1\lvert u_n - c \rvert \leq \frac{k^n}{1 - k}\,\lvert u_1 - u_0 \rvert \qquad \text{et} \qquad \lvert u_n - c \rvert \leq \frac{k}{1 - k}\,\lvert u_n - u_{n-1} \rvert

La première, dite a priori, permet de fixer à l'avance le nombre d'itérations nécessaires ; la seconde, dite a posteriori, transforme l'écart entre deux termes consécutifs en majoration de l'erreur réelle. C'est le critère d'arrêt utilisé dans les algorithmes.

Piège classique

Croire qu'une suite définie par un+1=f(un)u_{n+1} = f(u_n) converge dès que ff possède un point fixe. Pour f(x)=2xf(x) = -2x, le point fixe est 00, et pourtant la suite issue de u0=1u_0 = 1 vaut 1,2,4,8,1, -2, 4, -8, \ldots : elle diverge, car f=2>1\lvert f' \rvert = 2 > 1. Le point fixe est alors dit répulsif. Seule la contraction garantit qu'il est attractif.