Je réserve ici
Maths expertes (option Tle) Terminale

Ressources · Terminale · Maths expertes (option Tle)

Suites de matrices

Suites de matrices colonnes Uₙ₊₁ = AUₙ + B : rappels sur les suites arithmético-géométriques, écriture matricielle d’évolutions couplées, cas homogène Uₙ = AⁿU₀, calcul de Aⁿ avec une matrice diagonale, état stable, expression de Uₙ et convergence.

À propos de cette page
Cinq problèmes corrigés dans l’esprit du bac : bassins d’une station d’épuration, modèle proie-prédateur dans un lac avec un seuil en Python, médicament diffusé entre le sang et les tissus, escargots en deux classes d’âge, suite récurrente d’ordre 2 ramenée à une suite de matrices.

5 problèmes, du plus simple au plus difficile. Écris tes réponses, puis « Vérifier ».

Résoudre un problème : la méthode

  1. Je traduis l’énoncé : une équation par grandeur, en comptant ce qui reste, ce qui arrive des autres compartiments et ce qui est ajouté ; puis j’écris A et B, et je contrôle avec U1.
  2. Je trie : B nulle, la suite est homogène et Un=AnU0 ; sinon je cherche l’état stable C en résolvant (I−A)C=B et je pose Vn=Un−C.
  3. Je calcule An avec la matrice P donnée : je vérifie AP=PD, je calcule P−1, puis P−1U0 (ou P−1V0) avant de multiplier par Dn et par P.
  4. Je contrôle les expressions obtenues avec n=0 et n=1, en les comparant aux premiers termes calculés directement.
  5. J’interprète la limite dans le contexte (équilibre, extinction, croissance, proportions) et je cherche un seuil à la calculatrice ou avec une boucle Python à affectation simultanée.
Problème 1 · Deux bassins d’une station d’épuration★★★Moyen
Une station d’épuration comporte deux bassins, A et B, dont on relève les volumes d’eau, en m³, chaque heure. D’une heure à la suivante, le bassin A conserve 60 % de son volume et en transfère 30 % vers le bassin B, le reste étant évacué après traitement ; le bassin B conserve 50 % de son volume et en transfère 20 % vers le bassin A, le reste étant évacué. De plus, chaque heure, 500 m³ d’eau usée arrivent dans le bassin A et 150 m³ dans le bassin B. Au départ, le bassin A contient 2 600 m³ et le bassin B 1 100 m³. On note an et bn les volumes après n heures, Un =(anbn), et on donne P =(121−3) et D =(0,8000,3).
  1. 1. Traduire l’énoncé par une relation Un+1=AUn+B, en justifiant chaque coefficient. Matrice A : ligne du haut : et ; ligne du bas : et
    Rédaction attendue. Le bassin A garde 60% de son eau, reçoit 20% de celle de B et 500 m³ : an+1=0,6an+0,2bn+500. Le bassin B reçoit 30% de A, garde 50% de son eau et reçoit 150 m³ : bn+1=0,3an+0,5bn+150. Donc A =(0,60,20,30,5) et B =(500150).
  2. 2. U1 : a1= et b1=
    a1 =0,6×2600+0,2×1100+500 =1560+220+500 =2280 et b1 =0,3×2600+0,5×1100+150 =780+550+150 =1480.
  3. 3. U2 : a2= et b2=
    a2=1368+296+500=2164 et b2=684+740+150=1574.
  4. 4. Démontrer que la relation admet un unique état stable C et le déterminer : et
    Rédaction attendue. C =(xy) vérifie (I−A)C=B avec I−A =(0,4−0,2−0,30,5), de déterminant 0,4×0,5−0,2×0,3 =0,14≠0 : I−A est inversible, l’état stable existe et il est unique. Le système 0,4x−0,2y=500, −0,3x+0,5y=150 devient 2x−y=2500 et −3x+5y=1500 ; avec y=2x−2500 : 7x=14000, x=2000 et y=1500.
  5. 5. Démontrer que la suite de terme général Vn=Un−C vérifie Vn+1=AVn, puis donner V0 : et
    Rédaction attendue. C=AC+B, donc Vn+1 =(AUn+B)−(AC+B) =A(Un−C) =AVn. La suite (Vn) est homogène : Vn=AnV0, avec V0 =U0−C =(600−400).
  6. 6. Vérifier que AP=PD, puis calculer P−1V0 : et
    AP =(0,6+0,21,2−0,60,3+0,50,6−1,5) =(0,80,60,8−0,9) =PD. det(P)=−3−2=−5, donc P−1 =15(321−1) et P−1V0 =15(1800−800600+400) =(200200).
  7. 7. En déduire l’expression de an sous la forme an =γ+α×0,8n+β×0,3n, avec α= et β=
    An=PDnP−1, donc Vn =P(200×0,8n200×0,3n) =(200×0,8n+400×0,3n200×0,8n−600×0,3n). Ainsi an =2000+200×0,8n+400×0,3n et bn =1500+200×0,8n−600×0,3n. Pour n=1 : 2000+160+120=2280 et 1500+160−180=1480.
  8. 8. Limite de (bn), en m³ :
    0,8n et 0,3n tendent vers 0, donc (an) tend vers 2000 et (bn) vers 1500. À long terme, les volumes se stabilisent à 2000 m³ et 1500 m³ : les arrivées compensent les transferts et les évacuations. Cet équilibre ne dépend pas des volumes de départ.
