Me contacter

Spécialité maths, Terminale • Exercices corrigés à Montréal

Exercices corrigés : les équations différentielles (Terminale spécialité)

Voici une série d'exercices corrigés de spécialité mathématiques sur les équations différentielles, au niveau de la classe de Terminale du programme français, tel qu'il est suivi au Lycée Marie de France et au Collège Stanislas à Montréal.

Deux réflexes se paient cher au baccalauréat : écrire y=Ceax+by=C\mathrm{e}^{ax}+b au lieu de y=Ceaxbay=C\mathrm{e}^{ax}-\frac{b}{a}, et confondre la limite avec une valeur atteinte. Les exercices 9 et 6 sont construits contre ces deux automatismes.

Rappel de cours

  • L'équation y=f(x)y'=f(x) se résout en cherchant une primitive : les solutions sont y=F(x)+Cy=F(x)+C, où FF est une primitive de ff sur l'intervalle considéré.
  • Les solutions de y=ayy'=ay sur R\mathbb{R} sont les fonctions xCeaxx\mapsto C\mathrm{e}^{ax}, avec CC réel quelconque.
  • Les solutions de y=ay+by'=ay+b (avec a0a\neq 0) sont les fonctions xCeaxbax\mapsto C\mathrm{e}^{ax}-\frac{b}{a}. La constante ba-\frac{b}{a} est la solution particulière constante, celle qui annule yy'.
  • Structure des solutions : si gg est une solution particulière de y=ay+by'=ay+b, toute solution s'écrit y=g+hy=g+hhh est solution de l'équation homogène y=ayy'=ay.
  • Pour toute condition initiale y(x0)=y0y(x_{0})=y_{0}, il existe une unique solution qui la vérifie : la condition détermine CC, et une seule.
  • Si a<0a<0, toute solution tend vers ba-\frac{b}{a} quand x+x\to+\infty : c'est un équilibre stable, approché mais jamais atteint, sauf si y0=bay_{0}=-\frac{b}{a}.
  • Si a>0a>0, toute solution autre que la constante ba-\frac{b}{a} s'en éloigne : l'équilibre est instable.
  • Méthode standard pour un modèle : poser z=yz=y-\ell\ell est l'équilibre, ce qui ramène à z=azz'=az.

Partie A : Les bases (/50)

Exercice 1 : Vérifier qu'une fonction est solution

Dans chaque cas, dites si la fonction proposée est solution de l'équation différentielle donnée sur R\mathbb{R}, en détaillant le calcul.

  • a) f(x)=3e2xf(x)=3\mathrm{e}^{2x} et l'équation y=2yy'=2y.
  • b) g(x)=5e0,4x+7g(x)=5\mathrm{e}^{-0{,}4x}+7 et l'équation y=0,4y+2,8y'=-0{,}4y+2{,}8.
  • c) h(x)=xexh(x)=x\mathrm{e}^{x} et l'équation y=y+exy'=y+\mathrm{e}^{x}.
  • d) k(x)=x2k(x)=x^{2} et l'équation y=2yy'=2y.
Voir la correction

a) f(x)=6e2xf'(x)=6\mathrm{e}^{2x} et 2f(x)=6e2x2f(x)=6\mathrm{e}^{2x}. Les deux expressions coïncident pour tout réel xx : ff est solution.

b) g(x)=5×(0,4)e0,4x=2e0,4xg'(x)=5\times(-0{,}4)\mathrm{e}^{-0{,}4x}=-2\mathrm{e}^{-0{,}4x}. D'autre part 0,4g(x)+2,8=2e0,4x2,8+2,8=2e0,4x-0{,}4g(x)+2{,}8=-2\mathrm{e}^{-0{,}4x}-2{,}8+2{,}8=-2\mathrm{e}^{-0{,}4x}. Égalité pour tout xx : gg est solution. On remarque que 7=2,80,47=-\frac{2{,}8}{-0{,}4}, c'est bien la solution particulière constante.

