Université Saad Dahlab Blida

Génie des Procédés 17 juin 2021

Méthodes Numériques
L2 S4 1 heure EMD

Ne rien prêter surtout l’effaceur.

Exercice 1 (13 pts)

On donne la fonction tabulée ci-dessous

i 0 1 2 3 4
x 0 ,5 1 1,5 2 2,5
y 1,0000 1,2500 1,6667 2,2000 2,8182

Effectuer les calculs à 10-4

1. Développer le tableau des différences finies. Quel est le degré du polynôme.
2. Donner un encadrement de P(1,39).
3. Calculer PN(1,39) par la première formule de Newton, en prenant toutes les différences
finies.
4. Déterminer yL(1,39) linéairement par la formule de Lagrange.
5. Discuter le résultat entre PN et PL. Qu’est-ce que vous suggérez (sans faire de calcul) ?
6. Déterminer la droite de régression y = a0+ a1x par la méthode des moindres carrés.

Exercice 2 (7pts)

a-Arrondir les nombres suivants à trois chiffres significatifs exacts. Déterminer l’erreur
absolue et indiquer la précision.

x x* ∆ ∆≤ 0,5.10-n précision
2,1514
0,16152
0,01204
-0,0015281

-392,85

b- Déterminer l’erreur relative du produit x.y ; si on donne x=369,7 , ∆x=0,05 et

y=0,0042131 ; ∆y= 0,5.10-7

1
Université Saad Dahlab Blida
Génie des Procédés 17 juin 2021
Méthodes Numériques
L2 S4 1 heure EMD

Corrigé

Méthodes d'interpolation de
Newton Tableau des différences finies
i x y Δy Δ2y Δ3y Δ4y
0,166
0 0,5000 1,0000 0,2500 7 -0,0500 0,0182
0,116
1 1,0000 1,2500 0,4167 7 -0,0318
0,084
2 1,5000 1,6667 0,5333 8
3 2,0000 2,2000 0,6182
4 2,5000 2,8182

On a un polynôme du quatrième degré

x0 x−xo
q=
h x h q
0,5 0,5 1,39 1,78

Encadrement de x : 1< x < 1,5 alors 1,25 < y(1,5) <1,6667

2. La première formule d’interpolation de Newton donne :

2 3 4
∆ y0 ∆ y0 ∆ y0
y= y 0+ q ∆ y 0 +q (q−1) +q (q−1)(q−2) +q (q−1)(q−2)(q−3)
2 6 24

y= 1,5635

PN (1,39) =1,5635

Cette valeur se trouve dans l’intervalle précédemment établi.

4. Interpolation linéaire par la formule de Lagrange.

i x y Formule d'interpolation linéaire de Lagrange

1 1,0000 1,2500 x PL(1,39)
2 1,5000 1,6667 1,3900 1,3417
.

(x−x 1)
y ( 1 ,39 ) =
( x −x2 )
(x ¿ ¿ 2−x 1) y 2+ ¿
(x ¿ ¿ 1−x 2) y 1 ¿

yL= 1,5750

2
Université Saad Dahlab Blida
Génie des Procédés 17 juin 2021
Méthodes Numériques
L2 S4 1 heure EMD

La méthode de Newton donne une valeur plus précise par ce qu’on a utilisé un polynôme du
quatrième degré. Cependant, l’interpolation linéaire par la formule de Lagrange est moins
précise pour cet exemple. Pour accroitre la précision de la méthode de Lagrange, il suffit
d’accroitre le nombre de points d’interpolation.

En d’autres mots, il faut augmenter le degré du polynôme, ce qui va engendrer un calcul

supplémentaire.

4. Détermination de la droite de régression par la méthode des moindres carrés

i x y x^2 xy
0 0,5000 1,0000 0,2500 0,5000
1 1,0000 1,2500 1,0000 1,2500
2 1,5000 1,6667 2,2500 2,5000
3 2,0000 2,2000 4,0000 4,4000
4 2,5000 2,8182 6,2500 7,0455
somme
5 7,5000 8,9348 13,7500 15,6955

Sous forme matricielle

( ) ( )
4 4
5 ∑ xi ∑ yi
4 4
0
( )
ao =
a1 4
0

∑ xi ∑ xi2 ∑ xi yi
0 0 0

(75, 5 7 ,5 ao
13 ,75 a1)( ) (
=
8,9348
15,6955 )
a 0=0,4111 a1=0,9173
Ex2

Le nombre exact x, et x* c’est le nombre arrondi à 3 chiffres significatifs exacts.

x x* ∆ ∆≤ .0, 5.10-n précision

2,1514 2,15 0,0014 0,510-2 10-2
0,16152 0,162 0,0005 0,5.10-3 10-3
0,01204 0,0120 0,00004 0,5.10-4 10-4
-0,0015281 -0,00153 0,0000019 0,5.10-5 10-5

0,5.100 1 unité
-392,85 -393 0,15
625,55 626 0,45 0,5.100 1 unité

