Skip to content

Exercice 1 : Calcul matriciel

Algèbre linéaire et Géométrie 3D : Calcul matriciel (Corrigé pas à pas)

Section titled “Algèbre linéaire et Géométrie 3D : Calcul matriciel (Corrigé pas à pas)”

Ce problème aborde un classique des concours de la filière PTSI : l’étude des puissances d’une matrice carrée d’ordre 33 à l’aide de polynômes annulateurs et de la formule du binôme de Newton.

Les compétences évaluées ici sont fondamentales pour la suite de votre cursus :

  1. La maîtrise du calcul matriciel et de la manipulation des symboles de sommation.
  2. L’exploitation rigoureuse d’une relation de récurrence matricielle.
  3. L’inversion de matrices par des méthodes polynomiales.

On considère la matrice AA à coefficients réels définie par : A=(111011001)A = \begin{pmatrix} 1 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 1 \end{pmatrix}

On cherche à déterminer l’expression générale de la puissance AnA^n pour tout entier naturel nn, puis à en déduire l’inverse de AA.


Partie 1 : Décomposition de la matrice et calculs préliminaires

Section titled “Partie 1 : Décomposition de la matrice et calculs préliminaires”

Question 1. Écrire la matrice AA sous la forme A=I3+NA = I_3 + N, où I3I_3 est la matrice identité d’ordre 33 et NN une matrice nilpotente dont on précisera l’ordre (c’est-à-dire le plus petit entier kk tel que Nk=0N^k = 0).

Corrigé : On pose : A=(100010001)+(011001000)=I3+NA = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} + \begin{pmatrix} 0 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} = I_3 + N avec N=(011001000)N = \begin{pmatrix} 0 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix}.

Calculons les puissances successives de NN : N2=(011001000)(011001000)=(001000000)N^2 = \begin{pmatrix} 0 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} 0 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} = \begin{pmatrix} 0 & 0 & 1 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}

Puis : N3=N2×N=(001000000)(011001000)=(000000000)=03,3N^3 = N^2 \times N = \begin{pmatrix} 0 & 0 & 1 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} 0 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} = \begin{pmatrix} 0 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix} = 0_{3,3}

La matrice NN est donc nilpotente d’ordre 33.


Question 2. Justifier que les matrices I3I_3 et NN commutent, puis utiliser la formule du binôme de Newton pour exprimer AnA^n pour tout n∈Nn \in \mathbb{N}.

Corrigé : La matrice identité I3I_3 commute avec toute matrice carrée du même ordre, en particulier avec NN, c’est-à-dire que I3N=NI3=NI_3 N = N I_3 = N.

Puisque I3I_3 et NN commutent, nous pouvons appliquer la formule du binôme de Newton pour tout n∈Nn \in \mathbb{N} : An=(I3+N)n=∑k=0n(nk)I3n−kNk=∑k=0n(nk)NkA^n = (I_3 + N)^n = \sum_{k=0}^{n} \binom{n}{k} I_3^{n-k} N^k = \sum_{k=0}^{n} \binom{n}{k} N^k

Puisque Nk=03,3N^k = 0_{3,3} pour tout k≥3k \ge 3, tous les termes de la somme d’indice k≥3k \ge 3 sont nuls.

  • Cas où n=0n = 0 : A0=I3A^0 = I_3 (résultat cohérent).
  • Cas où n=1n = 1 : A1=∑k=01(1k)Nk=I3+N=AA^1 = \sum_{k=0}^{1} \binom{1}{k} N^k = I_3 + N = A.
  • Cas où n=2n = 2 : A2=I3+2N+N2A^2 = I_3 + 2N + N^2.
  • Cas où n≥3n \ge 3 : An=(n0)I3+(n1)N+(n2)N2A^n = \binom{n}{0} I_3 + \binom{n}{1} N + \binom{n}{2} N^2

Explicitons cette expression en remplaçant les coefficients binomiaux pour n≥3n \ge 3 :

  • (n0)=1\binom{n}{0} = 1
  • (n1)=n\binom{n}{1} = n
  • (n2)=n(n−1)2\binom{n}{2} = \frac{n(n-1)}{2}

Ainsi : An=I3+nN+n(n−1)2N2A^n = I_3 + n N + \frac{n(n-1)}{2} N^2

Vérifions si cette formule reste valable pour n=0,1,2n = 0, 1, 2 :

  • Pour n=0n=0 : I3+0+0=I3I_3 + 0 + 0 = I_3 (vrai).
  • Pour n=1n=1 : I3+N+0=I3+N=AI_3 + N + 0 = I_3 + N = A (vrai).
  • Pour n=2n=2 : I3+2N+N2I_3 + 2N + N^2 (vrai).

La formule est donc valable pour tout n∈Nn \in \mathbb{N}.

Question 3. Expliciter sous forme matricielle les coefficients de AnA^n pour n∈Nn \in \mathbb{N}.

Corrigé : Calculons explicitement la matrice AnA^n : An=(100010001)+n(011001000)+n(n−1)2(001000000)A^n = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} + n \begin{pmatrix} 0 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} + \frac{n(n-1)}{2} \begin{pmatrix} 0 & 0 & 1 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}

