Algèbre linéaire 201-NYC / Maths 105 • Complément québécois de Première et cégep à Montréal

Exercices corrigés : diagonalisation et valeurs propres (201-NYC)

Voici la série d'exercices corrigés d'algèbre linéaire sur la diagonalisation, les valeurs propres et les vecteurs propres. La partie A installe la mécanique : la définition Av=λvA\vec{v}=\lambda\vec{v} et ce qu'elle signifie géométriquement, le polynôme caractéristique det(AλI)=0\det(A-\lambda I)=0, la détermination des sous-espaces propres par résolution d'un système homogène, et la construction de PP et DD vérifiant A=PDP1A=PDP^{-1}. La partie B monte au niveau examen : le critère exact de diagonalisabilité avec une matrice qui échoue au test, le calcul de AnA^{n} pour de grandes puissances, un problème complet où deux populations couplées se découplent dans la base propre, la diagonalisation orthogonale des matrices symétriques, le théorème de Cayley-Hamilton pour calculer inverse et puissances, et les valeurs propres des matrices AkA^{k}, A1A^{-1} et A+cIA+cI.

Cette série s'adresse aux étudiants de cégep et aux élèves du complément québécois de mathématiques au Lycée Marie de France et au Collège Stanislas. Selon le cégep, la diagonalisation figure en fin du cours 201-NYC, dans un cours d'algèbre linéaire ultérieur (souvent numéroté 201-105) ou en première année universitaire : vérifiez le plan de cours de votre enseignant avant de travailler cette série.

L'idée directrice mérite d'être énoncée d'emblée, car elle donne son sens à tous les calculs : une matrice quelconque mélange les coordonnées, mais dans la BONNE base, la même transformation se réduit à des dilatations le long d'axes indépendants. Diagonaliser, c'est trouver cette base. Tout le reste, polynôme caractéristique compris, n'est que la procédure pour y parvenir.

Série autocorrigée Tape tes réponses sous chaque question : la page te dit juste ou faux avant d'ouvrir la correction. Avec un compte, chaque bonne réponse du premier coup rapporte des points.

Ce chapitre fait partie de Algèbre linéaire et géométrie vectorielle, 201-NYC
Avant de commencer Fiche de révision : les pièges et la méthode de ce chapitre

Avant ce chapitre

Ces notions sont supposées acquises ici. Si le premier exercice résiste, le blocage vient presque toujours de l'une d'elles, pas du chapitre lui-même.

Remonter plus loin : la chaîne complète (5 chapitres) ↓

Le chemin de remédiation, du plus ancien au plus proche. Un élève qui reprend ce chapitre de zéro le reprend dans cet ordre.

  1. 1Les systèmes d'équationsSecondaire 4 SN4
  2. 2Combinaisons linéaires, dépendance et bases
  3. 3Matrices, systèmes et déterminants
  4. 4Systèmes linéaires : Gauss et Gauss-Jordan
  5. 5Déterminants, cofacteurs et règle de Cramer

Rappel de cours

  • Définition : λ\lambda est une VALEUR PROPRE de AA s'il existe un vecteur NON NUL v\vec{v} tel que Av=λvA\vec{v}=\lambda\vec{v}. Ce v\vec{v} est un VECTEUR PROPRE associé à λ\lambda. Géométriquement, AA ne fait qu'ÉTIRER v\vec{v} sans changer sa direction.
  • POLYNÔME CARACTÉRISTIQUE : les valeurs propres sont les racines de det(AλI)=0\det(A-\lambda I)=0. Pour une matrice 2×22\times 2, il s'écrit λ2tr(A)λ+det(A)=0\lambda^{2}-\operatorname{tr}(A)\lambda+\det(A)=0, ce qui donne un contrôle immédiat.
  • SOUS-ESPACE PROPRE associé à λ\lambda : Eλ=ker(AλI)E_{\lambda}=\ker(A-\lambda I), obtenu en résolvant le système homogène (AλI)v=0(A-\lambda I)\vec{v}=\vec{0}. C'est toujours un sous-espace, de dimension au moins 1.
  • Multiplicité ALGÉBRIQUE de λ\lambda : son ordre comme racine du polynôme caractéristique. Multiplicité GÉOMÉTRIQUE : dimEλ\dim E_{\lambda}. On a toujours géométrique \leq algébrique.
  • CRITÈRE DE DIAGONALISABILITÉ : AA (de taille nn) est diagonalisable si et seulement si la somme des dimensions des sous-espaces propres vaut nn, c'est-à-dire si pour CHAQUE valeur propre, multiplicité géométrique = multiplicité algébrique.
  • Cas favorable : si AA possède nn valeurs propres DISTINCTES, elle est automatiquement diagonalisable. Toute matrice SYMÉTRIQUE réelle est diagonalisable, avec des vecteurs propres orthogonaux.
  • Diagonalisation : A=PDP1A=PDP^{-1}, où les COLONNES de PP sont les vecteurs propres et où DD est diagonale avec les valeurs propres dans le MÊME ORDRE. Conséquence : An=PDnP1A^{n}=PD^{n}P^{-1}, et DnD^{n} s'obtient en élevant chaque terme diagonal à la puissance nn.

Partie A : Les bases (/50)

Exercice 1 : Valeurs propres et vecteurs propres : la définition

Un vecteur propre est un vecteur que la matrice ne fait qu'étirer, sans le faire tourner. Tout part de cette image géométrique.

  • a) Vérifiez que v=(11)\vec{v}=\begin{pmatrix}1\\1\end{pmatrix} est un vecteur propre de A=(3113)A=\begin{pmatrix}3&1\\1&3\end{pmatrix} et donnez la valeur propre associée.
  • b) Vérifiez de même pour w=(11)\vec{w}=\begin{pmatrix}1\\-1\end{pmatrix}. Que remarquez-vous sur les directions de v\vec{v} et w\vec{w} ?
  • c) Le vecteur (10)\begin{pmatrix}1\\0\end{pmatrix} est-il un vecteur propre de AA ? Justifiez par le calcul.
  • d) Pourquoi exige-t-on v0\vec{v}\neq\vec{0} dans la définition ? Que se passerait-il sans cette condition ?

Tape tes réponses, la page te dit juste ou faux 0/5

a)
b)
c)
d)
Voir la correction

Réponses

  • a) Av=4vA\vec v=4\vec v : λ=4\lambda=4
  • b) Aw=2wA\vec w=2\vec w ; vw\vec v\perp\vec w
  • c) Non : (3; 1)(3;\ 1) n'est pas colinéaire
  • d) 0\vec 0 rendrait tout λ\lambda propre

a) On calcule Av=(3113)(11)=(3+11+3)=(44)=4(11)=4vA\vec{v}=\begin{pmatrix}3&1\\1&3\end{pmatrix}\begin{pmatrix}1\\1\end{pmatrix}=\begin{pmatrix}3+1\\1+3\end{pmatrix}=\begin{pmatrix}4\\4\end{pmatrix}=4\begin{pmatrix}1\\1\end{pmatrix}=4\vec{v}. La relation Av=λvA\vec{v}=\lambda\vec{v} est vérifiée avec λ=4\lambda=4 : v\vec{v} est bien un vecteur propre, de valeur propre 4. Géométriquement, la matrice multiplie par 4 tout vecteur dirigé selon la diagonale y=xy=x, sans le faire dévier.

b) Aw=(3113)(11)=(3113)=(22)=2wA\vec{w}=\begin{pmatrix}3&1\\1&3\end{pmatrix}\begin{pmatrix}1\\-1\end{pmatrix}=\begin{pmatrix}3-1\\1-3\end{pmatrix}=\begin{pmatrix}2\\-2\end{pmatrix}=2\vec{w}. Donc w\vec{w} est un vecteur propre de valeur propre 2. Remarque sur les directions : vw=(1)(1)+(1)(1)=0\vec{v}\cdot\vec{w}=(1)(1)+(1)(-1)=0, les deux vecteurs propres sont ORTHOGONAUX. Ce n'est pas un hasard : AA est symétrique (a12=a21=1a_{12}=a_{21}=1), et les vecteurs propres d'une matrice symétrique associés à des valeurs propres distinctes sont toujours orthogonaux. L'exercice 8 y revient.

c) A(10)=(31)A\begin{pmatrix}1\\0\end{pmatrix}=\begin{pmatrix}3\\1\end{pmatrix}. Pour que ce soit un vecteur propre, il faudrait (31)=λ(10)=(λ0)\begin{pmatrix}3\\1\end{pmatrix}=\lambda\begin{pmatrix}1\\0\end{pmatrix}=\begin{pmatrix}\lambda\\0\end{pmatrix}, ce qui exige simultanément λ=3\lambda=3 et 1=01=0 : impossible. Donc (10)\begin{pmatrix}1\\0\end{pmatrix} n'est PAS un vecteur propre. Géométriquement, AA fait tourner ce vecteur en plus de l'allonger, puisque l'image (31)\begin{pmatrix}3\\1\end{pmatrix} n'est pas colinéaire au vecteur de départ. Le test pratique est donc simple : calculer AuA\vec{u} et vérifier si le résultat est PROPORTIONNEL à u\vec{u}.

d) Sans la condition v0\vec{v}\neq\vec{0}, la définition perdrait tout contenu : A0=0=λ0A\vec{0}=\vec{0}=\lambda\vec{0} est vraie pour TOUT scalaire λ\lambda, si bien que n'importe quel nombre serait valeur propre de n'importe quelle matrice. La notion ne distinguerait plus rien. C'est aussi pourquoi les sous-espaces propres se définissent comme des NOYAUX : Eλ=ker(AλI)E_{\lambda}=\ker(A-\lambda I) contient nécessairement 0\vec{0} (tout sous-espace le contient), et la condition « λ\lambda est valeur propre » se traduit exactement par « ce noyau contient AUTRE CHOSE que 0\vec{0} », c'est-à-dire dimEλ1\dim E_{\lambda}\geq 1, ou encore det(AλI)=0\det(A-\lambda I)=0. On voit ici pourquoi le déterminant apparaît : il s'annule précisément quand le système homogène admet des solutions non triviales.

