Aller au contenu
Télécharger l'application

Accueil › ECG › Maths appliquées

ECG2 · Maths appliquées · semaine 30

Équations et systèmes différentiels en Python

Illustration informatique du cours sur les équations et les systèmes différentiels. Représentations graphiques de trajectoires et approche qualitative sur des exemples simples. La résolution numérique d'une équation différentielle n'est pas au programme.

Sommaire

Le cours

Résumé en cinq points

  1. L'objectif est d'illustrer graphiquement le cours, pas de résoudre numériquement, ce qui est hors programme.
  2. ==On calcule d'abord l'expression exacte des solutions, puis on l'évalue sur une grille de temps créée par linspace.==
  3. L'équilibre s'obtient en annulant le second membre, et le signe des valeurs propres décide de la stabilité.
  4. Faire varier la condition initiale ne change pas la limite pour les solutions convergeant vers un même équilibre stable ; quand l'équilibre est instable, elle peut changer le comportement.
  5. Un paramètre de vitesse de réaction modifie la rapidité de convergence ; dans le modèle étudié, il ne change pas la position de l'équilibre.
Pourquoi ce chapitre

Le cours de mathématiques du troisième semestre traite les systèmes différentiels linéaires à coefficients constants, de la forme X′=AXX'=AX, et leur résolution quand AA est diagonalisable en dimension 22 ou 33.

RemarqueObjectif

Le comportement des trajectoires est gouverné par le signe des valeurs propres de AA. On trace des trajectoires, on fait varier les paramètres, on fait varier les conditions initiales, et on observe la convergence ou la divergence.

MéthodeCompétences attendues

  1. Tracer une solution dont on connaît l'expression.
  2. Interpréter le rôle d'un paramètre.
  3. Reconnaître visuellement une convergence vers un équilibre.

AttentionRésolution numérique

Aucun algorithme de résolution numérique n'est exigé.

Classiques et HP

Méthodes et approfondissements du chapitre · le statut est précisé pour ECG2.

Classiques et HP réservés

Prolongements et démonstrations de Maths appliquées.

Essayer 7 jours gratuitement Sans engagement · cours, exercices et fiches complets

Fin de l'aperçu gratuit

Encore 13 parties de cours et 9 corrigés.

Essayer 7 jours gratuitement Sans engagement · cours, exercices et fiches complets

Méthodes et exercices

Exercices classiques corrigés

Les exercices se suivent : garde leur ordre. Difficulté et durée indicatives.

Exercice 1. Stock de capital ★

Application du cours · 15 min · Niveau khôlle · Équation scalaire linéaire

Un stock KK reçoit un investissement constant de 8 par unité de temps et se déprécie au taux 0,40{,}4, donc K′=8−0,4K.K' = 8 - 0{,}4K. Donner l'équilibre et l'expression de KK.

Écrire un programme qui trace KK pour K0∈{0,10,20,30,40}K_0 \in \{0, 10, 20, 30, 40\} avec l'équilibre, puis commenter.

Voir l'indice

Technique résolution exacte puis boucle sur les conditions initiales.

Voir le corrigé offert
  1. L'équilibre annule 8 − 0,4K, donc K* = 20. 2. K(t) = 20 + (K₀ − 20) e−0,4t. 3. La constante de temps vaut 2,5. Sur [0, 15], l'écart est divisé par e⁶, soit environ 400. 4. Programme.
import numpy as np
import matplotlib.pyplot as plt
t = np.linspace(0, 15, 300)
for K0 in [0, 10, 20, 30, 40]:
    K = 20 + (K0 - 20) * np.exp(-0.4 * t)
    plt.plot(t, K, label="K0 = " + str(K0))
plt.plot(t, 20 + 0 * t, "k--")
plt.xlabel("temps")
plt.ylabel("stock")
plt.legend()
plt.show()
  1. Commentaire. Toutes les courbes tendent vers 20, en croissant si K₀ < 20 et en décroissant si K₀ > 20. La courbe issue de 20 est constante. L'équilibre est stable, la limite ne dépend pas de K₀.

piège : Écrire math.exp au lieu de np.exp. La fonction du module math refuse un tableau.

Exercice 2. Paramètres d'un marché ★★

Classiques de khôlle et concours · 25 min · Niveau khôlle · Équation scalaire linéaire

On note D(p)=α−βpD(p) = \alpha - \beta p et S(p)=γ+δp,S(p) = \gamma + \delta p, avec α>γ\alpha > \gamma et β,δ,k\beta, \delta, k strictement positifs.

Le prix suit p′=k(D(p)−S(p)).p' = k(D(p) - S(p)). Déterminer l'équilibre.

