DERNIÈRE IMPRESSION LE 22 mai 2016 à 11:38

Matrices et suites

Table des matières

1 Matrice 2
1.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Opération sur les matrice . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2.1 Addition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2.2 Multiplication par un scalaire (réel) . . . . . . . . . . . . . . 4
1.2.3 Transposition d’une matrice . . . . . . . . . . . . . . . . . . . 4
1.2.4 Produit de deux matrices . . . . . . . . . . . . . . . . . . . . 4
1.3 Écriture matricielle d’un système linéaire . . . . . . . . . . . . . . . 6
1.4 Inversion d’une matrice . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4.2 Condition pour qu’une matrice d’ordre 2 soit inversible . . 7
1.5 Puissance n-ième d’une matrice carrée d’ordre 2 ou 3 . . . . . . . . 8
1.6 Diagonalisation d’une matrice d’ordre 2 . . . . . . . . . . . . . . . . 9

2 Étude de suite à l’aide de matrice 11

2.1 Un premier exemple : un système fermé . . . . . . . . . . . . . . . . 11
2.1.1 Le problème . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.1.2 Résolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.2 Étude d’une suite du type : Xn+1 = Xn M + B . . . . . . . . . . . . . . 13
2.2.1 Le problème . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2.2 Résolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

3 Marche aléatoire 16
3.1 Marche aléatoire simple sur un segment . . . . . . . . . . . . . . . . 17
3.2 Marche aléatoire aux sommets d’un tétraèdre . . . . . . . . . . . . . 17
3.3 Un retour en arrière est-il possible ? . . . . . . . . . . . . . . . . . . . 18

4 Traitement de l’image 19
4.1 Numériser les images . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.2 Opérations sur les images . . . . . . . . . . . . . . . . . . . . . . . . 19

PAUL MILAN 1 TERMINALE S SPÉ

 TABLE DES MATIÈRES

Introduction
Le but de ce chapitre est de résoudre quelques problèmes liés à des variables
discrète par l’intermédiaire d’un nouvel outil que constitue les matrices. Il s’agit
de mettre en évidence la pertinence d’introduire des matrices pour résoudre quel-
ques problèmes concrets. Bien que l’introduction des matrices dans le nouveau
programme doit se faire "naturellement", il me semble préférable d’introduire
directement les matrices, en donnant des exemples concrets, sans théorie exces-
sives, afin ensuite de traiter quelques exemples cités dans le programme.

1 Matrice
1.1 Définition

Définition 1 : Une matrice M(m × n) est un tableau de nombres possèdant

m lignes et n colonnes. On écrit alors :
 
a11 a12 . . . a1n
 a21 a22 . . . a2n 
M= ... ...

... ...
am1 am2 . . . amn

Les nombres aij sont les éléments ou coefficients de la matrice M. aij est situé
à l’intersection
  de la ie ligne et de la je colonne. On note parfois la matrice M
par aij .

 
1 2 0
Exemple : Soit A matrice (2 × 3) définie par :
4 3 −1
On a par exemple les coefficients a21 = 4 et a13 = 0
Application : Voici les productions (en milliers) de deux usines de cycles appar-
tenant à une même enseigne pour le premier semestre de l’année 2010 :

VTT adultes Vélos enfants VTC BMX Vélos de course

Usine 1 12,99 13,20 5,58 1,53 1,95
Usine 2 4,62 4,98 2,16 0,51 0,78

On peut alors créer une matrice production P dont les lignes correspondent aux
usines et les colonnes aux différents type de cycles. La matrice P est alors une
matrice (2 × 5) :
 
12, 99 13, 20 5, 58 1, 53 1, 95
P=
4, 62 4, 98 2, 16 0, 51 0, 78

Remarque : Quelques matrices particulières

• Si m = 1, la matrice M est appelée matrice ou vecteur ligne, par exemple :


M= 1 5 8

PAUL MILAN 2 TERMINALE S SPÉ

 1. MATRICE

• Si n = 1 , la matrice M est appelée matrice ou vecteur colonne, par exemple :

 
1
M =  3
−4
• Si m = n, la matrice M est appelée matrice carrée d’ordre m. Par exemple la
matrice carrée d’ordre 2 :  
4 5
M=
3 −2
• Une matrice carrée est symétrique si et seulement si aij = a ji ∀i 6= j. Par
exemple la matrice symétrique d’ordre 2 :
 
4 −1
M=
−1 4
• On définit la matrice unité Im d’ordre m par la matrice carrée d’ordre m qui
possède que des "1" sur sa diagonale et des "0" ailleurs. Par exemple la matrice
unité d’ordre 3 :  
1 0 0
I 3 = 0 1 0
0 0 1
• On définit une matrice diagonale d’ordre m par la matrice carrée d’ordre m qui
ne possède des éléments non nuls que sur sa diagonale. Par exemple :
 
2 0 0
D = 0 1 0
0 0 −3
• On définit une matrice triangulaire d’ordre m par une matrice carrée d’ordre m
qui possède un triangle composé uniquement de "0". Si la diagonale est compo-
sée de "0", on dit alors que la matrice est strictement triangulaire. Par exemple :
   
1 4 5 0 4 5
T = 0 2 7 T = 0 0 7
0 0 6 0 0 0
Matrice triangulaire supérieure Matrice strictement triangulaire

1.2 Opération sur les matrice

1.2.1 Addition