-11234567-2-1123456xAx = 2xyAyx est VECTEUR PROPRE : Ax reste sur la même droitey ne l'est pas : sa direction change

Exercice 2 : Le polynôme caractéristique

Chercher les valeurs propres revient à chercher pour quels λ\lambda le système (AλI)v=0(A-\lambda I)\vec{v}=\vec{0} admet une solution non triviale. Le déterminant tranche.

  • a) Expliquez pourquoi la condition « λ\lambda est valeur propre de AA » équivaut à det(AλI)=0\det(A-\lambda I)=0.
  • b) Déterminez les valeurs propres de A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix}.
  • c) Vérifiez vos résultats avec les deux relations λ1+λ2=tr(A)\lambda_{1}+\lambda_{2}=\operatorname{tr}(A) et λ1λ2=det(A)\lambda_{1}\lambda_{2}=\det(A).
  • d) Déterminez les valeurs propres de B=(200130451)B=\begin{pmatrix}2&0&0\\1&3&0\\4&5&-1\end{pmatrix}. Quelle particularité rend le calcul immédiat ?

Tape tes réponses, la page te dit juste ou faux 0/6

b)
c)
d)
Voir la correction

Réponses

  • a) Solutions non triviales : det(AλI)=0\det(A-\lambda I)=0
  • b) λ=2\lambda=2 et λ=5\lambda=5
  • c) 2+5=72+5=7, 2×5=102\times 5=10
  • d) Triangulaire : 22, 33, 1-1

a) Par définition, λ\lambda est valeur propre s'il existe v0\vec{v}\neq\vec{0} avec Av=λvA\vec{v}=\lambda\vec{v}, ce qui s'écrit Avλv=0A\vec{v}-\lambda\vec{v}=\vec{0}, soit (AλI)v=0(A-\lambda I)\vec{v}=\vec{0} (l'insertion de II est nécessaire pour que la soustraction ait un sens matriciel). La question devient donc : le système homogène (AλI)v=0(A-\lambda I)\vec{v}=\vec{0} admet-il une solution NON TRIVIALE ? Or un système homogène carré admet des solutions non triviales si et seulement si sa matrice n'est pas inversible, c'est-à-dire si et seulement si son déterminant est nul. D'où l'équivalence λ\lambda valeur propre     det(AλI)=0\iff\det(A-\lambda I)=0 ✓. Tout le chapitre repose sur cette traduction, et elle réutilise directement le critère d'inversibilité vu au chapitre des déterminants.

b) AλI=(4λ213λ)A-\lambda I=\begin{pmatrix}4-\lambda&2\\1&3-\lambda\end{pmatrix}, de déterminant (4λ)(3λ)2=124λ3λ+λ22=λ27λ+10(4-\lambda)(3-\lambda)-2=12-4\lambda-3\lambda+\lambda^{2}-2=\lambda^{2}-7\lambda+10. En annulant : λ27λ+10=0\lambda^{2}-7\lambda+10=0, soit (λ2)(λ5)=0(\lambda-2)(\lambda-5)=0. Les valeurs propres sont donc λ1=2\lambda_{1}=2 et λ2=5\lambda_{2}=5.

c) Trace : tr(A)=4+3=7\operatorname{tr}(A)=4+3=7, et λ1+λ2=2+5=7\lambda_{1}+\lambda_{2}=2+5=7 ✓. Déterminant : det(A)=4×32×1=10\det(A)=4\times 3-2\times 1=10, et λ1λ2=2×5=10\lambda_{1}\lambda_{2}=2\times 5=10 ✓. Ces deux relations sont la lecture des coefficients du polynôme caractéristique, qui s'écrit toujours λ2tr(A)λ+det(A)\lambda^{2}-\operatorname{tr}(A)\lambda+\det(A) pour une matrice 2×22\times 2. Elles fournissent une vérification en cinq secondes, et permettent même de deviner les valeurs propres quand elles sont entières : ici, chercher deux nombres de somme 7 et de produit 10 donne immédiatement 2 et 5. C'est le contrôle à faire systématiquement à l'examen.

d) La matrice BB est TRIANGULAIRE INFÉRIEURE (tous les termes au-dessus de la diagonale sont nuls). Or le déterminant d'une matrice triangulaire est le produit de ses termes diagonaux, donc det(BλI)=(2λ)(3λ)(1λ)\det(B-\lambda I)=(2-\lambda)(3-\lambda)(-1-\lambda), qui s'annule pour λ=2\lambda=2, λ=3\lambda=3 et λ=1\lambda=-1. Les valeurs propres d'une matrice triangulaire sont donc exactement ses termes DIAGONAUX, sans aucun calcul. La particularité vaut aussi pour les matrices triangulaires supérieures et, cas limite, pour les matrices diagonales. C'est cohérent avec l'objectif du chapitre : diagonaliser une matrice, c'est la ramener à une forme où ses valeurs propres se lisent directement sur la diagonale.

Exercice 3 : Les sous-espaces propres

Une fois les valeurs propres connues, chaque λ\lambda donne un système homogène à résoudre. Le sous-espace obtenu est le sous-espace propre, et c'est sa DIMENSION qui décidera de tout au prochain exercice.

  • a) Pour A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix} (valeurs propres 2 et 5), déterminez E2=ker(A2I)E_{2}=\ker(A-2I) et donnez-en une base.
  • b) Déterminez de même E5E_{5} et donnez-en une base.
  • c) Vérifiez directement que vos deux vecteurs propres satisfont Av=λvA\vec{v}=\lambda\vec{v}. Forment-ils une base de R2\mathbb{R}^{2} ?
  • d) Pourquoi un vecteur propre n'est-il jamais unique ? Que peut-on dire de tous les vecteurs propres associés à une même valeur propre ?

Tape tes réponses, la page te dit juste ou faux 0/5

a)
b)
c)
d)
Voir la correction

Réponses

  • a) E2=vect{(1; 1)}E_2=\text{vect}\{(-1;\ 1)\}
  • b) E5=vect{(2; 1)}E_5=\text{vect}\{(2;\ 1)\}
  • c) Déterminant 3-3 : base propre
  • d) Toute la direction est propre

a) A2I=(2211)A-2I=\begin{pmatrix}2&2\\1&1\end{pmatrix}. Le système (A2I)v=0(A-2I)\vec{v}=\vec{0} s'écrit 2x+2y=02x+2y=0 et x+y=0x+y=0 : les deux équations sont proportionnelles et se réduisent à x=yx=-y. En posant y=ty=t : v=(t,t)=t(1,1)\vec{v}=(-t,t)=t(-1,1). Donc E2=vect{(11)}E_{2}=\text{vect}\left\{\begin{pmatrix}-1\\1\end{pmatrix}\right\}, une droite, et dimE2=1\dim E_{2}=1. Contrôle de cohérence utile : le système DOIT être dégénéré (deux équations proportionnelles), sinon la seule solution serait 0\vec{0} et λ=2\lambda=2 ne serait pas valeur propre. Si l'échelonnement ne dégénère pas, c'est qu'il y a une erreur dans le polynôme caractéristique.

b) A5I=(1212)A-5I=\begin{pmatrix}-1&2\\1&-2\end{pmatrix}. Le système donne x+2y=0-x+2y=0, soit x=2yx=2y. En posant y=ty=t : v=(2t,t)=t(2,1)\vec{v}=(2t,t)=t(2,1). Donc E5=vect{(21)}E_{5}=\text{vect}\left\{\begin{pmatrix}2\\1\end{pmatrix}\right\} et dimE5=1\dim E_{5}=1. Là encore les deux lignes sont proportionnelles ✓.

c) Vérification pour λ=2\lambda=2 : A(11)=(4+21+3)=(22)=2(11)A\begin{pmatrix}-1\\1\end{pmatrix}=\begin{pmatrix}-4+2\\-1+3\end{pmatrix}=\begin{pmatrix}-2\\2\end{pmatrix}=2\begin{pmatrix}-1\\1\end{pmatrix} ✓. Pour λ=5\lambda=5 : A(21)=(8+22+3)=(105)=5(21)A\begin{pmatrix}2\\1\end{pmatrix}=\begin{pmatrix}8+2\\2+3\end{pmatrix}=\begin{pmatrix}10\\5\end{pmatrix}=5\begin{pmatrix}2\\1\end{pmatrix} ✓. Forment-ils une base de R2\mathbb{R}^{2} ? Ils sont deux dans un espace de dimension 2, il suffit donc de vérifier l'indépendance : le déterminant 1211=12=30\begin{vmatrix}-1&2\\1&1\end{vmatrix}=-1-2=-3\neq 0, ils sont indépendants et forment bien une BASE. C'est ce fait qui rendra AA diagonalisable : il existe une base de R2\mathbb{R}^{2} entièrement composée de vecteurs propres.

d) Un vecteur propre n'est jamais unique parce que tout multiple non nul d'un vecteur propre est encore un vecteur propre pour la même valeur propre : si Av=λvA\vec{v}=\lambda\vec{v}, alors A(cv)=cAv=cλv=λ(cv)A(c\vec{v})=cA\vec{v}=c\lambda\vec{v}=\lambda(c\vec{v}) pour tout c0c\neq 0 ✓. Géométriquement, c'est la DIRECTION qui est propre, pas un vecteur particulier de cette direction. Plus généralement, l'ensemble des vecteurs propres associés à λ\lambda, auquel on adjoint 0\vec{0}, forme exactement le sous-espace Eλ=ker(AλI)E_{\lambda}=\ker(A-\lambda I) : il est donc fermé pour l'addition et la multiplication scalaire. Conséquence pratique pour l'examen : deux étudiants peuvent donner (11)\begin{pmatrix}-1\\1\end{pmatrix} et (33)\begin{pmatrix}3\\-3\end{pmatrix} comme vecteur propre et avoir tous deux raison. Seule la direction compte, et c'est pourquoi la réponse attendue est une BASE du sous-espace propre plutôt qu'un vecteur imposé.

