Analyse Numérique - 2022-2023: Exercice 1

Université d’Angers
Département de Mathématiques
L3 M, MA et DL
Analyse numérique – 2022-2023

Feuille d’exercices no 3
Exercice 1. On considère l’équation différentielle y 0 = 2ty avec la condition initiale y(0) = 1.

1) Déterminer la solution exacte, calculer y(1) et esquisser la courbe intégrale t 7→ y(t).
2) Résoudre numériquement (sur papier) ce problème de Cauchy en utilisant la méthode d’Euler
et puis celle du point milieu pour calculer des valeurs approchées pour y(1). Effectuer chaque
calcul deux fois ; d’abord en prenant n = 1 et puis n = 2. (Pour n = 2 on veut approcher la
solution exacte en utilisant deux sous-intervalles et en obtenant les valeurs successives y1 ≈ y(.5)
et y2 ≈ y(1)).
3) Représenter sur le même graphe les lignes polygonales correspondant aux solutions numériques.
Comparer graphiquement les erreurs locales pour les deux méthodes.
Exercice 2 (python). On considère l’équation différentielle y 0 = y 2 avec la condition initiale

y(0) = −1.
1) Utiliser les méthodes d’Euler et du point milieu pour résoudre numériquement ce problème de
Cauchy sur l’intervalle [0, 1]. Développez vos propres fonctions python !
2) (Re)écrire les algorithmes de telle manière qu’on puisse les utiliser pour résoudre d’autres
problèmes de Cauchy.
3) Tracer dans le même système de coordonnées, pour chacune des méthodes, le graphe de l’erreur
|y(1) − yn | en fonction de n. Des coordonnées logarithmiques pourraient être utiles.
Exercice 3. On considère les équations différentielles suivantes :
(i) y 0 = t sin y (ii) y 0 = t cos y.
Faire correspondre chacun des champs directionnels ci-dessous à une des équations. Justifier les
choix.
y y
t t
Exercice 4. On considère l’équation de l’oscillateur harmonique amorti
y 00 = −y − 2ky 0 avec y(0) = 1 et y 0 (0) = 0,
où k est un paramètre réel positif ou nul.

1) Écrire la solution du problème en fonction de k. Il faudrait distinguer les cas k < 1, k = 1 et
k > 1.
2) (python) Donner la solution numérique dans les cas k = 1/3, k = 1 et k = 1.5 sur l’intervalle
[0, 5] en utilisant la méthode du point milieu (développée dans l’exercice 2) avec n = 300 pas.
1
26 février, 2023
Comparer les solutions numériques avec les solutions exactes en les dessinant dans un même système
de coordonnées.
3) Utiliser la fonction odeint de la bibliothèque scipy.integrate pour résoudre les problèmes du
point précédent. Introduire dans les comparaisons graphiques les nouvelles solutions numériques.
Exercice 5 (python — suite de l’exercice 4). Dessiner le champ de vecteurs (bidimensionnel) dans
l’espace des phases du problème linéaire associé à l’équation y 00 = −y − 2ky 0 . Faire la liaison entre
l’aspect du champ de vecteurs dans le voisinage de l’origine et les solutions explicites. (On pourrait
utiliser la fonction quiver.)
Exercice 6. On considère l’équation non linéaire x00 = −x + x2 .

1) Écrire l’équation d’ordre 1 en (x, y) associée.
2) Déterminer, dans l’espace des phases, les points stationnaires. Ce sont les orbites les plus
simples ou les zéros du champ de vecteurs !
3) Étudier la linéarisation du champ de vecteurs dans le voisinage des points stationnaires
Éventuellement indiquer les trajectoires du problème linéarisé autour des points stationnaires
4) En déduire les autres orbites dans l’espace des phases. Esquisser-les ! On rappelle que
• deux orbites ne peuvent s’intersecter sauf en un point stationnaire
• une orbite périodique contient en son intérieur au moins un point stationnaire
• si une orbite entre dans une région bornée qu’elle ne quitte plus par la suite, alors la région
contient au moins un point fixe.
5) (python) Vérifier l’esquisse du point précédent en représentant le champ de vecteurs en util-
isant la fonction quiver.
6) Classifier les orbites non stationnaires (celles qui se rassemblent forment une classe).
Exercice 7. On considère l’équation différentielle


x 0 = 1
2 x − 12 y − 1
30 x(x2 + 5y 2 )3/4
y 0 = 1
2 x + 21 y − 1
25 y(x2 + y 2 )3/4 .
1) Calculer J la jacobienne en (0, 0) du champ de vecteurs associé. Le champ de vecteurs associé

à la jacobienne représente la linéarisation du champ initial dans un voisinage de l’origine.
2) Calculer det(J) et tr(J). En inspectant les cas des champs de vecteurs vus lors de l’étude de
l’oscillateur harmonique, quelle sera l’allure des trajectoires dans le voisinage de l’origine ?
3) (python) En utilisant odeint, résoudre numériquement le problème de Cauchy associé à
plusieurs conditions initiales sur l’intervalle [0, 50]. On pourrait considérer (x(0), y(0)) = (−0.15, 0.1)
ou (−0.15, 0). Dessiner les trajectoires et vérifier vos interprétations concernant le comportement
autour de l’origine. Dessiner aussi la solution si (x(0), y(0)) = (−6, −6). Utiliser la fonction quiver
pour représenter dans le même espace des phases le champ de vecteurs.
4) Reprendre l’exercice pour l’équation z 00 + 21 z 0 [(z 0 )2 − 1] + z = 0.
Exercice 8.
1) Identifier (éventuellement en explicitant l’algorithme) le tableau numérique qui caractérise la
méthode de Heun.
2) En suivant la preuve pour l’ordre de la méthode du point milieu (de Runge), étudier l’ordre
de la méthode de Heun pour le problème de Cauchy
y 0 (t) = (1 + t2 ) sin y(t) et y(0) = 1.
Indication : Il faut faire la comparaison avec le développement de Taylor de la solution exacte à
l’ordre 3.

Analyse Numérique - 2022-2023: Exercice 1

Transféré par

Droits d'auteur :

Formats disponibles

Analyse Numérique - 2022-2023: Exercice 1

Transféré par

Informations du document

Description originale:

Titre original

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

Analyse Numérique - 2022-2023: Exercice 1

Transféré par

Droits d'auteur :

Formats disponibles

Université d’Angers

Analyse numérique – 2022-2023

Exercice 1. On considère l’équation différentielle y 0 = 2ty avec la condition initiale y(0) = 1.

Exercice 2 (python). On considère l’équation différentielle y 0 = y 2 avec la condition initiale

Exercice 3. On considère les équations différentielles suivantes :

(i) y 0 = t sin y (ii) y 0 = t cos y.

Exercice 4. On considère l’équation de l’oscillateur harmonique amorti

y 00 = −y − 2ky 0 avec y(0) = 1 et y 0 (0) = 0,

où k est un paramètre réel positif ou nul.

Exercice 6. On considère l’équation non linéaire x00 = −x + x2 .

Exercice 7. On considère l’équation différentielle

1) Calculer J la jacobienne en (0, 0) du champ de vecteurs associé. Le champ de vecteurs associé

Vous aimerez peut-être aussi