Suite de matrices colonnes

icône de pdf
Signaler

I. Propriétés

On considère une suite (Un)(U_n) de matrices colonnes telle que
Un+1=AUn+BU_{n+1} = A U_n + B pour tout entier naturel nn.

\circ\quad S’il existe une matrice XX telle que AX+B=XAX + B = X, alors la suite (Vn)(V_n) définie par
Vn=UnXV_n = U_n - X vérifie : Vn+1=AVnV_{n+1} = A V_n

Dans ce cas, on a : Un=An(U0X)+XU_n = A^n (U_0 - X) + X pour tout entier nn.

\circ\quad Si (Un)(U_n) est une suite convergente, alors elle converge vers une matrice UU vérifiant
AU+B=UAU + B = U.

Démonstration :

Soit une suite (Un)(U_n) de matrices colonnes telle que
Un+1=AUn+BU_{n+1} = A U_n + B pour tout entier nn.

\circ\quad S’il existe une matrice XX telle que AX+B=XAX + B = X, on définit la suite (Vn)(V_n) par :
Vn=UnXV_n = U_n - X

Alors : Vn+1=Un+1X=AUn+B(AX+B)=AUnAX=A(UnX)=AVnV_{n+1} = U_{n+1} - X = A U_n + B - (AX + B) = A U_n - AX = A (U_n - X) = A V_n

La suite (Vn)(V_n) est donc une suite géométrique de raison AA, donc : Vn=AnV0V_n = A^n V_0 pour tout entier nn.

Or V0=U0XV_0 = U_0 - X, donc : Vn=An(U0X)V_n = A^n (U_0 - X) pour tout entier nn

Et comme Vn=UnXV_n = U_n - X, alors : Un=Vn+X=An(U0X)+XU_n = V_n + X = A^n (U_0 - X) + X pour tout entier nn.

\circ\quad Si (Un)(U_n) est une suite convergente, soit UU sa limite. Alors :
limn+Un+1=limn+Un=U\displaystyle\lim_{n \to +\infty} U_{n+1} = \lim_{n \to +\infty} U_n = U

D'une part :
limn+Un+1=limn+(AUn+B)=A(limn+Un)+B=AU+B\displaystyle\lim_{n \to +\infty} U_{n+1} = \lim_{n \to +\infty} (A U_n + B) = A \left( \lim_{n \to +\infty} U_n \right) + B = AU + B

D'autre part : limn+Un+1=U\displaystyle\lim_{n \to +\infty} U_{n+1} = U

Donc UU vérifie bien : AU+B=UAU + B = U

II. Un exercice rédigé

Dans un entrepôt logistique, deux zones de stockage, notées A et B, contiennent chacune 100 palettes de produits. Un logiciel ajuste automatiquement les quantités stockées chaque jour à 8h, selon les règles suivantes :

\checkmark 50 % des palettes en zone A partent définitivement sur camions pour livraison,

\checkmarkpuis 75 % sont transférées de B vers A ;

\checkmarkà chaque cycle, 200 palettes supplémentaires sont livrées en zone A, et 300 en zone B.

On souhaite savoir la taille des zones de stockage indispensable à gérer ce flux, afin que tout puisse toujours être réceptionné.

On note ana_n et bnb_n les quantités de palettes (en centaines) dans les zones A et B au jour nn.
On modélise ces données à l’aide de la suite de matrices colonnes Un=(anbn)U_n = \begin{pmatrix} a_n \\ b_n \end{pmatrix}, où ana_n et bnb_n représentent le nombre de palettes (en centaines) le jour nn. On sait donc que U0=(11)U_0 = \begin{pmatrix} 1 \\ 1 \end{pmatrix}.

On définit les matrices : M=(0,50,7500,25)M = \begin{pmatrix} 0{,}5 \quad 0{,}75 \\ 0 \quad 0{,}25 \end{pmatrix} et C=(23)C = \begin{pmatrix} 2 \\ 3 \end{pmatrix}

\circ\quad 1. Justifier que, pour tout entier naturel nn, on a : Un+1=MUn+CU_{n+1} = M U_n + C

On considère la matrice P=(1 031)P = \begin{pmatrix} 1 \quad ~0 \\ 3 \quad -1 \end{pmatrix}

\circ\quad 2a. Calculer P2P^2. En déduire que PP est inversible et préciser son inverse.

\circ\quad 2b. Montrer que la matrice MM est diagonalisable, en déterminant une matrice diagonale DD telle que M=PDP1M = P D P^{-1}.

\circ\quad 2c. En déduire que, pour tout entier naturel nn, Mn=PDnP1M^n = P D^n P^{-1}.

\circ\quad 2d. Montrer que Mn=(0,5n3×0,5n00,25n)M^n = \begin{pmatrix} 0{,}5^n \quad 3 \times 0{,}5^n \\ 0 \quad 0{,}25^n \end{pmatrix}