Exercice 4 : Diagonaliser : construire P et D

Diagonaliser, c'est réécrire la matrice dans la base de ses vecteurs propres. La seule difficulté technique est de respecter l'ORDRE des colonnes.

  • a) Pour A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix}, écrivez PP et DD telles que A=PDP1A=PDP^{-1}.
  • b) Calculez P1P^{-1} et vérifiez explicitement que PDP1=APDP^{-1}=A.
  • c) Que se passe-t-il si l'on échange les deux colonnes de PP sans toucher à DD ? Illustrez le problème.
  • d) Interprétez géométriquement l'égalité A=PDP1A=PDP^{-1} : que fait chacune des trois matrices, lues de droite à gauche ?

Tape tes réponses, la page te dit juste ou faux 0/7

a)
b)
c)
d)
Voir la correction

Réponses

  • a) P=(1211)P=\begin{pmatrix}-1&2\\1&1\end{pmatrix}, D=diag(2;5)D=\text{diag}(2;5)
  • b) P1=(13231313)P^{-1}=\begin{pmatrix}-\frac13&\frac23\\\frac13&\frac13\end{pmatrix} ; PDP1=APDP^{-1}=A
  • c) PDP1=(3214)AP'DP'^{-1}=\begin{pmatrix}3&-2\\-1&4\end{pmatrix}\neq A
  • d) Changer de base, dilater, revenir

a) On place les vecteurs propres en COLONNES de PP et les valeurs propres correspondantes sur la diagonale de DD, dans le même ordre. Avec v2=(11)\vec{v}_{2}=\begin{pmatrix}-1\\1\end{pmatrix} pour λ=2\lambda=2 et v5=(21)\vec{v}_{5}=\begin{pmatrix}2\\1\end{pmatrix} pour λ=5\lambda=5 : P=(1211)P=\begin{pmatrix}-1&2\\1&1\end{pmatrix} et D=(2005)D=\begin{pmatrix}2&0\\0&5\end{pmatrix}.

b) detP=(1)(1)(2)(1)=3\det P=(-1)(1)-(2)(1)=-3, donc P1=13(1211)=(13231313)P^{-1}=\frac{1}{-3}\begin{pmatrix}1&-2\\-1&-1\end{pmatrix}=\begin{pmatrix}-\frac{1}{3}&\frac{2}{3}\\\frac{1}{3}&\frac{1}{3}\end{pmatrix}. Vérification : PD=(1211)(2005)=(21025)PD=\begin{pmatrix}-1&2\\1&1\end{pmatrix}\begin{pmatrix}2&0\\0&5\end{pmatrix}=\begin{pmatrix}-2&10\\2&5\end{pmatrix}. Puis PDP1=(21025)(13231313)=(23+10343+10323+5343+53)=(4213)=APDP^{-1}=\begin{pmatrix}-2&10\\2&5\end{pmatrix}\begin{pmatrix}-\frac{1}{3}&\frac{2}{3}\\\frac{1}{3}&\frac{1}{3}\end{pmatrix}=\begin{pmatrix}\frac{2}{3}+\frac{10}{3}&-\frac{4}{3}+\frac{10}{3}\\-\frac{2}{3}+\frac{5}{3}&\frac{4}{3}+\frac{5}{3}\end{pmatrix}=\begin{pmatrix}4&2\\1&3\end{pmatrix}=A ✓.

c) Si l'on échange les colonnes de PP pour obtenir P=(2111)P'=\begin{pmatrix}2&-1\\1&1\end{pmatrix} tout en gardant D=(2005)D=\begin{pmatrix}2&0\\0&5\end{pmatrix}, l'égalité tombe : on affirmerait alors que (21)\begin{pmatrix}2\\1\end{pmatrix} est associé à la valeur propre 2, alors qu'il l'est à 5. Le calcul le confirme, PDP1=(3214)AP'DP'^{-1}=\begin{pmatrix}3&-2\\-1&4\end{pmatrix}\neq A. La règle est donc stricte : la ii-ème colonne de PP doit être un vecteur propre associé au ii-ème terme diagonal de DD. En revanche, échanger simultanément les colonnes de PP ET les termes de DD est parfaitement licite et donne une autre diagonalisation, tout aussi valable : la diagonalisation n'est donc pas unique, seul l'APPARIEMENT compte.

d) Lue de droite à gauche sur un vecteur x\vec{x}, l'égalité Ax=PDP1xA\vec{x}=PDP^{-1}\vec{x} se décompose en trois étapes. D'abord P1P^{-1} CHANGE DE BASE : il exprime x\vec{x} dans la base des vecteurs propres, c'est-à-dire qu'il donne les coordonnées de x\vec{x} le long des deux directions privilégiées. Ensuite DD agit dans cette base, et son action est la plus simple possible : elle multiplie la première coordonnée par 2 et la seconde par 5, sans aucun mélange, puisque DD est diagonale. Enfin PP REVIENT à la base canonique. La transformation compliquée AA est donc une simple double dilatation, à condition de la regarder dans le bon repère. C'est tout le sens du chapitre : AA ne mélange les coordonnées que parce qu'on l'observe dans une base mal choisie, et diagonaliser revient à choisir le repère où la transformation devient transparente.

Exercice 5 : Valeurs propres complexes : une matrice sans direction propre réelle

Le polynôme caractéristique d'une matrice réelle peut avoir des racines complexes. La matrice n'a alors aucun vecteur propre réel, ce qui a une lecture géométrique très nette.

  • a) Déterminez les valeurs propres de la rotation A=(0110)A=\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} (rotation de 9090^{\circ}).
  • b) Expliquez géométriquement pourquoi cette matrice ne possède aucun vecteur propre réel.
  • c) Déterminez les valeurs propres de B=(2112)B=\begin{pmatrix} 2 & -1 \\ 1 & 2 \end{pmatrix}.
  • d) Vérifiez sur BB que les relations λ1+λ2=tr(B)\lambda_{1}+\lambda_{2}=\operatorname{tr}(B) et λ1λ2=det(B)\lambda_{1}\lambda_{2}=\det(B) restent valables avec des valeurs propres complexes conjuguées.

Tape tes réponses, la page te dit juste ou faux 0/6

a)
b)
c)
d)
Voir la correction

Réponses

  • a) λ=±i\lambda=\pm i
  • b) Rotation : aucune direction conservée
  • c) 2±i2\pm i
  • d) Somme 4, produit 5

a) Le polynôme caractéristique est det(AλI)=λ11λ=λ2+1\det(A-\lambda I)=\begin{vmatrix} -\lambda & -1 \\ 1 & -\lambda \end{vmatrix}=\lambda^{2}+1. En annulant : λ2=1\lambda^{2}=-1, donc λ=i\lambda=i et λ=i\lambda=-i. Les deux valeurs propres sont imaginaires pures, conjuguées l'une de l'autre. On retrouve d'ailleurs tr(A)=0=i+(i)\operatorname{tr}(A)=0=i+(-i) et det(A)=1=i(i)\det(A)=1=i\cdot(-i).

b) Un vecteur propre réel serait un vecteur réel non nul dont la direction est préservée par AA. Mais AA est une rotation de 9090^{\circ} : elle fait tourner TOUT vecteur du plan d'un quart de tour, donc aucun vecteur réel non nul ne conserve sa direction. Il ne peut donc exister aucun vecteur propre réel, ce qui est cohérent avec l'absence de valeur propre réelle. Les directions propres n'apparaissent qu'en passant au plan complexe.

c) det(BλI)=2λ112λ=(2λ)2+1=λ24λ+5\det(B-\lambda I)=\begin{vmatrix} 2-\lambda & -1 \\ 1 & 2-\lambda \end{vmatrix}=(2-\lambda)^{2}+1=\lambda^{2}-4\lambda+5. Le discriminant vaut 1620=4<016-20=-4<0, donc les racines sont complexes : λ=4±42=4±2i2=2±i\lambda=\dfrac{4\pm\sqrt{-4}}{2}=\dfrac{4\pm 2i}{2}=2\pm i. Les valeurs propres sont 2+i2+i et 2i2-i.

d) Somme : (2+i)+(2i)=4=tr(B)(2+i)+(2-i)=4=\operatorname{tr}(B) ✓. Produit : (2+i)(2i)=4i2=4+1=5=det(B)(2+i)(2-i)=4-i^{2}=4+1=5=\det(B) ✓. Les deux relations, démontrées pour des valeurs propres réelles, restent vraies pour des valeurs propres complexes, car elles ne sont que la lecture des coefficients du polynôme caractéristique λ2tr(B)λ+det(B)\lambda^{2}-\operatorname{tr}(B)\lambda+\det(B). Remarquez que la partie réelle commune, 22, est la demi-trace : c'est général pour une paire conjuguée, et cela donne un contrôle rapide. La matrice BB agit comme une similitude, rotation composée d'une homothétie, ce que confirmera le lien avec les nombres complexes de module 5\sqrt{5}.

-4-3-2-11234-4-3-2-11234rotation de 50°AUCUN vecteur ne garde sa direction :les valeurs propres sont COMPLEXES

Partie B : Niveau examen (/50)

Exercice 6 : Le critère de diagonalisabilité et la matrice qui échoue

Toutes les matrices ne sont pas diagonalisables, et l'examen teste presque toujours ce point. Le critère se ramène à une comparaison de deux multiplicités.

  • a) Énoncez le critère exact de diagonalisabilité, en distinguant multiplicité algébrique et multiplicité géométrique.
  • b) Montrez que A=(3103)A=\begin{pmatrix}3&1\\0&3\end{pmatrix} n'est PAS diagonalisable, en calculant les deux multiplicités.
  • c) Montrez au contraire que B=(3003)B=\begin{pmatrix}3&0\\0&3\end{pmatrix} est diagonalisable, bien qu'elle ait aussi une valeur propre double. Que change la différence entre AA et BB ?
  • d) Démontrez que si une matrice n×nn\times n possède nn valeurs propres DISTINCTES, alors elle est diagonalisable. Utilisez ce résultat pour conclure sans calcul sur (5702)\begin{pmatrix}5&7\\0&-2\end{pmatrix}.

