TD5 An 2015

Université Claude Bernard, Lyon 1 Licence Sciences & Technologies
43, boulevard du 11 novembre 1918 Spécialité : Mathématiques

69622 Villeurbanne cedex, France Analyse numérique L3- Automne 2015
Série d’exercices no 5/6

Interpolation polynomiale
Exercice 1. Formule des Différences Divisées (Un Classique)

Nous supposons que f : [a, b] ! R est une fonction n + 1 fois continûment différentiable. La
formule de Newton qui consiste à écrire le polynôme Pn aux points x0 ,..., xn sous la forme
Pn (x) = a0 + a1 (x x0 ) + ... + an (x x0 )...(x xn 1 ),
permet de construire le polynôme Pn à l’aide d’une récurrence. En effet,

n
Y1
Pn (x) = Pn 1 (x) + an (x xk ).
k=0
Autrement dit, connaissant Pn 1 , il suffit de calculer an pour connaître Pn .
a) Montrer que le polynôme d’interpolation de Lagrange de la fonction f aux points distincts

(xi )1in est donné par
Xn i 1
Y
Pn (x) = f [x0 , ..., xi ] (x xk ),
i=0 k=0
où f [.] désigne les différences divisées de f définies par

8
< f [xi ] = f (xi ),
1 pour tout i = 0, ..., n.
: f [x0 , ..., xk ] = (f [x1 , ..., xi ] f [x0 , ..., xi 1 ]),
xk x0
Montrer ensuite que f [x0 , ..., xn ] est invariant par permutations.
b) Montrer qu’il existe ⇠ 2 [a, b] tel que
f (n) (⇠)
f [x0 , ..., xn ] = .
n!
c) Montrer que
Mn+1
|P n(x) f (x)|  |⇡n (x)|,
(n + 1)!
1
où
i 1
Y
Mn+1 = max |f (n+1) (x)|, et ⇡n (x) = (x xi ).
axb
i=0
N.B. : Remarquons bien ici que l’estimation n’est pas forcément quelque chose de petit (voir Phé-
nomène de Runge).
Application.
Trouver l’interpolation de Lagrange de la fonction x ! f (x) = sin(⇡x/2) aux points x0 = 0,
x1 = 1 et x2 = 2. Puis à l’aide des questions précédentes établir une estimation d’erreur.
Exercice 2. Convergence de l’interpolatio de Lagrange Soit Ln le polynôme d’interpolation de

Lagrange de la fonction
1
f (x) = , 1  x  1,
x ↵
aux n + 1 points distincts x0 , ..., xn de l’intervalle [ 1, 1].
1. Calculer les dérivées successives de la fonction f .
2. Montrer que si ↵ > 3, et si les n + 1 points x0 , ..., xn sont équidistants, nous avons alors
lim kf Ln k1 = 0.
n!+1
3. Considérons toujours la fonction f

1
f (x) = , 1  x  1,
x ↵
aux n + 1 points distincts x0 , ..., xn équidistants de l’intervalle [ 1, 1]. Dans la pratique
nous n’agissons pas du tout comme ce qui précède. Nous préférons utiliser des polynômes
de degré peu élevé sur chaque petit intervalle [xi , xi+1 ].
Écrire l’approximation de Lagrange de degré 1, fn de f sur chaque intervalle [xi , xi+1 ],
i = 0, ..., n 1
4. Montrer que si ↵ 6= [ 1, 1], nous avons
c
kf L n k1 
n2
et donc que fn converge uniformément vers f lorsque n tend vers l’infini.
Exercice 3. Interpolation Polynomiale de Hermite

Soient x0 , ..., xn , n + 1 points distincts de l’intervalle [a, b], (a, b 2 R, a < b) et f de classe
C ([a, b], R) une fonction dont on connaît les valeurs et celles de sa dérivée en ces (n + 1) points
distincts.
Nous cherchons un polynôme Hn de degré minimal tel que
Hn (xi ) = f (xi ) et Hn0 (xi ) = f 0 (xi ), i = 0, ..., n.

2
Nous rappelons que les fonctions de base de l’interpolation de Lagrange, c’est à dire les polynômes
de degré n tels que li (xj ) = ij pour i, j = 0, ..., n sont donnés pour tout i = 0, ..., n par
Yn
x xj
li (x) = , pour tout x 2 R.
j=0
xi xj
j6=i
Nous allons montrer le résultat suivant :

“Le polynôme Hn s’écrit
n
X n
X
Hn (x) = f (xi )hi (x) + f 0 (xi )h̃i (x)
i=0 i=0
avec
hi (x) = (1 2)li0 (xi )(x xi ))li2 (x), et h̃i (x) = (x xi )li2 (x).
De plus, si f 2 C 2(n+1) ([a, b], R)
kf (2(n+1)) k1 n
|f (x) Hn (x)|  ⇧i=0 (x xi )2 .”
(2n + 2)!
1. Montrer que pour i, j = 0, ..., n
hi (xj ) = i,j , h0i (xj ) = 0,
et
h̃i (xj ) = 0, h̃0i (xj ) = i,j .
2. En déduire qu’il existe un unique polynôme Hn de degré 2n + 1 vérifiant les conditions

requises.
3. En déduire une majoration de l’erreur |f (x) Hn (x)| .
Exercice 4. Moindre carrés discrets

