See discussions, stats, and author profiles for this publication at: [Link]
net/publication/342424107
Interpolation de Lagrange avec Python - Master de Recherche - Math
Fondamentale (MR2MF)
Presentation · June 2020
DOI: 10.13140/RG.2.2.13554.53442
CITATIONS READS
0 15,270
1 author:
Mohamed Mabrouk
Military Academy, Tunisia
20 PUBLICATIONS 34 CITATIONS
SEE PROFILE
All content following this page was uploaded by Mohamed Mabrouk on 24 June 2020.
The user has requested enhancement of the downloaded file.
TD - Interpolation de Lagrange avec Python
2ème MASTER RECHERCHE MATH FONDAMENTALE
A. U. 2019/2020
Enseignant: MABROUK Mohamed
1. Introduction théorique.
le terme latin « interpolare » signifie « refaire, donner une nouvelle forme ».
Les polynômes de Lagrange au cœur de l’interpolation.
Les polynômes de Lagrange offrent un cadre idéal pour l'interpolation. Soit x0, x1,…, xn, n+1 réels
(points) distincts. Il existe n+1 polynômes Li de degrés n pour i = 0 à n définis par:
Ou bien :
Avec :
Interpoler la fonction f par un polynôme p de degré ≤ n, aux n+1 points x0, x1,…, xn,
revient à trouver le polynôme p tel que :
Si un tel polynôme existe, il s’écrit de manière unique :
En prenant x = xj , on a alors :
Par conséquent, le polynôme de Lagrange p(x) interpolant une fonction f s’écrit sous
la forme suivante :
Exemple : f(x) = x2 , n=2 → 3 points (x0=0, x1=1, x2=2).
Selon :
2. Simulation Numérique.
2.1. Données.
• numpy.poly1d(argument) : renvoie une fonction polynomiale. On peut donc l'utiliser afin de manipuler
le symbole x (polynôme) dans les Li et le P.
Exemple d’utilisation :
numpy.poly1d([ 1 , 0 ]) renvoie le polynôme 1x1+0x0 = x
numpy.poly1d([ 3 , 2 , 6 ]) renvoie le polynôme 3x2+2x1+6x0 = 3x2+2x+6
• [Link](p,x) : évalue un polynôme p (défini par numpy.poly1d) pour une valeur spécifique de x.
• lagrange([xi],[f(xi)]) est fonction du module « interpolate » de la bibliothèque « scipy » qui permet de
renvoyer le polynôme de Lagrange interpolant f aux points xi.
Pour nôtre TD on prendra comme exemple de fonction à interpoler : f(x) = 1/(1+x2) sur [-5,5].
2.2. Travail demandé :
1. Définir et représenter graphiquement la fonction f dans l’intervalle [-5,5].
2. Ecrire une fonction Lagrange(f, a, b, n) qui calcule le polynôme de dégrée n interpolant
la fonction f en n+1 points d’abscisse xi équidistants dans l’intervalle [a,b].
3. Tester votre fonction « Lagrange (f,a,b,n)» pour n=3, 5 et 10 et afficher, pour chacun des
cas, le polynôme de Lagrange correspondant.
4. Représenter, sur le même graphique, la courbe de f et les polynômes de Lagrange trouvés
dans la question 3. Constatation ?
Réponses
1. Définir et représenter graphiquement la fonction f dans l’intervalle [-5,5].
2. Ecrire une fonction Lagrange(f, a, b, n) qui calcule le polynôme de dégrée n interpolant la
fonction f en n+1 points d’abscisse xi équidistants dans l’intervalle [a,b].
3. Tester votre fonction « Lagrange (f,a,b,n)» pour n=3, 5 et 10 et afficher, pour chacun
des cas, le polynôme de Lagrange correspondant.
4. Représenter, sur le même graphique, la courbe de f et les polynômes de Lagrange
trouvés dans la question 3. Constatation ?
Un phénomène d'oscillation
apparaît aux bords de la
courbes de f(x) et s'amplifie
lorsque n augmente
==> Phénomène de Runge ou
Effet des bords.
View publication stats