Vous avez dit bizarre … comme c'est bizarre …

Suites numériques instables

Parmi les surprises ou bizarreries que l'on peut observer lors de l'utilisation d'une calculatrice ou d'un tableur, il y celles très connues des suites instables. Il s'agit de suites dont le calcul formel démontre qu'elles sont convergentes, mais dont les valeurs affichées par un tableur "dérapent" complètement à partir d'un certain rang, alors qu'avant ce rang fatal, tout avait l'air de se passer conformément au résultat théorique. Les exemples sont légions (dans une prochaine édition, je proposerai des exemples donnant lieu à des activités tout à fait consistantes du niveau de terminale S).

Voici un exemple très simple mais qui est vraiment très surprenant, évoqué par un collègue, Émile Sinturel, et que je vous propose d'analyser et de généraliser :

 

Soit la suite définie par :

On constate qu'on est bien dans un cas d'instabilité puisqu'au lieu de rester constante, la suite affichée a tout l'air de tendre vers -∞.

On a l'habitude de voir ce phénomène lorsque les valeurs en jeu ne sont pas des nombres décimaux, ce qui provoque une très faible erreur d'arrondi, laquelle est amplifiée par la réitération due à la récurrence, et finit par donner des valeurs qui n'ont plus rien à voir. Mais ce qui est très troublant ici, c'est que 0,3 est un bête nombre décimal à 1 chiffre après la virgule. En principe, pas de quoi affoler un simple tableur … C'est ce qu'on croit … et ce n’est pas ce qu'on voit !

 

Pour mener des investigations, généralisons un peu. Si la suite précédente est constante, c'est parce que u(n+1)=11un-3 =11un-10u0.
Du coup u1=11u0-10u0=uo. Une démonstration par récurrence permet d'en déduire le résultat.

 

Cas général : Soit la suite définie par

Tous les termes sont égaux à a.

 

Si a=3, on a vu que la suite "dérape" (pour n=14 , l'ordre de grandeur de l'erreur est de 0,01 ). Qu'en est-il pour d'autres valeurs de a ?

 

Pour prospecter, c'est là que GeoGebra va prendre toute son ampleur :

  • Créer un curseur : variable de 0 à 1 , pas de 0,05 .
  • Dans le tableur, générer les rangs dans la ligne 1, A2: , B2: = 11A2-10$A2 à recopier vers la droite.
  • Créer le nuage de points (Sélection de la plage, clic droit, Créer/Liste de points).

Il ne reste plus qu'à faire varier le curseur et observer…

On constate qu'il apparait 3 sortes de suites, celles qui semblent tendre vers +∞, celles qui semblent tendre vers -∞ , et celle qui, comme cela devrait être toujours le cas, sont constantes. Donc elles ne dérapent pas toutes ! Mais quelles sont les valeurs de qui sont bien sages, et quelles sont les rebelles ? Le résultat est spectaculaire : sur [0 ; 1], parmi les valeurs de que l'on peut tester avec le curseur, ce sont les valeurs 0, 0,25 ,  0,5,  0,75 et 1 qui renvoient bien une suite constante. Pour les autres valeurs, ça dérape, et systématiquement au rang 14.

 

Mais alors, qu'ont de particulier ces 5 nombres ? Il s'écrivent respectivement 20, 2-2, 2-1, 3x2-2 et 21 . Bon sang, mais c'est bien sûr un ordinateur, ça représente les nombres en binaire ! Comment écrit-on ces nombres en binaire ? Révision … 

En décimal

En binaire
 0
0,25 = 0x2 + 0x2-1 + 1x2-2 0,01
0,5 = 0x2 + 1x2-1 0,1
0,75 = 0x2 + 1x2-1 + 1x2-2 0,11
1 1

 

Et comment s'écrivent les autres, les rebelles ? Quelques exemples :

En décimal En binaire
 0,3 = 2-2 + 2-5 + 2-6 + 2-9  + ...  0,010011001 ... 
 0,4 = 2-2 + 2-3 + 2-6 + 2-7  + ...  0,0110011 ...
 0,7 = 2-1 + 2-3 + 2-4 + 2-7  + ...

 0,1011001 ... 

 

Les valeurs du terme initial qui font déraper la suite sont celles dont l'écriture à virgule en binaire est illimitée ! Les nombres non "binaux" en quelque sorte ! Les valeurs de ces nombres sont donc inévitablement arrondies. Pour autant, cela n'explique toujours pas le dérapage. En effet, arrondi ou pas, quel que soit a, 11a-10a=a et donc on devrait bien avoir u1=u0 et les suivants aussi …

Mais souvenons-nous : un outil de calcul dispose d'un nombre de chiffres fixé pour écrire les nombres. En conséquence, la précision de l'arrondi de a et de 10a n'est pas la même (pour a<1, la partie entière de 10a consomme 1 chiffre alors que celle de a n'en consomme pas). D'où le décalage fatal !

 

Reprenons l'exemple initial pour aller plus loin en faire une activité pour les élèves :

Pour le tableur 3 ≈ 10 × 0,3. Pour comprendre pourquoi la suite diverge, étudions la suite 

Il s'agit d'une suite arithmético-géométrique, objet d'exercices en lycée. Posons vn=un-0,3. On démontre que (vn) est une suite géométrique de 1er terme v0=ε et de raison q=11.

 

On en déduit : ∀n∈N, v= v× q= ε × 11n,

d'où u= v+ 0,3 = ε × 11+ 0,3.

Selon le signe de ε, la limite de (un) est +∞ ou -∞. Les valeurs obtenues sur tableur montrent que c'est -∞, c’est-à-dire que ε<0.

Enfin, 1113<1014<1114. Or nous savons que les calculs sur tableur se font avec une erreur relative de 10-14. C'est ce que l'on observe avec le tableur

 

Rémy Coste

page suivante

Retour sommaire
Retour sommaire

Les chantiers de pédagogie mathématique n°160 mars 2014

La Régionale Île-de-France APMEP, 26 rue Duméril, 75013 PARIS

Nantes 2017

Les inscriptions sont possibles jusqu’au dimanche 15 octobre (tarif préférentiel jusqu'au samedi 16 septembre)

 

Le groupe M.:A.T.H.

Le groupe M.:A.T.H. (Mathématiques : Approche par des Textes Historiques) poursuit ses séances de lecture de textes historiques.

lire la suite

Nouvelle annonce

Des postes pour la formation des professeurs des écoles en mathématiques à l'ESPE de l'académie de Versailles sont à pourvoir au 1er septembre 2017.

Voir l'annonce

Notre bulletin Les Chantiers

Le N°173 des Chantiers pédagogique de mathématique est en ligne

Lire le N°173

juin 2017

Concours 2017

"Maths et langages"

 

Voir les résultats du concours

 

Réunion du comité

Prochaine réunion 

septembre 2017

Le coin des annonces

Vous pouvez déposer vos annonces sur le site.

 

Repères de l'IREM

La revue Repères IREM, l'une des trois revues du réseau des IREM, est disponible en consultation à l'IREM de Paris

Algorithmique au bac S 2012

Dans le n°154 des Chantiers de Pédagogie Mathématique de Septembre 2012 , un article de Dominique Baroux et de Cécile Prouteau à propos des exercices d’algorithmique du bac S 2012

 

Dernière mise à jour

 

samedi 24 juin 2017