Équations différentielles
Équation différentielle du 1er ordre linéaire à coefficients constants
Forme générale : \(y' + ay = b\), avec \(a, b \in \mathbb{R}\), \(a \neq 0\).
La solution générale est la somme d'une solution particulière \(y_p\) et de la solution générale de l'équation homogène associée \(y' + ay = 0\) :
Équation homogène \(y' + ay = 0\)
Les solutions sont exactement les fonctions :
Trouver une solution particulière de \(y' + ay = b\)
- Si \(b\) est une constante, chercher \(y_p\) constante : \(y_p = \dfrac{b}{a}\).
- Si \(b = P(x)\,e^{\lambda x}\), chercher \(y_p = Q(x)\,e^{\lambda x}\) avec \(\deg Q = \deg P\) (ou \(\deg P + 1\) si \(\lambda = -a\)).
Solution générale — forme complète
Pour \(y' + ay = b\) (constantes) :
Équation du 2nd ordre \(y'' + py' + qy = 0\)
On résout l'équation caractéristique \(r^2 + pr + q = 0\), de discriminant \(\Delta = p^2 - 4q\) :
Problème de Cauchy
Données une équation différentielle et une condition initiale \(y(x_0) = y_0\), on détermine la valeur de la constante \(C\) en substituant.
Preuve : solutions de \(y' + ay = 0\)
Soit \(y\) une solution. Posons \(z(x) = y(x)\,e^{ax}\). Alors :
Donc \(z\) est constante : \(z(x) = C\), soit \(y(x) = Ce^{-ax}\). Réciproquement, toute fonction de cette forme est bien solution. \(\square\)
Preuve : unicité de la solution au problème de Cauchy
Soient \(y_1\) et \(y_2\) deux solutions de \(y' + ay = b\) vérifiant \(y(x_0) = y_0\). Alors \(w = y_1 - y_2\) vérifie \(w' + aw = 0\) et \(w(x_0) = 0\). Donc \(w(x) = Ce^{-ax}\) avec \(Ce^{-ax_0} = 0\), d'où \(C = 0\) et \(y_1 = y_2\). \(\square\)
Preuve : structure de l'ensemble des solutions
L'ensemble des solutions de \(y' + ay = b\) est un sous-espace affine de dimension 1 : si \(y_p\) est une solution particulière, toute solution s'écrit \(y = y_p + y_h\) où \(y_h\) est solution de l'équation homogène. Preuve : \((y - y_p)' + a(y - y_p) = (y' + ay) - (y_p' + ay_p) = b - b = 0\). \(\square\)
Résoudre \(y' + ay = b\) — méthode complète
- Résoudre l'homogène \(y' + ay = 0\) : solution \(y_h = Ce^{-ax}\).
- Trouver une solution particulière \(y_p\) (ici \(y_p = b/a\) si \(b\) est constante).
- Solution générale : \(y = Ce^{-ax} + \dfrac{b}{a}\).
- Condition initiale : substituer \(y(x_0) = y_0\) pour trouver \(C\).
Variation de la constante (méthode générale)
Pour \(y' + a(x)y = b(x)\), on pose \(y = C(x)\,e^{-A(x)}\) où \(A\) est une primitive de \(a\). On trouve \(C'(x) = b(x)\,e^{A(x)}\), d'où \(C(x)\) par intégration.
Modèle exponentiel \(y' = ky\)
Solution : \(y(t) = y_0\,e^{kt}\) où \(y_0 = y(0)\). Si \(k > 0\) : croissance exponentielle. Si \(k < 0\) : décroissance (demi-vie \(T_{1/2} = -\ln 2 / k\)).
Modèle logistique \(y' = ky(L - y)\)
Solution explicite obtenue par séparation des variables. La population tend vers \(L\) (capacité limite). Utile pour modéliser une croissance limitée.
La solution de \(y' + ay = 0\) est \(Ce^{-ax}\) avec \(C \in \mathbb{R}\) quelconque. Écrire uniquement \(e^{-ax}\) est incomplet — c'est une solution particulière, pas la solution générale.
Pour \(y' - 2y = 4\), la solution particulière constante vérifie \(-2y_p = 4\), donc \(y_p = -2\). Erreur classique : oublier de la calculer et écrire \(y = Ce^{2x}\) seulement.
Pour \(y' + 3y = 0\), la solution est \(Ce^{-3x}\) (signe moins). Pour \(y' - 3y = 0\), c'est \(Ce^{+3x}\). Toujours lire le signe devant \(ay\) dans l'équation normalisée.
Déterminer \(C\) en fin de résolution uniquement, une fois la solution générale complète écrite. Substituer trop tôt conduit à des erreurs de calcul.
Python — méthode d'Euler
Résolution numérique de \(y' = f(x, y)\), \(y(x_0) = y_0\), avec pas \(h\) :
def euler(f, x0, y0, h, n):
"""Résout y' = f(x, y) par la méthode d'Euler."""
x, y = x0, y0
for _ in range(n):
y = y + h * f(x, y)
x = x + h
return x, y
# Exemple : y' = -2y, y(0) = 1 → solution exacte y = e^{-2x}
f = lambda x, y: -2 * y
xf, yf = euler(f, x0=0, y0=1, h=0.01, n=100)
print(f"y(1) ≈ {yf:.5f}") # ≈ 0.13534
import math
print(f"Exact : {math.exp(-2):.5f}") # ≈ 0.13534
Python — résolution symbolique avec sympy
import sympy as sp
x = sp.Symbol('x')
y = sp.Function('y')
# Résoudre y' + 2y = 4
eq = sp.Eq(y(x).diff(x) + 2*y(x), 4)
sol = sp.dsolve(eq)
print(sol) # y(x) = C1*exp(-2*x) + 2
# Avec condition initiale y(0) = 3
sol_ci = sp.dsolve(eq, ics={y(0): 3})
print(sol_ci) # y(x) = exp(-2*x) + 2
Vérifier sa solution en la re-dérivant et en substituant dans l'équation. Si \(y = Ce^{-ax} + b/a\), alors \(y' = -aCe^{-ax}\) et \(y' + ay = -aCe^{-ax} + aCe^{-ax} + b = b\). ✓
Éq. diff. & Suites — méthode d'Euler discrète
La suite \(u_{n+1} = u_n + h\cdot f(t_n, u_n)\) est la version discrète de \(y' = f(t, y)\). Étudier la convergence de la suite revient à analyser la précision de l'approximation numérique.
Modèle de refroidissement de Newton
La température \(T(t)\) d'un corps vérifie \(T' = -k(T - T_{\text{ext}})\). Solution :
Problème type BAC : déterminer \(k\) à partir de deux mesures de température, puis prédire le moment où \(T\) atteint une valeur cible.
Éq. diff. & Probabilités — modèle de population
En temps continu, la fraction \(p(t)\) vérifiant une propriété suit souvent \(p' = \lambda(1 - p)\), de solution \(p(t) = 1 - (1-p_0)e^{-\lambda t}\). Limite : \(p(t) \to 1\) quand \(t \to +\infty\).
-
1. La solution générale de \(y' - 3y = 0\) est :
\(y' = 3y \Rightarrow y = Ce^{3x}\). L'équation est \(y' - 3y = 0\), donc \(a = -3\) dans \(y' + ay = 0\) : solution \(Ce^{-(-3)x} = Ce^{3x}\).
-
2. Quelle est la solution de \(y' + 2y = 6\) vérifiant \(y(0) = 4\) ?
Solution générale : \(y = Ce^{-2x} + 3\). En \(x=0\) : \(C + 3 = 4\), donc \(C = 1\). Solution : \(y = e^{-2x} + 3\).
-
3. Une quantité radioactive vérifie \(N' = -0{,}1 N\) avec \(N(0) = 500\). La demi-vie est :
\(N(T) = N_0/2 \Rightarrow e^{-0.1 T} = 1/2 \Rightarrow T = \ln 2 / 0.1 \approx 6.93\). (Notons que \(10\ln 2 = \ln 2/0.1\), donc A et D sont identiques — A est la bonne formulation.)
-
4. L'équation caractéristique de \(y'' - 5y' + 6y = 0\) a pour racines :
\(r^2 - 5r + 6 = 0 \Rightarrow (r-2)(r-3) = 0\). Racines : \(r = 2\) et \(r = 3\). Solution générale : \(y = C_1 e^{2x} + C_2 e^{3x}\).