Définition 2 : L’addition ou la soustraction de deux matrices de même di-

mension A et B est égale à la matrice C dont chaque coefficient est obtenu en
additionnant ou soustrayant chaque coefficient de la matrice A au coefficient
correspondant de la matrice B. Par exemple :
     
1 2 0 5 2 3 6 4 3
+ =
4 3 −1 1 3 4 5 6 3

Remarque : L’addition de deux matrices ne posent donc aucun problème.

PAUL MILAN 3 TERMINALE S SPÉ

 TABLE DES MATIÈRES

1.2.2 Multiplication par un scalaire (réel)

Définition 3 : Le produit de la matrice A par un réel λ, est égal à la matrice

B dont chaque coefficient est obtenu en multipliant chaque coefficient de la
matrice A par λ. Par exemple :
   
1 2 0 2 4 0
2· =
4 3 −1 8 6 −2

Remarque : Cette opération ne pose donc aucun problème. Ces deux opérations
sont identiques à celles utilisées par les vecteurs. Les matrices et les vecteurs ont
donc une même structure appelée : espace vectoriel sur R.

1.2.3 Transposition d’une matrice

Définition 4 : La transposée t M d’une matrice M(m × n) est la matrice

(n × m) obtenue en échangeant les lignes et les colonnes de la matrice M. Par
exemple  
  1 4
1 2 0
Si M = alors t M = 2 3
4 3 −1
0 −1

Remarque :
• La transposée d’un vecteur colonne est un vecteur ligne
• Si M est une matrice carrée symétrique, alors : t M = M

1.2.4 Produit de deux matrices

Définition 5 : Le produit d’un vecteur ligne par un vecteur colonne est

égal à la somme des produits de chaque coefficient du vecteur ligne avec le
coefficient correspondant du vecteur colonne. Par exemple :
 

5
4 3 −1 × 2 = 4 × 5 + 3 × 2 − 1 × 3 = 23

3

Remarque : Cette opération correspond au produit scalaire de deux vecteurs.

On généralise cette opération à deux matrices quelconques A et B pourvu que

le nombre de colonnes de la matrice A correspondent au nombre de lignes de la
matrice B.

PAUL MILAN 4 TERMINALE S SPÉ

 1. MATRICE

Définition 6 : Le produit de la matrice A(m × n) par la matrice B(n × p)

est égal à la matrice C(m × p) dont chaque coefficient cij est égal au produit
scalaire de la ligne i de la matrice A par la colonne j de la matrice B. Par
exemple :
 
  5 1  
1 2 0 1 × 5 + 2 × 2 + 0 × 3 1 × 1 + 2 × 3 + 0 × 5
× 2 3 =
4 3 −1 4 × 5 + 3 × 2−1 × 3 4 × 1 + 3 × 3−1 × 5
3 5
 
9 7
=
23 8

Remarque : Le produit de deux matrices est :

• associatif : A (B × C) = (A × B) C = ABC
• distributif par rapport à l’addition : A (B + C) = AB + AC
• non commutatif : AB 6= BA en général.
Exemple : Une association de consommateurs compare les prix de cinq pro-
duits p1 , p2 , p3 , p4 , p5 distincts dans trois magasins différents. Les observations
fournissent les données suivantes :

Produit p1 Produit p2 Produit p3 Produit p4 Produit p5

magasin 1 1 5 2 3 4
magasin 2 1,1 4,7 1,8 3,1 3,8
magasin 3 0,9 5,1 1,9 3,2 4

On peut stocker les prix des produits sous la forme d’une matrice P (3 × 5).
 
1 5 2 3 4
P = 1, 1 4, 7 1, 8 3, 1 3, 8
0, 9 5, 1 1, 9 3, 2 4

Pour comparer la dépense d’une ménagère selon les magasins, on considère un «

panier » indiquant pour chaque produit la quantité achetée.
On appelle q1 , q2 , q3 , q4 et q5 , les quantités correspondant aux 5 produits. par
exemple 2, 1, 3, 3, 2
Le panier d’une ménagère peut être représenté par un vecteur colonne Q (1 × 5) :
   
q1 2
 q2  1
   
Q=  q3  = 3
  
 q4  3
q5 2

Soit Π1 , Π2 et Π3 les prix du panier de la ménagère dans chacun des trois maga-
sins. On note alors Π (1 × 3) le vecteur colonne correspondant à ces trois prix :
 
Π1
Π = Π2 

Π3

PAUL MILAN 5 TERMINALE S SPÉ

 TABLE DES MATIÈRES

On peut donc traduire, le prix du panier de la ménagère dans chacun des trois
magasin par l’égalité matricielle suivante :

Π = P×Q

En remplaçant par les données de notre exemple, on a :

 
    2
Π1 1 5 2 3 4 1
 
Π2  = 1, 1 4, 7 1, 8 3, 1 3, 8 × 3
 
Π3 0, 9 5, 1 1, 9 3, 2 4 3
2

Ce qui donne :
     
Π1 1×2+5×1+2×3+3×3+4×2 30
Π2  = 1, 1 × 2 + 4, 7 × 1 + 1, 8 × 3 + 3, 1 × 3 + 3, 8 × 2 = 29, 2
Π3 0, 9 × 2 + 5, 1 × 1 + 1, 9 × 3 + 3, 2 × 3 + 4 × 2 30, 2

1.3 Écriture matricielle d’un système linéaire

