Interpolation linéaire

Un article de Wikipédia, l'encyclopédie libre.
Aller à : navigation, rechercher
Les points rouges correspondent aux points (xk,yk), et la courbe bleue représente la fonction d'interpolation, composée de segments de droite.

L’interpolation linéaire est la méthode la plus simple pour estimer la valeur prise par une fonction continue entre deux points déterminés (interpolation). Elle consiste à utiliser pour cela la fonction affine (de la forme f(x) = m.x +b) passant par les deux points déterminés. Cette technique était d'un emploi systématique lorsque l'on ne disposait que de tables numériques pour le calcul avec les fonctions transcendantes : les tables comportaient d'ailleurs à cet effet en marge les « différences tabulaires », auxiliaire de calcul servant à l'interpolation linéaire.

Enfin l'interpolation linéaire est la base de la technique de quadrature numérique par la méthode des trapèzes.

Principe[modifier | modifier le code]

Supposons que l'on connaisse les valeurs prises par une fonction f en deux points x_a et x_b :

  • f(x_a) =y_a et
  • f(x_b) =y_b.

La méthode consiste à approcher la fonction f par la fonction affine \bar{f} telle que \bar{f}(x_a) =y_a et \bar{f}(x_b) =y_b ; cette fonction a pour équation (trois formulations équivalentes) :

 \bar{f}(x) = \frac{y_a-y_b}{x_a-x_b} x + \frac{x_a \cdot y_b-x_b\cdot y_a}{x_a-x_b}

que l'on peut aussi écrire (formule de Taylor-Young au premier ordre) :

 \bar{f}(x) = y_a + (x-x_a) \frac{y_b-y_a}{x_b-x_a}

ou encore :

 \bar{f}(x) = \frac{x_b-x}{x_b-x_a} y_a + \frac{x-x_a}{x_b-x_a} y_b

Cette dernière formule correspond à la moyenne pondérée.

Exemple[modifier | modifier le code]

Exemple 1[modifier | modifier le code]

Par exemple, si nous souhaitons déterminer f(2,5) alors que l'on connaît les valeurs de f(2) = 0,9093 et f(3) = 0,1411, cette méthode consiste à prendre la moyenne des deux valeurs sachant que 2,5 est le milieu des deux points. On obtient par conséquent f(2,5) \approx \frac{0,9093+0,1411}{2}=0,5252.

Exemple 2[modifier | modifier le code]

Soit à calculer sin 0,712 84
Les tables de Laborde donnent :

x sin x Δ (différence
tabulaire)
0,712 0,653 349 2
756 7
0,713 0,654 105 9
756 1
0,714 0,654 862 0

Dans ce cas :

x_a = 0,712 et x_b = 0,713, car 0,712 < 0,71284 < 0,713
y_a = 0,653 349 2
y_b - y_a = \Delta = 756 7 \times 10^{-4}
et enfin x - x_a = 0,712 84 - 0,712 = 84 \times 10^{-5}

On a donc : \sin 0,712 84 \approx 0,653 349 2 + \frac{756 7 \times 10^{-7} \times 84 \times 10^{-5}}{10^{-3}} = 0,653 984 8.

En pratique, on pose les calculs de la façons suivante :

sin 0,712 \approx 0,653 349 2
Δ × . . . . 84 . . . . . 635 628
sin 0,712 84 \approx 0,653 984 8

Le fabricant des tables explique qu'en procédant ainsi, l'erreur absolue est inférieure à 2 × 10-7. La même méthode est expliquée dans les tables de logarithmes de Bouvart et Ratinet.

Précision[modifier | modifier le code]

Cette méthode est rapide et aisée mais sa précision dépend beaucoup de l'écart  |x_b-x_a|.

Le calcul différentiel et intégral permet, moyennant certaines hypothèses, de calculer l’erreur d'interpolation au pire de cas. En particulier, si la fonction f que l’on souhaite interpoler est deux fois continûment dérivable et si x \in [x_a;x_b] alors l’erreur d’interpolation est donnée par :

 |f(x)-\bar{f}(x)| \le C(x_b-x_a)^2 \quad\quad C = \frac18 \max_{y\in[x_a,x_b]} |f''(y)|.

En d’autres termes, le majorant de l’erreur est proportionnel au carré de la distance entre les nœuds. D'autres méthodes d'interpolation permettent, en prenant appui sur trois ou davantage de points (au lieu de deux), d'espérer une diminution l'erreur d'interpolation : par exemple la méthode des différences divisées de Newton, ou l’interpolation polynomiale. Attention tout de même au phénomène de Runge, phénomène où l'erreur d'interpolation explose quand on augmente le nombre de points.

L’interpolation linéaire peut être utilisée pour la recherche des zéros d'une équation (Méthode de la fausse position) ou le calcul numérique d'intégrales (méthode des trapèzes).

Bibliographie[modifier | modifier le code]

  • Jean Laborde (préf. de Jean Kuntzmann), Tables numériques de fonctions élémentaires, éd. Dunod, 1961 (rééd. 1968, 1976) - 149 pages
  • André Delachet, Les logarithmes, Presses Universitaires de France, coll. « Que sais-je ? »,‎ 1960 (réimpr. 1964, 1968, 1972, 1977), 128 p., « Le calcul logarithmique », p. 47-80