c) h(x)=ex+xex=(1+x)exh'(x)=\mathrm{e}^{x}+x\mathrm{e}^{x}=(1+x)\mathrm{e}^{x}, et h(x)+ex=xex+ex=(1+x)exh(x)+\mathrm{e}^{x}=x\mathrm{e}^{x}+\mathrm{e}^{x}=(1+x)\mathrm{e}^{x}. Égalité pour tout xx : hh est solution. Cette équation n'est pas du type y=ay+by'=ay+b, le second membre dépend de xx.

d) k(x)=2xk'(x)=2x et 2k(x)=2x22k(x)=2x^{2}. L'égalité 2x=2x22x=2x^{2} n'a lieu que pour x=0x=0 et x=1x=1. Une fonction n'est solution que si l'égalité vaut sur tout l'intervalle : kk n'est donc pas solution. Vérifier en un seul point ne prouve rien.

Exercice 2 : Résoudre en cherchant une primitive

Résolvez chaque équation différentielle sur l'intervalle indiqué, puis déterminez l'unique solution vérifiant la condition donnée.

  • a) y=3x24x+1y'=3x^{2}-4x+1 sur R\mathbb{R}, avec y(0)=2y(0)=2.
  • b) y=1xy'=\frac{1}{x} sur ]0 ; +[]0\ ;\ +\infty[, avec y(1)=0y(1)=0.
  • c) y=e2xy'=\mathrm{e}^{-2x} sur R\mathbb{R}, avec y(0)=1y(0)=1.
  • d) y=2xx2+1y'=\frac{2x}{x^{2}+1} sur R\mathbb{R}, avec y(0)=5y(0)=5.
Voir la correction

a) Les solutions sont les primitives : y(x)=x32x2+x+Cy(x)=x^{3}-2x^{2}+x+C. La condition y(0)=2y(0)=2 donne C=2C=2, d'où y(x)=x32x2+x+2y(x)=x^{3}-2x^{2}+x+2.

b) y(x)=lnx+Cy(x)=\ln x+C sur ]0 ; +[]0\ ;\ +\infty[. Avec y(1)=ln1+C=C=0y(1)=\ln 1+C=C=0, on obtient y(x)=lnxy(x)=\ln x.

c) Une primitive de e2x\mathrm{e}^{-2x} est 12e2x-\frac{1}{2}\mathrm{e}^{-2x}, donc y(x)=12e2x+Cy(x)=-\frac{1}{2}\mathrm{e}^{-2x}+C. La condition y(0)=12+C=1y(0)=-\frac{1}{2}+C=1 donne C=32C=\frac{3}{2}, d'où y(x)=3212e2xy(x)=\frac{3}{2}-\frac{1}{2}\mathrm{e}^{-2x}. Le facteur 12-\frac{1}{2} vient de la dérivée de 2x-2x : l'oublier est l'erreur la plus fréquente.

d) Le numérateur est la dérivée du dénominateur, qui est strictement positif : y(x)=ln(x2+1)+Cy(x)=\ln(x^{2}+1)+C. Avec y(0)=ln1+C=C=5y(0)=\ln 1+C=C=5, on obtient y(x)=ln(x2+1)+5y(x)=\ln(x^{2}+1)+5.

Exercice 3 : L'équation homogène

On rappelle que les solutions de y=ayy'=ay sur R\mathbb{R} sont les fonctions xCeaxx\mapsto C\mathrm{e}^{ax}.

  • a) Résolvez y=3yy'=3y, puis donnez la solution vérifiant y(0)=2y(0)=-2.
  • b) Résolvez y+5y=0y'+5y=0 avec y(0)=2y(0)=2.
  • c) Une population vérifie P=0,25PP'=0{,}25P avec P(0)=100P(0)=100, le temps étant en années. Au bout de combien d'années la population a-t-elle doublé ?
  • d) Une masse de matière radioactive vérifie m=0,1mm'=-0{,}1m, le temps étant en jours. Calculez sa demi-vie, puis vérifiez qu'elle ne dépend pas de la masse initiale.
Voir la correction

a) Les solutions sont y(x)=Ce3xy(x)=C\mathrm{e}^{3x}. La condition y(0)=C=2y(0)=C=-2 donne y(x)=2e3xy(x)=-2\mathrm{e}^{3x}.