Définition 7 : Soit le système (S) linéaire (n × n) suivant :



 a11 x1 + a12 x2 + · · · + a1n xn = b1

a x +a x +···+a x = b
21 1 22 2 2n n 2

 ...


an1 x1 + an2 x2 + · · · + ann xn = bn
     
a11 a12 ... a1n x1 b1
 a21 a22 ... a2n   x2   b2 
On pose M = 
. . .
, X=  et B= 
... ... . . . . . . . . .
an1 an2 ... ann xn bn
On a alors l’écriture matricielle du système (S) est :

MX = B

(
2x − 3y = 5
Exemple : Soit le système suivant :
5x − 4y = 1

Son écriture matricielle est donc :

    
2 −3 x 5
=
5 −4 y 1

PAUL MILAN 6 TERMINALE S SPÉ

 1. MATRICE

1.4 Inversion d’une matrice

1.4.1 Définition

Définition 8 : Une matrice carrée M est dites inversible (ou régulière) si,
et seulement si, il existe une matrice carrée, appelée matrice inverse et notée
M−1 , telle que :
M × M −1 = M −1 × M = I
Si M−1 n’existe pas, on dit que la matrice M est singulière

 
4 3
Exemple : Soit la matrice A carrée d’ordre 2, définie par : .
2 1
 
−0, 5 1, 5
Montrons que la matrice B définie par : est la matrice inverse de
1 −2
la matrice A.

       
4 3 −0, 5 1, 5 −2 + 3 6 − 6 1 0
A×B = × = =
2 1 1 −2 −1 + 1 3 − 2 0 1
       
−0, 5 1, 5 4 3 −2 + 3 −1, 5 + 1, 5 1 0
B×A = × = =
1 −2 2 1 4−4 3−2 0 1

1.4.2 Condition pour qu’une matrice d’ordre 2 soit inversible

Définition 9 : Soit M une matrice carrée d’ordre 2, on appelle déterminant

de la matrice M, noté det(M), le nombre réel tel que :
   
a b a b
M= alors det(M) =    = ad − bc
c d c d

Exemple : Pour la matrice A précédente, on a :

 
4 3
det(A) =   = 4 × 1 − 3 × 2 = −2
2 1

Théorème 1 : Une matrice carrée d’ordre deux est inversible si et seulement

si son déterminant est différent de 0.

M−1 existe ⇔ det(M) 6= 0

On a alors :
   
a b −1 1 d −b
M= et det(M) 6= 0 alors M =
c d det(M) −c a

PAUL MILAN 7 TERMINALE S SPÉ

 TABLE DES MATIÈRES

 
a b
Démonstration : On se donne une matrice A = et on cherche alors
  c d
x y
une matrice B = telle que :
z t

A × B = B × A = I2

La première égalité se traduit par :

     
a b x y 1 0
× =
c d z t 0 1

On obtient alors les deux systèmes suivants :

( (
ax + bz = 1 ay + bt = 0
et
cx + dz = 0 cy + dt = 1

Ces systèmes admettent des solutions si leur déterminant est non nul. La condi-
tion est donc :  
a b
 c d 6= 0 ⇔ det(A) 6= 0
 

Par substitution, on obtient les solutions suivantes :

d −c −b a
x= z= y=
t=
ad − bc ad − bc ad − bc ad − bc
 
1 d −b
On obtient alors la matrice B suivante :
ad − bc −c a
On vérifie ensuite que l’égalité suivante est vérifiée : B × A = I2
 
4 3
Exemple : Déterminer la matrice inverse de la matrice A =
2 1
On calcule : det(A) = −2 (calculé plus haut). Comme le déterminant de la ma-
trice A est non nul, la matrice A admet une matrice inverse A−1 telle que :
   
−1 1 1 −3 −0, 5 1, 5
A = =
−2 −2 4 1 −2

On retrouve la matrice B de l’exemple de la définition

1.5 Puissance n-ième d’une matrice carrée d’ordre 2 ou 3

Définition 10 : On appelle puissance n-ième d’une matrice carrée M, la

matrice, notée Mn telle que :

Mn = M
| ×M×
{z · · · × M}
n fois

PAUL MILAN 8 TERMINALE S SPÉ

 1. MATRICE

 
1 2
Exemple : On donne A = . Calculer A2 et A3
3 0
     
1 2 1 2 7 2
A2 = × =
3 0 3 0 3 6
     
7 2 1 2 13 14
A3 = × =
3 6 3 0 21 6
 
1 2 3
Exemple : On donne B = 3 2 1. A l’aide de votre calculatrice, calculer B2 ,
2 3 1
3
B ainsi que la matrice inverse (valeurs approchée décimale).
Pour la TI 82, on sélectionne la touche ≥. Il faut d’abord éditer la matrice en
donnant la dimension (ici 3 × 3) puis rentrer les coefficients de la matrice en vali-
dant à chaque coefficient. On quitte, puis on sélectionne la matrice et on l’élève à
la puissance 2 puis 3. Pour trouver la matrice inverse, on utilise la touche x-1. On
trouve alors :
     
13 15 8 74 80 62 −0.086 −0.083 0.417
B2 = 11 13 12 B3 = 74 84 58 B−1 =  0.583 −0.417 0.083 
13 13 10 71 82 62 −0.333 0.666 −0.333
Remarque : Les problèmes rencontrés font intervenir des puissances de matrices.
On peut faire les calculs avec une calculatrice pour des matrices de petite taille et
des puissances raisonnables, on peut faire les calculs à l’aide d’un logiciel pour
des puissances explicites, à l’aide d’un logiciel de calcul formel pour obtenir, dans
les bons cas, des formules « closes », donnant l’expression des coefficients en
fonction de l’exposant, mais il est plus difficile d’obtenir, dans le cas général, ce
que nous avons appelé des limites.
Certaines matrices sont très bien adaptées pour calculer la puissance n-ième.
C’est le cas particulièrement des matrices diagonales. En effet, pour trouver la
puissance n-ième d’une matrice diagonale, il suffit simplement d’élever à la puis-
sance n les coefficients de la diagonale les autres coefficients restant nuls.
   n 
a 0 0 a 0 0
D = 0 b 0 alors Dn =  0 bn 0 
0 0 c 0 0 cn
On peut ainsi être amener à diagonaliser une matrice.

1.6 Diagonalisation d’une matrice d’ordre 2

Définition 11 : On dit qu’une matrice carrée d’ordre 2 est diagonalisable

s’il existe une matrice carrée P d’ordre 2 inversible et les réels α et β tels que :
 
α 0
A=P P −1
0 β

On a alors : 
αn 0
n
A =P P −1
0 βn

PAUL MILAN 9 TERMINALE S SPÉ

 TABLE DES MATIÈRES

Théorème 2 : Une matrice carrée A d’ordre 2 est diagonalisable si, et

seulement si, il existe deux réels α et β (non nécessairement distincts) et deux
vecteurs colonnes non colinéaires V et W telles que :

AV = αV et AW = βW

Les réels α et β (s’ils existent) s’appellent les valeurs propres de la matrice A.

Les vecteurs colonnes associés s’appellent alors les vecteurs propres.

 
1 2
Exemple : On pose A = .
2 4
Diagonaliser A puis déterminer An en fonction de n.
   
a c
On pose V = et W =
b d
On doit donc avoir AV = αV et AW = βW
ce qui se traduit par les systèmes suivants :
( (
a + 2b = αa c + 2d = βc
et
2a + 4b = αb 2c + 4d = βd

qui peuvent se mettre sous la forme :

( (
a(1 − α) + 2b = 0 c(1 − β) + 2d = 0
et
2a + b(4 − β) = 0 2c + d(4 − β) = 0

Si ces systèmes admettent une solution unique, cette solution ne peut-être que
(0, 0). Les vecteurs colonnes seront égaux au vecteur nul donc colinéaires, ce qui
est impossible. Donc ces systèmes doivent avoir une droite solution, ce qui n’est
possible que si le déterminant de ces systèmes est nul. Soit λ, le nombre réel
qui représente soit α soit β. On doit avoir un déterminant nul pour le système
générique suivant : (
a(1 − λ) + 2b = 0
2a + b(4 − λ) = 0
 
1 − λ 2 

 2 = 0 ⇔ (1 − λ)(4 − λ) − 4 = 0 ⇔ λ2 − 5λ = 0
4 − λ
Cette équation du second degré, admet deux solutions distinctes 0 et 5 qui sont
les valeurs de α et β pour diagonaliser la matrice A.
Remarque : si cette équation n’avait pas admis deux racines distinctes, la matrice
A n’aurait pas été diagonalisable.
Il reste à déterminer des vecteurs propres associés aux valeurs de α et β. Comme
les systèmes n’ont pas de solution unique, on cherche des coordonnées simples :

α = 0 donne a = −2b si b = 1 alors a = −2

β = 5 donne 2c = d si c = 1 alors d = 2

PAUL MILAN 10 TERMINALE S SPÉ

 2. ÉTUDE DE SUITE À L’AIDE DE MATRICE

On obtient donc les deux vecteurs colonnes suivants :

   
−2 1
V= et W =
1 2
On en déduit alors :
!
− 52 1
   
−2 1 1 2 −1 5
P= et sa matrice inverse P −1 = − =
1 2 5 −1 −2 1 2
5 5

On obtient alors :
!
− 52 1
 n    n   
n α 0 −1 −2 1 0 0 5 5n −1 2 × 5n −1
A =P P = =
0 βn 1 2 0 5n 1 2 2 × 5n −1 4 × 5n −1
5 5

Remarque : La diagonalisation est donc une opération complexe qui est utili-
sable que si les conditions s’y prêtent !

2 Étude de suite à l’aide de matrice

2.1 Un premier exemple : un système fermé
2.1.1 Le problème
On conserve dans une enceinte une population d’êtres unicellulaires qui ne
peuvent se trouver que dans deux états physiologiques désignés par A et B. On
désigne par an et bn les effectifs - exprimés en milliers d’individus - des deux
sous-populations (correspondant à chacun des deux états A et B) à l’instant n. Des
observations menées sur une assez longue période permettent d’estimer que 95%
des unicellulaires se trouvant à l’instant n dans l’état A n’ont pas changé d’état
à l’instant n + 1, non plus que 80% de ceux se trouvant à l’instant n dans l’état
B. L’effectif total s’élève à 500 000 individus. Cet effectif reste constant durant le
temps.
1) Traduire, avec des données, le système donnant an+1 et bn+1 en fonction de an
et bn
2) La population à l’instant 0 satisfait a0 = 375. A l’aide d’un tableur, faire le
calcul des effectifs an et bn pour les valeur de n de 0 à 50.
Peut-on faire une conjecture sur le comportement des suite ( an ) et bn ) ?
Effectuer de nouveaux essais en prenant d’autres valeurs initiales. Quelle nou-
velle conjecture peut-on faire sur le comportement des suite ( an ) et (bn ) ?
3) Traduire le système d’équation à l’aide de la notation matricielle. En déduire
l’écriture matricielle donnant an et bn en fonction des conditions initiales a0 et
b0 .
 
0, 95 0, 2
4) On donne pour M = , grâce à une diagonalisation, on obtient :
0, 05 0, 8
 
n 1 4 + 0, 75n 4 − 4 × 0, 75n
M =
5 1 − 0, 75n 1 + 4 × 0, 75n
a) Exprimer an et bn en fonction de n et de a0 .
b) Conclure

PAUL MILAN 11 TERMINALE S SPÉ

 TABLE DES MATIÈRES

2.1.2 Résolution
1) Si 95% des cellules dans l’état A reste dans cet état à l’instant n + 1 alors 5%
passe à l’état B. De même si 80% des cellules dans l’état B reste dans cet état à
l’instant n + 1 alors 20% passe à l’état A. On obtient donc le système suivant :
(
an+1 = 0, 95an + 0, 2bn
bn+1 = 0, 05an + 0, 8bn

2) Si la population à l’instant initial satisfait a0 = 375 alors b0 = 500 − 375 =