Tape tes réponses, la page te dit juste ou faux 0/7

a)
b)
c)
d)
Voir la correction

Réponses

  • a) Géométrique = algébrique pour chaque λ\lambda
  • b) Multiplicités 2 et 1 : non
  • c) E3=R2E_3=\mathbb{R}^2 : oui
  • d) Valeurs propres distinctes : oui

a) Multiplicité ALGÉBRIQUE de λ\lambda : son ordre de multiplicité comme racine du polynôme caractéristique. Multiplicité GÉOMÉTRIQUE de λ\lambda : la dimension du sous-espace propre Eλ=ker(AλI)E_{\lambda}=\ker(A-\lambda I), autrement dit le nombre de vecteurs propres indépendants associés à λ\lambda. On a toujours 11\leq géométrique \leq algébrique. CRITÈRE : AA de taille nn est diagonalisable si et seulement si, pour CHAQUE valeur propre, la multiplicité géométrique égale la multiplicité algébrique; de façon équivalente, si et seulement si la somme des dimensions des sous-espaces propres vaut nn, c'est-à-dire s'il existe une base de l'espace formée de vecteurs propres.

b) AA est triangulaire, donc ses valeurs propres sont ses termes diagonaux : λ=3\lambda=3, avec multiplicité ALGÉBRIQUE 2 (le polynôme caractéristique est (3λ)2(3-\lambda)^{2}). Multiplicité géométrique : A3I=(0100)A-3I=\begin{pmatrix}0&1\\0&0\end{pmatrix}, et le système (A3I)v=0(A-3I)\vec{v}=\vec{0} donne la seule équation y=0y=0, avec xx libre. Donc E3=vect{(10)}E_{3}=\text{vect}\left\{\begin{pmatrix}1\\0\end{pmatrix}\right\} et dimE3=1\dim E_{3}=1. Comparaison : géométrique =1<2==1<2= algébrique. Le critère échoue, AA n'est PAS diagonalisable ✓. Concrètement, il n'existe qu'UNE direction propre dans un espace de dimension 2 : impossible de construire une base de vecteurs propres, et donc impossible de former une matrice PP inversible.

c) B=3IB=3I est également triangulaire de valeur propre 3 avec multiplicité algébrique 2. Mais ici B3I=(0000)B-3I=\begin{pmatrix}0&0\\0&0\end{pmatrix}, la matrice NULLE : le système (B3I)v=0(B-3I)\vec{v}=\vec{0} n'impose AUCUNE condition, toute solution convient, donc E3=R2E_{3}=\mathbb{R}^{2} et dimE3=2\dim E_{3}=2. Géométrique =2==2= algébrique : le critère est satisfait et BB est diagonalisable. Elle l'est même déjà, puisqu'elle est diagonale, avec P=IP=I. Ce qui change entre AA et BB est le terme hors diagonale : B=3IB=3I dilate uniformément toutes les directions, si bien que TOUT vecteur non nul est propre; AA, elle, ajoute un CISAILLEMENT qui ne laisse invariante qu'une seule direction. Deux matrices avec le même polynôme caractéristique peuvent donc avoir des comportements opposés : le polynôme caractéristique seul ne décide jamais de la diagonalisabilité, et c'est la raison d'être de la multiplicité géométrique.

d) Démonstration. Soient λ1,,λn\lambda_{1},\ldots,\lambda_{n} les nn valeurs propres distinctes, et v1,,vn\vec{v}_{1},\ldots,\vec{v}_{n} des vecteurs propres associés. Le point clef est que des vecteurs propres associés à des valeurs propres DISTINCTES sont linéairement indépendants. Voyons-le pour deux : supposons c1v1+c2v2=0c_{1}\vec{v}_{1}+c_{2}\vec{v}_{2}=\vec{0}. En appliquant AA : c1λ1v1+c2λ2v2=0c_{1}\lambda_{1}\vec{v}_{1}+c_{2}\lambda_{2}\vec{v}_{2}=\vec{0}. En multipliant la première relation par λ1\lambda_{1} et en soustrayant : c2(λ2λ1)v2=0c_{2}(\lambda_{2}-\lambda_{1})\vec{v}_{2}=\vec{0}. Comme λ2λ1\lambda_{2}\neq\lambda_{1} et v20\vec{v}_{2}\neq\vec{0}, il vient c2=0c_{2}=0, puis c1=0c_{1}=0 ✓ (le cas général se traite de même par récurrence). On dispose donc de nn vecteurs indépendants dans un espace de dimension nn : ils en forment une base, entièrement composée de vecteurs propres, et AA est diagonalisable ✓. Application : (5702)\begin{pmatrix}5&7\\0&-2\end{pmatrix} est triangulaire, ses valeurs propres sont 5 et 2-2, qui sont DISTINCTES et au nombre de 2 dans un espace de dimension 2. Elle est donc diagonalisable, sans avoir eu à calculer le moindre vecteur propre. Attention à ne pas retourner l'implication : la condition est SUFFISANTE, pas nécessaire, comme le montre B=3IB=3I du point c), diagonalisable avec une valeur propre double.

Exercice 7 : Problème : puissances d'une matrice et deux populations couplées

La diagonalisation n'est pas un exercice formel : elle permet de calculer AnA^{n} pour nn grand, et donc de prédire l'évolution à long terme d'un système couplé. C'est l'application qui justifie tout le chapitre.

  • a) Pour A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix} diagonalisée en P=(1211)P=\begin{pmatrix}-1&2\\1&1\end{pmatrix} et D=(2005)D=\begin{pmatrix}2&0\\0&5\end{pmatrix}, démontrez la formule An=PDnP1A^{n}=PD^{n}P^{-1}.
  • b) Calculez A3A^{3} par cette formule, puis vérifiez par multiplication directe.
  • c) Deux populations évoluent chaque année selon {xk+1=4xk+2ykyk+1=xk+3yk\begin{cases}x_{k+1}=4x_{k}+2y_{k}\\y_{k+1}=x_{k}+3y_{k}\end{cases}, avec x0=3x_{0}=3 et y0=1y_{0}=1. Exprimez (xkyk)\begin{pmatrix}x_{k}\\y_{k}\end{pmatrix} en fonction de kk.
  • d) Vers quelle proportion xkyk\frac{x_{k}}{y_{k}} tend le système quand kk\to\infty ? Interprétez ce résultat en termes de valeurs propres.

Tape tes réponses, la page te dit juste ou faux 0/6

b)
c)
d)
Voir la correction

Réponses

  • a) An=PDnP1A^n=PD^nP^{-1}
  • b) A3=(86783947)A^3=\begin{pmatrix}86&78\\39&47\end{pmatrix}
  • c) xk=2k+85k3x_k=\frac{2^k+8\cdot 5^k}{3}, yk=2k+45k3y_k=\frac{-2^k+4\cdot 5^k}{3}
  • d) Rapport limite 2

a) Démonstration par récurrence, ou plus directement par télescopage. Pour n=2n=2 : A2=(PDP1)(PDP1)=PD(P1P)DP1=PDIDP1=PD2P1A^{2}=\left(PDP^{-1}\right)\left(PDP^{-1}\right)=PD\left(P^{-1}P\right)DP^{-1}=PDIDP^{-1}=PD^{2}P^{-1} ✓, où l'on a utilisé P1P=IP^{-1}P=I. Le même mécanisme se répète : dans An=(PDP1)nA^{n}=\left(PDP^{-1}\right)^{n}, tous les facteurs intérieurs P1PP^{-1}P se simplifient deux à deux, et il ne reste que PDnP1PD^{n}P^{-1}. L'intérêt est décisif : DnD^{n} se calcule sans effort, puisqu'une matrice diagonale élevée à la puissance nn se réduit à (λ1n00λ2n)\begin{pmatrix}\lambda_{1}^{n}&0\\0&\lambda_{2}^{n}\end{pmatrix}. On remplace donc nn multiplications matricielles par deux multiplications et deux exponentiations de nombres.

b) D3=(800125)D^{3}=\begin{pmatrix}8&0\\0&125\end{pmatrix}. Puis PD3=(1211)(800125)=(82508125)PD^{3}=\begin{pmatrix}-1&2\\1&1\end{pmatrix}\begin{pmatrix}8&0\\0&125\end{pmatrix}=\begin{pmatrix}-8&250\\8&125\end{pmatrix}, et A3=PD3P1=(82508125)(13231313)=(83+2503163+250383+1253163+1253)=(86783947)A^{3}=PD^{3}P^{-1}=\begin{pmatrix}-8&250\\8&125\end{pmatrix}\begin{pmatrix}-\frac{1}{3}&\frac{2}{3}\\\frac{1}{3}&\frac{1}{3}\end{pmatrix}=\begin{pmatrix}\frac{8}{3}+\frac{250}{3}&-\frac{16}{3}+\frac{250}{3}\\-\frac{8}{3}+\frac{125}{3}&\frac{16}{3}+\frac{125}{3}\end{pmatrix}=\begin{pmatrix}86&78\\39&47\end{pmatrix}. Vérification directe : A2=(4213)2=(1814711)A^{2}=\begin{pmatrix}4&2\\1&3\end{pmatrix}^{2}=\begin{pmatrix}18&14\\7&11\end{pmatrix}, puis A3=A2A=(1814711)(4213)=(72+1436+4228+1114+33)=(86783947)A^{3}=A^{2}A=\begin{pmatrix}18&14\\7&11\end{pmatrix}\begin{pmatrix}4&2\\1&3\end{pmatrix}=\begin{pmatrix}72+14&36+42\\28+11&14+33\end{pmatrix}=\begin{pmatrix}86&78\\39&47\end{pmatrix} ✓.

