Corrigé ES2 CP2 2022

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

Annaba

E STI 2ème année CP-ST Analyse numérique

Corrigé de l’Examen Semestriel N◦ 2 (Juin 2022)

EXERCICE 1. [10 pts]


1. Calculer la spline cubique naturelle passant par les points (0, 1), (1, 12 ), et (2, 51 ).

1
2. Ces points sont ceux de la fonction f (t ) = 1+t 2.

π
R3
a. Donner une approximation de f (t )d t
0

b. Estimer l’erreur.

Corrigé

1. Puisqu’on dispose de 3 points, alors n = 2 et on doit donc calculer deux polynômes de degré 3 :
½
P 1 (t ) sur [0, 1]
P (t ) =
P 2 (t ) sur [1, 2]

Afin de déterminer ces deux polunômes, on doit calculer les cœfficients f 0" , f 1" et f 2" .

D’abord, comme il s’agit de la spline cubique naturelle, alors on a f 0" = f 2" = 0 et f 1" vérifie donc l’équa-
tion :
2 f 1" = 6 f [0, 1, 2]

1 pt

La table des différrence divisées est données par :

ti f (t i ) f [t i , t i +1 ] f [t i , t i +1 , t i +2 ]
0 1
− 12
1 1
1 2 10
3
− 10
1
2 5

1 pt
1 3
donc f [0, 1, 2] = 10 et f 1" = 10 .

1 pt

1 Février 2022
Annaba
E STI 2ème année CP-ST Analyse numérique

Maintenant, et puisque h = 1 on obtient :

Dans [0, 1] :
1 3 9 1 11
P 1 (t ) = t − (t − 1) + t = t 3 − t + 1
20 20 20 20

2 pt

Dans [1, 2] :
µ ¶
1 3 1 1 1 1 1 7
P 2 (t ) = − (t − 2) − − (t − 2) + (t − 1) = − (t − 2)3 − t +
20 2 20 5 20 4 10

2 pt

1
2. Les nœuds sont ceux de f (t ) = 1+t 2 :

a. Puisque P 1 (t ) est une interpolation de f sur [0, π3 ] ⊂ [0, 1] :


π π
Z3 Z1 Z3 ¸1 ¸ π3
t 4 11 2 (t − 2)4 t 2
· ·
7
f (t )d t ∼ P 1 (t )d t + P 2 (t )d t = − t +t + − − + t =
80 40 0 80 8 10 1
0 0 0

2 pt
¯π π ¯
¯R3 1
R 3
R ¯ ¯ π ¯
b. L’erreur est donnée par Er = ¯ f (t )d t − P 1 (t )d t + P 2 (t )d t ¯ = ¯[arctan(t )]03 − 1, 8857¯ = ...
¯ ¯ ¯ ¯
¯0 0 1 ¯

1 pt

2 Février 2022
Annaba
E STI 2ème année CP-ST Analyse numérique

2. On considère le problème suivant :


y 0 = t y,
½
t ∈ [0, 3]
y(0) = 2
a. On a ∀t ∈ [0, 3], ∀y 1 , y 2
¯ ¯ ¯ ¯ ¯ ¯
¯ f (t , y 1 ) − f (t , y 2 )¯ = t ¯ y 1 − y 2 ¯ ≤ 3 ¯ y 1 − y 2 ¯

0.5 pt
donc f est Lipschitzienne, alors le problème admet une solution unique.

0.5 pt
b. D’abord, on pose t n = 1 et puisque t n = t 0 + nh = n2 , alors n = 2 et l’approximation de y(1) est
donc y 2 .

0.5 pt
Maintenant, passons au calcul.
Le schéma est donné par :
( h ¡ ³ ¢´i
y n+1 = y n + h6 f t n , y n + 5 f t n + 3h 3h
¢ ¡
5 , y n + 5 f t ,
n ny
y0 donnée
en partant de y 0 = 2 on obtient
y 1 = 2, 25

0.5 pt
y 2 = 3, 20625

0.5 pt
Estimation de l’erreur :

La solution exacte du problème est donnée par :


t2
y(t ) = 2e 2

1 pt
1
Alors, la solution exacte est y(1) = 2e = 3, 2974 l’erreur est donc donnée par :
2

E = ¯ y(1) − y 1 ¯ = |3, 2974 − 3, 20625| = 0, 09115 = 9, 115 × 10−2


¯ ¯

1 pt

3 Février 2022

Vous aimerez peut-être aussi