ChronoMath, une chronologie des MATHÉMATIQUES
à l'usage des professeurs de mathématiques, des étudiants et des élèves des lycées & collèges

Arrondis des ordinateurs, erreurs de mesure    arrondi d'une somme , d'un produit , d'une racine carrée        » La notion élémentaire d'arrondi et de troncature | numération binaire, codage et arrondis | exemple d'accumulation d'erreurs d'arrondi

On aborde ici l'aspect mathématique des erreurs d'arrondi engendrées par les ordinateurs (ou les calculatrices), ce qui revient également à étudier les conséquences des erreurs de mesures en physique ou en astronomie.

Soit x un nombre entré ou calculé. Suivant sa valeur, l'ordinateur va l'arrondir par excès ou par défaut. La valeur machine de x sera alors :

Plus précisément :

Il suit que les erreurs absolues et relatives commises sur x sont respectivement :

 ➔  Nous allons étudier trois cas démontrant la difficulté d'apprécier les erreurs d'arrondi d'un calculateur électronique.

Soit z = f(x,y) le résultat d'une opération élémentaire : un produit, une somme, etc. En notant f '_x(x,y) et f '_y(x,y) les dérivées partielles de f par rapport à x et y, on a :

soit, sensiblement, en passant aux valeurs absolues :

A chaque mise en mémoire d'une donnée ou d'un résultat intermédiaire z, l'ordinateur (resp. l'observateur dans le cas d'une mesure physique) peut commettre une erreur e_s dépendant de sa capacité de stockage (resp. de la précision de ses calculs) induisant une erreur complémentaire ze_s : par exemple, le produit de deux entiers n'engendrant pas d'erreur de stockage peut en produire une sur le résultat. L'erreur absolue commise sur z vérifiera alors :

et l'erreur relative sera :

On a ici f '_x(x,y) = y et f '_y(x,y) = x. On obtient alors :

On voit là que l'erreur est la somme des erreurs relatives, ce qui est parfaitement naturel et acceptable. Notons que si x et y sont des données entières ne dépassant pas la capacité de traitement de la machine, on a e_x = e_y = 0 et la seule erreur possible est l'erreur e_s de stockage du produit. Cette erreur peut se produire si x et y sont grands : perte de chiffres significatifs par arrondi.

Exemple :

100009 * 100009 = 10001800081  : le calcul avec certaines calculatrices fournit 1.000180008E+10, soit 10001800080 : l'erreur absolue sur ce produit semble être 1 et l'erreur relative de 10^-10. Pour le calcul de 100009 * 100009 - 10001800081, la machine répond 1 : il faut toujours distinguer l'affichage et le résultat effectivement en mémoire. Sur ces machines, les ennuis commencent véritablement dès que le résultat possède plus de 10 ou 12 chiffres significatifs.

On a ici f '_x(x,y) = f '_y(x,y) = 1. On obtient :

L'erreur dépend là, non seulement des valeurs de x et y, mais aussi du résultat théorique de l'opération. Si x + y est petit, deux nombres quasiment opposés par exemple, ou bien, dans le cas x - y, si x est sensiblement égal à y, l'erreur peut être considérable.

 !   On voit donc avec quelle méfiance il faut regarder un résultat affiché par l'ordinateur ou un résultat d'opérations sur des mesures entachées d'erreurs.

 ➔   Donnons l'exemple d'une "bête" résolution de l'équation du second degré : x² - 10x - 1/4 = 0. Les solutions sont :

et considérons le petit programme JavaScript suivant calculant les solutions x1 et x2 de l'équation donnée au moyen des formules usuelles rappelées ci-dessus.

On les affiche à l'écran. On calcule le produit des racines (en principe -0,25) et on vérifie si les solutions calculées par l'ordinateur annulent effectivement le trinôme.

La seconde méthode évite l'addition de nombres "presque opposés" par la formule bien connue x'.x" = c/a et on affiche les mêmes tests que dans la première méthode.

Avec x et y positifs : on a ici f(x,y) = √(xy) , donc  f '_x(x,y) = y/2√(xy) et f '_y(x,y) = x/2√(xy)

Dans un problème d'arithmétique, si x et y sont donnés entiers, et ne dépassant pas la capacité de l'ordinateur, leur produit peut provoquer un dépassement de capacité.

Supposons cependant ici x < 100 et y < 100 : les erreurs commises sur x, y et xy sont alors nulles. S'il s'agit de vérifier que xy est un carré parfait, en admettant que l'algorithme de recherche de la racine carrée de la machine soit efficace (i.e. dont la précision est supérieure à la capacité de stockage), alors on peut être assuré que xy est un carré parfait si l'on a (en JavaScript) :

En effet, si xy est entier, sa racine carrée sera calculée comme étant de la forme :

de sorte que la valeur machine sera dans tous les cas l'entier n (la barre symbolisant la limite de la capacité de stockage).