c) Le système s'écrit uk+1=Auk\vec{u}_{k+1}=A\vec{u}_{k} avec uk=(xkyk)\vec{u}_{k}=\begin{pmatrix}x_{k}\\y_{k}\end{pmatrix}, donc uk=Aku0\vec{u}_{k}=A^{k}\vec{u}_{0}. La méthode efficace consiste à décomposer u0\vec{u}_{0} dans la base propre : cherchons aa et bb tels que (31)=a(11)+b(21)\begin{pmatrix}3\\1\end{pmatrix}=a\begin{pmatrix}-1\\1\end{pmatrix}+b\begin{pmatrix}2\\1\end{pmatrix}. Cela donne a+2b=3-a+2b=3 et a+b=1a+b=1; en additionnant, 3b=43b=4, donc b=43b=\frac{4}{3} et a=143=13a=1-\frac{4}{3}=-\frac{1}{3}. Comme AkA^{k} agit sur chaque vecteur propre en le multipliant par λk\lambda^{k} : uk=132k(11)+435k(21)\vec{u}_{k}=-\frac{1}{3}\cdot 2^{k}\begin{pmatrix}-1\\1\end{pmatrix}+\frac{4}{3}\cdot 5^{k}\begin{pmatrix}2\\1\end{pmatrix}, soit xk=2k3+85k3x_{k}=\frac{2^{k}}{3}+\frac{8\cdot 5^{k}}{3} et yk=2k3+45k3y_{k}=-\frac{2^{k}}{3}+\frac{4\cdot 5^{k}}{3}. Contrôle en k=0k=0 : x0=13+83=3x_{0}=\frac{1}{3}+\frac{8}{3}=3 ✓ et y0=13+43=1y_{0}=-\frac{1}{3}+\frac{4}{3}=1 ✓. Contrôle en k=1k=1 : x1=23+403=14x_{1}=\frac{2}{3}+\frac{40}{3}=14 et directement 4(3)+2(1)=144(3)+2(1)=14 ✓.

d) On calcule xkyk=2k3+85k32k3+45k3=2k+85k2k+45k\frac{x_{k}}{y_{k}}=\frac{\frac{2^{k}}{3}+\frac{8\cdot 5^{k}}{3}}{-\frac{2^{k}}{3}+\frac{4\cdot 5^{k}}{3}}=\frac{2^{k}+8\cdot 5^{k}}{-2^{k}+4\cdot 5^{k}}. En divisant haut et bas par 5k5^{k} : (25)k+8(25)k+4\frac{\left(\frac{2}{5}\right)^{k}+8}{-\left(\frac{2}{5}\right)^{k}+4}. Comme 25<1\frac{2}{5}<1, le terme (25)k\left(\frac{2}{5}\right)^{k} tend vers 0, et la limite vaut 84=2\frac{8}{4}=2. La proportion tend donc vers 2, c'est-à-dire exactement le rapport des composantes du vecteur propre (21)\begin{pmatrix}2\\1\end{pmatrix}. Interprétation : la valeur propre DOMINANTE est celle de plus grand module, ici λ=5\lambda=5, et à long terme sa contribution écrase celle de λ=2\lambda=2 puisque 5k5^{k} croît infiniment plus vite que 2k2^{k}. Le système s'aligne donc sur la direction propre dominante, quelle que soit la condition initiale (pourvu que sa composante bb sur cette direction soit non nulle). C'est le principe qui gouverne les chaînes de Markov, les modèles de populations structurées et l'algorithme PageRank de Google : la valeur propre dominante fixe le TAUX de croissance (ici, la population est multipliée par 5 à chaque période à long terme) et son vecteur propre fixe la RÉPARTITION d'équilibre.

2468101224681012état initialdirection propre DOMINANTEquelle que soit la répartition de départ,le système s'aligne sur cette direction

Exercice 8 : Matrices symétriques : vecteurs propres orthogonaux et diagonalisation orthogonale

Les matrices symétriques réelles occupent une place à part : elles sont toujours diagonalisables, et par une matrice de passage ORTHOGONALE. On l'illustre puis on le démontre.

  • a) Pour A=(3113)A=\begin{pmatrix} 3 & 1 \\ 1 & 3 \end{pmatrix}, déterminez les valeurs propres et un vecteur propre pour chacune, et vérifiez qu'ils sont orthogonaux.
  • b) Normalisez les deux vecteurs propres et formez la matrice PP dont ils sont les colonnes. Vérifiez que PP est orthogonale, c'est-à-dire que PTP=IP^{T}P=I.
  • c) Vérifiez la diagonalisation orthogonale A=PDPTA=PDP^{T}, en utilisant P1=PTP^{-1}=P^{T}.
  • d) Démontrez, pour une matrice symétrique quelconque, que deux vecteurs propres associés à des valeurs propres distinctes sont orthogonaux.

Tape tes réponses, la page te dit juste ou faux 0/3

a)
b)
Voir la correction

Réponses

  • a) λ=4\lambda=4 : (1;1)(1;1) ; λ=2\lambda=2 : (1;1)(1;-1)
  • b) P=12(1111)P=\frac{1}{\sqrt2}\begin{pmatrix}1&1\\1&-1\end{pmatrix}, PTP=IP^TP=I
  • c) PDPT=APDP^T=A
  • d) (λμ)uv=0(\lambda-\mu)\,\vec u\cdot\vec v=0

a) Le polynôme caractéristique est λ26λ+8=(λ4)(λ2)\lambda^{2}-6\lambda+8=(\lambda-4)(\lambda-2), donc les valeurs propres sont λ1=4\lambda_{1}=4 et λ2=2\lambda_{2}=2. Pour λ=4\lambda=4 : A4I=(1111)A-4I=\begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix} donne x=yx=y, d'où v4=(11)\vec{v}_{4}=\begin{pmatrix} 1 \\ 1 \end{pmatrix}. Pour λ=2\lambda=2 : A2I=(1111)A-2I=\begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix} donne x=yx=-y, d'où v2=(11)\vec{v}_{2}=\begin{pmatrix} 1 \\ -1 \end{pmatrix}. Produit scalaire : v4v2=(1)(1)+(1)(1)=0\vec{v}_{4}\cdot\vec{v}_{2}=(1)(1)+(1)(-1)=0, les vecteurs propres sont orthogonaux.

b) Chaque vecteur a pour norme 12+12=2\sqrt{1^{2}+1^{2}}=\sqrt{2}, donc les vecteurs normalisés sont 12(11)\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 1 \end{pmatrix} et 12(11)\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ -1 \end{pmatrix}. On pose P=12(1111)P=\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix} et D=(4002)D=\begin{pmatrix} 4 & 0 \\ 0 & 2 \end{pmatrix}. Vérification de l'orthogonalité : PTP=12(1111)(1111)=12(2002)=IP^{T}P=\dfrac{1}{2}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}=\dfrac{1}{2}\begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}=I. Les colonnes de PP forment donc une base ORTHONORMÉE, ce qui est le propre d'une matrice orthogonale.

c) Comme PP est orthogonale, P1=PTP^{-1}=P^{T}, et la diagonalisation A=PDP1A=PDP^{-1} s'écrit A=PDPTA=PDP^{T}. On vérifie : PD=12(1111)(4002)=12(4242)PD=\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 4 & 0 \\ 0 & 2 \end{pmatrix}=\dfrac{1}{\sqrt{2}}\begin{pmatrix} 4 & 2 \\ 4 & -2 \end{pmatrix}, puis PDPT=12(4242)(1111)=12(6226)=(3113)=APDP^{T}=\dfrac{1}{2}\begin{pmatrix} 4 & 2 \\ 4 & -2 \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}=\dfrac{1}{2}\begin{pmatrix} 6 & 2 \\ 2 & 6 \end{pmatrix}=\begin{pmatrix} 3 & 1 \\ 1 & 3 \end{pmatrix}=A ✓. L'avantage de la diagonalisation orthogonale est qu'elle évite le calcul d'un inverse : la transposée, immédiate, en tient lieu.

d) Soit AA symétrique (AT=AA^{T}=A), et u\vec{u}, v\vec{v} des vecteurs propres pour des valeurs propres distinctes λμ\lambda\neq\mu : Au=λuA\vec{u}=\lambda\vec{u} et Av=μvA\vec{v}=\mu\vec{v}. On calcule le produit scalaire (Au)v(A\vec{u})\cdot\vec{v} de deux façons. D'une part, (Au)v=λuv(A\vec{u})\cdot\vec{v}=\lambda\vec{u}\cdot\vec{v}. D'autre part, en utilisant la symétrie, (Au)v=(Au)Tv=uTATv=uTAv=u(Av)=μuv(A\vec{u})\cdot\vec{v}=(A\vec{u})^{T}\vec{v}=\vec{u}^{T}A^{T}\vec{v}=\vec{u}^{T}A\vec{v}=\vec{u}\cdot(A\vec{v})=\mu\,\vec{u}\cdot\vec{v}. En égalant : λuv=μuv\lambda\,\vec{u}\cdot\vec{v}=\mu\,\vec{u}\cdot\vec{v}, soit (λμ)uv=0(\lambda-\mu)\,\vec{u}\cdot\vec{v}=0. Comme λμ\lambda\neq\mu, on conclut uv=0\vec{u}\cdot\vec{v}=0 : les deux vecteurs propres sont orthogonaux. C'est ce résultat, joint au fait qu'une matrice symétrique réelle a toujours toutes ses valeurs propres réelles, qui garantit la diagonalisation orthogonale : le théorème spectral.

-4-3-2-11234-4-3-2-11234v₁v₂Av₁ = λ₁v₁une matrice SYMÉTRIQUE a toujours des vecteurspropres ORTHOGONAUX : ce n'est pas un hasard

Exercice 9 : Le théorème de Cayley-Hamilton : inverse et puissances sans diagonaliser