125. On rentre ces données dans les cellules B2 et C2. On écrit les formules
suivantes dans les cellules B3 et C3 puis on recopie ces formules vers le bas :

B3 : = 0.95 ∗ B2 + 0.2 ∗ C2 et C3 : = 0.05 ∗ B2 + 0.8 ∗ C2

On obtient alors le tableau suivant :

A B C
1 n an bn
2 0 375 125
3 1 381,25 118,75
4 2 385,937 5 114,062 5
5 3 389,453 125 110,546 875
6 4 392,089 843 8 107,910 156 3
7 5 394,067 382 8 105,932 617 2
... ... ... ...
... ... ... ...
35 33 399,998 116 5 100,001 883 5
36 34 399,998 587 4 100,001 412 6
37 35 399,998 940 5 100,001 059 5
38 36 399,999 205 4 100,000 794 6

Le système se stabilise et le nombre de cellules dans l’état A semble converger

vers 400 000 et dans l’état B vers 100 000.
Inversons par exemple les valeurs de départ. Prenons a0 = 125 et b0 = 375.
On obtient alors le tableau suivant :

A B C
1 n an bn
2 0 125 375
3 1 193,75 306,25
4 2 245,312 5 254,687 5
5 3 283,984 375 216,015 625
6 4 312,988 281 3 187,011 718 8
7 5 334,741 210 9 165,258 789 1
... ... ... ...
... ... ... ...
35 33 399,979 281 7 100,020 718 3
36 34 399,984 461 3 100,015 538 7
37 35 399,988 345 9 100,011 654 1
38 36 399,991 259 5 100,008 740 5

