Correction de la série d’exercice N 2 : Interpolation polynômiale 2eme Année
Univercité sétif 1 - 2019/2020
Exercice 1
p
Soit f (x) = 3
x
1- Le polynôme d’interpolation
a/ De Lagrange
xi 1 8 27
yi = f (xi ) 1 2 3
Le polynôme sera de degré 2 avec :
Q
2
(x xj )
X
2
j=0; j6=i
P2L (x) = yi Li (x) et Li (x) = ,
Q
2
i=0 (xi xj )
j=0; j6=i
(x x1 )(x x2 ) (x 8)(x 27) 2
L0 (x) = (x0 x1 )(x0 x2 )
= (1 8)(1 27)
= x 35x+216
182
:
(x x0 )(x x2 ) (x 1):(x 27) 2
L1 (x) = (x1 x0 )(x1 x2 )
= (8 1):(8 27)
= x +28x 133
27
:
(x x0 )(x x1 ) (x 1):(x 8) 2
L2 (x) = (x2 x0 )(x2 x1 )
= (27 1)(27 8)
= x 494 9x+8
:
On trouve P2L (x) = y0 L0 (x) + y1 L1 (x) + y2 L2 (x) ; ce qui donne :
1 1 1 35 28 9 216 27 8
P2L (x) = + x2 + + x+ + :
182 133 494 182 133 494 182 133 494
6 2 43 1434
= x + x+ :
1729 247 1729
b/ Newton en fonction des di¤érences divisées (D:D)
1
8
>
> X 2
>
> N
P2 (x) = i i (x) telle que
>
>
>
> i=0
<
On a i (x) = 1
>
>
>
>
> i (x) = (x x0 ) = (x 1)
>
>
>
:
i (x) = (x x0 ) (x x1 ) = (x 1) (x 8) = x2 9x + 8:
On forme le tableau des (D:D) :
xi 0 1 2
x0 = 0 f (x0 ) = 1
2 1 1
8 1
= 7
1 1
x1 = 2 f (x2 ) = 2 19 7 6
= - 1729
27 1
3 2 1
27 8
= 19
x2 = 6 f (x3 ) = 3
P2N (x) = 1 + 17 (x 1) 6
1729
(x2 9x + 8) ;
Alors on trouve
6 2 43 1434
P2N (x) = x + x+ :
1729 247 1729
p
3
C/ L’erreur de l’interpolation commise pour la valeur 20:
a/ L’erreur exacte
Eexact = jf (20) P (20)j = j2; 7144 2; 9230j = 0; 2086:
b/L’erreur théorique
M (3) Q
2
ET heo = (20 xi ) avec M3 = max f (3) (x) ;
(3)! i=0 [x0 ;x2 ]=[1;27]
Alors
p 1
f (x) = 3
x = x3 ;
0 2
f (x) = 13 x 3 ;
00 5
2
f (x) = 9
x 3 ;
(3) 8
10
f (x) = 27
x 3
11 (3)
80
f (4) (x) = 81
x 3 f (4) (x) <0 alors f (x) est décroissante sur [1; 27] ;
2
donc M3 = f (3) (1) = 10
27
:
Q
2
(20 xi ) = j(20 x0 ) (20 x1 ) (20 x2 )j = j(20 1) (20 8) (20 27)j = 1596;
i=0
10
alors ET heo = 27
6
1596 = 98; 51:
toujour Eexact < ET heo .
Exercice 2
On a f (x) = x cos ( x)
1
xi x0 = 0 x1 = 2
x2 = 1
1
yi = f (xi ) 1 2
2
Le polynôme d’interpolation sera de degré 2
1-Le polynôme d’interpolation Newton en fonction des di¤érences …nies ([Link]P )
On construit le tableau des ([Link]P )
xi 40 yi = yi 4yi 42 yi
0 1
1 3
2
+1= 2
1 1 3 3
2 2 2 2
= 0
1 3
2 2
= 2
1 2
On a h = xi+1 xi = 21 ; et d’après la formule de la fonction de Newton on obtient :
4y0 4 2 y0
P2N (x) = y0 + h:1!
(x x0 ) + h2 :2!
(x x0 ) (x x1 )
3 3
1
= 1+ 1
2
:1!
(x 0) + 2
1 2
(x 0) x 2
2 ( )
2
2!
alors P2N (x) = 3x 1:
1
3- Calcul de l’erreur exacte et théorique point x = 4
:
1
on résoudre l’équation sin 10
= sin 2
x on trouve x = 5
donc Eext = f 14 P2 1
4
= j 0; 4571 + 0; 25j = 0; 2071:
M3
Q
2
1
et Etheo = (3)! 4
xi ;
i=0
3
avec M3 = max jf 3 (x)j
x2[0;1]
f (x) = x cos ( x)
0
f (x) = 1 + sin ( x)
f " (x) = 2
cos ( x)
000 3
f (x) = sin ( x) :
1
M3 = max f 3 2
= 3
x2[0;1]
3 1 1 1 1
alors Et = 3! 4
0 4 2 4
1 = [Link]
Exercice 3
1-Le polynome de Lagrange
xi 0 1 2 3
yi = f (xi ) 2 1 2 3
Le polynôme sera de degré 2 avec :
Q
3
(x xj )
X
3
j=0; j6=i
P3L (x) = yi Li (x) et Li (x) = ,
Q
3
i=0 (xi xj )
j=0; j6=i
(x x1 )(x x2 )(x x3 ) 3 2
L0 (x) = (x0 x1 )(x0 x2 )(x x3 )
= (x 1)(x 2)(x 3)
(0 1)(0 2)(0 3)
= x 6x +11x 6
6
:
(x x0 )(x x2 )(x x3 ) 3 2
L1 (x) = (x1 x0 )(x1 x2 )(x x3 )
= x:(x 2)(x 3)
1:(1 2)(1 3)
= x 5x2 +6x :
(x x0 )(x x1 )(x x1 ) 3 2
L2 (x) = (x2 x0 )(x2 x1 )(x x1 )
= x:(x 1)(x 3)
2(2 1)(2 3)
= x 4x2 +3x
(x x0 )(x x1 )(x x2) 3 2
L3 (x) = (x3 x0 )(x3 x1 )(x3 x2 )
= x:(x 1)(x 2)
3(3 1)(3 2)
= x 3x6 +2x
On trouve P3L (x) = y0 L0 (x) + y1 L1 (x) + y2 L2 (x) + y3 L3 (x) ; ce qui donne :
1 3 8
P3L (x) = 3
x + 2x2 3
x + 2:
1-Le polynôme d’interpolation Newton en fonction des di¤érences …nies ([Link]P )
On construit le tableau des ([Link]P )
4
xi yi 4yi 42 yi 43 yi
0 2
1–2= -1
1 1 1+1= 2
3-2=1 0-2= -2
2 2 1-1=0
3-2=1
3 3
On a h = xi+1 xi = 3 2 = 1; et d’après la formule de la fonction de Newton on
obtient :
4y0 4 2 y0 4 3 y0
P3N (x) = y0 + h:1!
(x x0 ) + h2 :2!
(x x0 ) (x x1 ) + h3 :3!
(x x0 ) (x x1 ) (x x2 )
1 1 4 3 y0
=2+ 1:1!
(x 0) + (1)2 2!
(x 0) (x 1) + 13 :3!
(x 0) (x 1) (x 2)
Alors
1 3 8
P3L (x) = 3
x + 2x2 3
x + 2:
