TP2 Interpollation

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 6

DEPARTEMENT SCIENCES ET TECHNOLOGIE

Module : TP Méthodes Numériques Enseignant : M. GASMI

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.

Par exemple, la fonction 𝑓 (𝑥) = 5 + 𝑥 2 − 3𝑥 3 est un polynôme de degré 3 et de coefficients


(5,0,1,-3)

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.

Les polynômes Li(x) ont la forme suivante :

𝑛
(𝑥 − 𝑥𝑖 )
𝐿𝑖 (𝑥) = ∏
(𝑥𝑖 − 𝑥𝑗 )
𝑗=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

Ce polynôme vérifie par construction f(xi)=yi.

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'interpolation de Newton respecte les même conditions que l'interpolation de Lagrange, à


savoir f(xi)=yi où pour tout point initial (xi,yi). Le polynôme sous-jacent est le même que celui
obtenu par la forme de Lagrange avec le but de calculer plus efficacement les valeurs du
polynôme interpolateur.

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

A partir des données du tableau, trouver la valeur de y à x=100m².


MODULE : TP METHODES NUMERIQUES ENSEIGNANT : M. GASMI

Solution

Exercice

Soit la fonction définie par 𝑓 (𝑥) = 3√𝑥


Ecrire sous Matlab un programme qui interpole la fonction f pour x=15 en utilisant la
méthode de newton.
Le polynôme d'interpolation construit sur (xi , f(xi)), i = 0 à 4, avec x0 = 0 , x1 = 1 , x2 = 8 ,
x3 = 27 , x4 = 64.

-Comparer les résultats

Vous aimerez peut-être aussi