PAUL MILAN 12 TERMINALE S SPÉ

 2. ÉTUDE DE SUITE À L’AIDE DE MATRICE

On peut donc faire la conjecture suivante : « quel que soit l’état initiale, le
système converge vers 400 000 cellules dans l’état A et 100 000 dans l’état B ».
3) La traduction du système avec les notations matricielle donne :
        n  
a n +1 0, 95 0, 2 an an 0, 95 0, 2 a0
= soit =
bn + 1 0, 05 0, 8 bn bn 0, 05 0, 8 b0
4) a) D’après les données on obtient :
    
an 1 4 + 0, 75n 4 − 4 × 0, 75n a0
= n n
bn 5 1 − 0, 75 1 + 4 × 0, 75 b0
En prenant b0 = 500 − a0 , on obtient :
4 0, 75n 4 n 4(0, 75n )
a n = a0 + a0 + 400 − a0 − 400(0, 75 ) + a0
5 5 5 5
= 0, 75n ( a0 − 400) + 400

1 0, 75n 1 4(0, 75n )

bn = a0 − a0 + 100 − a0 + 400(0, 75n ) − a0
5 5 5 5
= 0, 75n (400 − a0 ) + 100

b) On a : lim 0, 75n = 0, car −1 < 0, 75 < 1

n→+∞
les suites ( an ) et (bn ) convergent, quelque soit l’état initial, vers les valeurs
respectives de 400 et 100 soit respectivement 400 000 et 100 000 cellules.

2.2 Étude d’une suite du type : Xn+1 = Xn M + B

2.2.1 Le problème
On estime que les patients admis dans un certain service d’un hôpital peuvent
se trouver dans l’un des 4 états suivants : 1. Soins réguliers, 2. Chirurgie, 3. Soins
intensifs, 4. Sortie.
Cette estimation est décrite par le tableau suivant, dans lequel sont indiquées les
probabilités de passage d’un des états à un autre dans un intervalle de 24 heures
(probabilités obtenues par modélisation des fréquences observées sur une longue
période).
Tableau de circulation des malades entre les services :

Soins réguliers Chirurgie Soins intensifs Sortie

Soins réguliers 0,6 0,2 0 0,2
Chirurgie 0,1 0 0,8 0,1
Soins intensifs 0,5 0 0,33 0,17
Sortie 0 0 0 0

