الانتقال إلى المحتوى
DaoudMathالرياضيات مع محمد بن داود

Équations différentielles - Cours et exercices corrigés | 2e Bac PC et SVT

KORASATI · MATHÉMATIQUES

Retrouvez le cours intégral ci-dessous et sa version A4 mise en page, avec toutes les figures et les exercices corrigés.

Ouvrir et télécharger le cours en PDF

Sur téléphone, utilisez le bouton ci-dessus pour agrandir le document. La lecture est accessible sans inscription.


Objectifs. Résoudre les équations y′=ay+by'=ay+b et y″+ay′+by=0y''+ay'+by=0 à coefficients constants, déterminer les constantes grâce aux conditions initiales et vérifier les solutions obtenues.

Comprendre ce que l’on cherche

Une équation différentielle relie une fonction inconnue à une ou plusieurs de ses dérivées. La lettre yy représente donc une fonction x↦y(x)x\mapsto y(x), et non un nombre. Résoudre l’équation sur un intervalle II, c’est déterminer toutes les fonctions suffisamment dérivables sur II qui vérifient l’égalité pour chaque x∈Ix\in I.

Exemple - Vérifier une solution

Pour y(x)=3e2xy(x)=3e^{2x}, on a y′(x)=6e2x=2y(x)y'(x)=6e^{2x}=2y(x). Cette fonction est une solution de y′=2yy'=2y sur ℝ\mathbb{R}. La fonction nulle est aussi solution. En revanche, y(x)=3exy(x)=3e^x ne l’est pas, car y′=yy'=y et yy n’est jamais nulle.

Une vérification n’est pas une résolution complète

Trouver une fonction qui convient ne prouve pas qu’on a obtenu toutes les solutions. Pour donner la famille complète, on emploie le théorème adapté, puis on détermine les constantes si des conditions supplémentaires sont imposées.

L’équation du premier ordre

L’équation homogène : y′=ayy'=ay

Dans tout ce cours, les lettres aa et bb désignent des constantes réelles.

Toutes les solutions de y′=ayy'=ay

Sur un intervalle II, les solutions sont exactement les fonctions y(x)=Ceax,C∈ℝ.y(x)=Ce^{ax},\qquad C\in\mathbb{R}. Si a=0a=0, on retrouve les fonctions constantes y(x)=Cy(x)=C.

Justification. Pour une solution yy, on pose z(x)=e−axy(x)z(x)=e^{-ax}y(x). Alors z′(x)=e−ax(y′(x)−ay(x))=0.z'(x)=e^{-ax}\bigl(y'(x)-ay(x)\bigr)=0. Ainsi zz est constante sur II, donc y=Ceaxy=Ce^{ax}. Réciproquement, la dérivée de CeaxCe^{ax} est aCeaxaCe^{ax} : toute fonction de cette forme convient.

Une valeur initiale fixe la constante

La condition y(x0)=y0y(x_0)=y_0 donne Ceax0=y0Ce^{ax_0}=y_0, donc C=y0e−ax0C=y_0e^{-ax_0}. L’unique solution qui vérifie cette condition est y(x)=y0ea(x−x0).\boxed{y(x)=y_0e^{a(x-x_0)}}. Par exemple, y′=3yy'=3y et y(1)=2y(1)=2 conduisent à y(x)=2e3(x−1)y(x)=2e^{3(x-1)}. Toujours vérifier la valeur au point x0x_0, qui n’est pas nécessairement 00.

L’équation affine : y′=ay+by'=ay+b

Le cas a≠0a\ne0

La fonction constante yp=−b/ay_p=-b/a est une solution particulière : yp′=0y_p'=0 et ayp+b=0ay_p+b=0. Si yy est une autre solution, alors z=y−yp=y+b/az=y-y_p=y+b/a vérifie z′=azz'=az.

Solution générale et condition initiale

Si a≠0a\ne0, toutes les solutions de y′=ay+by'=ay+b sont y(x)=Ceax−ba,C∈ℝ.y(x)=Ce^{ax}-\frac ba,\qquad C\in\mathbb{R}. Avec y(x0)=y0y(x_0)=y_0, on obtient y(x)=−ba+(y0+ba)ea(x−x0).y(x)=-\frac ba+\left(y_0+\frac ba\right)e^{a(x-x_0)}.

Exemple - y′=−2y+6y'=-2y+6, avec y(0)=1y(0)=1