b) L'équation s'écrit y=5yy'=-5y, donc y(x)=Ce5xy(x)=C\mathrm{e}^{-5x}, et y(0)=C=2y(0)=C=2 : y(x)=2e5xy(x)=2\mathrm{e}^{-5x}. Attention au signe, le +5y+5y passe de l'autre côté.

c) P(t)=100e0,25tP(t)=100\mathrm{e}^{0{,}25t}. On cherche tt tel que e0,25t=2\mathrm{e}^{0{,}25t}=2, soit t=ln20,252,77t=\frac{\ln 2}{0{,}25}\approx 2{,}77 années, c'est-à-dire environ 2 ans et 9 mois.

d) m(t)=m0e0,1tm(t)=m_{0}\mathrm{e}^{-0{,}1t}. La demi-vie TT vérifie m0e0,1T=m02m_{0}\mathrm{e}^{-0{,}1T}=\frac{m_{0}}{2}, et m0>0m_{0}>0 se simplifie : e0,1T=12\mathrm{e}^{-0{,}1T}=\frac{1}{2}, donc T=ln20,16,93T=\frac{\ln 2}{0{,}1}\approx 6{,}93 jours. La masse initiale disparaît du calcul, c'est pourquoi la demi-vie caractérise le noyau et non l'échantillon.

Exercice 4 : L'équation avec second membre constant

Pour chaque équation, déterminez d'abord la solution particulière constante, puis l'ensemble des solutions, puis celle qui vérifie la condition initiale.

  • a) y=2y+6y'=2y+6.
  • b) y=0,5y+4y'=-0{,}5y+4 avec y(0)=2y(0)=2. Donnez la limite de la solution en ++\infty.
  • c) Reprenez l'équation de b) avec y(0)=12y(0)=12, et comparez les deux comportements.
  • d) Démontrez que si y0=bay_{0}=-\frac{b}{a}, alors la solution est constante.
Voir la correction

a) La solution constante y=ky=k vérifie 0=2k+60=2k+6, soit k=3k=-3. Les solutions sont y(x)=Ce2x3y(x)=C\mathrm{e}^{2x}-3.

b) La constante vérifie 0=0,5k+40=-0{,}5k+4, soit k=8k=8. Les solutions sont y(x)=Ce0,5x+8y(x)=C\mathrm{e}^{-0{,}5x}+8. La condition y(0)=C+8=2y(0)=C+8=2 donne C=6C=-6, d'où y(x)=6e0,5x+8y(x)=-6\mathrm{e}^{-0{,}5x}+8. Comme 0,5<0-0{,}5<0, e0,5x0\mathrm{e}^{-0{,}5x}\to 0 et limx+y(x)=8\lim_{x\to+\infty}y(x)=8, en croissant depuis 2.

c) Ici C=128=4C=12-8=4 et y(x)=4e0,5x+8y(x)=4\mathrm{e}^{-0{,}5x}+8. La limite est encore 8, mais la solution décroît depuis 12. Les deux courbes encadrent la droite d'équation y=8y=8 et s'en approchent par des côtés opposés sans jamais l'atteindre : c'est un équilibre stable.

d) Si yy est la fonction constante égale à ba-\frac{b}{a}, alors y=0y'=0 et ay+b=a×(ba)+b=b+b=0ay+b=a\times\left(-\frac{b}{a}\right)+b=-b+b=0. L'équation est vérifiée, donc cette constante est solution, et elle vaut bien y0y_{0} en 00. Par unicité de la solution vérifiant une condition initiale donnée, c'est la seule : la solution est constante.

Exercice 5 : Condition initiale décalée et structure des solutions

La condition initiale n'est pas toujours donnée en 00, et l'ensemble des solutions a une structure qu'il faut savoir démontrer.

  • a) Déterminez la solution de y=3yy'=3y vérifiant y(2)=1y(2)=1.
  • b) Déterminez la solution de y=2y+6y'=-2y+6 vérifiant y(1)=4y(1)=4.
  • c) Soient ff et gg deux solutions de y=ay+by'=ay+b sur R\mathbb{R}. Démontrez que fgf-g est solution de y=ayy'=ay.
  • d) Déduisez-en que toute solution de y=ay+by'=ay+b s'écrit CeaxbaC\mathrm{e}^{ax}-\frac{b}{a}, en prenant pour gg la solution constante.