Les informations chiffrées précédentes peuvent être stockées sous la forme d’une
matrice M (4 × 4) appelée matrice de transition :
 
0, 6 0, 2 0 0, 2
0, 1 0 0, 8 0, 1 
M= 0, 5 0 0, 33 0, 17


0 0 0 0

PAUL MILAN 13 TERMINALE S SPÉ

 TABLE DES MATIÈRES

On pourrait, à partir du tableau trace le graphe probabiliste suivant :

Sortie

0, 2 0, 17
0.1

0, 2
0.8
0, 6 Soins réguliers Chirurgie Soins intensifs 0, 33
0, 1
0, 5

Supposons qu’un certain jour n, la distribution

des patients suivant les quatre
états possibles s’écrive Xn = 12 5 6 3 . Le lendemain n + 1, la nouvelle dis-
tribution sera Xn+1 tel que
X n +1 = X n × M
Ce qui donne :  
0, 6 0, 2 0 0, 2

0, 1 0 0, 8 0, 1 

X n +1 = 12 5 6 3 × 
0, 5 0 0, 33 0, 17
 = 10, 7 2, 4 6 3, 9

0 0 0 0
Supposons qu’au jour 0, dix patients soient admis en soins réguliers
et qu’il n’y ait
aucun patient en cours de traitement. On note X0 = 10 0 0 0 la répartition
des malades le jour 0 et Xn la répartition des malades au n-ième jour, n entier
positif.
Supposons également que 10 patients soient admis chaque jour.
1) En utilisant la notation matricielle et votre calculatrice, déterminer la réparti-
tion des patients les jours 1 et 2 soit X1 et X2 .
2) Exprimer Xn+1 en fonction de Xn .
3) A l’aide d’un tableur, déterminer les répartitions pour les jours de 1 à 35. De
même à l’aide d’une calculatrice proposer un programme permettant à l’aide
de variables matricielles de donner la matrice répartition Xn pour les valeurs
suivante de n : n = 15, n = 30 et n = 50
Que constatez-vous ?
4) On admet que cette suite de matrices converge vers une répartition X. Dé-
terminer X à l’aide d’un calcul matriciel puis à l’aide de la résolution d’un
système et retrouver le résultat de la question précédente.

2.2.2 Résolution
1) On obtient :
 
0, 6 0, 2 0 0, 2

0, 1 0 0, 8 0, 1 


X1 = 10 0 0 0 0, 5 0 0, 33 0, 17 + 10 0 0 0 = 16 2 0 2


0 0 0 0
 
0, 6 0, 2 0 0, 2

0, 1 0 0, 8 0, 1 


X2 = 16 2 0 2 0, 5 0 0, 33 0, 17 + 10 0 0 0 = 19, 8 3, 2 1, 6 3, 4


0 0 0 0

PAUL MILAN 14 TERMINALE S SPÉ

 2. ÉTUDE DE SUITE À L’AIDE DE MATRICE

2) On obtient, en posant B = 10 0 0 0

X n +1 = X n M + B

3) On rentre la matrice X0 dans les cellules B2, C2, D2 et E2. On écrit les formules
suivantes dans les cellules B3, C3, D3 et E3 puis on recopie ces formules vers
le bas :
B3 : = 0, 6 ∗ B2 + 0, 1 ∗ C2 + 0, 5 ∗ D2 + 10

C3 : = 0, 2 ∗ B2 D3 : = 0, 8 ∗ C2 + 0, 33 ∗ D2

E3 : = 0, 2 ∗ B2 + 0, 1 ∗ C2 + 0, 17 ∗ D2
On obtient alors le tableau suivant :

A B C D E
1 n Soins régulier Chirurgie Soins intensifs Sortie
2 0 10 0 0 0
3 1 16 2 0 2
4 2 19,8 3,2 1,6 3,4
5 3 23 3,96 3,088 4,552
6 4 25,74 4,6 4,18704 5,52096
7 5 27,997 52 5,148 5,061 723 2 6,319 796 8
8 6 29,844 173 6 5,599 504 5,788 768 656 6,974 796 944
9 7 31,360 838 89 5,968 834 72 6,389 896 856 7,512 875 792
10 8 32,608 335 23 6,272 167 778 6,883 733 739 7,955 333 715
11 9 33,634 084 79 6,521 667 047 7,289 366 356 8,319 118 560
12 10 34,477 300 75 6,726 816 957 7,622 824 535 8,618 175 943
... ... ... ... ... ...
... ... ... ... ... ...
34 32 38,321 088 06 7,661 952 278 9,143 102 974 9,981 437 875
35 33 38,330 399 55 7,664 217 612 9,146 785 804 9,984 740 345
36 34 38,338 054 39 7,666 079 91 9,149 813 404 9,987 455 258
37 35 38,344 347 33 7,667 610 879 9,152 302 351 9,989 687 148

On appellera [A] la matrice M, la matrice [B] Variables : [A], [B], [C] matrices
la matrice B et la matrice [C] la matrice Xn . On N et I entiers naturels
peut alors proposer le programme ci-contre : Entrées et initialisation
Remarque : On aura auparavant rentré les Lire N
matrices [A] et [B] et dimensionné la matrice [B] → [C]
[C] en 1 × 4. Traitement
On trouve alors : pour I de 1 à N faire
X15 = 36, 911 7, 319 8, 585 9, 481

[C] × [A] + [B] → [C]

fin
X30 = 38, 296 7, 656 9, 133 9, 973

Sorties : Afficher [C]
X50 = 38, 372 7, 674 9, 163 9, 999

Ces résultats montrent que la situation tend à se stabiliser.

4) Si la suite de matrice converge, alors elle converge vers X tel que :

X = XM + B

Résolution matricielle : I désigne la matrice unité d’ordre 4.

• On isole la matrice X : X − XM = B

PAUL MILAN 15 TERMINALE S SPÉ

 TABLE DES MATIÈRES

• On factorise par X : X(I − M) = B

• On multiplie par la matrice inverse : X = B(I − M)−1
Un application numérique avec la calculatrice donne :


X = 38, 373 7, 675 9, 164 10, 000

Résolution par un système : Permet de déterminer les valeurs exactes

Si on pose : X = a b c d , on a donc :
 
0, 6 0, 2 0 0, 2


0, 1 0 0, 8 0, 1 

a b c d = a b c d   + 10 0 0 0
0, 5 0 0, 33 0, 17
0 0 0 0

On obtient alors le système suivant :

 

 0, 6a + 0, 1b + 0, 5c + 10 = a 
 − 0, 4a + 0, 1b + 0, 5c = −10 (1)

 0, 2a = b 
 a = 5b (2)
⇔

 0, 8b + 0, 33c = c 
 0, 8b − 0, 67c = 0 (3)

 

0, 2a + 0, 1b + 0, 17c = d d = 0, 2a + 0, 1b + 0, 17c (4)

