L3 : MAP, MATHINFO M325 Calcul Scientifique II

Université Paris–Sud Centre d’Orsay

AN-EDP Année 2013 − 2014

c J.-B.A.K. <[email protected]>

Devoir : durée 1 h
Attention :
— Vous devez créer un répertoire préfixé par votre nom et suffixé par M325, dans lequel se trouverons vos scripts.
— Aux termes de l’épreuve vous devez transmettre le répertoire zippé à l’adresse qui vous sera indiquée.
— Les réponses doivent figurer sous forme de commentaires dans vos scripts.

Dans les exercices qui vont suivre, on considère les schémas numérique pour la résolution de l’équation différentielle

x′ (t) = f (t, x(t)), t ∈]t0 , t0 + T [, x(t0 ) = x0 . (1)

r
Où f ∈ C ([t0 , t0 + T ]; R), (l’entier r, ainsi que les réels t0 et T seront précisés si nécessaire dans chaque exercice).
On se donne alors un entier N fixé, et on pose h = T /N et tn = t0 + nh (0 ≤ n ≤ N ). On désigne par xn une valeur
approchée de x(tn ) pour tout n ∈ {0, . . . N }.

Exercice-1 : Analyse de précision du schéma

Dans cet exercice, on suppose le problème (1) résolu par trois schémas distincts (S1), (S2), et (S3) tous à un pas.

On ne dispose d’aucune information sur ces schémas, si ce n’est les courbes de convergence données ici par l’évolution,
en fonction du pas de discrétisation (h), de l’erreur (Eh ) à l’instant final entre la solution exacte et la solution approchée.
Ces courbes sont fournies pour deux types d’échelles dans les Figure 1 et Figure 2. La Figure 2 contient en plus certaines
droites dont les pentes sont indiquées dans la légende.

A la lueur de ces graphiques, et pour le problème considéré, classer ces trois suivant chacune des priorités suivantes :
1. Meilleur ordre de convergence (on précisera l’ordre de convergence potentiel de chaque schéma).
2. Meilleure qualité de l’erreur à pas h fixé.
3. Plus petite constante de consistance (sous l’hypothèse des constantes de stabilité identiques).
Conseiller alors un utilisateur sur le choix d’un de ces schémas.
0.14 −2
S1
S2
S3
−4
0.12

−6
0.1

−8
0.08
log( Eh )

S1
Eh

−10 S2
S3
0.06 pente 1
−12 pente 2
pente 3
0.04
−14

0.02
−16

0 −18
0 0.05 0.1 0.15 0.2 0.25 −4.5 −4 −3.5 −3 −2.5 −2 −1.5 −1
h log( h )

F IGURE 1 – Erreur Eh à l’échelle normale. F IGURE 2 – Erreur Eh à l’échelle logarithmique.

Exercice-2 : Erreur locale

On considère ici le problème (1) avec les paramètres suivants :

1
f (t, y) = 4t y 2 , t0 = 0, T = 2, x(t0 ) = 1. (2)

Pour un pas h et une donnée initiale x0 = 1, on s’intéresse à l’approximation x1 de x(h), produite par trois schémas
(S1),(S2) et (S3), donnée pas :
(
x0 = 1,
 
x0 = 1, x0 = 1,
(S1) (S2) (S3)
 
x1 = x0 + hf (t0 , x0 ) x1 = x0 + hf t0 + h2 , x0 + h2 f (t0 , x0 ) x1 = 1 − 0.3h + 2h2 + h4

Q-2-1 : Donner en fonction de h, l’expression de x1 fournie par les schémas (S1) et (S2).

i i
Q-2-2 : Représenter sur le même graphique les courbes log(h) 7→ log(|E1h |), i = 1, 2, 3, où pour i = 1, 2, 3, E1h
désigne l’erreur générée pas le schéma (Si) sur le calcul de x1 . On admettra que la solution exacte du problème (1),
pour les données (2) est x(t) = (1 + t2 )2 .

Q-2-3 : Commenter les observations et dire si celles-ci étaient prévisibles ?

Exercice-3 : Applications

On considère encore l’équation différentielle (1) avec les expressions (2).

Pour sa discrétisation, on considère le schéma numérique suivant :


xn+2 = xn + 2hfn+1 , n = 0, . . . , N − 2,
(S) (3)
x0 = 1; x1 (à déterminer).
T
où on a posé fj = f (tj , xj ), tj = jh, j = 0, . . . , N, avec h = N.

Q-3-1 : Ecrire une fonction Matlab, monSchemaDeuxPas qui calcule la suite des solutions produites par le
schéma.

function [t,x] = monSchemaDeuxPas(f, t0, tf, x0, x1, h)

% ENTREE :
% f -> second-membre de l’edo ;
% t0 -> instant initial
% tf -> temps final;
% x0 -> donnée initiale;
% x1 -> approximation de x(t0+h);
% h-> pas;
% SORTIE :
% t-> liste des instants générés [t0,t1,...,tN]
% x-> liste des solutions générées aux instants t

Q-3-2 : Pour h = 0.1, représenter sur un même graphique les solutions obtenues pour les valeurs de x1 obtenues
par les schémas (S1), (S2) et (S3).

Q-3-3 : En faisant varier le pas h comme h = 2−i avec i = 2, 3, 4, 5, 6, et pour chaque valeur de x1 obtenue par
les schéma (S1 ), (S2) et (S3), on a calculé l’erreur entre la solution exacte et la solution fournie par le schéma (S) à
l’instant final : Eh = |x(tN ) − xN |. Les courbes d’évolution de Eh en fonction de h obtenues ont été représentées à
l’échelle logarithmique sur une même figure. On a obtenu la Figure 3.

1. Commenter ce graphique (donner les ordre de convergence) et dire si ces résultats étaient prévisibles.
On rappelle que pour le schéma (S) une estimation d’erreur est donnée par :
 
max |x(tn ) − xn | ≤ K1 max{|x(t0 ) − x0 |, |x(t1 ) − x1 |} + K2 h2 ,
0≤n≤N

où K1 et K2 sont des constantes positives indépendantes de h.

2. Fournir un script Matlab qui réalise cette expérience et comparer le graphique généré à celui de la Figure 3.
 1
S1
S2
S3
0 pente 1
pente 2
pente 3

−1

−2
log( Eh )

−3

−4

−5

−6
−4.5 −4 −3.5 −3 −2.5 −2 −1.5 −1
log( h )

F IGURE 3 – Courbes à l’échelle logarithme, de l’évolution de l’erreur Eh = |x(tN ) − xN | entre la solution exacte et
la solution produite par le schéma (S), pour différentes valeurs de x1 . Sur la légende, Si, i = 1, 2, 3 désigne la courbe
associée à la valeur de x1 issue du schéma (Si). On a aussi joint des droites de pentes respectives 1, 2 et 3. Les courbes
sont obtenues pour les données (2) avec h = 2−k , k = 2, 3, 4, 5, 6.