Voir la correction

a) Les solutions sont y(x)=Ce3xy(x)=C\mathrm{e}^{3x}. La condition y(2)=Ce6=1y(2)=C\mathrm{e}^{6}=1 donne C=e6C=\mathrm{e}^{-6}, d'où y(x)=e6e3x=e3x6=e3(x2)y(x)=\mathrm{e}^{-6}\mathrm{e}^{3x}=\mathrm{e}^{3x-6}=\mathrm{e}^{3(x-2)}. Écrire directement e3(x2)\mathrm{e}^{3(x-2)} est plus lisible et se vérifie d'un coup d'œil en x=2x=2.

b) La constante vérifie 0=2k+60=-2k+6, soit k=3k=3, donc y(x)=Ce2x+3y(x)=C\mathrm{e}^{-2x}+3. La condition y(1)=Ce2+3=4y(1)=C\mathrm{e}^{-2}+3=4 donne Ce2=1C\mathrm{e}^{-2}=1, soit C=e2C=\mathrm{e}^{2}, d'où y(x)=e2e2x+3=e22x+3y(x)=\mathrm{e}^{2}\mathrm{e}^{-2x}+3=\mathrm{e}^{2-2x}+3.

c) Posons u=fgu=f-g. Alors u=fg=(af+b)(ag+b)=a(fg)=auu'=f'-g'=(af+b)-(ag+b)=a(f-g)=au. Donc uu est solution de l'équation homogène y=ayy'=ay. Le second membre constant bb s'élimine par différence, c'est tout le mécanisme.

d) Prenons pour gg la fonction constante égale à ba-\frac{b}{a}, qui est solution d'après l'exercice 4. Si ff est une solution quelconque, alors fgf-g est solution de y=ayy'=ay d'après c), donc il existe CC tel que f(x)(ba)=Ceaxf(x)-\left(-\frac{b}{a}\right)=C\mathrm{e}^{ax}, c'est-à-dire f(x)=Ceaxbaf(x)=C\mathrm{e}^{ax}-\frac{b}{a}. Réciproquement toute fonction de cette forme est solution, par le calcul direct de l'exercice 1.

Partie B : Niveau baccalauréat (/50)

Exercice 6 : La loi de refroidissement de Newton

Un café servi à 85 °C refroidit dans une pièce maintenue à 20 °C. On note θ(t)\theta(t) sa température en degrés Celsius après tt minutes.

La loi de Newton affirme que la vitesse de refroidissement est proportionnelle à l'écart avec la température ambiante : θ(t)=k(θ(t)20)\theta'(t)=-k\bigl(\theta(t)-20\bigr), où k>0k>0.

Une mesure donne θ(10)=65\theta(10)=65.

  • a) On pose z(t)=θ(t)20z(t)=\theta(t)-20. Montrez que zz vérifie z=kzz'=-kz, et exprimez θ(t)\theta(t) en fonction de kk.
  • b) Déterminez kk à partir de la mesure à 10 minutes. Donnez une valeur arrondie à 10410^{-4}.
  • c) Au bout de combien de temps le café atteint-il 40 °C ?
  • d) Quelle est la limite de θ\theta en ++\infty ? Le café atteint-il un jour la température de la pièce ?
Voir la correction

a) z=θz'=\theta' car 20 est constante, donc z=k(θ20)=kzz'=-k(\theta-20)=-kz. Les solutions sont z(t)=Cektz(t)=C\mathrm{e}^{-kt}, d'où θ(t)=20+Cekt\theta(t)=20+C\mathrm{e}^{-kt}. La condition θ(0)=85\theta(0)=85 donne C=65C=65, donc θ(t)=20+65ekt\theta(t)=20+65\mathrm{e}^{-kt}.

b) θ(10)=20+65e10k=65\theta(10)=20+65\mathrm{e}^{-10k}=65, donc e10k=4565=913\mathrm{e}^{-10k}=\frac{45}{65}=\frac{9}{13}, d'où k=110ln1390,0368k=\frac{1}{10}\ln\frac{13}{9}\approx 0{,}0368 par minute.