De (2) a = 5b, on remplace dans les équations (1) et (3), on obtient alors le
système suivant :
( (
−1, 9b + 0, 5c = −10 −19b + 5c = −100
⇔
0, 8b − 0, 67c = 0 80b − 67c = 0

On obtient alors les valeurs :

6 700 8 000
b= ≃ 7,674 684 994 et c= ≃ 9,163 802 978
873 873
On en déduit à l’aides des équation (2) et (4) :
33 500
a= ≃ 38,373 424 97 et d = 10
873
Le système se stabilise donc par la répartition :
 
33 500 6 700 8 000
X= 10
873 873 873

3 Marche aléatoire
On s’intéresse au comportement à long terme d’une marche aléatoire. Il s’agit de
calculer les probabilités pour le héros d’une marche aléatoire dans un réseau de
se trouver après n pas en tel ou tel sommet (ou nœud) du réseau.

PAUL MILAN 16 TERMINALE S SPÉ

 3. MARCHE ALÉATOIRE

3.1 Marche aléatoire simple sur un segment

Le personnage se déplace d’un sommet à l’autre du graphe ci-dessous. S’il est en
A ou en B, il ne peut aller qu’en P, s’il est en P, il peut aller en A ou en B avec des
probabilités que nous considérons comme identiques.

Ab
P
b
Bb

On peut représenter la situation par une matrice M (dite de transition) qui in-
dique les probabilités de passage d’un sommet à un autre. La matrice M ci-
dessous représente la marche dans le réseau (A,P,B).

A P B
 
A 0 1 0
M = P  21 0 1
2
B 0 1 0

Les coefficients figurant sur chaque ligne donnent les probabilités de passage du
sommet qui donne son nom à la ligne à celui qui donne son nom à la colonne. La
diagonale ne contient de ce fait que des 0.
Si maintenant on s’intéresse aux probabilités de passer d’un sommet à un autre
en deux pas, cela revient à calculer les coefficients de la matrice M2 . On peut
généraliser ce résultat, les probabilités de passer d’un sommet à un autre en n pas
correspondent aux coefficient de la matrice Mn .
On obtient alors :
  1   
0 1 0 0 21
2
0 1 0
M =  21 0 21  M2 =  0 1 0  M3 =  21 0 12 
0 1 0 1 1 0 1 0
2 0 2

On constate que M3 = M
On a alors : M2k+1 = M et M2k = M ∀k ∈ N
On peut interpréter ce résultat (pour la première ligne) : par exemple, partant du
sommet A, le personnage est sûrement en P après un nombre impair de pas, en B
1
ou en A avec des probabilités après un nombre pair de pas.
2

3.2 Marche aléatoire aux sommets d’un tétraèdre

À la différence de la situation précé- A

dente, dans la marche aux sommets
d’un triangle comme dans la marche
aux sommets d’un tétraèdre, on peut
passer, à chaque étape, de tout sommet
donné à tout autre sommet donné.

On numérote les sommet A, B, C et D B

respectivement 1, 2, 3 et 4 D

C
PAUL MILAN 17 TERMINALE S SPÉ
 TABLE DES MATIÈRES

Dans l’hypothèse d’équiprobabilité, la matrice de transition M (donnant les pro-

babilités de passage d’un sommet à un autre) s’écrit :
 1 1 1
0 3 3 3
1 1 1
3 0 3 3
M=
1 1 1

3 3 0 3
1 1 1
3 3 3 0

Si on calcule les coefficients approchés des matrices M5 et M10 qui correspondent

à des chaînes de longueurs respectives de 5 et 10, on trouve :
   
0, 2469 0, 2510 0, 2510 0, 2510 0, 2500 0, 2500 0, 2500 0, 2500
0, 2510 0, 2469 0, 2510 0, 2510
 M10 = 0, 2500 0, 2500 0, 2500 0, 2500
 
M5 =  0, 2510 0, 2510 0, 2469 0, 2510 0, 2500 0, 2500 0, 2500 0, 2500
0, 2510 0, 2510 0, 2510 0, 2469 0, 2500 0, 2500 0, 2500 0, 2500

Les coefficients semblent se stabiliser tous autour de la valeur 0,25. Ce qui veut
dire que les différentes probabilités s’estompent rapidement. Les probabilités d’al-
ler d’un sommet quelconque à un autre sont très voisines de 0, 25 dès que n est
supérieur à 10.
Pour un résultat littéral, on peut montrer par récurrence que :
 
un vn vn vn
 vn un vn vn 
Mn =  vn vn un vn 


vn vn vn un

où les termes généraux des suites (un ) et (vn ) sont :

" #
1 n −1
     
1 1 1 n
un = 1− − et vn = 1− −
4 3 4 3

1
Ces deux suites convergent vers . On retrouve alors le résultat conjecturé.
4

3.3 Un retour en arrière est-il possible ?

Ayant quitté un sommet du tétraèdre, au bout de combien de pas aléatoires le
personnage peut-il compter y revenir ?
Soit X la variable aléatoire donnant, pour chaque marche, ce nombre de pas.
On a :
1 2 1
P ( X = 1) = 0 , P ( X = 2) = et P( X = 3) = ×
3 3 3
En effet, pour que le personnage soit en A, par exemple, après n pas sans y avoir
été dans aucune de ses positions précédentes, il est nécessaire qu’à chacun de
ses déplacements précédents il soit passe d’un sommet qui n’était pas A à un
autre qui n’était pas A non plus, choisissant donc l’un de deux sommets sur trois
possibles.
On peut alors vérifier par récurrence que :

PAUL MILAN 18 TERMINALE S SPÉ

 4. TRAITEMENT DE L’IMAGE

  n −2
2 1
P( X = n) = × pour n > 2.
3 3
2
Ce sont les termes d’une suite géométrique (à partir du rang 2) de raison et de
3
1
premier terme . On peut alors calculer la somme des (n − 1) termes à partir du
3
rang 2 :
  n −1
2
n 1−   n −1
1 3 2
∑ P( X = k) = 3 × = 1−
2 3
k =2 1−
3
Comme le terme au rang 1 est nul, en passant à la limite, on obtient :
n
lim∑ P( X = n) = 1
n→+∞
k =1

Ce qui veut dire que l’on est sur que l’on reviendra en A pourvu que l’on soit
patient !

4 Traitement de l’image
4.1 Numériser les images
On a extrait l’image ci-contre d’une photographie
d’Alan Turing disputant une course de 3 miles en
1946. Cette photographie a été reproduite sur un site
web consacré à l’un des « inventeurs » de l’infor-
matique dont l’adresse est donnée ci-dessous. Elle a
donc été « numérisée », c’est-à-dire transformée en
une suite de 0 et de 1. Le rectangle est décomposé en
un certain nombre de petits carrés, et à chacun de ces
carrés a été attribué un nombre qui représente une
nuance de gris. La finesse de la décomposition (le
nombre de carrés) est la définition de l’image. La dé-
finition de cette image particulière n’est pas bonne :
on devine les pixels (mot fabriqué avec les débuts des mots anglais picture ele-
ment).
Toute image n’utilisant que le noir et le blanc peut ainsi être représentée par un
tableau contenant autant de cases que l’image contient de pixels, chacune de ces
cases étant occupée par 0 ou 1. L’image est donc représentée par une matrice dont
tous les éléments sont 0 ou 1.
L’adresse du site consacré à Turing est :
http://www.turing.org.uk/turing/scrapbook/run.html

4.2 Opérations sur les images

Pour trouver le négatif B de l’image A, on remplace 1 par 0 et 0 par 1 dans la
matrice A

PAUL MILAN 19 TERMINALE S SPÉ

 TABLE DES MATIÈRES

   
1 0 0 1 1 0 0 1 1 0 0 1
1 1 0 0 0 1 0 0 1 1 1 0
   
1
A= 0 0 0 1 0 0
B= 1 1 1 0 1
1 0 0 0 0 0 0 1 1 1 1 1
0 1 1 0 1 0 1 0 0 1 0 1

On peut également coder des images en nuances de gris en attribuant à chaque

pixel un nombre compris entre 0 et 1, proche de 1 si la case est gris foncé, proche
de 0 si elle est gris clair. On peut également définir l’image négatif de l’image de
départ en lui associant la matrice dont les éléments sont les compléments à 1 des
éléments de la matrice de départ.

Les deux images ci-dessus sont le négatif l’une de l’autre. D’autres critères peuvent
être enregistrés dans les éléments de la matrice associée à une image, la lumino-
sité par exemple. Une multiplication de tous les éléments de la matrice représen-
tant la luminosité par un même facteur modifie la luminosité de l’ensemble.
Si deux images ont le même format et la même définition (associées aux ma-
trices A et B), il est possible de leur faire correspondre leur somme, associée à la
somme des matrices qui les définissent, en convenant qu’un coefficient supérieur
à 1 donne un pixel de couleur noire. On peut aussi leur faire correspondre leur
différence, avec cette fois la convention que tout pixel associé à un nombre négatif
est blanc, ou restituer l’image positive |A − B| en particulier pour différentier les
images et faire apparaître la trame des contours, horizontaux, verticaux, obliques.

PAUL MILAN 20 TERMINALE S SPÉ