Ici a=−2a=-2 et b=6b=6. La solution constante est −b/a=3-b/a=3. La famille générale s’écrit y(x)=3+Ce−2xy(x)=3+Ce^{-2x}. La condition initiale donne 1=3+C1=3+C, donc y(x)=3−2e−2x.y(x)=3-2e^{-2x}. Sa dérivée vaut 4e−2x4e^{-2x} et −2y+6=4e−2x-2y+6=4e^{-2x} ; l’équation et la condition initiale sont bien satisfaites. Pour x≥0x\ge0, la fonction croît de 11 vers 33 sans atteindre 33.

Le cas a=0a=0

L’équation devient y′=by'=b. On cherche donc les primitives d’une constante : y(x)=bx+C.y(x)=bx+C. Avec y(x0)=y0y(x_0)=y_0, la solution est y(x)=y0+b(x−x0)y(x)=y_0+b(x-x_0). La formule contenant b/ab/a ne doit jamais être utilisée lorsque a=0a=0.

Interpréter une solution dans un modèle

Une température TT peut, dans un modèle simple, vérifier T′(t)=−k(T(t)−Ta),k>0,T'(t)=-k\bigl(T(t)-T_a\bigr),\qquad k>0, où TaT_a est la température ambiante supposée constante. Avec T(0)=T0T(0)=T_0, T(t)=Ta+(T0−Ta)e−kt.T(t)=T_a+(T_0-T_a)e^{-kt}. Le signe de T0−TaT_0-T_a indique si la température descend ou monte vers TaT_a. La variable tt est limitée au domaine du modèle, généralement t≥0t\ge0 ; kk et tt doivent utiliser des unités compatibles.

Prudence avec les signes

Pour y′+ky=cy'+ky=c, commencer par écrire y′=−ky+cy'=-ky+c. Le coefficient aa est −k-k. Un signe erroné transforme une décroissance attendue en croissance exponentielle.

Le second ordre à coefficients constants

On étudie maintenant y″+ay′+by=0y''+ay'+by=0. Une solution doit être deux fois dérivable. On associe à l’équation différentielle l’équation caractéristique r2+ar+b=0,Δ=a2−4b.r^2+ar+b=0,\qquad \Delta=a^2-4b. L’essai y=erxy=e^{rx} explique ce polynôme : y″+ay′+by=(r2+ar+b)erxy''+ay'+by=(r^2+ar+b)e^{rx}. Le théorème suivant, admis ici, décrit toutes les solutions réelles.

Trois cas pour la solution générale

Si Δ>0\Delta>0 : les racines réelles distinctes sont r1,r2r_1,r_2 et y(x)=Aer1x+Ber2x.y(x)=Ae^{r_1x}+Be^{r_2x}. Si Δ=0\Delta=0 : la racine double est r=−a/2r=-a/2 et y(x)=(A+Bx)erx.y(x)=(A+Bx)e^{rx}. Si Δ<0\Delta<0 : on pose α=−a/2\alpha=-a/2 et β=−Δ/2>0\beta=\sqrt{-\Delta}/2>0, puis y(x)=eαx(Acos⁡(βx)+Bsin⁡(βx)).y(x)=e^{\alpha x}\bigl(A\cos(\beta x)+B\sin(\beta x)\bigr). Dans chaque cas, AA et BB sont des constantes réelles indépendantes.

Exemple - Deux racines distinctes

Pour y″−3y′+2y=0y''-3y'+2y=0, on résout r2−3r+2=(r−1)(r−2)=0r^2-3r+2=(r-1)(r-2)=0. Les solutions sont y=Aex+Be2xy=Ae^x+Be^{2x}. Avec y(0)=1y(0)=1 et y′(0)=0y'(0)=0, on a A+B=1,A+2B=0.A+B=1,\qquad A+2B=0. En soustrayant, B=−1B=-1, puis A=2A=2. La solution cherchée est y(x)=2ex−e2xy(x)=2e^x-e^{2x}.

Exemple - Une racine double

Pour y″−4y′+4y=0y''-4y'+4y=0, la racine double est 22, donc y=(A+Bx)e2xy=(A+Bx)e^{2x}. La dérivée est y′=(B+2A+2Bx)e2x.y'=\bigl(B+2A+2Bx\bigr)e^{2x}. Les conditions y(0)=1y(0)=1, y′(0)=3y'(0)=3 donnent A=1A=1 et B+2A=3B+2A=3, d’où B=1B=1 et y=(1+x)e2xy=(1+x)e^{2x}.

Ne pas perdre une constante

Quand la racine est double, écrire Aerx+BerxAe^{rx}+Be^{rx} ne convient pas : cela se réduit à une seule constante. Le facteur xx dans (A+Bx)erx(A+Bx)e^{rx} est indispensable.

Oscillations et méthode de résolution

Une équation avec Δ<0\Delta<0