Nous rappelons tout d’abord le théorème de la projection sur un sous-espace vectoriel :
“Soient F un sous-espace vectoriel de Rn+1 et y 2 Rn+1 . Alors il existe un unique v ⇤ 2 F tel que
kv ⇤ yk = min kv yk.
v2F
De plus, v ⇤ est la projection orthogonale de y sur l’espace F : v ⇤ = PF y est telle que
(v ⇤ y, v) = 0, v 2 F .”
Objectif de l’exercice : soient '1 (x) = 1 et '2 (x) = x et '3 (x) = x2 . Nous recherchons le
polynôme de degré 2 qui approche le mieux le nuage de points (xi , yi )1i4 suivant :
( 1, 1), (0, 0), (1, 1), (2, 2).
Autrement dit, nous souhaitons trouver '⇤ 2 V := vect('1 , '2 , '3 ) telle que
3
P
4 P
4
|'⇤ (xi ) yi |2 = min | ⇤
(xi ) yi | 2
i=1 2V i=1
1. Tracer le nuage de points.

P
3
2. Écrire le polynôme '(x) = uj 'j (x) et expliciter le problème sous la forme :
j=1
“trouver v ⇤ 2 F solution de
kv ⇤ yk = min kv yk, où y = (y1 , ..., y4 )T , et F = {v 2 R4 , v = Bu, u 2 R3 }.”
v2F
3. Montrer que ce problème admet une solution unique v ⇤ 2 F .

4. Montrer que si v ⇤ 2 F est solution, alors il existe un unique u⇤ 2 R3 solution de
B T Bu = B T Y .
5. Montrer que B est de rang 3. En déduire alors que B T B est définie positive.
6. Expliciter la solution du problème.
Exercice 5. Polynôme de Chebychev

Soit n 2 N, nous définissons le polynôme de Chebychev de première espèce par
Tn (x) = cos(n arccos(x)), x 2 [ 1, 1].
1. Montrer que les fonctions Tn satisfont la formule de récurrence
⇢
T0 (x) = 1, T1 (x) = x,
Tn+1 (x) = 2xTn (x) Tn 1 (x).
2. Montrer ensuite que les polynômes Tn (x) sont orthogonaux par rapport à la fonction poids
(1 x2 ) 1/2 ,
8
Z 1 < ⇡, si n = m = 0,
dx
Tn (x Tm (x)) = ⇡/2, si n = m 6= 0,
1 (1 x2 ) 1/2 :
0, si n 6= m.
3. Montrer que Tn (x) est un polynôme de degré n dont le coefficient de xn est 2n 1 .
k⇡
4. Nous posons tn (x) = 21 n Tn (x), yk = cos( ), k = 0, ..., n, calculer tn (yk ).
n
5. Soient x1 , ..., xn , n points quelconques de [ 1, 1]. Nous posons wn (x) = (x x1 )...(x xn ).
Supposons par l’absurde que kwn k1 < ktn k1 . Montrer alors que
⇢ n
t (yk ) wn (yk ) > 0, si k est pair,
tn (yk ) wn (yk ) < 0, si k est impair.
6. En déduire que kwn k1 ktn k1 .
7. Application : Soit Ln 1 le polynôme de Lagrange de la fonction f définie pour tout x > 2
par f (x) = ln(x + 2) aux points racines du polynôme de Chebychev Tn . Déterminer n tel
que
max | ln(x + 2) Ln 1 (x)|  2 10
.
1x1

TD5 An 2015

Transféré par

Droits d'auteur :

Formats disponibles

TD5 An 2015

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

TD5 An 2015

Transféré par

Droits d'auteur :

Formats disponibles

Université Claude Bernard, Lyon 1 Licence Sciences & Technologies

43, boulevard du 11 novembre 1918 Spécialité : Mathématiques

Série d’exercices no 5/6

Exercice 1. Formule des Différences Divisées (Un Classique)

Pn (x) = a0 + a1 (x x0 ) + ... + an (x x0 )...(x xn 1 ),

permet de construire le polynôme Pn à l’aide d’une récurrence. En effet,

Autrement dit, connaissant Pn 1 , il suffit de calculer an pour connaître Pn .

a) Montrer que le polynôme d’interpolation de Lagrange de la fonction f aux points distincts

où f [.] désigne les différences divisées de f définies par

b) Montrer qu’il existe ⇠ 2 [a, b] tel que

Exercice 2. Convergence de l’interpolatio de Lagrange Soit Ln le polynôme d’interpolation de

3. Considérons toujours la fonction f

Exercice 3. Interpolation Polynomiale de Hermite

Hn (xi ) = f (xi ) et Hn0 (xi ) = f 0 (xi ), i = 0, ..., n.

Nous allons montrer le résultat suivant :

De plus, si f 2 C 2(n+1) ([a, b], R)

2. En déduire qu’il existe un unique polynôme Hn de degré 2n + 1 vérifiant les conditions

Exercice 4. Moindre carrés discrets

De plus, v ⇤ est la projection orthogonale de y sur l’espace F : v ⇤ = PF y est telle que

( 1, 1), (0, 0), (1, 1), (2, 2).

1. Tracer le nuage de points.

3. Montrer que ce problème admet une solution unique v ⇤ 2 F .

Exercice 5. Polynôme de Chebychev

Vous aimerez peut-être aussi