\circ\quad 3a. Montrer que la matrice X=(104)X = \begin{pmatrix} 10 \\ 4 \end{pmatrix} vérifie X=MX+CX = M X + C

\circ\quad 3b. Pour tout entier naturel nn, on définit Vn=UnXV_n = U_n - X. Montrer que Vn+1=MVnV_{n+1} = M V_n

\circ\quad 3c. En déduire que Vn=MnV0V_n = M^n V_0, où V0=U0XV_0 = U_0 - X

\circ\quad 3d. Montrer que Un=(18×0,5n+9×0,25n+103×0,25n+4)U_n = \begin{pmatrix} -18 \times 0{,}5^n +9\times 0{,}25^n+ 10 \\ -3 \times 0{,}25^n + 4 \end{pmatrix}

\circ\quad 4a. Montrer que la suite (bn)(b_n) est croissante et majorée. En déduire sa limite.

\circ\quad 4b. Déterminer la limite de la suite (an)(a_n)

\circ Conclure

Solution : principe de la démonstration

\circ\quad 1. Justifier que, pour tout entier naturel nn, on a Un+1=MUn+CU_{n+1} = M U_n + C

On modélise la situation :

\checkmarkÀ chaque cycle, 50 % des palettes de A partent (il en reste donc 50 %, soit 0,5an0{,}5 a_n).

\checkmark75 % des palettes de B sont transférées vers A (soit 0,75bn0{,}75 b_n).

\checkmark200 palettes (soit 22 centaines) sont livrées dans A, et 300 palettes (soit 33 centaines) dans B.

\checkmark25 % des palettes de B restent (soit 0,25bn0{,}25 b_n).

Ainsi, on a :

an+1=0,5an+0,75bn+2a_{n+1} = 0{,}5 a_n + 0{,}75 b_n + 2
bn+1=0,25bn+3b_{n+1} = 0{,}25 b_n + 3

Ce qui s’écrit sous forme matricielle :

Un+1=MUn+CU_{n+1} = M U_n + C
avec M=(0,50,7500,25)M = \begin{pmatrix} 0{,}5 \quad 0{,}75 \\ 0 \quad 0{,}25 \end{pmatrix} et C=(23)C = \begin{pmatrix} 2 \\ 3 \end{pmatrix}

\circ\quad 2a. Calculer P2P^2. En déduire que PP est inversible et préciser son inverse.

On a : P=(1031)P = \begin{pmatrix} 1 \quad 0 \\ 3 \quad -1 \end{pmatrix}

Alors : P2=(1031)(1031)=(1001)=I2P^2 = \begin{pmatrix} 1 \quad 0 \\ 3 \quad -1 \end{pmatrix} \begin{pmatrix} 1 \quad 0 \\ 3 \quad -1 \end{pmatrix} = \begin{pmatrix} 1 \quad 0 \\ 0 \quad 1 \end{pmatrix} = I_2

Donc PP est inversible et P1=PP^{-1} = P

\circ\quad 2b. Montrer que la matrice MM est diagonalisable, en déterminant une matrice diagonale DD telle que M=PDP1M = P D P^{-1}.

Pour que M=PDP1M = P D P^{-1}, il suffit que PMP=PPDP1PPMP=PPDP^{-1}P soit :

D=PMPD=PMP. Après calculs, on trouve que la matrice diagonale DD est :
D=(0,5000,25)D = \begin{pmatrix} 0{,}5 \quad 0 \\ 0 \quad 0{,}25 \end{pmatrix}

\circ\quad 2c. En déduire que, pour tout entier naturel nn, Mn=PDnP1M^n = P D^n P^{-1}

Ce résultat se démontre par récurrence, on montre que si M=PDP1M = P D P^{-1}, alors Mn=PDnP1M^n = P D^n P^{-1} pour tout nNn \in \mathbb{N}.

\circ\quad 2d. Montrer que Mn=(0,5n3×0,5n00,25n)M^n = \begin{pmatrix} 0{,}5^n \quad 3 \times 0{,}5^n \\ 0 \quad 0{,}25^n \end{pmatrix}

On a Dn=(0,5n000,25n)D^n = \begin{pmatrix} 0{,}5^n \quad\quad 0 \\ 0 \quad\quad 0{,}25^n \end{pmatrix}

On calcule Mn=PDnP1=PDnPM^n = P D^n P^{-1} = P D^n P (car P1=PP^{-1} = P)

Ce produit donne :

Mn=(0,5n3×0,5n3×0,25n00,25n)M^n = \begin{pmatrix} 0{,}5^n \quad\quad 3 \times 0{,}5^n -3\times 0,25^n\\ 0 \quad\quad\quad\quad\quad\quad\quad 0{,}25^n \end{pmatrix}