Pour y″+2y′+5y=0y''+2y'+5y=0, le discriminant est Δ=4−20=−16\Delta=4-20=-16. Ainsi α=−1\alpha=-1 et β=2\beta=2 : y=e−x(Acos⁡2x+Bsin⁡2x).y=e^{-x}\bigl(A\cos2x+B\sin2x\bigr). On dérive avant d’utiliser les conditions initiales : y′=e−x((−A+2B)cos⁡2x+(−2A−B)sin⁡2x).y'=e^{-x}\bigl((-A+2B)\cos2x+(-2A-B)\sin2x\bigr). Avec y(0)=1y(0)=1 et y′(0)=0y'(0)=0, on obtient A=1A=1 et −A+2B=0-A+2B=0, donc B=1/2B=1/2. Le facteur exponentiel décroissant amortit les oscillations quand x→+∞x\to+\infty.

Le cas y″+ω2y=0y''+\omega^2y=0

Pour ω>0\omega>0, les solutions sont Acos⁡(ωx)+Bsin⁡(ωx)A\cos(\omega x)+B\sin(\omega x). Avec y(0)=y0y(0)=y_0 et y′(0)=v0y'(0)=v_0, y(x)=y0cos⁡(ωx)+v0ωsin⁡(ωx).y(x)=y_0\cos(\omega x)+\frac{v_0}{\omega}\sin(\omega x). La dérivée initiale est une information distincte de la valeur initiale. Deux conditions indépendantes fixent les deux constantes.

Figure disponible dans la version PDF ci-dessus.

Le dessin compare une oscillation d’amplitude constante à une oscillation amortie. Les deux courbes correspondent à des équations différentes.

Une procédure à conserver

  1. Mettre l’équation sous une forme reconnue et relever les coefficients avec leurs signes. Au second ordre, diviser d’abord par le coefficient constant non nul de y″y'' s’il n’est pas égal à 11.

  2. Donner la famille générale, avec toutes les constantes nécessaires.

  3. Calculer la dérivée si une condition porte sur y′y' ; remplacer ensuite xx par la valeur demandée.

  4. Résoudre le système sur les constantes. Vérifier l’équation et les conditions initiales, puis préciser le domaine pertinent.

Exercices corrigés : premier ordre

Exercice 1 - Une condition en un point différent de zéro

Résoudre y′=2yy'=2y sur ℝ\mathbb{R} avec y(1)=3y(1)=3. Calculer y(0)y(0) et déterminer la limite en −∞-\infty.

Correction. Les solutions générales sont y=Ce2xy=Ce^{2x}. Comme Ce2=3Ce^2=3, on obtient C=3e−2C=3e^{-2}, donc y(x)=3e2(x−1)y(x)=3e^{2(x-1)}. Ainsi y(0)=3e−2y(0)=3e^{-2}. Lorsque x→−∞x\to-\infty, 2(x−1)→−∞2(x-1)\to-\infty, donc y(x)→0y(x)\to0. La dérivée vaut 6e2(x−1)=2y(x)6e^{2(x-1)}=2y(x) et y(1)=3y(1)=3, ce qui vérifie les deux données.

Exercice 2 - Équation affine et équilibre

Résoudre y′+3y=6y'+3y=6 avec y(0)=5y(0)=5. Étudier les variations sur [0;+∞[[0;+\infty[ et déterminer quand y(x)≤2,3y(x)\le 2{,}3.

Correction. On écrit y′=−3y+6y'=-3y+6. L’équilibre est 22 et la solution générale est y=2+Ce−3xy=2+Ce^{-3x}. La condition initiale donne C=3C=3 : y(x)=2+3e−3x,y′(x)=−9e−3x<0.y(x)=2+3e^{-3x},\qquad y'(x)=-9e^{-3x}<0. La solution décroît strictement de 55 vers 22. L’inégalité demandée équivaut à 3e−3x≤0,33e^{-3x}\le0{,}3, soit e−3x≤1/10e^{-3x}\le1/10. La croissance de ln\ln donne −3x≤−ln⁡10-3x\le-\ln10, donc x≥ln⁡103≈0,768.x\ge\frac{\ln10}{3}\approx0{,}768. Le signe change quand on divise par −3-3. Enfin y′+3y=6y'+3y=6 et y(0)=5y(0)=5 se vérifient directement.

Exercice 3 - Refroidissement dans un modèle

Une température, en degrés Celsius, vérifie T′=−0,1(T−20)T'=-0{,}1(T-20), où tt est en minutes et t≥0t\ge0. On donne T(0)=80T(0)=80. Déterminer T(t)T(t), puis l’instant où T(t)=35T(t)=35.