3
Université Saad Dahlab Blida
Génie des Procédés 17 juin 2021
Méthodes Numériques
L2 S4 1 heure EMD

b-erreur relative d’un produit

z=x . y

∆z ∆x ∆ y
= +
z x y

Erreur relative du produit

∆x ∆y ∆z ∆z
x y ∆x ∆y %
x y z z
369,5 0,0042131 0,05 5.10-8 1,3510-4 1,1910-5 1,4710-4 0,015%

4
Université Saad Dahlab Blida
Génie des Procédés 17 juin 2021
Méthodes Numériques
L2 S4 1 heure EMD

Exercice 3

Tableau 1. Table des différences finies.

I x Y Δy Δ2y Δ3y Δ4y

5
Université Saad Dahlab Blida
Génie des Procédés 17 juin 2021
Méthodes Numériques
L2 S4 1 heure EMD

0 3,5 144,2500 6,4480 0,1980 0,0030 0,0000

1 3,55 150,6980 6,6460 0,2010 0,0030
2 3,6 157,3440 6,8470 0,2040
3 3,65 164,1910 7,0510
4 3,7 171,2420

Δ3y=0,0030 et Δ4y = 0,0000 ce qui implique que le polynôme est du troisième degré.

1ere formule de Newton

1.h x q=(x- P3, 1er

xo)/h formule de
Newton
0,05 3,52 0,4 P3(3,52) 146,8056320

3,68 3,6 P3(3,68) 168,3969280

P1 3,52 P1(3,52) 146,8292000

P2 3,52 P2(3,52) 146,8054400

h =0,05, x0= q=(x-xo)/h

3,5 ; x= 3,52
0,4
P1(3,52) 146,8292000

P2(3,52) 146,8054400

P3(3,52) 146,8056320

2 eme formule de Newton

X q=(x- P(x)
xn)/h
3,52 -3,60 P3(3,52) 146,8056320
3,68 -0,40 P3(3,68) 168,3969280
3,68 -0,40 P1(3,68) 168,4216000
3,68 -0,40 P2(3,68) 168,3971200

Calcul d’erreurs de la première formule de Newton

6
Université Saad Dahlab Blida
Génie des Procédés 17 juin 2021
Méthodes Numériques
L2 S4 1 heure EMD

erreurs absolues erreurs Nombre de chiffres

x=3,52 relatives significatifs exact s
abs(P3(x)- 2,36E-02 1,61E-04 4
P1(x)) 0,2410-1

abs(P2(x)- 2,38E-02 1,62E-04 4

P1(x)) 0,2410-1
abs(P3(x)- 1,92E-04 1,31E-06 6
P2(x)) 0,1910-3

Abs(P3(x)-y) 0 0 Exact
1,61E-04 =1,6110-4=0,1610-3

A titre illustratif on va comparer l’augmentation de la précision en augmentant le degré du

polynôme.

Calcul de l’ordre de grandeur de la précision x=3,52

erreurs absolues %
D1=abs(P3(x)- 2,36E-02 D1/D2 1,01
P1(x))
D2=abs(P2(x)- 2,38E-02 D2/D3 123,75
P1(x))
D3=abs(P3(x)- 1,92E-04
P2(x))
abs(P3(x)-f(x)) 0,000000
0

Conclusion

On remarque qu’en augmentant le degré du polynôme la précision sur la valeur approchée

augmente. C’est ce qu’il faut dire aux étudiants.

2 formule de Newton
X=3,68
erreurs absolues erreurs relatives Nombre de chiffres
significatifs exacts
abs(P3(x)- 2,47E-02 1,47E- 4
P1(x)) 04<0,510-3
abs(P2(x)- 2,45E-02 1,45E- 4
P1(x)) 04<0,510-3
abs(P3(x)- 1,92E-04 1,14E- 6
P2(x)) 06<0,510-5
abs(P3-y) 0,000000 0 Exact
0

7
Université Saad Dahlab Blida
Génie des Procédés 17 juin 2021
Méthodes Numériques
L2 S4 1 heure EMD

Calcul de l’ordre de grandeur de la précision , x=3,68

erreurs absolues %
D1=abs(P3(x)- D1/D2
P1(x)) 0,024672 1,0078431
D2=abs(P2(x)- D2/D3
P1(x)) 0,02448 127,5000000
D3=abs(P3(x)-
P2(x)) 0,000192
abs(P3(x)-f(x)) 21,59129
6

Le degré polynomial ascendant améliore la

précision de la valeur approchée.
Question
Est-ce que l’interpolation linéaire est acceptable en
pratique.
Réponse
Oui dans une certaine mesure.

y=f(x)
175.0000
170.0000
f(x) = 4.00000000005 x³ − 3.00000000052 x² + 1.00000000187 x + 5.99999999774
165.0000
160.0000
155.0000 y
150.0000
y

Polynomial (y)
145.0000
140.0000
135.0000
130.0000
3.45 3.5 3.55 3.6 3.65 3.7 3.75
x