TP2 Interpollation
TP2 Interpollation
TP2 Interpollation
TP II : Interpolation polynômiale
Les méthodes de l’interpolation polynômiale montrèrent qu'il existe une fonction interpolatrice
capable de relier un nombre arbitraire de points en un seul tenant. Donc il doit peut-être exister
un polynôme f(x) dont la courbe représentative passe par tous les points initiaux.
Pour rappel, un polynôme est une fonction usuellement présentée sous la forme d'une somme
de monômes 𝑓(𝑥) = 𝜃0 + 𝜃1 . 𝑥1 + 𝜃2 . 𝑥2 + ⋯ + +𝜃𝑛−1 . 𝑥𝑛−1 + +𝜃𝑛 . 𝑥𝑛 tel
que 𝜃0 , 𝜃1 , 𝜃2 , … . 𝜃𝑛−1 , 𝜃𝑛 sont les coefficients de ce polynôme et où la plus grande puissance
de x correspond au degré de ce polynôme. À noter qu'ici les coefficients seront classés par ordre
de puissance croissante par souci de cohérence avec les algorithmes.
METHODE DE LAGRANGE
L'idée est de réécrire le polynôme f(x) de degré n comme une combinaison linéaire de
polynômes Li(x) de degré n possédant les propriétés suivantes :
1. Les polynômes Li(x) forment une base dans l'espace vectoriel des polynômes de degré
inférieur ou égal à n.
1 𝑠𝑖 𝑖 = 𝑗
2. 𝐿𝑖 (𝑥𝑗 ) = {
0 𝑠𝑖 𝑖 ≠ 𝑗
MODULE : TP METHODES NUMERIQUES ENSEIGNANT : M. GASMI
Ces polynômes Li(x) sont communément nommés les polynômes d'interpolation de Lagrange.
𝑛
(𝑥 − 𝑥𝑖 )
𝐿𝑖 (𝑥) = ∏
(𝑥𝑖 − 𝑥𝑗 )
𝑗=0,𝑗≠𝑖
Il ne reste plus qu'à faire passer chaque polynôme Li(x) par l'ordonnée du point initial
correspondant et à en faire la somme pour obtenir finalement le polynôme interpolateur f(x) de
degré n :
𝑃𝑛 (𝑥) = ∑ 𝑦𝑖 . 𝐿𝑖 (𝑥)
𝑖=0
Algorithme
L’algorithme suivant est une traduction directe de la formule d’interpolation de
Lagrange.
X:=xi
Y:=Yi
S:=0
N:=length(X)
Pour :=1 a n faire
L :=1
Pour j :=1 a n faire
Si i<>j alors
L := L*((x –x[j])/(x[i] –x[j]))
Fin-Si
Fin-Pour
S :=S+L*Y[i]
Fin-Pour
Afficher (S)
MODULE : TP METHODES NUMERIQUES ENSEIGNANT : M. GASMI
Exemple
Soit par exemple une expérience où on enregistre la vitesse d’un véhicule en fonction de la
distance, les résultats sont donnés dans le tableau suivant
A partir des données du tableau, trouver la valeur de y à x = 550
Solution
MODULE : TP METHODES NUMERIQUES ENSEIGNANT : M. GASMI
Exercice
Soit une fonction 𝑓(𝑥)=𝑒^𝑥 définie sur l’intervalle [3.50,3.70] avec :
Ecrire un programme qui:
Déterminer le pas d’interpolation pour un nombre ‘n’ donné de sous intervalles (n =4);
Remplit un tableau avec les coordonnées des points d’appui ;
Interpole la fonction f(x) pour x=3.62 en utilisant la méthode de Lagrange.
METHODE DE NEWTON
L’idée est de se donner une autre base de telle sorte que l'ajout d'un nouveau point initial ne
nous obligera pas de recalculer les polynômes de cette base, mais simplement d'adjoindre un
nouveau polynôme à la base déjà existante.
En généralisant, ajoutons un ième point (xi,yi) à un ensemble de i−1 points déjà existants.
Pour un ensemble de i−1 points et son polynôme interpolateur fi−1(x), il est donc possible
d'ajouter un nouveau point et d'exprimer le nouveau polynôme interpolateur fi(x) en fonction
fi−1(x). Les produits des termes faisant apparaître les racines Ni(x)=(x−x0)(x−x1)…(x−xi−1)
sont appelés polynômes de Newton et forment une base.
𝑖−1
𝑁𝑖 (𝑥) = ∏(𝑥𝑖 − 𝑥𝑗 )
𝑗=1
Le polynôme interpolateur f(x) peut ainsi s'écrire comme une combinaison linéaire de ces
polynômes de Newton :
𝑓 (𝑥) = ∑ 𝛼𝑖 . 𝑁𝑖 (𝑥)
𝑖=1
MODULE : TP METHODES NUMERIQUES ENSEIGNANT : M. GASMI
Il ne reste à présent qu'à expliciter les coefficients αi que nous nommerons coefficients de
Newton. Ils peuvent s'exprimer sous la forme de différences divisées, usuellement notées […],
tel que αi=[y0,…,yi].
Algorithme
X:=xi
Y:=Yi
n:=length(X)-1
Pour i:=1 a n faire
F[i,0] :=Y[i]
Fin-pour
Pour i :=1 a n faire
Pour j :=0 a n-i faire
F[i,j]:= (𝑭[𝒊 − 𝟏, 𝒋 + 𝟏] − 𝑭[𝒊 − 𝟏, 𝒋])/(𝒙[𝒊 + 𝒋] − 𝒙[𝒋])
Fin-pour
Fin-pour
P:=1
S:=F[0,0]
Pour i :=0 a n-1 faire
P:=P*(x-x(j))
S:=S+F[i+1,0]
Fin-Pour
Afficher (S)
Exemple
Soit par exemple les résultats d’une étude économique où on a enregistré les prix des
appartements en fonction de leurs surfaces, les résultats sont donnés dans le tableau suivant
Solution
Exercice