Correction. La température d’équilibre est 2020. Les solutions sont T=20+Ce−0,1tT=20+Ce^{-0{,}1t}. Comme C=80−20=60C=80-20=60, T(t)=20+60e−0,1t.T(t)=20+60e^{-0{,}1t}. L’équation T(t)=35T(t)=35 donne e−0,1t=15/60=1/4e^{-0{,}1t}=15/60=1/4, puis −0,1t=−ln⁡4-0{,}1t=-\ln4. Ainsi t=10ln⁡4≈13,86t=10\ln4\approx13{,}86 minutes. Cette solution est unique car T′(t)=−6e−0,1t<0T'(t)=-6e^{-0{,}1t}<0. La limite est 2020, et la température reste strictement supérieure à 2020 pour tout temps fini. Ces résultats sont ceux du modèle supposé.

Cas limite à savoir traiter

Pour y′=0y+4y'=0y+4, y(2)=1y(2)=1, la solution est simplement y(x)=4x−7y(x)=4x-7. On ne remplace pas aa par 00 dans une formule qui divise par aa.

Exercices corrigés : second ordre

Exercice 4 - Deux racines réelles

Résoudre y″+y′−2y=0y''+y'-2y=0 avec y(0)=3y(0)=3 et y′(0)=0y'(0)=0.

Correction. Le polynôme caractéristique est r2+r−2=(r−1)(r+2)r^2+r-2=(r-1)(r+2). Ainsi y=Aex+Be−2xy=Ae^x+Be^{-2x} et y′=Aex−2Be−2xy'=Ae^x-2Be^{-2x}. Les conditions donnent A+B=3,A−2B=0,A+B=3,\qquad A-2B=0, d’où B=1B=1 et A=2A=2. Donc y(x)=2ex+e−2xy(x)=2e^x+e^{-2x}. On vérifie y″+y′−2y=(2ex+4e−2x)+(2ex−2e−2x)−2(2ex+e−2x)=0.y''+y'-2y=(2e^x+4e^{-2x})+(2e^x-2e^{-2x})-2(2e^x+e^{-2x})=0. Les valeurs initiales sont 33 et 00, comme demandé.

Exercice 5 - Racine double

Résoudre y″−6y′+9y=0y''-6y'+9y=0 avec y(0)=2y(0)=2 et y′(0)=5y'(0)=5.

Correction. Le polynôme est (r−3)2(r-3)^2. La solution générale est y=(A+Bx)e3xy=(A+Bx)e^{3x} et sa dérivée vaut (B+3A+3Bx)e3x(B+3A+3Bx)e^{3x}. On a A=2A=2 et B+3A=5B+3A=5, donc B=−1B=-1 : y(x)=(2−x)e3x,y′(x)=(5−3x)e3x,y″(x)=(12−9x)e3x.y(x)=(2-x)e^{3x},\quad y'(x)=(5-3x)e^{3x},\quad y''(x)=(12-9x)e^{3x}. Alors y″−6y′+9y=[12−9x−6(5−3x)+9(2−x)]e3x=0y''-6y'+9y=[12-9x-6(5-3x)+9(2-x)]e^{3x}=0. Les conditions initiales sont également vérifiées.

Exercice 6 - Une oscillation et une dérivée initiale

Résoudre y″+4y=0y''+4y=0 avec y(0)=1y(0)=1, y′(0)=−2y'(0)=-2. Donner une écriture avec un seul cosinus et déterminer la première racine positive.

Correction. On a y=Acos⁡2x+Bsin⁡2xy=A\cos2x+B\sin2x et y′=−2Asin⁡2x+2Bcos⁡2xy'=-2A\sin2x+2B\cos2x. Les conditions donnent A=1A=1 et 2B=−22B=-2, donc y(x)=cos⁡2x−sin⁡2x=2cos⁡(2x+π4).y(x)=\cos2x-\sin2x=\sqrt2\cos\left(2x+\frac\pi4\right). Les zéros vérifient 2x+π/4=π/2+kπ2x+\pi/4=\pi/2+k\pi, soit x=π/8+kπ/2x=\pi/8+k\pi/2 avec k∈ℤk\in\mathbb{Z}. La première racine positive est π/8\pi/8. La dérivée seconde est −4cos⁡2x+4sin⁡2x=−4y-4\cos2x+4\sin2x=-4y, ce qui vérifie l’équation.

Bilan

Premier ordre : une constante, puis une valeur initiale. Second ordre : deux constantes, souvent une valeur et une dérivée initiales. Le discriminant concerne le polynôme caractéristique ; les fonctions obtenues doivent ensuite satisfaire les données du problème.