Méthode de Ruffini-Horner - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Introduction

Connue sous le nom de méthode de Horner, règle de Ruffini ou algorithme de Ruffini-Horner, cette méthode se décline sur plusieurs niveaux. Elle permet de calculer la valeur d'un polynôme en \scriptstyle x_0 . Elle présente un algorithme simple effectuant la division euclidienne d'un polynôme par \scriptstyle X -  x_0 . Mais elle offre aussi une méthode de changement de variable X=\scriptstyle x_0 + Y dans un polynôme. C'est sous cette forme qu'elle est utilisée pour déterminer une valeur approchée d'une racine d'un polynôme.

Histoire

La méthode de Ruffini-Horner de recherche d'une valeur approchée de racine d'un polynôme est publiée à quelques années d'intervalle par Paolo Ruffini (1804-1807-1813) et par William George Horner (1819-1845 posthume) mais il semble bien que Horner n'ait pas eu connaissance des travaux de Ruffini. La méthode de Horner est ensuite popularisée par les mathématiciens De Morgan et J.R. Young. Dans leurs premières publications, ces deux auteurs utilisent des méthodes de dérivations pour effectuer le changement de variable X = x0 + Y. Par la suite, ils présentent des versions ne faisant appel qu'à des techniques algébriques. La méthode de Ruffini-Horner est difficilement exploitable si le polynôme possède deux racines trop proches. Ruffini n'évoque pas ce problème mais Horner propose une procédure spéciale pour ce cas-là.

En tant que technique de changement de variable, on retrouve des algorithmes analogues, en Chine, pour l'extraction de racine n-ième, dans les Neuf Chapitres (263 après J.C)et dans l'œuvre de Al Samaw'al (XIIe siècle). Mais il semble bien que Sharaf al-Dīn al-Tūsī (XIIe siècle) soit le premier à l'utiliser dans le cas général d'une équation de degré 3.

Quotient d'un polynôme par X - x0

Cette même méthode permet aussi d'obtenir la division d'un polynôme par \scriptstyle X-x_0 . Soit \scriptstyle P = a_nX^n + a_{n-1}X^{n-1} + ... + a_0 .

La division euclidienne de P par \scriptstyle X-x_0 donne

P = (X-x_0)Q + P(x_0)\,

où Q est un polynôme de degré n - 1.

Si on écrit \scriptstyle Q = q_{n-1}X^{n - 1} + q_{n-2}X^{n-2}+ ...+ q_1X + q_0 et si on identifie les coefficients de même degré dans les deux membres, on obtient :

q_{n-1} = a_n \,
q_{k-1} - q_kx_0=a_k \, pour tout k tel que 0 < k < n

Soit encore

q_{k-1}=q_kx_0+a_k\, pour tout k tel que 0 < k < n


Les n valeurs de la suite q calculées ici sont précisément les n valeurs successives calculées dans le paragraphe précédent pour évaluer \scriptstyle P(x_0) . La mémorisation de ces valeurs successives donne donc les coefficients du polynôme quotient, la dernière valeur étant celle du reste.

Application pratique : Division euclidienne de \scriptstyle 4X^3 - 7X^2 + 3X - 5 par \scriptstyle X-2

Il suffit de reprendre le tableau précédemment construit et de lire dans les cases de la seconde ligne les coefficients de Q.

Coefficients de P 4 − 7 3 − 5
Coeficients de Q 4 8 − 7 = 1 2 + 3 = 5 Reste = 10 − 5 = 5

Donc

4X^3 - 7X^2 + 3X - 5 = (X-2)(4X^2 + X + 5) + 5\,

Valeur d'un polynôme en un point

Soit \scriptstyle P =a_nX^n + a_{n-1}X^{n-1} + ... + a_0 un polynôme et \scriptstyle x_0 un nombre. Le calcul de \scriptstyle P(x_0) = a_nx_0^n + a_{n-1}x_0^{n-1} + ... + a_0 laisse à penser qu'il faut calculer chacune des puissances de \scriptstyle x_0 , multiplier celle-ci par son coefficient ak puis faire la somme de ce que l'on a trouvé.

Si on calcule \scriptstyle x_0^n en multipliant successivement \scriptstyle x_0 par lui-même, le nombre de produits nécessaire est alors de n + (n − 1) + ... + 2 + 1 = n(n + 1) / 2, quantité qui croît comme le carré du degré du polynôme.

On peut améliorer la vitesse du calcul de \scriptstyle x_0^n par une méthode d'exponentiation rapide, permettant de réduire le temps du calcul de \scriptstyle P(x_0) à une quantité qui croît comme nln(n).

La méthode de Horner consiste à améliorer encore ce résultat en effectuant le calcul comme suit :

P(x_0) = ((...((a_nx_0 + a_{n-1})x_0 + a_{n-2})x_0 + ... ) + a_1)x_0 + a_0\,

Le nombre de produits est alors réduit à n, de sorte que le temps de calcul d'une fonction polynomiale en un point a est seulement proportionnel au degré du polynôme. La méthode consiste donc à multiplier le premier coefficient par \scriptstyle x_0 et à lui ajouter le second coefficient. On multiplie alors le nombre obtenu par \scriptstyle x_0 et on lui ajoute le troisième coefficient, etc. Elle s'organise très bien à l'aide d'un tableau dans lequel chaque case de la seconde ligne est obtenue en multipliant le coefficient de la case de gauche par \scriptstyle x_0 et en lui ajoutant le coefficient de la case du dessus.

Coefficients de P an an - 1 an - 2 ... a1 a0
Facteur x0 an anx0 + an - 1 (anx0 + an - 1)x0 + an-2 ... q0 P(x0)=q0x0 + a0

Exemple pratique : Calcul de \scriptstyle 4X^3 - 7X^2 + 3X - 5 pour \scriptstyle X=2

Coefficients de P 4 − 7 3 − 5
Facteur 2 4 8 − 7 = 1 2 + 3 = 5 P(2) = 10 − 5 = 5

Cette méthode permet aussi d'effectuer une conversion rapide d'un nombre écrit en base \scriptstyle x_0 en écriture en base 10. En effet, si un nombre s'écrit, en base \scriptstyle x_0 , \scriptstyle \overline{a_na_{n-1}\cdots a_0} , ce nombre vaut \scriptstyle a_nx_0^n+a_{n-1}x_0^{n-1}+\cdots+ a_0 .

Exemple pratique : écriture en base 10 du nombre hexadécimal DA78

Coefficients 13 10 7 8
Facteur 16 13 13 × 16 + 10 = 218 218 ×16+ 7 = 3495 DA78 = 3495 × 16 + 8 = 55928
Page générée en 0.185 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales | Partenaire: HD-Numérique
Version anglaise | Version allemande | Version espagnole | Version portugaise