Toute matrice annule son propre polynôme caractéristique. Pour une matrice 2×22\times 2, cela fournit une relation A2=tr(A)Adet(A)IA^{2}=\operatorname{tr}(A)\,A-\det(A)\,I qui permet de calculer inverse et puissances sans passer par les vecteurs propres.

  • a) Énoncez le théorème de Cayley-Hamilton pour une matrice 2×22\times 2, puis vérifiez-le sur A=(4213)A=\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}.
  • b) Déduisez-en une expression de A1A^{-1} en fonction de AA et II, et calculez-la.
  • c) Déduisez-en que A3A^{3} s'écrit comme combinaison linéaire de AA et II, et donnez sa valeur.
  • d) Expliquez pourquoi toute puissance AnA^{n} et l'inverse A1A^{-1} s'expriment comme combinaison de AA et II, et comparez cette méthode à la diagonalisation.

Tape tes réponses, la page te dit juste ou faux 0/5

b)
c)
d)
Voir la correction

Réponses

  • a) A27A+10I=0A^2-7A+10I=0
  • b) A1=110(3214)A^{-1}=\frac{1}{10}\begin{pmatrix}3&-2\\-1&4\end{pmatrix}
  • c) A3=39A70IA^3=39A-70I
  • d) Récurrence universelle contre formule close

a) Théorème de Cayley-Hamilton : une matrice annule son polynôme caractéristique. Pour une matrice 2×22\times 2, celui-ci est λ2tr(A)λ+det(A)\lambda^{2}-\operatorname{tr}(A)\lambda+\det(A), donc A2tr(A)A+det(A)I=0A^{2}-\operatorname{tr}(A)\,A+\det(A)\,I=0. Ici tr(A)=7\operatorname{tr}(A)=7 et det(A)=10\det(A)=10, la relation attendue est A27A+10I=0A^{2}-7A+10I=0. Vérification : A2=(1814711)A^{2}=\begin{pmatrix} 18 & 14 \\ 7 & 11 \end{pmatrix}, 7A=(2814721)7A=\begin{pmatrix} 28 & 14 \\ 7 & 21 \end{pmatrix}, donc A27A+10I=(1828+101414771121+10)=(0000)A^{2}-7A+10I=\begin{pmatrix} 18-28+10 & 14-14 \\ 7-7 & 11-21+10 \end{pmatrix}=\begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix} ✓.

b) On isole le terme constant dans A27A+10I=0A^{2}-7A+10I=0 : 10I=7AA2=A(7IA)10I=7A-A^{2}=A(7I-A), donc A110(7IA)=IA\cdot\dfrac{1}{10}(7I-A)=I, ce qui identifie l'inverse : A1=110(7IA)=110((7007)(4213))=110(3214)A^{-1}=\dfrac{1}{10}(7I-A)=\dfrac{1}{10}\left(\begin{pmatrix} 7 & 0 \\ 0 & 7 \end{pmatrix}-\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}\right)=\dfrac{1}{10}\begin{pmatrix} 3 & -2 \\ -1 & 4 \end{pmatrix}. Vérification : AA1=110(4213)(3214)=110(100010)=IA\,A^{-1}=\dfrac{1}{10}\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}\begin{pmatrix} 3 & -2 \\ -1 & 4 \end{pmatrix}=\dfrac{1}{10}\begin{pmatrix} 10 & 0 \\ 0 & 10 \end{pmatrix}=I ✓. On retrouve d'ailleurs la formule classique de l'inverse 2×22\times 2, mais dérivée du théorème.

c) De A2=7A10IA^{2}=7A-10I, on multiplie par AA : A3=7A210A=7(7A10I)10A=49A70I10A=39A70IA^{3}=7A^{2}-10A=7(7A-10I)-10A=49A-70I-10A=39A-70I. Numériquement : A3=39(4213)70I=(15670783911770)=(86783947)A^{3}=39\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}-70I=\begin{pmatrix} 156-70 & 78 \\ 39 & 117-70 \end{pmatrix}=\begin{pmatrix} 86 & 78 \\ 39 & 47 \end{pmatrix}. C'est exactement le A3A^{3} obtenu par PD3P1PD^{3}P^{-1} dans le problème précédent : deux méthodes, un seul résultat.

d) La relation A2=7A10IA^{2}=7A-10I permet de remplacer tout A2A^{2} par une combinaison de AA et II. Par récurrence, chaque puissance AnA^{n} se ramène de proche en proche à une telle combinaison αnA+βnI\alpha_{n}A+\beta_{n}I : dès qu'un A2A^{2} apparaît, on le réduit. L'inverse aussi, comme on l'a vu en b). Comparaison avec la diagonalisation : la méthode An=PDnP1A^{n}=PD^{n}P^{-1} donne une formule CLOSE (dépendant explicitement de 2n2^{n} et 5n5^{n}) et fonctionne même pour des exposants symboliques, mais elle exige de diagonaliser, donc de trouver vecteurs propres et inverse de PP. Cayley-Hamilton, lui, ne demande que la trace et le déterminant, s'applique même aux matrices NON diagonalisables, mais calcule les puissances de proche en proche plutôt que par une formule close. Les deux approches sont complémentaires : formule explicite d'un côté, récurrence universelle de l'autre.

Exercice 10 : Démonstrations : valeurs propres de l'inverse, des puissances et des matrices translatées

Connaître les valeurs propres de AA donne gratuitement celles de nombreuses matrices dérivées, avec les MÊMES vecteurs propres. Ces résultats se démontrent en une ligne et évitent des calculs répétés.

  • a) Démontrez que si λ\lambda est une valeur propre de AA de vecteur propre v\vec{v}, alors λk\lambda^{k} est valeur propre de AkA^{k}, de même vecteur propre.
  • b) Démontrez que si AA est inversible et λ\lambda une valeur propre, alors λ0\lambda\neq 0 et 1λ\dfrac{1}{\lambda} est valeur propre de A1A^{-1}, de même vecteur propre.
  • c) Démontrez que λ+c\lambda+c est valeur propre de A+cIA+cI, de même vecteur propre, pour tout scalaire cc.
  • d) Application à A=(4213)A=\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}, de valeurs propres 22 et 55 : donnez sans nouveau calcul les valeurs propres de A3A^{3}, de A1A^{-1} et de A+2IA+2I, et vérifiez la relation de trace sur A+2IA+2I.

Tape tes réponses, la page te dit juste ou faux 0/4

d)
Voir la correction

Réponses

  • a) Akv=λkvA^k\vec v=\lambda^k\vec v
  • b) A1v=1λvA^{-1}\vec v=\frac1\lambda\vec v
  • c) (A+cI)v=(λ+c)v(A+cI)\vec v=(\lambda+c)\vec v
  • d) 88, 125125 ; 12\frac12, 15\frac15 ; 44, 77

a) On part de Av=λvA\vec{v}=\lambda\vec{v} et on applique AA de façon répétée. A2v=A(Av)=A(λv)=λAv=λ2vA^{2}\vec{v}=A(A\vec{v})=A(\lambda\vec{v})=\lambda A\vec{v}=\lambda^{2}\vec{v}. En itérant, ou par récurrence, Akv=λkvA^{k}\vec{v}=\lambda^{k}\vec{v}. Donc λk\lambda^{k} est valeur propre de AkA^{k}, avec le même vecteur propre v\vec{v} : élever la matrice à une puissance élève les valeurs propres à cette puissance, sans toucher aux directions propres.

b) Si AA est inversible, aucune valeur propre n'est nulle : λ=0\lambda=0 signifierait Av=0A\vec{v}=\vec{0} pour un v0\vec{v}\neq\vec{0}, donc un noyau non trivial, ce qui contredit l'inversibilité. On peut donc écrire 1λ\dfrac{1}{\lambda}. En appliquant A1A^{-1} à Av=λvA\vec{v}=\lambda\vec{v} : v=A1(λv)=λA1v\vec{v}=A^{-1}(\lambda\vec{v})=\lambda A^{-1}\vec{v}, d'où A1v=1λvA^{-1}\vec{v}=\dfrac{1}{\lambda}\vec{v}. Donc 1λ\dfrac{1}{\lambda} est valeur propre de A1A^{-1}, de même vecteur propre.

c) (A+cI)v=Av+cIv=λv+cv=(λ+c)v(A+cI)\vec{v}=A\vec{v}+cI\vec{v}=\lambda\vec{v}+c\vec{v}=(\lambda+c)\vec{v}. Donc λ+c\lambda+c est valeur propre de A+cIA+cI, avec le même vecteur propre. Ajouter cIcI translate toutes les valeurs propres de cc sans modifier les directions propres, puisque cIcI agit de la même façon sur tous les vecteurs.

d) Les valeurs propres de AA sont 22 et 55. Par a), celles de A3A^{3} sont 23=82^{3}=8 et 53=1255^{3}=125. Par b), celles de A1A^{-1} sont 12\dfrac{1}{2} et 15\dfrac{1}{5}. Par c), celles de A+2IA+2I sont 2+2=42+2=4 et 5+2=75+2=7. Vérification de la relation de trace sur A+2I=(6215)A+2I=\begin{pmatrix} 6 & 2 \\ 1 & 5 \end{pmatrix} : sa trace vaut 6+5=116+5=11, et la somme des valeurs propres 4+7=114+7=11 ✓. Tous ces résultats s'obtiennent sans reprendre le moindre polynôme caractéristique : c'est l'intérêt pratique de ces trois propriétés, qui reviennent constamment dans les problèmes de suites et de systèmes dynamiques. Notons enfin que AA et ATA^{T} ont exactement les mêmes valeurs propres, car det(ATλI)=det((AλI)T)=det(AλI)\det(A^{T}-\lambda I)=\det((A-\lambda I)^{T})=\det(A-\lambda I), même si leurs vecteurs propres, eux, diffèrent en général.

Partie C : les classiques (/50)

Exercice 11 : Une matrice 3×3 à valeur propre double

