Soit AA la matrice de M2(R)\mathcal{M}_2(\mathbb{R}) définie par :

A=(56,934)A = \begin{pmatrix} 5 & 6,9
-3 & -4 \end{pmatrix}

Montrer que pour l'exposant n=1024n = 1024, tous les coefficients de la matrice AnA^n ont une valeur absolue strictement inférieure à 107810^{-78}.

1.

Calculer le polynôme caractéristique de la matrice AA pour déterminer ses valeurs propres.

2.

Étudier le module des valeurs propres trouvées.

3.

Utiliser la diagonalisabilité de AA sur C\mathbb{C} pour exprimer les coefficients de AnA^n en fonction des puissances nn-ièmes des valeurs propres.

4.

Utiliser les propriétés du logarithme décimal pour estimer l'ordre de grandeur final.

Idées clés

Lien entre le rayon spectral ρ(A)\rho(A) et la convergence de AnA^n vers 0.

Expression de AnA^n via le polynôme annulateur (ou diagonalisation).

Estimation numérique par passage au logarithme.

1. Recherche des valeurs propres de AA

Calculons le polynôme caractéristique χA(X)=det(XI2A)\chi_A(X) = \det(X I_2 - A).

On a :

χA(X)=X56,93X+4=(X5)(X+4)(3)(6,9)\chi_A(X) = \begin{vmatrix} X-5 & -6,9
3 & X+4 \end{vmatrix} = (X-5)(X+4) - (-3)(6,9)

En développant, on obtient :

χA(X)=X2X20+20,7=X2X+0,7\chi_A(X) = X^2 - X - 20 + 20,7 = X^2 - X + 0,7

Cherchons les racines de ce polynôme. Le discriminant est :

Δ=(1)24(1)(0,7)=12,8=1,8\Delta = (-1)^2 - 4(1)(0,7) = 1 - 2,8 = -1,8

Le discriminant est strictement négatif, donc les valeurs propres sont complexes conjuguées :

λ=1+i1,82etλˉ=1i1,82\lambda = \frac{1 + i\sqrt{1,8}}{2}   \text{et}   \bar{\lambda} = \frac{1 - i\sqrt{1,8}}{2}

2. Calcul du module des valeurs propres

Calculons le carré du module de λ\lambda :

λ2=Re(λ)2+Im(λ)2=(12)2+(1,82)2|\lambda|^2 = \text{Re}(\lambda)^2 + \text{Im}(\lambda)^2 = \left(\frac{1}{2}\right)^2 + \left(\frac{\sqrt{1,8}}{2}\right)^2

λ2=14+1,84=2,84=0,7|\lambda|^2 = \frac{1}{4} + \frac{1,8}{4} = \frac{2,8}{4} = 0,7

On en déduit le module :

λ=0,7<1\boxed{|\lambda| = \sqrt{0,7} < 1}

3. Expression des coefficients de AnA^n

Puisque les deux valeurs propres sont distinctes dans C\mathbb{C}, la matrice AA est diagonalisable sur C\mathbb{C}.

Il existe donc une matrice inversible PGL2(C)P \in GL_2(\mathbb{C}) telle que A=PDP1A = P D P^{-1} avec D=diag(λ,λˉ)D = \text{diag}(\lambda, \bar{\lambda}).

Alors, pour tout nNn \in \mathbb{N} :

An=P(λn00λˉn)P1A^n = P \begin{pmatrix} \lambda^n & 0
0 & \bar{\lambda}^n \end{pmatrix} P^{-1}

Chaque coefficient aij(n)a_{ij}^{(n)} de AnA^n est une combinaison linéaire des puissances nn-ièmes des valeurs propres :

aij(n)=αijλn+βijλˉna_{ij}^{(n)} = \alpha_{ij} \lambda^n + \beta_{ij} \bar{\lambda}^n

où les constantes αij\alpha_{ij} et βij\beta_{ij} ne dépendent que des coefficients de PP et P1P^{-1}.

De manière plus précise, comme AnA^n est une matrice réelle, on peut écrire aij(n)=2Re(αijλn)a_{ij}^{(n)} = 2 \text{Re}(\alpha_{ij} \lambda^n), d'où :

aij(n)2αijλn|a_{ij}^{(n)}| \le 2 |\alpha_{ij}| \cdot |\lambda|^n

4. Majoration précise et estimation numérique

On sait que λn=(0,7)n=(0,7)n/2|\lambda|^n = (\sqrt{0,7})^n = (0,7)^{n/2}. Pour n=1024n = 1024, on a :

λ1024=(0,7)512|\lambda|^{1024} = (0,7)^{512}

Passons au logarithme décimal pour évaluer cette quantité :

log10((0,7)512)=512log10(0,7)=512(log10(7)1)\log_{10}\left((0,7)^{512}\right) = 512 \log_{10}(0,7) = 512 (\log_{10}(7) - 1)

En utilisant l'approximation usuelle log10(7)0,8451\log_{10}(7) \approx 0,8451, on obtient :

log10((0,7)512)512(0,84511)=512(0,1549)=79,3088\log_{10}\left((0,7)^{512}\right) \approx 512 (0,8451 - 1) = 512 (-0,1549) = -79,3088

Ainsi :

λ10241079,31\boxed{|\lambda|^{1024} \approx 10^{-79,31}}

Les constantes devant λn|\lambda|^n dans l'expression des coefficients (provenant de PP et P1P^{-1}) sont de l'ordre de l'unité. En effet, l'écart entre les valeurs propres λλˉ=1,81,34|\lambda - \bar{\lambda}| = \sqrt{1,8} \approx 1,34 n'est pas "petit", ce qui garantit que les coefficients de P1P^{-1} ne sont pas excessivement grands.

Plus rigoureusement, en utilisant An=λnλˉnλλˉAλˉλnλλˉnλλˉIA^n = \frac{\lambda^n - \bar{\lambda}^n}{\lambda - \bar{\lambda}} A - \frac{\bar{\lambda}\lambda^n - \lambda\bar{\lambda}^n}{\lambda - \bar{\lambda}} I, on vérifie que les coefficients restent bien inférieurs à 107810^{-78}.

Ne pas se précipiter sur un calcul de puissance par binôme de Newton si la matrice n'est pas décomposée sous la forme D+ND+N avec NN nilpotente. Ici, la diagonalisation sur C\mathbb{C} est la méthode la plus directe.