Problème 7 : Suite Récurrente Non Linéaire
Énoncé du problème
Soit u0 = 2 et la suite définie par :
un 1 cos un
un+1 = + +
2 un n2
Étudier la convergence et trouver un développement asymptotique.
Solution complète
Étude de la suite sans perturbation
Considérons d’abord la suite sans le terme de perturbation :
un 1
un+1 = f (un ) = +
2 un
Point fixe. Le point fixe u∗ vérifie :
u∗ 1 u∗ 1 √
u∗ = + ∗ ⇒ = ∗ ⇒ (u∗ )2 = 2 ⇒ u∗ = 2
2 u 2 u
Stabilité. La dérivée de f est :
1 1
f ′ (u) = − 2
2 u
√
Au point fixe u∗ = 2:
√ 1 1
f ′ ( 2) = − = 0
2 2
La dérivée nulle indique une convergence très rapide (au moins quadratique) pour la suite non perturbée.
Introduction de la perturbation
√
Posons un = 2 + vn où vn est une petite perturbation.
√
Développement de Taylor. Développons chaque terme autour de 2 :
√ √
un 2 + vn 2 vn
= = +
2 2 2 2
vn2 vn3
1 1 1 1 1 vn 4
=√ =√ · v = √ 1− √ + − √ + O(vn )
un 2 + vn 2 1 + √n2 2 2 2 2 2
√ √ √
cos un cos( 2 + vn ) cos 2 cos vn − sin 2 sin vn
2
= 2
=
n n n2
√ 2
vn
√ 3
cos 2(1 − + O(vn4 )) − sin 2(vn − v6n + O(vn5 ))
2
=
n2
√ √ √ 3
cos 2 sin 2 cos 2 2 vn
= − v n − vn + O
n2 n2 2n2 n2
Linéarisation
En substituant dans la relation de récurrence :
√ !
√ vn2 vn3
2 vn 1 vn 4
2 + vn+1 = + + √ − + √ − + O(vn )
2 2 2 2 2 2 4
√ √ √ 3 !
cos 2 sin 2 cos 2 2 vn
+ 2
− 2
vn − 2
vn + O
n n 2n n2
√
2
√
Or + √1 = 2, donc :
2 2
√
√ √ v
n vn vn2 vn3 cos 2
2 + vn+1 = 2+ − + √ − +
2 2 2 2 4 n2
√ √ 3
sin 2 cos 2 2 4 vn
− 2
vn − 2
vn + O(vn ) + O
n 2n n2
Ainsi :
√ √ √
vn2 vn3 cos 2 sin 2
3
cos 2 2 4 vn
vn+1 = √ − + 2
− 2
vn − 2
vn + O(vn ) + O
2 2 4 n n 2n n2
Analyse de la convergence
L’équation dominante est :
1
vn+1 = O(vn2 ) +O
n2
1
Ordre de grandeur. Supposons que vn = O . Alors : n
1 1 1
vn+1 = O + O = O
n2 n2 n2
Ce qui est cohérent avec vn = O n1 .
En fait, on peut montrer par récurrence que |vn | ≤ Cn pour une certaine constante C > 0.
Développement asymptotique
Premier terme
a 1
Posons vn = n
+o n
. Alors :
√
a2
a 1 cos 2 1
+o = √ + + o
n+1 n 2 2n2 n2 n2
1
Développons n+1
:
1 1 1 1
= − 2 +O
n+1 n n n3
Donc :
√
a2
a a 1 cos 2 1
− 2 +o = √ + 2
+o
n n n 2 2n2 n n2
En identifiant les termes d’ordre n1 , on obtient a = 0.
Deuxième terme
1 b
+ o n12 :
Ainsi, vn = o n
. Essayons alors vn = n2
√
b 1 cos 2 1
2
+o 2
= 2
+o
(n + 1) n n n2
1 1 2 1
Or (n+1)2
= n2
− n3
+O n4
, donc :
√
b 2b 1 cos 2 1
2
− 3 +o = +o
n n n2 n2 n2
√
D’où b = cos 2.
Troisième terme
√
Pour le terme suivant, posons vn = cosn2 2 + nc3 + o n13 .
Alors :
√ √ √
cos 2 c 1 cos 2 2 cos 2 c 1
vn+1 = 2
+ 3
+o 3
= 2
− 3
+ 3 +o
(n + 1) (n + 1) n n n n n3
Et le membre droit :
√ √ √ √ √ √
vn2 cos 2 cos2 2 sin 2 cos 2
cos 2 sin 2 1 1
√ + − vn + o = + √ − + o
2 2 n2 n2 n3 n2 2 2n4 n4 n3
1
En comparant les termes d’ordre n3
:
√ √
−2 cos 2 + c = 0 ⇒ c = 2 cos 2
Conclusion et théorème
La suite (un ) définie par u0 = 2 et
un 1 cos un
un+1 = + +
2 un n2
√
converge vers 2 et admet le développement asymptotique :
√ √
√
cos 2 2 cos 2 1
un = 2 + 2
+ 3
+o
n n n3
Démonstration. La convergence est assurée par :
√ √
1. Le terme principal f (u) = a un point fixe attractif en 2 avec f ′ ( 2) = 0
u
2
+ 1
u
2. La perturbation cosn2un est de l’ordre de O n12
√
3. L’analyse de stabilité montre que vn = un − 2 satisfait vn+1 = O(vn2 ) + O n12
4. Le développement asymptotique est obtenu par la méthode des développements limités successifs
Illustrations
Illustration 1 : Comportement qualitatif
un
un
asymptotique √
2
Illustration 2 : Erreur et développement asymptotique
Erreur et développement asymptotique
0.6
vn théorique
√
cos 2
n2
0.5 O(1/n) (comparaison)
0.4
2
√
0.3
un −
0.2
0.1
0
2 4 6 8 10 12 14 16 18 20
n
Illustration 3 : Convergence rapide due à f’(2)=0
Convergence rapide (échelle logarithmique)
1024
Convergence quadratique typique
Notre suite O(1/n2 )
10−28
10−80
2|
10−132
√
|un −
10−184
10−236
10−288
10−340
0 1 2 3 4 5 6 7 8 9 10
n
Interprétation physique
La suite étudiée combine deux phénomènes :
u
1. Convergence rapide due au schéma itératif u 7→ 2
+ u1 qui est une méthode de type Newton pour
résoudre u2 = 2
cos un
2. Perturbation externe n2
qui ralentit la convergence mais reste contrôlable
Le développement asymptotique montre que l’erreur décroît comme O(1/n2 ), ce qui est relativement
rapide et démontre la robustesse du schéma itératif principal face aux petites perturbations.