Équations différentielles

Bon nombre des lois les plus importantes de la nature, du mouvement des planètes à la croissance des populations, ne sont pas décrites par de simples formules mais par des équations qui relient une fonction à ses taux de variation. Celles-ci sont connues sous le nom d'équations différentielles. Elles constituent le langage de la science et de l'ingénierie modernes.
Dans ce chapitre, nous apprendrons ce que sont les équations différentielles et nous explorerons plusieurs techniques fondamentales pour trouver leurs solutions. Nous commencerons par des méthodes analytiques permettant d'obtenir des solutions exactes, telles que la séparation des variables. Nous étudierons ensuite des méthodes numériques, comme la méthode d'Euler, qui nous permettent d'approximer des solutions lorsque les solutions exactes sont inaccessibles.

Une équation différentielle est une équation qui contient une fonction inconnue et une ou plusieurs de ses dérivées.

• L'ordre d'une équation différentielle est l'ordre de la plus haute dérivée qu'elle contient.
• Une solution générale est une famille de fonctions qui satisfait à l'équation, impliquant généralement une ou plusieurs constantes arbitraires.
• Une condition initiale est une valeur spécifiée de la fonction ou de ses dérivées en un point particulier.
• Une solution particulière est une solution spécifique correspondant à une condition initiale donnée. Elle peut être obtenue en utilisant une condition initiale pour déterminer les valeurs des constantes arbitraires dans la solution générale.

Une pomme est lâchée sans vitesse initiale d'une hauteur de 10 mètres. Sa position verticale, \(y(t)\), est régie par l'équation différentielle du second ordre :$$\frac{d^2 y}{dt^2}=-g$$où \(g\) est la constante de l'accélération gravitationnelle.

1. Énoncer les conditions initiales pour la position \(y(0)\) et la vitesse \(y'(0)\).
2. Vérifier que la solution générale de cette équation est \(y(t) = -\frac{1}{2}gt^2 + At + B\).
3. Utiliser les conditions initiales pour trouver la solution particulière décrivant le mouvement de la pomme.

Alors que les techniques algébriques nous donnent des formules pour les solutions, une approche géométrique peut nous donner une intuition puissante sur le comportement des solutions. Une équation différentielle de la forme \(\dfrac{dy}{dx} = f(x, y)\) est une machine qui nous donne la pente d'une courbe solution en n'importe quel point \((x,y)\) du plan.
Un champ de pentes (ou champ de directions) est une représentation graphique de cette information. À chaque point d'une grille, nous dessinons un petit segment de droite avec la pente donnée par l'équation différentielle. Le champ de pentes résultant agit comme un ensemble de ``courants'' dans une rivière. Toute solution particulière de l'équation différentielle est une courbe qui ``suit'' ces courants, toujours tangente aux marqueurs de pente qu'elle traverse.

Pour esquisser un champ de pentes pour \(\dfrac{dy}{dx} = f(x, y)\) :

1. Choisir des points de la grille : Sélectionner un ensemble représentatif de points à coordonnées entières \((x, y)\).
2. Calculer les pentes : Pour chaque point, calculer la valeur de la pente \(m = f(x, y)\). Organiser ces valeurs dans un tableau.
3. Dessiner les segments : À chaque point de la grille, dessiner un court segment de droite avec la pente calculée.

Pour l'équation différentielle \(\dfrac{dy}{dx} = x-y\) :

1. Esquisser le champ de pentes pour les coordonnées entières où \(-1 \le x \le 2\) et \(-1 \le y \le 2\).
2. Sur votre esquisse, dessiner la courbe de la solution particulière qui passe par le point \((0,1)\).

Lorsque la dérivée d'une fonction ne dépend que de \(x\), on peut trouver la solution générale en réarrangeant l'équation et en intégrant directement.$$\begin{aligned}\frac{dy}{dx} &= f(x)\\ dy &= f(x)\, dx \\ \int dy &= \int f(x) \, dx \\ y &= \int f(x) \, dx + C\end{aligned}$$

Trouver la solution générale de \(\dfrac{dy}{dx} = 3x^2 \) et la solution particulière avec la condition initiale \(y(0)=4\).

Une équation différentielle à variables séparables est une équation qui peut être exprimée sous la forme \(\dfrac{dy}{dx} = g(x)h(y)\). Ces équations peuvent être résolues en séparant algébriquement les variables avant d'intégrer.

Cette méthode s'applique lorsque l'équation peut être réarrangée de sorte que tous les termes en \(y\) soient d'un côté et tous les termes en \(x\) de l'autre.

1. Séparer : Réarranger l'équation sous la forme \(\dfrac{1}{h(y)}\, dy = g(x)\, dx\).
2. Intégrer : Intégrer les deux membres de l'équation par rapport à leurs variables respectives :$$ \int \dfrac{1}{h(y)} \, dy = \int g(x) \, dx $$
3. Isoler \(y\) : Si possible, résoudre l'équation résultante pour \(y\) afin d'obtenir une solution explicite.

Résoudre \(\dfrac{dy}{dx} = -2xy^2\) avec la condition initiale \(y(1)=1\).

Lorsqu'une équation différentielle ne peut pas être résolue analytiquement, on peut approximer sa solution à l'aide de méthodes numériques. La méthode d'Euler est la plus simple, utilisant l'approximation linéaire pour avancer le long de la courbe solution.
La méthode commence par une condition initiale \((x_0, y_0)\) et procède par étapes. Définissons un petit pas constant \(h\). La coordonnée \(x\) de chaque point suivant est trouvée en ajoutant ce pas : \(x_{n+1} = x_n + h\).
Le cœur de la méthode consiste à approximer la dérivée \(\dfrac{dy}{dx}\) en utilisant la pente du segment de droite reliant deux points consécutifs, \((x_n, y_n)\) et \((x_{n+1}, y_{n+1})\) :$$ \dfrac{dy}{dx} \text{ en } (x_n, y_n) \approx \dfrac{y_{n+1} - y_n}{x_{n+1} - x_n} = \dfrac{y_{n+1} - y_n}{h}. $$Puisque l'équation différentielle nous donne la valeur exacte de la dérivée, \(\dfrac{dy}{dx} = f(x_n, y_n)\), nous pouvons les égaler :$$ \dfrac{y_{n+1} - y_n}{h}\approx f(x_n, y_n). $$En réarrangeant cette formule pour trouver la valeur suivante de \(y\), \(y_{n+1}\), on obtient l'étape itérative :$$ y_{n+1} \approx y_n + h \cdot f(x_n, y_n). $$La collection de segments de droite créée par cette procédure itérative forme une approximation polygonale de la courbe solution réelle.

Pour approximer la solution de \(\dfrac{dy}{dx} = f(x,y)\) avec la condition initiale \((x_0, y_0)\) et un pas \(h\), les coordonnées du point suivant \((x_{n+1}, y_{n+1})\) sont trouvées à partir du point précédent \((x_n, y_n)\) en utilisant les formules itératives :$$\begin{cases}x_{n+1} = x_n + h, \\ y_{n+1} = y_n + h \cdot f(x_n, y_n).\end{cases}$$