- Interpolation newtonienne
-
En analyse numérique, l'interpolation newtonienne, du nom d'Isaac Newton, est une méthode d'interpolation polynomiale permettant d'obtenir le polynôme de Lagrange comme combinaison linéaire de polynômes de la base newtonienne.
Contrairement à l'interpolation de Hermite par exemple, cette méthode ne diffère de l'interpolation lagrangienne que par la façon dont le polynôme est calculé, le polynôme d'interpolation qui en résulte est le même. Pour cette raison on parle aussi plutôt de la forme de Newton du polynôme de Lagrange.
Sommaire
Définition
Étant donné un ensemble de k+1 points
- (les xj tous distincts 2 à 2), l'interpolation polynomiale dans une base de Newton est une combinaison linéaire de polynômes appartenant à cette base
avec les polynômes de Newton définis de la manière suivante
et les coefficients comme ceci
où
est la notation pour les différences divisées.
Par conséquent, le polynôme d'interpolation peut être écrit ainsi
Idée principale
Résoudre un problème d'interpolation conduit à un problème d'algèbre linéaire où nous devons inverser une matrice. En utilisant une méthode standard d'interpolation polynomiale, on obtient une matrice du type matrice de Vandermonde qui est très difficile à inverser. En choisissant une base newtonienne de polynômes, on obtient une matrice triangulaire inférieure qui s'inverse beaucoup plus facilement.
Pour k+1 points, on construit la base de Newton ainsi
En utilisant ces polynômes, nous devons inverser
pour résoudre le problème d'interpolation polynomiale.
Cette matrice peut être inversée successivement en résolvant
On commence par j = 0 qui nous donne a0 puis pour j = 1, le calcul de a0 nous permet de déduire a1. Et ainsi de suite jusqu'à j = k.
Application
Comme le montre la définition des différences divisées, des points supplémentaires peuvent être ajoutés pour créer un nouveau polynôme d'interpolation sans recalculer les coefficients. De plus, si un point est modifié, il est inutile de recalculer l'ensemble des coefficients. Autre avantage, si les xi sont équi-répartis, le calcul des différences divisées devient nettement plus rapide. Par conséquent, l'interpolation polynomiale dans une base de Newton est priviligiée par rapport à une interpolation lagrangienne pour des raisons pratiques.
Lien externe
Interpolation polynômiale (sic) de type Newton et différences divisées sur math-linux.com
Wikimedia Foundation. 2010.