Problème 2 · Proies et prédateurs dans un lac★★★Moyen
Dans un lac, on suit deux populations de poissons, en milliers d’individus : une espèce herbivore et une espèce prédatrice. D’une année sur l’autre, un modèle « proie-prédateur » linéaire donne un+1=1,1un−0,2vn et vn+1=0,1un+0,8vn, où un est l’effectif des herbivores et vn celui des prédateurs après n années, avec u0=900 et v0=400. On pose Un =(unvn), P =(2111) et D =(1000,9).
  1. 1. Matrice A telle que Un+1=AUn : ligne du haut : et ; ligne du bas : et
    On lit le système ligne par ligne : A =(1,1−0,20,10,8). Le coefficient −0,2 traduit l’effet des prédateurs sur les herbivores.
  2. 2. U2 : u2= et v2=
    u1=990−80=910 et v1=90+320=410, puis u2=1001−82=919 et v2=91+328=419.
  3. 3. Vérifier que AP=PD, puis calculer l’inverse de P. Coefficients de P−1, ligne du haut : et ; ligne du bas : et
    AP =(2,2−0,21,1−0,20,2+0,80,1+0,8) =(20,910,9) et PD =(20,910,9). det(P)=2−1=1, donc P−1 =(1−1−12).
  4. 4. Justifier que Un=PDnP−1U0, puis calculer P−1U0 : et
    Rédaction attendue. La suite est homogène, donc Un=AnU0 d’après le cours. De AP=PD on tire A=PDP−1 (on multiplie à droite par P−1), puis An=PDnP−1 par récurrence, les facteurs P−1P se simplifiant. Donc Un=PDnP−1U0, avec P−1U0 =(900−400−900+800) =(500−100).
  5. 5. un=1000−k×0,9n avec k=
    DnP−1U0 =(500−100×0,9n) car 1n=1, puis Un =(1000−100×0,9n500−100×0,9n). Pour n=2 : 1000−81=919 et 500−81=419.
  6. 6. Démontrer que la suite (un) est strictement croissante. un+1−un=c×0,9n avec c=
    Rédaction attendue. un+1−un =−100×0,9n+1+100×0,9n =100×0,9n(1−0,9) =10×0,9n. Or 0,9n>0 pour tout n, donc un+1−un>0 : (un) est strictement croissante.
  7. 7. Limite de (vn), en milliers :
    −1<0,9<1, donc 0,9n tend vers 0 : (un) tend vers 1000 et (vn) vers 500. À long terme, environ 1000000 herbivores et 500000 prédateurs : deux herbivores pour un prédateur. C’est le coefficient 1 de D qui donne cet équilibre non nul.
  8. 8. Recopier et compléter la fonction ci-dessous pour qu’elle renvoie le plus petit entier n tel que un>990.
    def seuil():
        u, v = ..., ...
        n = 0
        while ...:
            u, v = ..., ...
            n = n + 1
        return n
    Valeur renvoyée :
    def seuil():
        u, v = 900, 400
        n = 0
        while u <= 990:
            u, v = (1.1*u - 0.2*v,
                    0.1*u + 0.8*v)
            n = n + 1
        return n
    L’affectation simultanée est indispensable : u et v doivent être calculés à partir des anciennes valeurs. La fonction renvoie 22. Par le calcul : un>990 équivaut à 100×0,9n<10, soit 0,9n<0,1, avec 0,921≈0,109 et 0,922≈0,098.