c) On résout 20+65ekt=4020+65\mathrm{e}^{-kt}=40, soit ekt=2065=413\mathrm{e}^{-kt}=\frac{20}{65}=\frac{4}{13}, donc t=1kln1341,17870,0367732,1t=\frac{1}{k}\ln\frac{13}{4}\approx\frac{1{,}1787}{0{,}03677}\approx 32{,}1 minutes. Le café passe donc environ 32 minutes au-dessus de 40 °C.

d) Comme k<0-k<0, ekt0\mathrm{e}^{-kt}\to 0 et limt+θ(t)=20\lim_{t\to+\infty}\theta(t)=20. Mais 65ekt65\mathrm{e}^{-kt} est strictement positif pour tout tt : la température du café reste strictement supérieure à 20 °C à chaque instant. Le modèle décrit une limite, pas une valeur atteinte. À 3 heures, l'écart résiduel vaut encore environ 0,09 °C, indétectable en pratique mais non nul dans le modèle.

Exercice 7 : Chute avec frottement fluide

Un parachutiste tombe sans vitesse initiale. Sa vitesse v(t)v(t), en mètres par seconde après tt secondes, vérifie

v(t)=9,80,2v(t)v'(t)=9{,}8-0{,}2\,v(t), avec v(0)=0v(0)=0.

Le terme 0,2v-0{,}2v modélise le frottement de l'air, proportionnel à la vitesse.

  • a) Déterminez la vitesse limite du parachutiste, c'est-à-dire la solution constante de l'équation.
  • b) Résolvez l'équation avec la condition initiale, et vérifiez que la vitesse limite est bien la limite de vv en ++\infty.
  • c) Calculez la vitesse au bout de 5 secondes, arrondie au dixième.
  • d) Au bout de combien de temps le parachutiste atteint-il 90 % de sa vitesse limite ?
Voir la correction

a) La solution constante v=v=\ell vérifie 0=9,80,20=9{,}8-0{,}2\ell, soit =9,80,2=49\ell=\frac{9{,}8}{0{,}2}=49 m/s, soit environ 176 km/h. À cette vitesse le frottement compense exactement le poids et l'accélération s'annule.

b) L'équation s'écrit v=0,2v+9,8v'=-0{,}2v+9{,}8, donc v(t)=Ce0,2t+49v(t)=C\mathrm{e}^{-0{,}2t}+49. La condition v(0)=C+49=0v(0)=C+49=0 donne C=49C=-49, d'où v(t)=49(1e0,2t)v(t)=49\left(1-\mathrm{e}^{-0{,}2t}\right). Comme e0,2t0\mathrm{e}^{-0{,}2t}\to 0, on a bien limt+v(t)=49\lim_{t\to+\infty}v(t)=49 m/s.

c) v(5)=49(1e1)49×0,632131,0v(5)=49\left(1-\mathrm{e}^{-1}\right)\approx 49\times 0{,}6321\approx 31{,}0 m/s. Après 5 secondes le parachutiste est déjà à 63 % de sa vitesse limite : le régime s'installe vite.

d) On résout 49(1e0,2t)=0,9×4949\left(1-\mathrm{e}^{-0{,}2t}\right)=0{,}9\times 49, soit e0,2t=0,1\mathrm{e}^{-0{,}2t}=0{,}1, donc t=ln100,211,5t=\frac{\ln 10}{0{,}2}\approx 11{,}5 secondes. Remarquez que ce temps ne dépend pas de 9,8 : seul le coefficient de frottement fixe la rapidité de l'installation du régime, la gravité ne fixe que la valeur limite.

Exercice 8 : La cuve de saumure

Une cuve contient 1 000 litres de solution salée renfermant 20 kg de sel. On y déverse en continu une saumure titrant 0,03 kg de sel par litre, au débit de 10 litres par minute, et on soutire le mélange, supposé homogène à chaque instant, au même débit de 10 litres par minute.

Le volume reste donc constant. On note q(t)q(t) la masse de sel présente dans la cuve, en kilogrammes, après tt minutes.

  • a) Justifiez que q(t)=0,30,01q(t)q'(t)=0{,}3-0{,}01\,q(t).
  • b) Résolvez cette équation avec la condition initiale.
  • c) Déterminez la limite de qq et interprétez-la en termes de concentration.
  • d) Au bout de combien de temps la cuve contient-elle 25 kg de sel ?
