Mathématiques de l’ingénieur 2012-2013
Contrôle. Durée 1 h 45.
Nom et prénom :
Exercice 1 .
1. En rappelant la définition du problème de Cauchy pour une équation différentielle ordinaire, donner :
a. Une condition suffisante, pour que ce problème admette une solution unique.
b. Un exemple vérifiant cette condition (Justifiez).
2. Montrer qu’on peut introduire au schéma numérique d’Euler explicite de deux manières différentes.
3. Rappeler la définition générale d’un Schéma numérique à un pas, puis montrer que le Schéma d’Euler
implicite coïncide avec cette définition.
4. Considérons le problème de Cauchy : Trouver y ∈ C 1 ([ 0,1] , ℝ ) solution de :
y '(t ) = y (t ) 2 + g (t ), t ∈ [ 0,1]
y (0) = 2
Où g (t ) = −3 + 2 t − (1 − t ) ( 2 − t ) et dont la solution exacte est : y (t ) = (1 − t )( 2 − t )
2 2
En utilisant la méthode d’Euler explicite avec y0 = y (0) , calculer :
c. y1 , y2 , y3 et y4 , pour h = 0.1 , b. y1 et y 2 pour h = 0.2 c. y1 pour h = 0.4 puis remplir le tableau :
h t=0.2 Erreur globale GE/h
GE= y (tn ) − yn
0.4 y1
0.2 y2 =
0.1 y4 =
Exacte y (0.4) =
Expliquer les résultats obtenus.
Exercice 2 .
Considérons le problème de Cauchy : Trouver y ∈ C ( I , ℝ ) solution de :
1
y '(t ) = f (t , y (t )), ∀ t ∈ I = [ t0 , t0 + T ]
(2.1)
y (t0 ) = y0
Où f est une fonction Lipschitzienne de rapport L .
1. Soit ( an ) n∈ℕ une suite de réels positifs telle que an +1 ≤ A + an + B, ∀n ∈ ℕ.
An − 1
a. Montrer que an ≤ a0 An + B , ∀n ∈ ℕ.
A −1
b. Posons A = 1 + K , avec K > 0. Montrer que (1 + K )n ≤ e nK , puis en déduire que :
e nK − 1
an ≤ a0 enK + B , ∀n ∈ ℕ.
K
2. Soit en = y (tn ) − yn l’erreur à l’instant t = tn , en appliquant la méthode d’Euler explicite au problème (2.1).
h2
a. Montrer que | en +1 |≤ (1 + h L) | en | + M , où M est une constante à déterminer.
2
b. En déduire que | en |≤ C h , où C est une constante à déterminer.
3. Le Schéma d’Euler modifiée pour l’approximation du problème de Cauchy (2.1) est donnée par :
1 1
yn +1 = yn + h f tn + h, yn + h f (tn , yn ) .
2 2
Montrer que l’erreur de troncature Tn = h 2 F ( f , f x , f y , f xx , f xy , f yx ) (tn , y (tn )) + O (h3 ) .
Où F ( f , f x , f y , f xx , f xy , f yx ) (tn , y (tn )) est une fonction à déterminer.
Pr : Amrani Bonne chance