Problème 3 · Diffusion d’un médicament entre deux compartiments★★★Difficile
Un patient reçoit une perfusion continue. On modélise l’organisme par deux compartiments, le sang et les tissus, et on relève les masses de principe actif, en mg, toutes les heures. D’une heure à la suivante, le sang conserve 70 % de sa masse et en transmet 20 % aux tissus, le reste étant éliminé ; les tissus conservent 80 % de leur masse et en restituent 10 % au sang, le reste étant éliminé. Chaque heure, la perfusion apporte 12 mg dans le sang. Au départ, l’organisme ne contient pas de principe actif. On note sn et tn les masses dans le sang et dans les tissus après n heures, Un =(sntn), et on donne P =(112−1) et D =(0,9000,6).
  1. 1. Déterminer les matrices A et B de la relation Un+1=AUn+B. Coefficients de A, ligne du haut : et ; ligne du bas : et
    Le sang reçoit 70% de lui-même, 10% des tissus et 12 mg : sn+1=0,7sn+0,1tn+12. Les tissus reçoivent 20% du sang et gardent 80% : tn+1=0,2sn+0,8tn. Donc A =(0,70,10,20,8) et B =(120).
  2. 2. Valeurs exactes de U3 : s3= et t3=
    U0 est nulle, donc U1 =B =(120). Puis s2=8,4+12=20,4 et t2=2,4. Enfin s3 =0,7×20,4+0,1×2,4+12 =14,28+0,24+12 =26,52 et t3=4,08+1,92=6.
  3. 3. Démontrer que C =(6060) est l’unique état stable. det(I−A)=
    Rédaction attendue. AC+B =(42+6+1212+48) =(6060) =C : C est un état stable. I−A =(0,3−0,1−0,20,2) a pour déterminant 0,3×0,2−0,1×0,2 =0,04≠0 : I−A est inversible, l’équation (I−A)X=B a une seule solution, et C est l’unique état stable.
  4. 4. Pour tout n, soit Vn=Un−C. Démontrer que Vn=AnV0, puis calculer P−1V0 : et
    Rédaction attendue. Vn+1 =(AUn+B)−(AC+B) =A(Un−C) =AVn : la suite (Vn) est homogène, donc Vn=AnV0, avec V0 =(−60−60). AP =(0,90,61,8−0,6) =PD et det(P)=−3, donc P−1 =13(112−1) et P−1V0 =13(−120−60) =(−40−20).
  5. 5. tn =60+c×0,9n+d×0,6n avec c= et d=
    Vn =PDnP−1V0 =P(−40×0,9n−20×0,6n) =(−40×0,9n−20×0,6n−80×0,9n+20×0,6n). Donc sn =60−40×0,9n−20×0,6n et tn =60−80×0,9n+20×0,6n. Pour n=2 : 60−32,4−7,2=20,4 et 60−64,8+7,2=2,4.
  6. 6. Limite commune de (sn) et de (tn), en mg :
    0,9n et 0,6n tendent vers 0 : les deux masses tendent vers 60 mg. À perfusion constante, l’organisme atteint un régime d’équilibre où l’apport compense exactement l’élimination.
  7. 7. À l’aide de la calculatrice, déterminer le nombre d’heures à partir duquel la masse dans le sang dépasse 55 mg :
    sn>55⟺40×0,9n+20×0,6n<5, et ce membre de gauche décroît. Pour n=19 il vaut environ 5,40 (s19≈54,6) ; pour n=20, environ 4,86 (s20≈55,1). La masse dans le sang dépasse 55 mg au bout de 20 heures. Le terme en 0,6n est déjà négligeable : c’est 0,9n qui fixe la durée de mise à l’équilibre.