Voir la correction

a) Le sel entrant apporte 0,03×10=0,30{,}03\times 10=0{,}3 kg par minute. Le sel sortant part à la concentration de la cuve, qui vaut q(t)1000\frac{q(t)}{1000} kg par litre, au débit de 10 L/min, soit 10q(t)1000=0,01q(t)\frac{10\,q(t)}{1000}=0{,}01\,q(t) kg par minute. La variation de la masse de sel est la différence entre ce qui entre et ce qui sort : q(t)=0,30,01q(t)q'(t)=0{,}3-0{,}01\,q(t).

b) La solution constante vérifie 0=0,30,01k0=0{,}3-0{,}01k, soit k=30k=30. Les solutions sont q(t)=Ce0,01t+30q(t)=C\mathrm{e}^{-0{,}01t}+30, et q(0)=C+30=20q(0)=C+30=20 donne C=10C=-10, d'où q(t)=3010e0,01tq(t)=30-10\mathrm{e}^{-0{,}01t}.

c) limt+q(t)=30\lim_{t\to+\infty}q(t)=30 kg. La concentration limite dans la cuve vaut alors 301000=0,03\frac{30}{1000}=0{,}03 kg par litre : c'est exactement la concentration de la saumure entrante. À long terme la cuve finit par titrer comme son alimentation, quelle que soit sa charge initiale en sel.

d) On résout 3010e0,01t=2530-10\mathrm{e}^{-0{,}01t}=25, soit e0,01t=0,5\mathrm{e}^{-0{,}01t}=0{,}5, donc t=ln20,0169,3t=\frac{\ln 2}{0{,}01}\approx 69{,}3 minutes. Autrement dit, l'écart à la valeur limite est divisé par deux toutes les 69 minutes environ, ce qui est la demi-vie de l'équation homogène associée.

Exercice 9 : Quatre erreurs à corriger

Chacune des affirmations suivantes est fausse. Corrigez-la et justifiez.

  • a) « Les solutions de y=2yy'=2y sont les fonctions x2exx\mapsto 2\mathrm{e}^{x}. »
  • b) « Les solutions de y=3y+12y'=-3y+12 sont les fonctions xCe3x+12x\mapsto C\mathrm{e}^{-3x}+12. »
  • c) « Puisqu'une condition initiale détermine la solution, deux conditions la déterminent encore mieux : cherchons la solution de y=2yy'=2y telle que y(0)=1y(0)=1 et y(1)=3y(1)=3. »
  • d) « Une solution de y=ay+by'=ay+b avec a<0a<0 finit par valoir exactement ba-\frac{b}{a}. »
Voir la correction

a) Faux : le coefficient aa va dans l'exposant, pas devant. Les solutions sont xCe2xx\mapsto C\mathrm{e}^{2x}, avec CC quelconque. La fonction proposée, f(x)=2exf(x)=2\mathrm{e}^{x}, vérifie f(x)=2ex=f(x)f'(x)=2\mathrm{e}^{x}=f(x) et non 2f(x)2f(x) : elle est solution de y=yy'=y.

b) Faux : la solution particulière constante n'est pas bb mais ba=123=4-\frac{b}{a}=-\frac{12}{-3}=4. Les solutions sont xCe3x+4x\mapsto C\mathrm{e}^{-3x}+4. Vérification : si y=Ce3x+4y=C\mathrm{e}^{-3x}+4 alors y=3Ce3x=3(y4)=3y+12y'=-3C\mathrm{e}^{-3x}=-3(y-4)=-3y+12. La fonction annoncée donnerait y=3y+36y'=-3y+36.

c) Faux : le problème est surdéterminé et n'a en général aucune solution. La condition y(0)=1y(0)=1 impose déjà C=1C=1, donc y(x)=e2xy(x)=\mathrm{e}^{2x} et y(1)=e27,39y(1)=\mathrm{e}^{2}\approx 7{,}39, ce qui contredit y(1)=3y(1)=3. Une équation du premier ordre n'admet qu'une seule constante d'intégration : une condition, et une seule, la fixe.

