Interpolation polynomiale

En mathématiques, en analyse numérique, l'interpolation polynomiale est une technique d'interpolation d'un ensemble de données ou d'une fonction par un polynôme. En d'autres termes, étant donné un ensemble de points (obtenu, par exemple, à la suite d'une expérience), on cherche un polynôme qui passe par tous ces points, p(x_i) = y_i, et éventuellement vérifie d'autres conditions, de degré si possible le plus bas.

Cependant, dans le cas de l'interpolation lagrangienne, par exemple, le choix des points d'interpolation est critique. L'interpolation en des points régulièrement espacés peut fort bien diverger même pour des fonctions très régulières (phénomène de Runge).

Définition

Les points rouges correspondent aux points

(x k, y k)

, et la courbe bleue représente le polynôme d'interpolation.

Dans la version la plus simple (interpolation lagrangienne), on impose simplement que le polynôme passe par tous les points donnés. Étant donné un ensemble de $n + 1$ points, i.e. couples $(x i, y i)$ (où les réels $x i$ sont distincts 2 à 2, les y_i pouvant être des réels, complexes ou éléments d'un espace vectoriel quelconque), on cherche à trouver un polynôme $p$ (à coefficients de la même nature que les y_i) de degré $n$ au plus, qui vérifie :

p(x_{i})=y_{i}{\mbox{ , }}i=0,\ldots ,n

.

Le théorème de l'unisolvance précise qu'il n'existe qu'un seul polynôme $p$ de degré inférieur ou égal à $n$ défini par un tel ensemble de $n + 1$ points.

L'interpolation d'Hermite consiste à chercher un polynôme qui non seulement prend les valeurs fixées aux abscisses données, mais dont également la dérivée, donc la pente de la courbe, prend une valeur imposée en chacun de ces points. Naturellement, il faut pour cela un polynôme de degré supérieur au polynôme de Lagrange. On peut aussi imposer encore la valeur des dérivées secondes, troisièmes, etc. en chaque point. La démarche de l'interpolation newtonienne utilisant les différences divisées est particulièrement adaptée pour construire ces polynômes.
La méthode des splines consiste à chercher des fonctions polynômiales par morceaux, c'est-à-dire sur chaque sous-intervalle [x_i-1,x_i], mais de plus bas degré (typiquement 3 pour les splines cubiques), en choisissant les coefficients pour obtenir une fonction continue et dérivable également aux points x_i .

Construction du polynôme d'interpolation de Lagrange

On voit aisément que la combinaison linéaire $p(x)=\sum y_{i}L_{i}(x)$ vérifie bien p(x_i) = y_i pour i = 0,..., n, si les polynômes (L_i) _{i = 0,...,n} vérifient L_i(x_j) = δ_ij = 1 si i = j, 0 sinon (voir symbole de Kronecker). Il est tout aussi évident que c'est bien le cas pour $L_{i}:=\prod _{j\neq i}{\frac {x-x_{j}}{x_{i}-x_{j}}}$ , où le produit porte sur tous les indices j dans { 0, ..., n } \ { i }. La propriété caractéristique L_i(x_j) = δ_ij implique immédiatement que la famille (L_i) est libre, donc une base de R_n[x], appelée la base de Lagrange (ou lagrangienne) relative à la famille (x_i) _{i = 0,...,n} .

Erreur d'interpolation

L'erreur d'interpolation lors de l'approximation d'une fonction $f$ , c'est-à-dire : lorsque $y i = f (x i)$ dans ce qui précède, est donnée par une formule de type Taylor-Young :

Si $f$ est $n + 1$ fois différentiable sur $I = [min(x 0,..., x n, x),max(x 0,..., x n, x)]$ alors

$\exists \xi \in I\quad f(x)-p_{n}(x)={\frac {f^{(n+1)}(\xi )}{(n+1)!}}\prod _{i=0}^{n}(x-x_{i}).$

L'existence d'un tel $ξ$ se démontre en appliquant de manière itérée le théorème de Rolle[1] :

Démonstration

Soit $x\in \mathbb {R}$ . Si $x$ est un point d'interpolation, $f (x) - p n (x) = 0$ et la formule est vérifiée.

Dans le reste de la démonstration, on suppose que $x$ n'est pas une abscisse d'interpolation. Introduisons une fonction auxiliaire $g$ :

g:\mathbb {R} \rightarrow \mathbb {R} :t\rightarrow g(t)=f(t)-p_{n}(t)-{\frac {f(x)-p_{n}(x)}{\prod _{i=0}^{n}(x-x_{i})}}\prod _{i=0}^{n}(t-x_{i}).

Cette fonction $g$ possède $n + 2$ racines distinctes :

x_{0},\dots ,x_{n},x.

Par application du théorème de Rolle, $g'$ , dérivée de $g$ , possède $n +1$ racines distinctes (toutes situées exactement entre deux racines successives de $g$ ). En appliquant encore $n$ fois le théorème de Rolle, on obtient que $\exists \xi \in I$ tel que $0=g^{(n+1)}(\xi )=f^{(n+1)}(\xi )-{\frac {f(x)-p_{n}(x)}{\prod _{i=0}^{n}(x-x_{i})}}(n+1)!$ (puisque la dérivée d'ordre $n +1$ de $p n$ est nulle).

En isolant $f (x) - p n (x)$ on obtient le résultat escompté :

f(x)-p_{n}(x)={\frac {f^{(n+1)}(\xi )}{(n+1)!}}\prod _{i=0}^{n}(x-x_{i})

Dans le cas particulier où $x i = x 0 + ih$ (points uniformément répartis), se produit en général une aggravation catastrophique de l'erreur d'interpolation, connue sous le nom de phénomène de Runge, lorsqu'on augmente le nombre de points pour un intervalle $[x 0, x n]$ donné.

Références

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Polynomial interpolation » (voir la liste des auteurs).

(en) Endre Süli et David F. Mayers, An Introduction to Numerical Analysis, Cambridge University Press, 2003, p. 183-184 sur Google Livres.

Voir aussi

Portail de l'analyse

Cet article est issu de Wikipedia. Le texte est sous licence Creative Commons - Attribution - Partage dans les Mêmes. Des conditions supplémentaires peuvent s'appliquer aux fichiers multimédias.

[1] (en) Endre Süli et David F. Mayers, An Introduction to Numerical Analysis, Cambridge University Press, 2003, p. 183-184 sur Google Livres.