On considère A=(211121112)A=\begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix}, que l'on peut écrire A=I+JA=I+J, où JJ est la matrice dont tous les coefficients valent 11.

  • a) Montrez que (1; 1; 1)(1;\ 1;\ 1) est un vecteur propre de AA et donnez sa valeur propre.
  • b) Montrez que Av=v+(x+y+z)(1; 1; 1)A\vec{v}=\vec{v}+(x+y+z)(1;\ 1;\ 1) pour v=(x; y; z)\vec{v}=(x;\ y;\ z). Déduisez-en une deuxième valeur propre et la dimension de son sous-espace propre.
  • c) Contrôlez avec la trace et le déterminant de AA.
  • d) Écrivez PP et DD telles que A=PDP1A=PDP^{-1}, et justifiez que PP est inversible.
  • e) Sachant que J2=3JJ^{2}=3J, montrez que An=I+4n13JA^{n}=I+\frac{4^{n}-1}{3}J et donnez les coefficients de A3A^{3}.

Tape tes réponses, la page te dit juste ou faux 0/8

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) λ=4\lambda=4
  • b) λ=1\lambda=1 sur le plan x+y+z=0x+y+z=0
  • c) Trace 6, déterminant 4
  • d) detP=3\det P=3 : diagonalisable
  • e) A3A^3 : 22 et 21

a) A(1; 1; 1)=(4; 4; 4)=4(1; 1; 1)A(1;\ 1;\ 1)=(4;\ 4;\ 4)=4(1;\ 1;\ 1) : valeur propre 44.

b) Jv=(x+y+z)(1; 1; 1)J\vec{v}=(x+y+z)(1;\ 1;\ 1), donc Av=v+JvA\vec{v}=\vec{v}+J\vec{v} donne la formule. Si x+y+z=0x+y+z=0, alors Av=vA\vec{v}=\vec{v} : tout vecteur non nul du plan x+y+z=0x+y+z=0 est propre pour la valeur propre 11. Ce plan est de dimension 22, de base (1; 1; 0)(1;\ -1;\ 0) et (1; 0; 1)(1;\ 0;\ -1).

c) La somme des valeurs propres, comptées avec leur multiplicité, vaut 4+1+1=6=trA4+1+1=6=\operatorname{tr}A. Leur produit vaut 4×1×1=44\times 1\times 1=4, et detA=2(41)1(21)+1(12)=611=4\det A=2(4-1)-1(2-1)+1(1-2)=6-1-1=4. Les deux contrôles concordent.

d) P=(111110101)P=\begin{pmatrix} 1 & 1 & 1 \\ 1 & -1 & 0 \\ 1 & 0 & -1 \end{pmatrix}, dont les colonnes sont (1; 1; 1)(1;\ 1;\ 1), (1; 1; 0)(1;\ -1;\ 0) et (1; 0; 1)(1;\ 0;\ -1), et D=diag(4; 1; 1)D=\text{diag}(4;\ 1;\ 1). detP=1(10)1(10)+1(0+1)=30\det P=1(1-0)-1(-1-0)+1(0+1)=3\neq 0 : les trois vecteurs propres sont indépendants. La valeur propre double a bien un sous-espace propre de dimension 22 : AA est diagonalisable.

e) Par récurrence : vrai pour n=1n=1, puisque 413=1\frac{4-1}{3}=1. Si An=I+cnJA^{n}=I+c_{n}J, alors An+1=(I+cnJ)(I+J)=I+(1+cn+3cn)JA^{n+1}=(I+c_{n}J)(I+J)=I+(1+c_{n}+3c_{n})J, et 1+4cn=1+4n+143=4n+1131+4c_{n}=1+\frac{4^{n+1}-4}{3}=\frac{4^{n+1}-1}{3}. Pour n=3n=3 : 633=21\frac{63}{3}=21, donc A3A^{3} a 2222 sur la diagonale et 2121 ailleurs.

Le fil : une valeur propre double ne bloque rien si son sous-espace propre est de dimension 2 ; écrire A=I+JA=I+J rend les valeurs propres presque évidentes.

Exercice 12 : Les valeurs propres lues sur la géométrie

On considère la projection orthogonale P=12(1111)P=\frac{1}{2}\begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix} sur la droite y=xy=x, la réflexion S=(0110)S=\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} d'axe y=xy=x, l'homothétie H=3IH=3I et le cisaillement C=(1201)C=\begin{pmatrix} 1 & 2 \\ 0 & 1 \end{pmatrix}.

  • a) Sans polynôme caractéristique, donnez les valeurs propres de PP et des vecteurs propres associés. Vérifiez par le calcul.
  • b) Même question pour SS.
  • c) Même question pour HH.
  • d) Déterminez les valeurs propres de CC et la dimension du sous-espace propre. CC est-elle diagonalisable ?
  • e) Exprimez PP en fonction de SS et de II, et retrouvez ainsi les valeurs propres de PP à partir de celles de SS.

Tape tes réponses, la page te dit juste ou faux 0/8

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) 1 et 0
  • b) 1 et 1-1
  • c) 3 double, tout vecteur
  • d) 1 double, dimension 1 : non
  • e) P=S+I2P=\frac{S+I}{2}

a) Un vecteur de la droite y=xy=x n'est pas modifié par la projection : (1; 1)(1;\ 1) est propre pour 11. Un vecteur perpendiculaire à la droite est écrasé sur l'origine : (1; 1)(1;\ -1) est propre pour 00. Calcul : P(1; 1)=(1; 1)P(1;\ 1)=(1;\ 1) et P(1; 1)=(0; 0)P(1;\ -1)=(0;\ 0).

b) La réflexion laisse fixe l'axe, (1; 1)(1;\ 1) est propre pour 11, et retourne la direction perpendiculaire, (1; 1)(1;\ -1) est propre pour 1-1. Calcul : S(1; 1)=(1; 1)=(1; 1)S(1;\ -1)=(-1;\ 1)=-(1;\ -1).

c) HH multiplie tout vecteur par 33 : 33 est l'unique valeur propre, double, et TOUT vecteur non nul est propre. Le sous-espace propre est le plan entier.

d) CC est triangulaire : valeur propre 11, double. CI=(0200)C-I=\begin{pmatrix} 0 & 2 \\ 0 & 0 \end{pmatrix} impose y=0y=0 : le sous-espace propre est l'axe des xx, de dimension 11. Géométriquement, seule la direction horizontale n'est pas inclinée par le cisaillement. Dimension 1<21<2 : CC n'est pas diagonalisable.

e) S+I2=12(1111)=P\frac{S+I}{2}=\frac{1}{2}\begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}=P. Les valeurs propres de SS, 11 et 1-1, deviennent celles de PP par λλ+12\lambda\mapsto\frac{\lambda+1}{2} : 11 et 00, avec les mêmes vecteurs propres, comme le prévoit l'exercice 10.

Le fil : une direction que la transformation ne fait qu'étirer est une direction propre ; la géométrie donne souvent les valeurs propres avant tout calcul.

Exercice 13 : Diagonalisable ou non selon un paramètre

Pour tout réel aa, on considère Aa=(1a02)A_{a}=\begin{pmatrix} 1 & a \\ 0 & 2 \end{pmatrix}, Ba=(2a02)B_{a}=\begin{pmatrix} 2 & a \\ 0 & 2 \end{pmatrix} et Ca=(11a1)C_{a}=\begin{pmatrix} 1 & 1 \\ a & 1 \end{pmatrix}.

  • a) Donnez les valeurs propres de AaA_{a}. Pour quelles valeurs de aa la matrice est-elle diagonalisable ?
  • b) Même question pour BaB_{a}.
  • c) Calculez les valeurs propres de CaC_{a} en fonction de aa, pour a0a\geq 0.
  • d) Discutez la diagonalisabilité de CaC_{a} dans R\mathbb{R} selon le signe de aa.
  • e) Pour a=4a=4, déterminez des vecteurs propres de C4C_{4}.

Tape tes réponses, la page te dit juste ou faux 0/7

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) Toujours diagonalisable
  • b) Seulement pour a=0a=0
  • c) 1±a1\pm\sqrt a
  • d) Oui si a>0a>0, non sinon
  • e) (1; 2)(1;\ 2) et (1; 2)(1;\ -2)

a) AaA_{a} est triangulaire : valeurs propres 11 et 22, distinctes quel que soit aa. Deux valeurs propres distinctes en dimension 2 suffisent : AaA_{a} est diagonalisable pour tout aa. Un vecteur propre pour 22 est (a; 1)(a;\ 1).

b) BaB_{a} a la seule valeur propre 22, double. Ba2I=(0a00)B_{a}-2I=\begin{pmatrix} 0 & a \\ 0 & 0 \end{pmatrix} : si a0a\neq 0, le sous-espace propre est la droite y=0y=0, de dimension 11, et BaB_{a} n'est pas diagonalisable ; si a=0a=0, B0=2IB_{0}=2I est déjà diagonale. Diagonalisable si et seulement si a=0a=0.

c) det(CaλI)=(1λ)2a=0\det(C_{a}-\lambda I)=(1-\lambda)^{2}-a=0 donne (1λ)2=a(1-\lambda)^{2}=a, donc λ=1±a\lambda=1\pm\sqrt{a}.

d) Si a>0a>0 : deux valeurs propres réelles distinctes, CaC_{a} est diagonalisable. Si a=0a=0 : valeur propre 11 double, C0I=(0100)C_{0}-I=\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} a un noyau de dimension 11, C0C_{0} n'est pas diagonalisable. Si a<0a<0 : valeurs propres complexes 1±ia1\pm i\sqrt{-a}, aucune direction propre réelle, CaC_{a} n'est pas diagonalisable dans R\mathbb{R}.

e) Valeurs propres 33 et 1-1. C43I=(2142)C_{4}-3I=\begin{pmatrix} -2 & 1 \\ 4 & -2 \end{pmatrix} donne y=2xy=2x : vecteur propre (1; 2)(1;\ 2). C4+I=(2142)C_{4}+I=\begin{pmatrix} 2 & 1 \\ 4 & 2 \end{pmatrix} donne y=2xy=-2x : vecteur propre (1; 2)(1;\ -2). Contrôle : C4(1; 2)=(3; 6)C_{4}(1;\ 2)=(3;\ 6).

