TD03 Et 04

Pôle Universitaire Léonard de Vinci
MESIMA350223 – Méthodes numériques
TD No.3 et 4 (durée 3h)

Interpolation polynômiale. Méthode de moindres carrés.
Exercice I (Application d’interpolation)

Soit f la fonction définie par f (x) = sin(πx) + cos(2πx). On considère 3 points x0 = 0,
x1 = 1/4 et x2 = 1/2.
n o
(1) Déterminer le polynôme P (x) d’interpolation lié aux 3 points (xi , f (xi ) .
i=0,1,2
(2) En déduire une valeur approchée de sin(π/5) + cos(2π/5).

(3) Donner une estimation des erreurs de la valeur approchée obtenue en question (2) par
rapport à sa valeur exacte, sans calculer cette dernière.
(4) Que ce que l’on peut faire si l’on veut obtenir une valeur approchée plus précise ?
Exercice II (Interpolation de Newton)

Etant donné 10 points suivant
(x0 , y0 ) = (0, 2), (x1 , y1 ) = (1, 6), (x2 , y2 ) = (2, 12), (x3 , y3 ) = (3, 20), (x4 , y4 ) = (4, 30),
(x5 , y5 ) = (5, 42), (x6 , y6 ) = (6, 56), (x7 , y7 ) = (7, 72), (x8 , y8 ) = (8, 90), (x9 , y9 ) = (9, 110).
(1) Déterminer le polynôme d’interpolation lié aux points (noté P (x)).

(2) Supposons qu’il y a une perturbation sur le point (x5 , y5 ) qui devient (5, 42+ε). On note
Pε (x) le polynôme d’interpolation associé au nuage des points perturbés. Déterminer
la valeur de Pε (10).
Exercice III
−1 0 1 2
Etant donné quatre points : , , , , on cherche des polynômes des
−1 0 1 8
moindres carrés.
(1) Déterminer la parabole (i.e. un polynôme de degré 2) des moindres carrés.
(2) Déterminer le polynôme de degré 3 des moindres carrés.
(3) Calculer le polynôme d’intrepolation de ces 4 points. Comparer à la question (2) : que
l’on peut conclure?
Exercice IV (Application de la méthode de moindres carrés)

Dix adolescents s’exercent à lancer le poids, du bras droit puis du bras gauche. Les résultats
(distances en mètres) obtenus sont les suivants :
1
Construire un modèle qui permet d’établir un lien linéaire entre les distances de lance par le
bras droit et par le bras gauche. Utiliser le modèle afin de répondre les questions suivantes.
(1) Un adolescent qui lance à 6.5m du bras droit peut espérer lancer à combien de mètres
avec le bras gauche ?
(2) Un adolescent lance le poids à 4.2m du bras gauche. Quelle sera sa performance avec
le bras droit ?
Exercice V (Interpolation en Python)

Soit (xi , yi )i=0,··· ,n un nuage des n + 1 points donnés. Le but est de
— construire le polynôme d’interpolation (de Langrange) associé à ces points ;
— appliquer votre programme dans différents cas,
— et savoir analyser les résultats obtenus.
(1) Construction du programme

• Analyser le programme suivant
import numpy as np
def Li(i,X):
x=np.poly1d([1, 0])
L=1
n=len(X)
for j in range(n):
if j != i:
L=L*(x-X[j])/(float(X[i]-X[j]))
return L
def Interpol_Lagrange(X,Y):
s=0
n=len(X)
for i in range(n):
s=s+Li(i,X)*Y[i]
return s
Les X et Y représentent quels donnés dans la formule mathématique ? Chaque boucle

dans ce programme correspond quelle opération de la formule mathématique ?
2
• Créer le programme principal qui
— initialise le nuage des points et les stocke dans les tableaux X et Y ;
— détermine le polynôme d’interpolation en appelant la fonction
Interpol Lagrange(X, Y );
— trace sur le même graphe les points saisis et la courbe du polynôme d’interpo-
lation dans un intervalle qui contient tous points saisis.
(2) Application et analyse du cas 1

On considère le nuage des points suivant :
(0, 1), (−1, −1), (1, 5), (3, 19), (4, 29).
• Déterminer le polynôme d’interpolation (noté P (x)) et tracer sa courbe.
• Changer légèrement la coordonnée d’un point (remplacer (0, 1) par (0, 1 + ε), par
exemple), puis déterminer le polynôme d’interpolation (noté Pε (x)). Tracer sa courbe
sur le même graphe. Observer comment la courbe de Pε (x) change en fonction de ε.
(3) Application et analyse du cas 2

Soit (xi , yi )i=0,··· ,n un nuage des n + 1 points donnés par xi = i et yi = i3 .
• Déterminer son polynôme d’interpolation et tracer sa courbe en utilisant le programme
réalisé en (1) et en supposant n = 5.
• Si on choisit n = 20 ? n = 50 ? Comment expliquer les résultats obtenus ?
(4) Pour ceux qui veulent aller plus loin ...

On a vu en cours que l’interpolation polynômiale peut se servir de l’approximation d’une fonc-
tion f . On va montrer ici que la stratégie du choix des points à interpoler a une conséquence
très importante sur la qualité de cette approximation.
On considère la fonction f définie par
1
f (x) = ,
1 + x2
et on voudrait approcher cette f sur l’intervalle [−5, 5] par des polynômes d’interpolation.
(4.A) D’abord, on travaille avec N abscisses d’interpolation équidistantes.
(a) On prend N = 3, i.e. on choisit 3 points équidistants x0 = −5, x1 = 0 et x2 = 5.

Utiliser votre programme afin de
— donner le polynôme d’interpolation dans ce cas ;
— tracer, dans la même figure, la courbe de la fonction f en couleur rouge et la courbe
du polynôme d’interpolation en bleu, tous les deux sur l’intervalle [−5, 5].
(b) Re-faire la question (a) en prenant N = 6 puis N = 11, où N désigne le nombre des
abscisses d’interpolation équidistantes.
(c) En comparant les trois cas étudiés, quelle est votre conclusion ?
3
(4.B) Maintenant, on travaille avec des points de Tchebycheff {xi }i=1,··· ,N définis par
2i − 1
xi = 5 cos π , i = 1, · · · , N.
2N
(d) Re-faire la question (a) en prenant respectivement N = 3, N = 6 et N = 11 avec ce
choix de xi .
(e) Quelle amélioration est apportée par le choix des points de Tchebycheff ?

TD03 Et 04

Transféré par

Droits d'auteur :

Formats disponibles

TD03 Et 04

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TD03 Et 04

Transféré par

Droits d'auteur :

Formats disponibles

Pôle Universitaire Léonard de Vinci

MESIMA350223 – Méthodes numériques

TD No.3 et 4 (durée 3h)

Exercice I (Application d’interpolation)

(2) En déduire une valeur approchée de sin(π/5) + cos(2π/5).

Exercice II (Interpolation de Newton)

(1) Déterminer le polynôme d’interpolation lié aux points (noté P (x)).

Exercice IV (Application de la méthode de moindres carrés)

Exercice V (Interpolation en Python)

(1) Construction du programme

Les X et Y représentent quels donnés dans la formule mathématique ? Chaque boucle

(2) Application et analyse du cas 1

(3) Application et analyse du cas 2

(4) Pour ceux qui veulent aller plus loin ...

(a) On prend N = 3, i.e. on choisit 3 points équidistants x0 = −5, x1 = 0 et x2 = 5.

Vous aimerez peut-être aussi