Problème 4 · Escargots répartis en deux classes d’âge★★★Difficile
Une population d’escargots est recensée chaque année et répartie en jeunes (moins d’un an) et adultes. Modèle retenu : chaque adulte donne naissance en moyenne à 3 jeunes par an ; la moitié des jeunes survit et devient adulte l’année suivante ; la moitié des adultes survit d’une année sur l’autre ; les jeunes ne se reproduisent pas. On note jn et an les nombres de jeunes et d’adultes après n années, avec j0=100 et a0=60, et Un =(jnan). On donne A =(030,50,5), P =(231−1) et D =(1,500−1).
  1. 1. Justifier que Un+1=AUn, puis donner U1 : j1= et a1=
    Rédaction attendue. Les jeunes de l’année n+1 sont les naissances de l’année, dues aux seuls adultes : jn+1=3an. Les adultes de l’année n+1 sont les jeunes devenus adultes et les adultes survivants : an+1=0,5jn+0,5an. Ce système s’écrit Un+1=AUn avec A =(030,50,5). Ainsi j1=180 et a1=50+30=80.
  2. 2. U3 : j3= et a3=
    U2 =(24090+40) =(240130), puis U3 =(390120+65) =(390185).
  3. 3. Vérifier que AP=PD. P−1=15N, coefficients de N : ligne du haut : et ; ligne du bas : et
    AP =(0+30−31+0,51,5−0,5) =(3−31,51) et PD =(3−31,51). det(P)=−2−3=−5, donc P−1 =−15(−1−3−12) =15(131−2).
  4. 4. jn =112×1,5n+k×(−1)n avec k=
    Un=AnU0=PDnP−1U0, et P−1U0 =15(100+180100−120) =(56−4). Donc Un =P(56×1,5n−4(−1)n) : jn=112×1,5n−12(−1)n et an=56×1,5n+4(−1)n. Pour n=2 : 112×2,25−12=240 et 126+4=130.
  5. 5. Démontrer que le quotient jnan converge. Limite :
    Rédaction attendue. En divisant numérateur et dénominateur par 1,5n, non nul : jnan =112−12qn56+4qn avec q=−11,5=−23. Comme −1<q<1, qn tend vers 0, et par quotient jnan tend vers 11256=2. À long terme, il y a deux jeunes pour un adulte, quelle que soit la répartition de départ ; la population, elle, est multipliée par environ 1,5 chaque année.
  6. 6. La fonction suivante renvoie le nombre d’années au bout duquel la population totale dépasse 10000 escargots.
    def annee():
        j, a = 100, 60
        n = 0
        while j + a <= 10000:
            j, a = (3*a,
                    0.5*j + 0.5*a)
            n = n + 1
        return n
    Valeur renvoyée :
    La population totale vaut jn+an =168×1,5n−8(−1)n : environ 9680 pour n=10 et 14540 pour n=11. Le seuil est franchi la onzième année, et la fonction renvoie 11.
Problème 5 · Une suite récurrente d’ordre 2★★★Difficile
On considère la suite (un) définie par u0=0, u1=1 et, pour tout entier naturel n, un+2=un+1+2un. Chaque terme dépend des deux précédents : on se ramène à une suite de matrices colonnes en posant Xn =(un+1un). On donne P =(211−1) et D =(200−1).
  1. 1. u6=
    u2=1+0=1, u3=1+2=3, u4=3+2=5, u5=5+6=11, u6=11+10=21.
  2. 2. Justifier que Xn+1=AXn et donner A : ligne du haut : et ; ligne du bas : et
    Rédaction attendue. Xn+1 =(un+2un+1) =(un+1+2unun+1). La première ligne s’écrit 1×un+1+2×un et la seconde 1×un+1+0×un : Xn+1=AXn avec A =(1210).
  3. 3. Vérifier que AP=PD, puis calculer la matrice 3P−1 : ligne du haut : et ; ligne du bas : et
    AP =(2+21−221) =(4−121) et PD =(4−121). det(P)=−2−1=−3, donc P−1 =−13(−1−1−12) =13(111−2) et 3P−1 =(111−2).
  4. 4. En déduire que un=2n−(−1)nk pour tout entier naturel n, avec k=
    Xn=AnX0=PDnP−1X0 avec X0 =(10) et P−1X0 =13(11). Alors Xn =13P(2n(−1)n) =13(2n+1+(−1)n2n−(−1)n). Le terme un est la seconde coordonnée : un=2n−(−1)n3. Contrôle : u6=64−13=21.
  5. 5. Démontrer que un+1+un=2n pour tout n. En déduire u8 à partir de u6, sans la formule de un :
    Rédaction attendue. un+1+un =2n+1−(−1)n+1+2n−(−1)n3 =3×2n+(−1)n−(−1)n3 =2n, car −(−1)n+1=(−1)n. Donc u7=26−u6=64−21=43, puis u8=27−u7=128−43=85. Contrôle : 256−13=85.
  6. 6. Limite de un+1un (pour n⩾1) :
    En divisant par 2n : un+1un =2−(−1)n+1×0,5n1−(−1)n×0,5n. Comme 0,5n tend vers 0 et que (−1)n reste entre −1 et 1, le quotient tend vers 2 : à long terme, chaque terme vaut environ le double du précédent.
  7. 7. Valeur renvoyée par terme(10), avec la fonction suivante :
    def terme(n):
        a, b = 0, 1
        for k in range(n):
            a, b = b, b + 2*a
        return a
    À chaque passage dans la boucle, le couple (a,b) passe de (uk,uk+1) à (uk+1,uk+2) : c’est exactement Xk+1=AXk. Après 10 passages, a vaut u10=1024−13=341.
Bloqué sur ce chapitre ?

Moi c'est Ben, j'ai créé ce site pour aider le plus d'élèves possible, partout en France. Je peux aussi reprendre avec toi ce qui coince en cours particulier : en visio, ou à domicile si tu es à Marseille.

Prof de maths à Marseille · Cours particuliers au lycée · Aide aux devoirs