d) Faux : c'est une limite, pas une valeur atteinte. La solution s'écrit y(x)=Ceaxbay(x)=C\mathrm{e}^{ax}-\frac{b}{a}, et CeaxC\mathrm{e}^{ax} ne s'annule jamais quand C0C\neq 0. Par exemple y(x)=6e0,5x+8y(x)=-6\mathrm{e}^{-0{,}5x}+8 vérifie y(x)<8y(x)<8 pour tout xx, tout en tendant vers 8. Le seul cas où la valeur est atteinte est C=0C=0, c'est-à-dire y0=bay_{0}=-\frac{b}{a}, et la solution est alors constante.

Exercice 10 : Problème de synthèse : un lac et son quota de pêche

Une population de poissons dans un lac est modélisée par la fonction PP, où P(t)P(t) est le nombre d'individus après tt années. En l'absence de pêche, la population croît au taux annuel de 5 %. Un quota autorise le prélèvement de 400 poissons par an, réparti régulièrement.

Le modèle retenu est donc P(t)=0,05P(t)400P'(t)=0{,}05\,P(t)-400.

  • a) Déterminez la population d'équilibre PP^{*}, c'est-à-dire la solution constante.
  • b) On pose z(t)=P(t)Pz(t)=P(t)-P^{*}. Montrez que z=0,05zz'=0{,}05z, puis exprimez P(t)P(t) en fonction de P(0)P(0).
  • c) Le lac compte initialement 9 000 poissons. Décrivez l'évolution et calculez la population après 20 ans.
  • d) Le lac compte initialement 7 500 poissons. Montrez que la population s'éteint et calculez l'année de l'extinction. Quel quota annuel maximal aurait été soutenable dans ce cas ?
Voir la correction

a) La solution constante vérifie 0=0,05P4000=0{,}05P^{*}-400, soit P=4000,05=8000P^{*}=\frac{400}{0{,}05}=8\,000 poissons. À cet effectif, la croissance naturelle produit exactement 400 poissons par an, que la pêche prélève.

b) z=Pz'=P' puisque PP^{*} est constante, donc z=0,05P400=0,05(P8000)=0,05zz'=0{,}05P-400=0{,}05\left(P-8\,000\right)=0{,}05z. Ainsi z(t)=z(0)e0,05tz(t)=z(0)\mathrm{e}^{0{,}05t}, et P(t)=8000+(P(0)8000)e0,05tP(t)=8\,000+\bigl(P(0)-8\,000\bigr)\mathrm{e}^{0{,}05t}.

c) Avec P(0)=9000P(0)=9\,000, on obtient P(t)=8000+1000e0,05tP(t)=8\,000+1\,000\mathrm{e}^{0{,}05t}, fonction croissante qui tend vers ++\infty. Après 20 ans, P(20)=8000+1000e10718P(20)=8\,000+1\,000\mathrm{e}\approx 10\,718 poissons.

d) Avec P(0)=7500P(0)=7\,500, on obtient P(t)=8000500e0,05tP(t)=8\,000-500\mathrm{e}^{0{,}05t}, fonction décroissante. Elle s'annule quand e0,05t=16\mathrm{e}^{0{,}05t}=16, soit t=ln160,0555,5t=\frac{\ln 16}{0{,}05}\approx 55{,}5 : la population s'éteint au cours de la 56e année. L'équilibre est instable, puisque a=0,05>0a=0{,}05>0 : de part et d'autre de 8 000, on s'en éloigne. Pour une population initiale de 7 500, un quota QQ est soutenable si l'équilibre Q0,05\frac{Q}{0{,}05} ne dépasse pas 7 500, soit Q375Q\leq 375 poissons par an. Vingt-cinq poissons de trop dans le quota suffisent à condamner le lac.

Voir aussi

Vous cherchez un tuteur de spécialité maths en Terminale à Montréal ?

Contactez-moi pour une première séance. On travaille sur des exercices calibrés sur le niveau réel des contrôles au Lycée Marie de France et au Collège Stanislas.

Site par Studio Squalli