Accueil🇫🇷Chercher

Interpolation d'Hermite

En analyse numérique, l'interpolation d'Hermite, nommée d'après le mathématicien Charles Hermite, est une extension de l'interpolation de Lagrange, qui consiste, pour une fonction dérivable donnée et un nombre fini de points donnés, à construire un polynôme qui est à la fois interpolateur (c'est-à-dire dont les valeurs aux points donnés coïncident avec celles de la fonction) et osculateur (c'est-à-dire dont les valeurs de la dérivée aux points donnés coïncident avec celles de la dérivée de la fonction). Cette méthode d'interpolation permet d'éviter les phénomènes de Runge dans l'interpolation numérique ou, plus simplement, de manipuler des polynômes ayant des propriétés proches de celles de la fonction interpolée.

Comparaison graphique entre interpolation lagrangienne (en rouge) et hermitienne (en bleu) de la fonction (en noir) en trois points équidistants -1, 1/2, 2.

Définition du problème à l'ordre 1

Soit f une fonction de classe C1 d'une variable définie sur un segment [a, b] et à valeurs réelles et soient n + 1 points (x0, x1, ... , xn) de [a, b] distincts deux à deux. L'objectif est de construire un polynôme P de degré minimal tel que :

.

Puisque l'on impose 2n + 2 valeurs pour déterminer le polynôme P, celui-ci sera donc de degré au plus 2n + 1.

Construction

Une méthode de construction de P consiste à prendre les carrés des polynômes de Lagrange associés aux points x0, x1, ... , xn :

,

de degré 2n et vérifiant :

.

Un polynôme P de la forme

satisfait donc les 2n + 2 conditions si et seulement si les polynômes Pi vérifient :

,

ce qui équivaut à :

.

La solution la plus simple est de choisir

et P est alors de degré au plus 2n + 1.

Unicité

L'unicité du polynôme interpolateur d'Hermite se montre de façon similaire à celle du polynôme interpolateur de Lagrange : soient deux polynômes P et R vérifiant les hypothèses voulues. On a donc deux polynômes de degré au plus 2n + 1 dont les valeurs et les dérivées coïncident en n + 1 points. Ainsi, PR est divisible par (Xx0)2(Xx1)2…(Xxn)2 qui est un polynôme de degré 2n + 2. Puisque PR est de degré au plus 2n + 1, il est forcément nul.

Majoration de l'erreur

L'erreur d'approximation causée par l'interpolation d'Hermite est donnée par le résultat suivant :

Si f est de classe C2n+2, alors pour tout x dans [a, b], il existe ξ dans ]a, b[ tel que

.

Ainsi, l'interpolation d'Hermite est d'un ordre très supérieur à celui de l'interpolation lagrangienne (d'ordre n + 1).

Extension aux ordres supérieurs

L'interpolation d'Hermite peut être étendue à l'interpolation des valeurs des dérivées supérieures, en cherchant, pour une fonction f de classe Cm sur [a, b], un polynôme interpolateur P vérifiant :

.

Le polynôme à construire est donc de degré minimal (n + 1)(m + 1) – 1. Une méthode pour le définir consiste à introduire les polynômes

où les Lk sont les polynômes de Lagrange définis précédemment et les Qk,l sont des polynômes de degré m – l tels que

et pour tout q tel que .

Ainsi, par construction, on a :

.

Le polynôme P recherché s'écrit alors :

.

Cette méthode apporte plus de régularité à l'interpolation. Elle reste cependant d'un faible intérêt pratique au regard des calculs qu'elle implique[1].

Applications

L'interpolation d'Hermite en deux points est la base des splines cubiques. Voir également l'article Spline cubique d'Hermite.

L'interpolation d'Hermite peut également être utilisée dans la résolution de problèmes aux limites non linéaires[2].

Notes et références

  1. « lumimath.univ-mrs.fr/~jlm/trav… »(Archive.orgWikiwixArchive.isGoogle • Que faire ?).
  2. (en) R. E. Grundy, « The application of Hermite interpolation to the analysis of non-linear diffusive initial-boundary value problems », IMA Journal of Applied Mathematics, vol. 70, no 6, , p. 814–838 (lire en ligne)

Lien externe

(en) Eric W. Weisstein, « Hermite's Interpolating Polynomial », sur MathWorld

Cet article est issu de wikipedia. Text licence: CC BY-SA 4.0, Des conditions supplémentaires peuvent s’appliquer aux fichiers multimédias.