Num_TD5

UFHB Année
2022-2023
Maths-Info
TD5 : Analyse Numérique des équations diérentielles
ordinaires
Exercice 1.
Soit le problème de Cauchy suivant :
(
y ′′ (t) + 2y ′ (t) = 2y(t), t ∈ [0, 1]
(1)
y(0) = 1 , y ′ (0) = 2.
1. Ecrire cette équation sous la forme d'un système diérentiel de deux équations d'ordre un.
2. On note h le pas numérique. On applique la méthode d'Euleur à ce système, écrire l'algorithme
correspondant
3. Même question pour la méthode de Runge Kutta d'ordre 2.
4. Même question pour la méthode de Runge Kutta d'ordre 4.
5. Donner la solution exacte de 1, ainsi que y(t) et y (t) pour t = 0, 2.
′
6. Comparer les résultats précédents à ceux numériquement obtenus grâce au schéma d'Euler
en prenant h = 0, 1.
Exercice 2.
Soit f : [a, b] × R → R uniformément L-Lipschitz par rapport la deuxième variable. Soit le
problème de Cauchy : (
y ′ (t) = f (t, y(t))
y(x) = ya .
Soit N > 0 , on pose h = et t = a + ih pour 0 ≤ i ≤ N .
b−a
i
1. considère la méthode d'Euler

N
yn+1 = yn + hf (tn , yn ).
2. Montrer que la m´ethode est consistante et stable.
3. On considère le schéma de Runge-Kutta 2 (Euler modié) suivant
h h
yn+1 = yn + hf (tn + , yn + f (tn , yn )).
2 2
consistance? ordre? convergence?
4. On considère le schéma de Taylor à un pas suivant
yn+1 = yn + hΦ(tn , yn , h)
avec Φ(t , y , h) = f (t , y ) + g (t , y ) et g (t , y ) = (t , y ) + f (t , y ) (t , y ).
h ∂f ∂f
Montrer que ce schéma est consistant et donner l'ordre de consistance. Montrer que si g
n n n n 2 2 n n 2 n n ∂t n n n n ∂y n n
et L -Lipschitz le schéma est consistant et convergent. Construire des schémas d'ordres

2
supérieurs.
2
1
5. Appliquer les schémas précédents à l'équation
(
y ′ (t) = y(t) sin(t), t ∈ [0, Π]
y(0) = 1 .
Exercice 3.
On considère l'équation diérentielle :
y ′ (t) = f (t, y(t)), y(0) = y0 ,
où f : R → R est une fonction de classe C et y 2

∈ R. On suppose qu'il existe une constante C > 0
telle que, pour tout t ∈ R et y ∈ R :
0
∂f
(t, y) ≤ C.
∂y
Pour résoudre numériquement cette équation diérentielle, on considère une famille de méthodes à
un pas dénie de la manière suivante : pour n ∈ N, t = nh où h est le pas de temps; on considère
des approximation y de y(t ) que l'on calcule par la formule de récurrence :
n
n n
yn+1 = yn + hΦ(tn , yn , h)
où Φ est une fonction de la forme :

Φ(t, z, h) = a1 f (t, z) + a2 f (t + p1 h, z + p2 hf (t, z)),
avec a ,a ,p ,p ∈ R.
1 2 1 2
1. Rappelez très brièvement pourquoi cette condition implique l'existence et l'unicité de la

solution y(.) de l'équation diérentielle qui est dénie sur [0; +∞[.
2. Donner des conditions sur les paramètres a ,a ,p ,p pour que cette méthode soit stable et
consistante.
1 2 1 2
3. En déduire pour quelles valeurs de a ,a ,p ,p , cette méthode est convergente.

1 2 1 2
4. Donner des conditions sur les paramètres a ,a ,p ,p pour que cette méthode soit d'ordre 2.
1 2 1 2
5. En déduire la forme générale des fonctions Φ(dénies comme ci-dessus) qui donnent une
méthode d'ordre 2. (On pourra montrer qu'il s'agit d'une famille à un paramètre et utiliser
par exemple le paramètre α = a .) 1
Exercice 4.
On étudie la méthode numérique (M) de résolution de l'équation diérentielle y (t) = f (t, y(t)) ′
dénie par
(
yn+1 = yn + hn Φ(tn , yn , hn )
Φ(t, y, h) = αf (t, y) + βf (t + h2 , y + h2 f (t, y)) + γf (t + h, y + hf (t, y))
où α, β, γ sont des réels compris entre 0 et 1.

1. Pour quelles valeurs du triplet (α, β, γ) retrouve-t-on
la méthode d'Euler?
la méthode du point milieu?
la méthode de Heun?
2
2. Dans cette question et la suivante, on supposera que la fonction f (t, y) est de classe C ∞
sur [t , t + T ] × R, et k− lipschitzienne en y. Pour quelles valeurs de (α, β, γ) la méthode

proposée est-elle stable?
0 0
3. Quelles relations doivent satisfaire (α, β, γ) pour que la méthode soit consistante? conver-
gente? d'ordre ≥ 1 ? d'ordre ≥ 2 ?
4. Donner des conditions sur les paramètres a ,a ,p ,p pour que cette méthode soit d'ordre 2.
1 2 1 2
5. La méthode (M) peut-elle être d'ordre supérieur?

Exercice 5.
On considère la méthode de Runge-Kutta dénie par le tableau
0 0 0 0
1/4 1 0 0 .
3/4 -9/20 6/5 0
1 1/9 1/3 5/9
On considère l'équation y = f (t, y) = t + y + 1.
′
1. Résoudre cette équation et calculer f (0, 0) pour tout n ∈ N.

[n]
2. Que peut-on dire de (0, 0, 0) ?

∂nΦ
3. Déterminer l'ordre de cette méthode.

∂hn

Num_TD5

Transféré par

Droits d'auteur :

Formats disponibles

Num_TD5

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Num_TD5

Transféré par

Droits d'auteur :

Formats disponibles

UFHB Année

1. considère la méthode d'Euler

et L -Lipschitz le schéma est consistant et convergent. Construire des schémas d'ordres

où f : R → R est une fonction de classe C et y 2

où Φ est une fonction de la forme :

1. Rappelez très brièvement pourquoi cette condition implique l'existence et l'unicité de la

3. En déduire pour quelles valeurs de a ,a ,p ,p , cette méthode est convergente.

où α, β, γ sont des réels compris entre 0 et 1.

sur [t , t + T ] × R, et k− lipschitzienne en y. Pour quelles valeurs de (α, β, γ) la méthode

5. La méthode (M) peut-elle être d'ordre supérieur?

1. Résoudre cette équation et calculer f (0, 0) pour tout n ∈ N.

2. Que peut-on dire de (0, 0, 0) ?

3. Déterminer l'ordre de cette méthode.

Vous aimerez peut-être aussi