Département de Mathématiques
Filière : m.a - Semestre : S5 Année Universitaire 2019/2020
Module : Analyse Numérique H. Allouche
Examen
Session ordinaire
NB : Il sera tenu compte de la clarté de la rédaction et de la propreté de la copie
Exercice :
La méthode d’Adams-Bashforth à trois (3) pas pour la résolution numérique du problème
de Cauchy
′
y (t) = F (t, y(t)) t ∈ I = [t0 , t0 + T ]
(23)
y(t0 ) = µ
est la suivante :
h
yn+1 = yn + (23fn − 16fn−1 + αfn−2 )
12
1. Que désigne-t-on par fi (ou c’est quoi fi ) ? est elle une méthode explicite ou implicite ?
comment se fait le démarrage du calcul ?
23 −16 α
2. Quelle est l’idée derrière le calcul des coefficients , et
12 12 12
3. Admettons que le problème à résoudre est
′
y (t) = 1 t ∈ I = [t0 , t0 + T ]
y(t0 ) = µ
k−1
X
et on veut que la méthode d’Adams-Bashforth à k pas yn+1 = yn + h bi,k fn−i
i=0
fournisse la solution exacte, quelle condition doivent vérifier les coefficients bi,k , i =
0, · · · k − 1 pour que cela se réalise ?
4. Déterminer la valeur du coefficient α.
5. On pose
1 1
εn = (y(tn+1 ) − y(tn ))− (23F (tn , y(tn )) − 16F (tn−1 , y(tn−1 )) + αF (tn−2 , y(tn−2 )))
h 12
(a) Qu’appelle-t-on εn ?
1 tn+1 ′
Z
(b) Sachant que εn peut être écrite sous la forme : εn = (y (t)−p(t))dt où p(t)
h tn
est un polynôme de degré ≤ 3. C’est quoi le polynôme p(t), écrire son expression.
124
(c) Sachant que si Pk est le polynôme d’interpolation d’une fonction f aux abscisses
x0 , x1 , · · · , xk et si f est de classe C k+1 sur J (un intervalle contenant les xi ) alors
l’erreur d’interpolation en un point x ∈ J est :
(x − x0 ) · · · (x − xk ) (k+1)
Ek (x) = f (x) − Pk (x) = f (ξ) où ξ ∈ J.
(k + 1)!
Montrer que si la solution y(t) est de classe C τ sur I alors il existe une constante
K telle que |εn | ≤ Kh3 et préciser la valeur ”minimale” de τ
(d) Si la méthode est stable que peut-on déduire ?
La méthode d’Adams-Moulton à 3 pas pour la résolution du problème (23) est
h
yn+1 = yn + (5fn+1 + 8fn − fn−1 )
12
6. C’est quoi fn+1 ? est-elle une méthode explicite ou implicite ?
7. Montrer que pour h suffisamment petit, la méthode est bien définie.
8. Expliquer comment se déroule le calcul de la solution approchée à l’aide de cette
méthode.
9. On considère la méthode de prédiction-correction où le schéma prédicteur est celui
d’Adams-Bashforth à 3 pas et le schéma correcteur est celui d’Adams-Moulton à 3
pas.
(l)
Comment la suite yn+1 que génère cette méthode est-elle définie ?
l∈ N
(0)
N.B : N’oublier pas de préciser la valeur initiale yn+1 de cette suite.
(l)
10. Montrer que pour h suffisamment petit, la suite yn+1 converge.
N
l∈
11. Compléter l’algorithme qui retourne la solution approchée du problème (23) par la
méthode de prédiction-correction que vous avez décrit dans les questions précédentes
125
fonction (y( ), t( ) : réel) ←− Pred Corr(F( , ) :fonction, Yi( ), a, b : réel, N : entier, eps : réel)
début
var h,f0,f1,f2,f3,x0,x1 : réel
n : entier
h ←− (b-a)/N
t(0) ←− a
t(1) ←− t(0)+h
t(2) ←− t(1)+h
y(0) ←− Yi(0)
y(1) ←− Yi(1)
y(2) ←− Yi(2)
f0 ←− F(t(0),y(0))
f1 ←− · · ·
f2 ←− · · ·
pour n ←− 3 à N faire
t(n) ←− t(n-1)+h
x0 ←− y(n-1)+h*(23*f2-16*f1+α*f0)/12
f3 ←− F(t(n),x0)
x1 ←− · · ·
tant que (abs(x1-x0)≥ eps) faire
x0 ←− · · ·
f3 ←− · · ·
x1 ←− · · ·
fin tant que
y(n) ←− x1
f0 ←− · · ·
f1 ←− · · ·
f2 ←− · · ·
fin pour
retourne (y,t)
fin
126