2024-2025
1A Génie Mécanique
EMI - UM5 de Rabat
TP4
Méthode d’interpolation
Nous avons vu dans le cours que le polynôme d’interpolation de Lagrange peut être calculé
par la méthode de Lagrange. Cette méthode consiste à trouver une base de polynômes,
qui sont les polynômes de Lagrange, dans laquelle on exprime le polynôme d’interpolation.
Étant donné une fonction f et un ensemble (xj )0≤j≤n de points d’interpolation, le polynôme
d’interpolation p s’exprime comme
n
X
p(x) = f (xj )Lj,n (x), ∀x ∈ R,
j=0
avec (Lj,n )0≤j≤n la base de polynômes de Lagrange donnée par
n
Y x − xk
Lj,n (x) = .
k=0
xj − xk
k̸=j
1. Implémentation de la méthode de Lagrange
Écrire une fonction Methode Lagrange qui implémente la méthode précédente avec
en arguments d’entrée un vecteur x interp de points d’interpolation et un vecteur x de
points d’abscisse. Cette fonction retournera la matrice L dont chaque colonne correspond
au vecteur ℓi = Li,n (x).
On obtient alors le vecteur y = p(x) grâce au produit matrice-vecteur y = Lz où
z = f (x interp). Le produit matrice-vecteur peut se faire grâce à la fonction [Link].
2. Tracé de la fonction Runge et du polynôme d’interpolation
Tracer ensuite sur une même figure la fonction Runge sur l’intervalle [−3, 3] ainsi que
son polynôme d’interpolation obtenu avec la méthode de Lagrange pour 5, 10 et 15 points
d’interpolation.
Méthode de Newton
La méthode de Newton, comme pour celle de Lagrange, permet d’exprimer le polynôme
d’interpolation de Lagrange dans une base de polynômes, qui sont les polynômes de
Newton. Étant donné une fonction f et un ensemble (xj )0≤j≤n de points d’interpolation,
le polynôme d’interpolation p s’exprime comme
n
X
p(x) = zj Nj,n (x), ∀x ∈ R,
j=0
1
avec (Nj,n )0≤j≤n la base de polynômes de Newton donnée par
1, pour j = 0
Nj,n (x) = Q
j−1 (x − x ), pour j ≥ 1
k=0 k
et (zj )0≤j≤n les coefficients qui sont solution du système triangulaire
1 0 0 ··· 0
1 (x1 − x0 ) 0 ··· 0 z0 y0
.. .. z1 y1
1 (x − x ) (x − x )(x − x )
2 0 2 0 2 1 . . .. = .. .
. . . ... . .
.. .. .. 0
Qn−1 zn yn
1 (xn − x0 ) (xn − x0 )(xn − x1 ) ··· j=0 (xn − xj )
1. Implémentation de la méthode de Newton
Écrire une fonction Methode Newton qui implémente la méthode précédente avec en
arguments d’entrée un vecteur x interp de points d’interpolation et un vecteur x de points
d’abscisse. Cette fonction retournera la matrice N dont chaque colonne correspond au
vecteur ni = Ni,n (x).
2. Calcul des coefficients de Newton
Écrire une fonction Coeff Newton qui permet de calculer les coefficients (zj )0≤j≤n en
résolvant le système triangulaire inférieur par un algorithme de descente. Cette fonction
prendra en entrée le vecteur x interp des points d’interpolation ainsi que la fonction f à
interpoler et rendra le vecteur z des coefficients.
3. Tracé de la fonction Runge et du polynôme d’interpolation
Tracer ensuite sur une même figure la fonction Runge sur l’intervalle [−3, 3] ainsi que
son polynôme d’interpolation obtenu avec la méthode de Newton pour 5, 10 et 15 points
d’interpolation.