4correction Interro 01-1

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

Université Batna 2 -- Faculté MI -- Département de Mathématiques

Interrogation(1) en Analyse Numérique 1 -- 2ème année Licence.


Mercredi 28/11/ 2018
Exercice :(06 pts) TD
𝜋
1) Calculer l’interpolant P2 de la fonction 𝑓(𝑥) = cos(𝑥) en les points 𝑥𝑖 = 2 𝑖, avec 𝑖 = 0,1,2.
𝜋
2) Soit P3 l’interpolant de la même fonction en les points 𝑥𝑖 = 𝑖, avec 𝑖 = 0,1,2,3.
2
a) Y’a-t-il une relation entre P2 et P3 ? Est-elle généralisable ?
b) Calculer P3.
𝜋
c) Soit Et(x) l’erreur théorique au point x. estimer Et(𝜋) et Et( ).
3
---------------------------------------------------------------------------------------------
Correction interrogation 01

1) On choisit la méthode de NEWTON car elle permet d’utiliser les calculs de la question 1) pour
répondre à la question 2).
Calculons l’interpolant de Newton 𝑝2 . On commence par construire le tableau des différences divisées :

i 𝑥𝑖 𝑓(𝑥𝑖 ) Différence d’ordre 1 Différence d’ordre 2


0 0 1
1 𝜋 0 2
2 𝑓[𝑥0 , 𝑥1 ] = −
𝜋
2 𝜋 -1 2 𝑓[𝑥0 , 𝑥1 , 𝑥2 ] =0
𝑓[𝑥1 , 𝑥2 ] = −
𝜋
(𝟎. 𝟓𝟎𝒑𝒕) (𝟎. 𝟓𝟎𝒑𝒕) + (𝟎. 𝟓𝟎𝒑𝒕) (𝟎. 𝟓𝟎𝒑𝒕)
On a alors, 𝑝2 (𝑥) = 𝑓(𝑥0 ) + 𝑓[𝑥0 , 𝑥1 ](𝑥 − 𝑥0 ) + 𝑓[𝑥0 , 𝑥1 , 𝑥2 ](𝑥 − 𝑥0 )(𝑥 − 𝑥1 ), (0. 5pt)
2
Donc 𝑝2 (𝑥) = 1 − 𝜋 𝑥. (𝟎. 𝟓𝒑𝒕)

2) a) Oui, il ya une relation entre P2 et P3 : 𝑝3 (𝑥) = 𝑝2 (𝑥) + 𝐷3 . 𝑁3 (𝑥) dans la formule de Newton
(𝟎. 𝟓𝒑𝒕) et elle est généralisable : 𝑝𝑛+1 (𝑥) = 𝑝𝑛 (𝑥) + 𝐷𝑛+1 . 𝑁𝑛+1 (𝑥) pour tout 𝑛 ≥ 0 (𝟎. 𝟓𝒑𝒕)

b) Pour calculer l’interpolant 𝑝3 , il suffit de calculer une différence divisée en plus, i.e. ajouter une
ligne au tableau précédant. Alors

i 𝑥𝑖 𝑓(𝑥𝑖 ) Différence d’ordre 1 Différence d’ordre 2 Différence d’ordre 3


0 0 1
1 𝜋 0 2
2 𝑓[𝑥0 , 𝑥1 ] = −
𝜋
2 𝜋 -1 2 𝑓[𝑥0 , 𝑥1 , 𝑥2 ] =0
𝑓[𝑥1 , 𝑥2 ] = −
𝜋
3 3𝜋 0 2 4 8
𝑓[𝑥2 , 𝑥3 ] = 𝑓[𝑥1 , 𝑥2 , 𝑥3 ] = 2 𝑓[𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 ] =
2 𝜋 𝜋 3𝜋 3
(𝟎. 𝟐𝟓𝒑𝒕) (𝟎. 𝟐𝟓𝒑𝒕) (𝟎. 𝟐𝟓𝒑𝒕)
𝑝3 (𝑥) = 𝑝2 (𝑥) + 𝑓[𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 ](𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )
2 4 8
Donc, 𝑝3 (𝑥) = 1 − 3𝜋 𝑥 − 𝜋2 𝑥 2 + 3𝜋3 𝑥 3 . (𝟎. 𝟐𝟓𝒑𝒕).
d) Et(𝜋) = 0 car 𝑥2 = 𝜋 est un neoud. (𝟎. 𝟓𝒑𝒕).
𝜋 𝜋 𝜋 𝜋 𝜋 1 2 𝜋 4 𝜋 8 𝜋 11
Et( ) = 𝑓 ( ) - 𝑝3 ( ) = cos (3) − 𝑝3 (3) =
3 3 3 2
- (1 −
3𝜋 3
− 𝜋2 ( 3 )2 + 3𝜋3 ( 3 )3) = 162 . (𝟎. 𝟓𝒑𝒕).

Vous aimerez peut-être aussi