An=(100010001)+(0nn00n000)+(00n(n−1)2000000)A^n = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} + \begin{pmatrix} 0 & n & n \\ 0 & 0 & n \\ 0 & 0 & 0 \end{pmatrix} + \begin{pmatrix} 0 & 0 & \frac{n(n-1)}{2} \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}

En sommant terme à terme : An=(1nn+n(n−1)201n001)A^n = \begin{pmatrix} 1 & n & n + \frac{n(n-1)}{2} \\ 0 & 1 & n \\ 0 & 0 & 1 \end{pmatrix}

Simplifions le coefficient de la première ligne, troisième colonne : n+n2−n2=2n+n2−n2=n2+n2=n(n+1)2n + \frac{n^2 - n}{2} = \frac{2n + n^2 - n}{2} = \frac{n^2 + n}{2} = \frac{n(n+1)}{2}

Finalement : An=(1nn(n+1)201n001)A^n = \begin{pmatrix} 1 & n & \frac{n(n+1)}{2} \\ 0 & 1 & n \\ 0 & 0 & 1 \end{pmatrix}


Partie 3 : Application à l’inversion de la matrice AA

Section titled “Partie 3 : Application à l’inversion de la matrice AAA”

Question 4. En utilisant la relation N=A−I3N = A - I_3 et le fait que N3=0N^3 = 0, trouver un polynôme PP tel que P(A)=0P(A) = 0. En déduire que AA est inversible et exprimer A−1A^{-1} en fonction de AA.

Corrigé : Puisque N=A−I3N = A - I_3 et que N3=03,3N^3 = 0_{3,3}, nous avons immédiatement : (A−I3)3=03,3(A - I_3)^3 = 0_{3,3}

Développons cette expression à l’aide des identités remarquables ou du binôme : A3−3A2+3A−I3=03,3A^3 - 3A^2 + 3A - I_3 = 0_{3,3}

Le polynôme P(X)=X3−3X2+3X−1P(X) = X^3 - 3X^2 + 3X - 1 est donc un polynôme annulateur de la matrice AA.

Pour montrer que AA est inversible et calculer son inverse, réécrivons l’égalité précédente en isolant la matrice identité I3I_3 : I3=3A−3A2+A3I_3 = 3A - 3A^2 + A^3 I3=A(3I3−3A+A2)I_3 = A(3I_3 - 3A + A^2)

Cette relation s’écrit sous la forme A×B=I3A \times B = I_3 avec B=3I3−3A+A2B = 3I_3 - 3A + A^2. Par conséquent, la matrice AA est inversible et son inverse est : A−1=A2−3A+3I3A^{-1} = A^2 - 3A + 3I_3

Vérifions ce résultat à l’aide de notre formule des puissances (avec n=−1n = -1 ou directement par le calcul polynomial) : Calculons A2A^2 : A2=(123012001)A^2 = \begin{pmatrix} 1 & 2 & 3 \\ 0 & 1 & 2 \\ 0 & 0 & 1 \end{pmatrix}

Alors : A−1=(123012001)−3(111011001)+3(100010001)A^{-1} = \begin{pmatrix} 1 & 2 & 3 \\ 0 & 1 & 2 \\ 0 & 0 & 1 \end{pmatrix} - 3\begin{pmatrix} 1 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 1 \end{pmatrix} + 3\begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}

A−1=(1−3+32−3+03−3+001−3+32−3+0001−3+3)=(1−1001−1001)A^{-1} = \begin{pmatrix} 1-3+3 & 2-3+0 & 3-3+0 \\ 0 & 1-3+3 & 2-3+0 \\ 0 & 0 & 1-3+3 \end{pmatrix} = \begin{pmatrix} 1 & -1 & 0 \\ 0 & 1 & -1 \\ 0 & 0 & 1 \end{pmatrix}


  1. Cohérence dimensionnelle et de type : Toutes les matrices manipulées sont carrées d’ordre 33 à coefficients réels, ce qui garantit la compatibilité des opérations (additions et multiplications).
  2. Cas particuliers :
    • Pour n=0n = 0, la formule générale donne A0=I3A^0 = I_3, ce qui est exact.
    • Pour n=1n = 1, A1=(111011001)=AA^1 = \begin{pmatrix} 1 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 1 \end{pmatrix} = A, ce qui valide la formule.
  3. Vérification de l’inverse : On peut tester le résultat en calculant le produit A×A−1A \times A^{-1} : (111011001)(1−1001−1001)=(1−1+1−1+101−1+1001)=I3\begin{pmatrix} 1 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} 1 & -1 & 0 \\ 0 & 1 & -1 \\ 0 & 0 & 1 \end{pmatrix} = \begin{pmatrix} 1 & -1+1 & -1+1 \\ 0 & 1 & -1+1 \\ 0 & 0 & 1 \end{pmatrix} = I_3 Le calcul est correct. L’utilisation d’un polynôme annulateur s’avère ici beaucoup plus rapide et élégante qu’un pivot de Gauss pour inverser la matrice.