\circ\quad 3a. Montrer que la matrice X=(104)X = \begin{pmatrix} 10 \\ 4 \end{pmatrix} vérifie X=MX+CX = M X + C

On calcule :

MX=(0,5×10+0,75×40×10+0,25×4)=(5+31)=(81)M X = \begin{pmatrix} 0{,}5 \times 10 + 0{,}75 \times 4 \\ 0 \times 10 + 0{,}25 \times 4 \end{pmatrix} = \begin{pmatrix} 5 + 3 \\ 1 \end{pmatrix} = \begin{pmatrix} 8 \\ 1 \end{pmatrix}

MX+C=(81)+(23)=(104)=XM X + C = \begin{pmatrix} 8 \\ 1 \end{pmatrix} + \begin{pmatrix} 2 \\ 3 \end{pmatrix} = \begin{pmatrix} 10 \\ 4 \end{pmatrix} = X

Donc X=MX+CX = M X + C

\circ\quad 3b. Pour tout entier naturel nn, on définit Vn=UnXV_n = U_n - X. Montrer que Vn+1=MVnV_{n+1} = M V_n

On a : Vn+1=Un+1X=MUn+CX=M(UnX)+(MX+CX)V_{n+1} = U_{n+1} - X = M U_n + C - X = M (U_n - X) + (M X + C - X)

Or MX+C=XM X + C = X, donc : Vn+1=MVnV_{n+1} = M V_n

On obtient bien une suite géométrique de matrices.

\circ\quad 3c. En déduire que Vn=MnV0V_n = M^n V_0, où V0=U0XV_0 = U_0 - X

V0=(11)(104)=(93)V_0 = \begin{pmatrix} 1 \\ 1 \end{pmatrix} - \begin{pmatrix} 10 \\ 4 \end{pmatrix} = \begin{pmatrix} -9 \\ -3 \end{pmatrix}

En appliquant la propriété du cours, on obtient : Vn=MnV0V_n = M^n V_0

\circ\quad 3d. Montrer que Un=(18×0,5n+9×0,25n+103×0,25n+4)U_n = \begin{pmatrix} -18 \times 0{,}5^n + 9 \times 0{,}25^n + 10 \\ -3 \times 0{,}25^n + 4 \end{pmatrix}

On a :

Vn=MnV0=(0,5n3×0,5n3×0,25n00,25n)(93)V_n = M^n V_0 = \begin{pmatrix} 0{,}5^n \quad\quad 3 \times 0{,}5^n -3\times 0,25^n\\ 0 \quad\quad\quad\quad\quad\quad\quad 0{,}25^n \end{pmatrix} \begin{pmatrix} -9 \\ -3 \end{pmatrix}

Produit :

aa : 9×0,5n9×0,5n+9×0,25n=18×0,5n+9×0,25n-9 \times 0{,}5^n - 9 \times 0{,}5^n +9\times 0,25^n= -18 \times 0{,}5^n+9\times 0,25^n
bb : 3×0,25n-3 \times 0{,}25^n

Donc :

Vn=(18×0,5n+9×0,25n3×0,25n)V_n = \begin{pmatrix} -18 \times 0{,}5^n +9\times 0,25^n\\ -3 \times 0{,}25^n \end{pmatrix}

Et donc :

Un=Vn+X=(18×0,5n+9×0,25n+103×0,25n+4)U_n = V_n + X = \begin{pmatrix} -18 \times 0{,}5^n +9\times 0,25^n+ 10 \\ -3 \times 0{,}25^n + 4 \end{pmatrix}

\circ\quad 4a. Montrer que la suite (bn)(b_n) est croissante et majorée. En déduire sa limite.

On a : bn=3×0,25n+4b_n = -3 \times 0{,}25^n + 4

La suite (0,25n)(0{,}25^n) est strictement décroissante, donc (3×0,25n)(-3 \times 0{,}25^n) est strictement croissante.

Donc (bn)(b_n) est croissante.

De plus, 0,25n00{,}25^n \to 0, donc bn4b_n \to 4

Comme bn<4b_n \lt 4, la suite est majorée.

Conclusion : (bn)(b_n) est croissante et majorée, donc converge vers 44

\circ\quad 4b. Déterminer la limite de la suite (an)(a_n)

On a : an=18×0,5n+9×0,25n+10a_n = -18 \times 0{,}5^n + 9 \times 0{,}25^n + 10

Les deux premiers termes tendent vers 0, donc : liman=10\lim a_n = 10

\circ Conclusion

À long terme, les quantités de palettes se stabilisent :

\checkmarken zone A : an10a_n \to 10, soit 10001000 palettes

\checkmark en zone B : bn4b_n \to 4, soit 400400 palettes

Les capacités minimales des zones de stockage doivent donc être d’au moins 10001000 palettes pour A et 400400 pour B.