Dire ce que change kk, puis ce que change une hausse de α\alpha.

Proposer un tracé.

Voir le début du corrigé
  1. p' = k(α − γ) − k(β + δ)p. 2. L'équilibre vaut p* = (α − γ)/(β + δ), strictement positif car α > γ. …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 3. Retrouver l'équation à partir d'un tracé ★★

Classiques de khôlle et concours · 20 min · Niveau khôlle · Identification de l'équation

Une courbe solution de y′=ay+by' = ay + b part de 10, décroît vers une asymptote d'ordonnée 4, et vaut 7 vers t=0,35t = 0{,}35.

Retrouver aa et bb.

Voir le début du corrigé
  1. L'asymptote donne y* = 4 = −b/a. …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 4. Une dette qui s'emballe ★★

Classiques de khôlle et concours · 20 min · Niveau khôlle · Équation scalaire linéaire

Une dette DD porte un intérêt continu de taux 0,50{,}5, et on rembourse 2 par unité de temps, donc D′=0,5D−2.D' = 0{,}5D - 2. Étudier l'équilibre.

Pour D0=2D_0 = 2, trouver l'instant où la dette s'annule.

Tracer pour D0∈{2,3,4,5,6}D_0 \in \{2, 3, 4, 5, 6\} sur [0,3][0, 3].

Voir le début du corrigé
  1. L'équilibre vaut 44, et D(t)=4+(D0−4)e0,5tD(t) = 4 + (D_0 - 4)e^{0{,}5t}.

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 5. Système stable et direction d'arrivée ★★

Classiques de khôlle et concours · 45 min · Niveau écrit Ecricome · Système linéaire et valeurs propres

Soit AA de lignes (−2,1)(-2, 1) et (2,−3)(2, -3).

Résoudre X′=AXX' = AX avec X(0)=(3,0)X(0) = (3, 0).

Tracer la trajectoire dans le plan de phase et commenter.

Voir le début du corrigé
  1. La trace vaut −5-5 et le déterminant 44. …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 6. Deux marchés qui échangent leurs clients ★★

Classiques de khôlle et concours · 35 min · Niveau écrit BCE, EDHEC ou emlyon · Système linéaire et valeurs propres

On pose x′=−0,3x+0,2yx' = -0{,}3x + 0{,}2y et y′=0,3x−0,2y,y' = 0{,}3x - 0{,}2y, où xx et yy sont les clientèles de deux enseignes.

Montrer que x+yx + y est constant, résoudre, et tracer xx pour plusieurs répartitions initiales de même total 100.

Voir le début du corrigé
  1. x′+y′=0x' + y' = 0, donc x+y=Nx + y = N est constant.

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 7. Sensibilité autour d'un point selle ★★

Classiques de khôlle et concours · 30 min · Niveau écrit BCE, EDHEC ou emlyon · Système linéaire et valeurs propres

Soit AA de lignes (1,2)(1, 2) et (2,1)(2, 1).

Déterminer les conditions initiales donnant une trajectoire bornée.

Pour X0=(1,−1+ε)X_0 = (1, -1 + \varepsilon), estimer l'instant où la trajectoire s'écarte nettement de l'origine.

Voir le début du corrigé
  1. Les valeurs propres sont 33 et −1-1, avec V1=(1,1)V_1 = (1, 1) et V2=(1,−1)V_2 = (1, -1).

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 8. Discussion selon un paramètre ★★★

Approfondissement · 60 min · Niveau écrit ESSEC ou HEC · Discussion selon un paramètre

Soit m réel et A(m) de lignes (m, 1) et (1, m). Discuter selon m le comportement des solutions de X' = A(m)X et illustrer par un programme.

Voir le début du corrigé
  1. V1=(1,1)V_1 = (1, 1) est propre pour m+1m + 1 et V2=(1,−1)V_2 = (1, -1) pour m−1m - 1, quel que soit mm.

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Exercice 9. Contrôler un calcul avec odeint ★

Application du cours · 15 min · Type oral avec Python · Résolution numérique avec odeint

L'énoncé fournit odeint.

Pour le système de l'exercice 5, écrire la fonction dérivée, obtenir la solution approchée, et mesurer l'écart maximal avec la formule exacte.

Voir le début du corrigé
  1. La fonction s'écrit def derivee(X, t): suivie de return [-2 * X[0] + X[1], 2 * X[0] - 3 * X[1]]. …

La suite du corrigé dans Sup et Spé+ : essai gratuit 7 jours.

Fiche de cours rédigée à partir du programme officiel de la voie EC, enseignement annuel d'informatique et algorithmique de seconde année, partie I « Programme du troisième semestre », paragraphe 2 « Équations et systèmes différentiels ».