Le fil : des valeurs propres distinctes garantissent la diagonalisation ; une valeur propre multiple exige de vérifier la dimension du sous-espace propre.

Exercice 14 : Problème : les modes propres de deux chariots à ressorts

Deux chariots de même masse mm glissent sans frottement entre deux murs. Trois ressorts identiques de raideur kk relient le mur de gauche au premier chariot, les deux chariots entre eux, et le second chariot au mur de droite. On note x1x_{1} et x2x_{2} les déplacements des chariots par rapport à l'équilibre, et l'on prend km=4 s2\frac{k}{m}=4\ \mathrm{s^{-2}}.

La deuxième loi de Newton donne x1=4(2x1+x2)x_{1}''=4(-2x_{1}+x_{2}) et x2=4(x12x2)x_{2}''=4(x_{1}-2x_{2}), soit X=KXX''=-KX. On admet que si Kv=ω2vK\vec{v}=\omega^{2}\vec{v}, alors X(t)=vcos(ωt)X(t)=\vec{v}\cos(\omega t) est une solution : un MODE PROPRE, où les deux chariots oscillent à la même pulsation ω\omega.

  • a) Écrivez la matrice KK. Quelle propriété a-t-elle ?
  • b) Calculez les valeurs propres de KK.
  • c) Déterminez un vecteur propre pour chacune, et décrivez le mouvement des chariots dans chaque mode.
  • d) Calculez les deux pulsations ω\omega et les périodes correspondantes, au centième.
  • e) On écarte le premier chariot de 2 cm2\ \mathrm{cm}, le second restant en place, et on lâche sans vitesse. Décomposez X(0)=(2; 0)X(0)=(2;\ 0) sur les deux modes et écrivez x1(t)x_{1}(t) et x2(t)x_{2}(t).

Tape tes réponses, la page te dit juste ou faux 0/10

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) K=(8448)K=\begin{pmatrix}8&-4\\-4&8\end{pmatrix} symétrique
  • b) 4 et 12
  • c) (1;1)(1;1) en phase, (1;1)(1;-1) en opposition
  • d) 22 et 3,463{,}46 rad/s ; 3,143{,}14 s et 1,811{,}81 s
  • e) x1=cos2t+cos12tx_1=\cos 2t+\cos\sqrt{12}t, x2=cos2tcos12tx_2=\cos 2t-\cos\sqrt{12}t

a) K=(8448)K=\begin{pmatrix} 8 & -4 \\ -4 & 8 \end{pmatrix} : elle est SYMÉTRIQUE, donc diagonalisable avec des vecteurs propres orthogonaux.

b) λ216λ+(6416)=λ216λ+48=(λ4)(λ12)\lambda^{2}-16\lambda+(64-16)=\lambda^{2}-16\lambda+48=(\lambda-4)(\lambda-12) : valeurs propres 44 et 1212.

c) Pour 44 : K4I=(4444)K-4I=\begin{pmatrix} 4 & -4 \\ -4 & 4 \end{pmatrix}, donc x1=x2x_{1}=x_{2}, vecteur (1; 1)(1;\ 1) : les chariots bougent ensemble, EN PHASE, et le ressort du milieu n'est jamais étiré. Pour 1212 : K12I=(4444)K-12I=\begin{pmatrix} -4 & -4 \\ -4 & -4 \end{pmatrix}, donc x1=x2x_{1}=-x_{2}, vecteur (1; 1)(1;\ -1) : les chariots bougent en OPPOSITION, symétriquement.

d) ω1=4=2 rad/s\omega_{1}=\sqrt{4}=2\ \mathrm{rad/s} et ω2=123,46 rad/s\omega_{2}=\sqrt{12}\approx 3{,}46\ \mathrm{rad/s}. Périodes T=2πωT=\frac{2\pi}{\omega} : π3,14 s\pi\approx 3{,}14\ \mathrm{s} et 2π121,81 s\frac{2\pi}{\sqrt{12}}\approx 1{,}81\ \mathrm{s}. Le mode en opposition est plus rapide, parce que le ressort du milieu y travaille aussi.

e) (2; 0)=1(1; 1)+1(1; 1)(2;\ 0)=1\cdot(1;\ 1)+1\cdot(1;\ -1). Chaque mode évolue à sa pulsation : x1(t)=cos(2t)+cos(12t)x_{1}(t)=\cos(2t)+\cos(\sqrt{12}\,t) et x2(t)=cos(2t)cos(12t)x_{2}(t)=\cos(2t)-\cos(\sqrt{12}\,t), en centimètres. Les deux pulsations n'étant pas dans un rapport entier simple, le mouvement n'est pas périodique : l'énergie passe d'un chariot à l'autre.

Le fil : diagonaliser un système couplé le découpe en modes indépendants ; les valeurs propres donnent les fréquences, les vecteurs propres les formes du mouvement.

Exercice 15 : Problème : les axes principaux d'une ellipse inclinée

On étudie la courbe d'équation 5x2+4xy+2y2=65x^{2}+4xy+2y^{2}=6. On pose X=(xy)X=\begin{pmatrix} x \\ y \end{pmatrix} et A=(5222)A=\begin{pmatrix} 5 & 2 \\ 2 & 2 \end{pmatrix}.

  • a) Vérifiez que XTAX=5x2+4xy+2y2X^{T}AX=5x^{2}+4xy+2y^{2}.
  • b) Calculez les valeurs propres de AA.
  • c) Déterminez un vecteur propre pour chacune et vérifiez qu'ils sont orthogonaux.
  • d) Dans le repère orthonormé porté par les vecteurs propres normalisés, avec les coordonnées xx' selon le vecteur de la plus grande valeur propre et yy' selon l'autre, l'équation devient 6x2+y2=66x'^{2}+y'^{2}=6. Justifiez, puis donnez les demi-axes de l'ellipse, au centième.
  • e) Donnez l'angle du grand axe avec l'axe des xx, au centième de degré, et un point de la courbe à distance maximale de l'origine.

Tape tes réponses, la page te dit juste ou faux 0/6

b)
c)
d)
e)
Voir la correction

Réponses

  • a) XTAX=5x2+4xy+2y2X^TAX=5x^2+4xy+2y^2
  • b) 6 et 1
  • c) (2; 1)(2;\ 1) et (1; 2)(1;\ -2)
  • d) Demi-axes 1 et 62,45\sqrt6\approx 2{,}45
  • e) 63,43-63{,}43^\circ

a) AX=(5x+2y; 2x+2y)AX=(5x+2y;\ 2x+2y) et XTAX=x(5x+2y)+y(2x+2y)=5x2+4xy+2y2X^{T}AX=x(5x+2y)+y(2x+2y)=5x^{2}+4xy+2y^{2}. Le coefficient 44 du terme croisé se partage en deux 22 hors diagonale, ce qui rend AA symétrique.

b) λ27λ+6=(λ1)(λ6)\lambda^{2}-7\lambda+6=(\lambda-1)(\lambda-6) : valeurs propres 66 et 11.

c) Pour 66 : A6I=(1224)A-6I=\begin{pmatrix} -1 & 2 \\ 2 & -4 \end{pmatrix}, donc x=2yx=2y : (2; 1)(2;\ 1). Pour 11 : AI=(4221)A-I=\begin{pmatrix} 4 & 2 \\ 2 & 1 \end{pmatrix}, donc y=2xy=-2x : (1; 2)(1;\ -2). Produit scalaire 22=02-2=0 : orthogonaux, comme pour toute matrice symétrique.

d) Avec PP orthogonale de colonnes 15(2; 1)\frac{1}{\sqrt{5}}(2;\ 1) et 15(1; 2)\frac{1}{\sqrt{5}}(1;\ -2), on pose X=PXX=PX' : XTAX=XT(PTAP)X=XTDX=6x2+y2X^{T}AX=X'^{T}(P^{T}AP)X'=X'^{T}DX'=6x'^{2}+y'^{2}. L'équation x2+y26=1x'^{2}+\frac{y'^{2}}{6}=1 est celle d'une ellipse de demi-axes 11, selon (2; 1)(2;\ 1), et 62,45\sqrt{6}\approx 2{,}45, selon (1; 2)(1;\ -2). Le terme croisé a disparu : c'est l'intérêt de la diagonalisation orthogonale.

e) Le grand axe est dirigé par (1; 2)(1;\ -2), d'angle arctan(2)63,43\arctan(-2)\approx -63{,}43^{\circ}. Le point 615(1; 2)\sqrt{6}\cdot\frac{1}{\sqrt{5}}(1;\ -2) a pour carrés de coordonnées 65\frac{6}{5} et 245\frac{24}{5} et pour produit 125-\frac{12}{5} : 565+4(125)+2245=6485+485=65\cdot\frac{6}{5}+4\cdot\left(-\frac{12}{5}\right)+2\cdot\frac{24}{5}=6-\frac{48}{5}+\frac{48}{5}=6. Il est sur la courbe, à la distance 6\sqrt{6} de l'origine.

Le fil : une forme quadratique se lit sur sa matrice symétrique ; ses vecteurs propres donnent les axes, ses valeurs propres les longueurs.

Chapitre précédent Transformations linéaires du plan Chapitre suivant Modèles matriciels : Markov, Leontief, Leslie

Ce chapitre resservira dans

Les chapitres qui le réclament en amont, plus tard dans l'année ou dans les années suivantes.

Voir aussi

Vous cherchez un tuteur en algèbre linéaire à Montréal ?

Contactez-moi pour une première séance. On travaille l'algèbre linéaire au niveau réel des évaluations, y compris la diagonalisation, le chapitre qui donne enfin un sens géométrique à tout le calcul matriciel de la session.

Site par Studio Squalli