Aller au contenu principal

Réduction des endomorphismes

Valeurs propres et vecteurs propres​

Ici EE désigne un espace vectoriel et f∈L(E)f \in \mathcal{L}(E).

Définition 1. Soit λ∈K\lambda \in \mathbb{K} ; on dit que λ\lambda est une valeur propre (notée vp) de ff s'il existe x∈Ex \in E non nul tel que f(x)=λxf(x) = \lambda x. On dit que x∈Ex \in E est vecteur propre de ff si x≠0x \neq 0 et s'il existe λ∈K\lambda \in \mathbb{K} tel que f(x)=λxf(x) = \lambda x. L'ensemble des valeurs propres de ff est le spectre de ff, noté sp(f)sp(f) ; c'est un sous-ensemble de K\mathbb{K}.

Définition 1 bis. Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}). De même : λ\lambda est valeur propre de AA si ∃X≠0\exists X \neq 0 tel que AX=λXAX = \lambda X, où X∈RnX \in \mathbb{R}^n, et XX est vecteur propre de AA si X≠0X \neq 0 et ∃λ∈K\exists \lambda \in \mathbb{K} tel que AX=λXAX = \lambda X.

Remarque. Soit f∈L(E)f \in \mathcal{L}(E), où dim⁡(E)=n\dim(E) = n. Soit BB une base de EE et A=matB(f)A = mat_B(f). Alors λ\lambda est valeur propre de ff donc λ\lambda est valeur propre de AA. De plus, xx est vecteur propre de ff donc XX est vecteur propre de AA, où X=matB(x)X = mat_B(x).

Proposition 2. Soit f∈L(E)f \in \mathcal{L}(E) et EE un espace vectoriel de dimension finie (n=dim⁡(E)n = \dim(E)). λ∈sp(f)⇔f−λId\lambda \in sp(f) \Leftrightarrow f - \lambda I_d est non injectif ⇔Ker(f−λId)≠{0}⇔rg(f−λId)<n\Leftrightarrow \text{Ker}(f - \lambda I_d) \neq \{0\} \Leftrightarrow rg(f - \lambda I_d) < n.

Proposition 2 bis. Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}). Alors λ∈sp(A)⇔Ker(A−λIn)≠{0}⇔A−λIn∉GLn(K)⇔rg(A−λIn)<n\lambda \in sp(A) \Leftrightarrow \text{Ker}(A - \lambda I_n) \neq \{0\} \Leftrightarrow A - \lambda I_n \notin GL_n(\mathbb{K}) \Leftrightarrow rg(A - \lambda I_n) < n. En particulier : A∈GLn(K)⇔0∉sp(A)A \in GL_n(\mathbb{K}) \Leftrightarrow 0 \notin sp(A).

Démonstration. Soit λ∈sp(f)\lambda \in sp(f). Donc ∃x≠0\exists x \neq 0 dans EE tel que f(x)=λxf(x) = \lambda x. Donc (f−λId)(x)=0(f - \lambda I_d)(x) = 0. Donc x∈Ker(f−λId)x \in \text{Ker}(f - \lambda I_d). Donc Ker(f−λId)≠{0}\text{Ker}(f - \lambda I_d) \neq \{0\} et f−λIdf - \lambda I_d est non injectif. La réciproque est vraie : si Ker(f−λId)≠{0}\text{Ker}(f - \lambda I_d) \neq \{0\}, ∃x≠0\exists x \neq 0 dans Ker(f−λId)\text{Ker}(f - \lambda I_d), c'est-à-dire {(f−λId)(x)=0x≠0⇔{f(x)=λxx≠0⇔λ∈sp(f)\begin{cases} (f - \lambda I_d)(x) = 0 \\ x \neq 0 \end{cases} \Leftrightarrow \begin{cases} f(x) = \lambda x \\ x \neq 0 \end{cases} \Leftrightarrow \lambda \in sp(f). Puisque dim⁡(E)<+∞\dim(E) < +\infty, f−λIdf - \lambda I_d non injectif, c'est-à-dire f−λIdf - \lambda I_d non surjectif, donc rg(f−λId)<nrg(f - \lambda I_d) < n. □\square

Définition - Proposition 3. Soit f∈L(E)f \in \mathcal{L}(E) et λ∈sp(f)\lambda \in sp(f). Si xx est tel que x≠0x \neq 0 et f(x)=λxf(x) = \lambda x, on dit que λ\lambda et xx sont valeur propre et vecteur propre associés. On appelle sous-espace propre associé à λ\lambda, noté parfois sep(λ)sep(\lambda), le sev Ker(f−λId)\text{Ker}(f - \lambda I_d). Il est formé des vecteurs propres de ff associés à λ\lambda, ainsi que de 0.

Démonstration. Soit x∈sep(λ)x \in sep(\lambda), alors (f−λId)(x)=0(f - \lambda I_d)(x) = 0, c'est-à-dire f(x)=λxf(x) = \lambda x : xx est vecteur propre associé à λ\lambda (ou bien x=0x = 0). □\square

Remarque. Soit f∈L(E)f \in \mathcal{L}(E) et λ∈sp(f)\lambda \in sp(f). Alors sep(λ)sep(\lambda) est un sev de EE.

Remarque. sep(λ)sep(\lambda) est un sev stable par ff. De plus, la restriction de ff à ce sev est l'homothétie de rapport λ\lambda :

f∣sep(λ)=λIdf_{|sep(\lambda)} = \lambda I_d

Rappel. Soit EE un espace vectoriel ; l'homothétie de rapport λ\lambda est l'application λId\lambda I_d. On la note hλh_\lambda ; on a ∀x∈E,h(x)=λx\forall x \in E, h(x) = \lambda x. Sa matrice dans n'importe quelle base est (λ…0⋱0…λ)=λIn\begin{pmatrix} \lambda & \dots & 0 \\ & \ddots & \\ 0 & \dots & \lambda \end{pmatrix} = \lambda I_n. On remarque que sp(h)={λ}sp(h) = \{\lambda\} et tout vecteur est vecteur propre associé.

Exemples. Soit u∈L(R2)u \in \mathcal{L}(\mathbb{R}^2) défini par u(e1)=2e1u(e_1) = 2e_1 et u(e2)=−e2u(e_2) = -e_2. Alors 2∈sp(u)2 \in sp(u) et e1∈sep(2)e_1 \in sep(2), −1∈sp(u)-1 \in sp(u) et e2∈sep(−1)e_2 \in sep(-1). La matrice de uu dans B=(e1,e2)B = (e_1, e_2) est (200−1)\begin{pmatrix} 2 & 0 \\ 0 & -1 \end{pmatrix}.

Soit A=(1…1⋮⋮1…1)A = \begin{pmatrix} 1 & \dots & 1 \\ \vdots & & \vdots \\ 1 & \dots & 1 \end{pmatrix}. Déterminer les valeurs propres de AA, et les vecteurs propres associés.

On résout l'équation aux valeurs propres AX=λXAX = \lambda X ; ici λ\lambda et XX sont inconnus. Notons X=(x1⋮xn)X = \begin{pmatrix} x_1 \\ \vdots \\ x_n \end{pmatrix}. Alors AX=λX⇔{x1+⋯+xn=λx1⋮x1+⋯+xn=λxnAX = \lambda X \Leftrightarrow \begin{cases} x_1 + \dots + x_n = \lambda x_1 \\ \vdots \\ x_1 + \dots + x_n = \lambda x_n \end{cases}. Donc λx1=⋯=λxn\lambda x_1 = \dots = \lambda x_n. Donc soit λ=0\lambda = 0, soit x1=⋯=xnx_1 = \dots = x_n.

  • Si λ=0\lambda = 0, alors x1+⋯+xn=0x_1 + \dots + x_n = 0.
  • Si x1=⋯=xnx_1 = \dots = x_n, alors nx1=λx1n x_1 = \lambda x_1 ; si x1=0x_1 = 0, alors X=0X = 0. Or on cherche des X≠0X \neq 0. Donc x1≠0x_1 \neq 0.

Conclusion : sp(A)={0,n}sp(A) = \{0, n\}. De plus, sep(0)={X∈Rn,x1+⋯+xn=0}sep(0) = \{X \in \mathbb{R}^n, x_1 + \dots + x_n = 0\}. C'est un hyperplan de Rn\mathbb{R}^n. Une base en est (1−10⋮0),(10−10⋮0),…,(10⋮0−1)\begin{pmatrix} 1 \\ -1 \\ 0 \\ \vdots \\ 0 \end{pmatrix}, \begin{pmatrix} 1 \\ 0 \\ -1 \\ 0 \\ \vdots \\ 0 \end{pmatrix}, \dots, \begin{pmatrix} 1 \\ 0 \\ \vdots \\ 0 \\ -1 \end{pmatrix}. On a aussi sep(n)={X∈Rn,x1=⋯=xn}sep(n) = \{X \in \mathbb{R}^n, x_1 = \dots = x_n\}. C'est une droite engendrée par (1⋮1)\begin{pmatrix} 1 \\ \vdots \\ 1 \end{pmatrix}.

Remarque. Soit ϕ:{Rn→RX↦∑i=1nxi\phi : \begin{cases} \mathbb{R}^n \to \mathbb{R} \\ X \mapsto \sum_{i=1}^{n} x_i \end{cases}. C'est une forme linéaire et donc sep(0)=Ker(ϕ)sep(0) = \text{Ker}(\phi) est bien un hyperplan. La base donnée doit être vérifiée (libre et génératrice).

Remarque. Donner les vecteurs propres de ff veut dire décrire les sepsep en donnant leur dimension ainsi qu'une base.

Proposition 4. Soit f∈L(E)f \in \mathcal{L}(E), et soient λ1,…,λN\lambda_1, \dots, \lambda_N des valeurs propres de ff. Alors les sepsep associés sont en somme directe.

Démonstration. (À savoir.) Par récurrence sur NN :

  • Initialisation. Si N=2N = 2. Supposons que x1∈sep(λ1)x_1 \in sep(\lambda_1) et x2∈sep(λ2)x_2 \in sep(\lambda_2) sont tels que x1+x2=0x_1 + x_2 = 0. On applique ff : f(x1+x2)=f(x1)+f(x2)=0=λ1x1+λ2x2f(x_1 + x_2) = f(x_1) + f(x_2) = 0 = \lambda_1 x_1 + \lambda_2 x_2, donc {x1+x2=0λ1(x1+x2)=0\begin{cases} x_1 + x_2 = 0 \\ \lambda_1 (x_1 + x_2) = 0 \end{cases}, puis (λ2−λ1)x2=0(\lambda_2 - \lambda_1) x_2 = 0. Or λ2≠λ1\lambda_2 \neq \lambda_1, donc x2=0x_2 = 0, puis x1=−x2=0x_1 = -x_2 = 0.
  • Hérédité. Supposons la propriété vraie au rang NN. Montrons qu'elle est vraie au rang N+1N + 1 : soient (xi)i=1,…,N+1(x_i)_{i=1,\dots,N+1} avec xi∈sep(λi)x_i \in sep(\lambda_i) tels que ∑i=1N+1xi=0\displaystyle\sum_{i=1}^{N+1} x_i = 0. Alors, de même, f(∑xi)=∑1N+1λixi=0f(\sum x_i) = \displaystyle\sum_{1}^{N+1} \lambda_i x_i = 0. On multiplie ∑1N+1xi\displaystyle\sum_{1}^{N+1} x_i par λN+1\lambda_{N+1}, et on le soustrait à ∑1N+1λixi=0\displaystyle\sum_{1}^{N+1} \lambda_i x_i = 0. On obtient ∑i=1N(λN+1−λi)xi=0\displaystyle\sum_{i=1}^{N} (\lambda_{N+1} - \lambda_i) x_i = 0. Or les sep(λi)sep(\lambda_i), pour i=1,…,Ni = 1, \dots, N, forment NN sepsep. Par hypothèse de récurrence, ils sont en somme directe. Donc ∑i=1N(λN+1−λi)xi=0⇒(λN+1−λi)xi=0,∀i≤N\displaystyle\sum_{i=1}^{N} (\lambda_{N+1} - \lambda_i) x_i = 0 \Rightarrow (\lambda_{N+1} - \lambda_i) x_i = 0, \forall i \leq N. La propriété est montrée au rang N+1N + 1.

□\square

Exemple (hors programme). Soit E=C∞(R)E = C^\infty(\mathbb{R}). Soit D:{E→Ef↦f′D : \begin{cases} E \to E \\ f \mapsto f' \end{cases}. Alors D(f+λg)=(f+λg)′=f′+λg′D(f + \lambda g) = (f + \lambda g)' = f' + \lambda g'. Donc D∈L(E)D \in \mathcal{L}(E). De plus λ∈sp(D)⇔∃f≠0,D(f)=λf⇔f′=λf⇔f(x)=Ceλx\lambda \in sp(D) \Leftrightarrow \exists f \neq 0, D(f) = \lambda f \Leftrightarrow f' = \lambda f \Leftrightarrow f(x) = C e^{\lambda x}, où C∈RC \in \mathbb{R}. Ainsi, tout λ∈R\lambda \in \mathbb{R} est vp et sp(D)=Rsp(D) = \mathbb{R}. Ici, dim⁡(E)=+∞\dim(E) = +\infty ; en appliquant P4, on retrouve l'exercice 2 de la feuille 1 (les vp forment une famille libre puisque les sepsep sont en somme directe).

Polynôme caractéristique​

Définition - Proposition 1. Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}). Alors la fonction {K→KX↦det(A−XIn)\begin{cases} \mathbb{K} \to \mathbb{K} \\ X \mapsto det(A - X I_n) \end{cases} est un polynôme, appelé "polynôme caractéristique". On le note souvent χA\chi_A. De même, si f∈L(E)f \in \mathcal{L}(E), on définit χf(X)=det(f−XID)\chi_f(X) = det(f - X I_D), le polynôme caractéristique de ff.

Remarque. Si AA est une matrice de ff, alors par définition du déterminant on a χf=χA\chi_f = \chi_A. On peut le vérifier : en effet si A∼BA \sim B, alors AA et BB représentent le même endomorphisme, et B=P−1APB = P^{-1} A P, où P∈GLn(K)P \in GL_n(\mathbb{K}). Alors χB(X)=det(B−XIn)=det(P−1AP−XP−1P)=det(P−1(A−XIn)P)=det(P−1)det(A−XIn)det(P)=det(A−XIn)=χA(X)\chi_B(X) = det(B - X I_n) = det(P^{-1} A P - X P^{-1} P) = det(P^{-1} (A - X I_n) P) = det(P^{-1}) det(A - X I_n) det(P) = det(A - X I_n) = \chi_A(X). Ainsi "2 matrices semblables ont le même polynôme caractéristique". La réciproque est fausse : soient A1=(0100)A_1 = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} et A2=(0000)A_2 = \begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}. Alors χA1(X)=det(A1−XI2)=∣−X10−X∣=X2\chi_{A_1}(X) = det(A_1 - X I_2) = \begin{vmatrix} -X & 1 \\ 0 & -X \end{vmatrix} = X^2 et χA2=∣−X00−X∣=X2\chi_{A_2} = \begin{vmatrix} -X & 0 \\ 0 & -X \end{vmatrix} = X^2. Donc χA1=χA2\chi_{A_1} = \chi_{A_2}. Mais A1A_1 n'est pas semblable à A2A_2.

Exemple. Soit hλ:{E→Ex↦λxh_\lambda : \begin{cases} E \to E \\ x \mapsto \lambda x \end{cases}. Alors dans n'importe quelle base BB, matB(hλ)=(λ0⋱0λ)mat_B(h_\lambda) = \begin{pmatrix} \lambda & & 0 \\ & \ddots & \\ 0 & & \lambda \end{pmatrix}. En effet, si B=(e1,…,en)B = (e_1, \dots, e_n), hλ(e1)=λe1=λe1+0.e2+⋯+0enh_\lambda(e_1) = \lambda e_1 = \lambda e_1 + 0.e_2 + \dots + 0 e_n. Donc la première colonne est (λ0⋮0)\begin{pmatrix} \lambda \\ 0 \\ \vdots \\ 0 \end{pmatrix}, et pareil pour les autres colonnes ; ainsi : χhλ(X)=det(hλ−XIn)=det(λIn−XIn)=det((λ−X)In)=(λ−X)ndet(In)=(λ−X)n\chi_{h_\lambda}(X) = det(h_\lambda - X I_n) = det(\lambda I_n - X I_n) = det((\lambda - X) I_n) = (\lambda - X)^n det(I_n) = (\lambda - X)^n.

Remarque. Dans DP1, on a admis que χf\chi_f est un polynôme. Cela se voit en développant : ∣a11−Xa12a1n⋮⋱⋮an1…ann−X∣\begin{vmatrix} a_{11} - X & a_{12} & a_{1n} \\ \vdots & \ddots & \vdots \\ a_{n1} & \dots & a_{nn} - X \end{vmatrix}.

Proposition 2. Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}). Alors χA\chi_A est de la forme :

χA(X)=(−1)n(Xn−tr(A)Xn−1+…⏟coefficients pas gentils+(−1)ndet(A))\chi_A(X) = (-1)^n \Big(X^n - tr(A) X^{n-1} + \underbrace{\dots}_{\text{coefficients pas gentils}} + (-1)^n det(A)\Big)

Démonstration. (Partielle.) On développe ∣a11−Xa21…a1n⋮⋮⋮⋮an1……ann−X∣\begin{vmatrix} a_{11} - X & a_{21} & \dots & a_{1n} \\ \vdots & \vdots & \vdots & \vdots \\ a_{n1} & \dots & \dots & a_{nn} - X \end{vmatrix} par rapport à la première colonne :

(a11−X)∣a22−X…a2n⋮⋱⋮an2…ann−X∣+∑i≥2ai1cofacteur(i,1)⏟deˊterminant de matrices du meˆme type de taille n−1=⋯=(a11−X)…(ann−X)+termes de degreˊ≤n−1(a_{11} - X) \begin{vmatrix} a_{22} - X & \dots & a_{2n} \\ \vdots & \ddots & \vdots \\ a_{n2} & \dots & a_{nn} - X \end{vmatrix} + \sum_{i \geq 2} a_{i1} \underbrace{\text{cofacteur}(i, 1)}_{\text{déterminant de matrices du même type de taille } n - 1} = \dots = (a_{11} - X) \dots (a_{nn} - X) + \text{termes de degré} \leq n - 1

=(−X)n+(−1)n−1(∑i≥1aii)Xn−1+⋯+= (-X)^n + (-1)^{n-1} \Big(\displaystyle\sum_{i \geq 1} a_{ii}\Big) X^{n-1} + \dots + termes de degré (dur...). De manière plus accessible : χA(0)=det(A−0.In)=det(A)\chi_A(0) = det(A - 0.I_n) = det(A). Ainsi, le terme constant de χA\chi_A est det(A)det(A). □\square

Démonstration. (À savoir.)

  • Si n=2n = 2, on obtient avec A=(abcd)A = \begin{pmatrix} a & b \\ c & d \end{pmatrix} que

    χA(X)=X2−(tr(A))X+det(A)=X2−(a+d)X+(ad−bc)\chi_A(X) = X^2 - (tr(A)) X + det(A) = X^2 - (a + d) X + (ad - bc)
  • deg(χA)=ndeg(\chi_A) = n, et χA(0)=det(A)\chi_A(0) = det(A).

□\square

Exercice. Retrouver par le calcul que, si n=2n = 2, χA(X)=X2−tr(A)X+det(A)\chi_A(X) = X^2 - tr(A) X + det(A). On a χA(X)=∣a−Xbcd−X∣=…\chi_A(X) = \begin{vmatrix} a - X & b \\ c & d - X \end{vmatrix} = \dots

Théorème III.3. Soit f∈L(E)f \in \mathcal{L}(E). Alors le spectre de ff est constitué de l'ensemble des racines de χf\chi_f. Autrement dit, λ∈sp(f)⇔χf(λ)=0\lambda \in sp(f) \Leftrightarrow \chi_f(\lambda) = 0. Ou encore sp(f)=χf−1({0})sp(f) = \chi_f^{-1}(\{0\}).

Démonstration. Soit λ∈sp(f)\lambda \in sp(f). Ceci équivaut à f−λIdf - \lambda Id non injectif ⇔det(f−λId)=0⇔χf(λ)=0\Leftrightarrow det(f - \lambda I_d) = 0 \Leftrightarrow \chi_f(\lambda) = 0. □\square

Exercice. Calculer les valeurs propres et les vecteurs propres de A=(1221)A = \begin{pmatrix} 1 & 2 \\ 2 & 1 \end{pmatrix}.

On utilise le T3 : on calcule χA\chi_A et on cherche ses racines : χA(X)=∣1−X221−X∣=(1−X)2−4=(X+1)(X−3)\chi_A(X) = \begin{vmatrix} 1 - X & 2 \\ 2 & 1 - X \end{vmatrix} = (1 - X)^2 - 4 = (X + 1)(X - 3). Donc les racines de χA\chi_A sont −1-1 et 3. Ainsi, puisque les valeurs propres sont les racines de χA\chi_A, on a que sp(A)={−1,3}sp(A) = \{-1, 3\}. On cherche les vecteurs propres associés, c'est-à-dire qu'on en déduit les sepsep :

  • pour λ=−1\lambda = -1, on résout AX=−X⇔{x+2y=−x2x+y=−y⇔x+y=0AX = -X \Leftrightarrow \begin{cases} x + 2y = -x \\ 2x + y = -y \end{cases} \Leftrightarrow x + y = 0. Ainsi, (x,y)∈Ker(A+ID)⇔x+y=0(x, y) \in \text{Ker}(A + I_D) \Leftrightarrow x + y = 0. Autrement dit : sep(−1)=Ker(A+Id)={(x,y)∈R,x=−y}sep(-1) = \text{Ker}(A + I_d) = \{(x, y) \in \mathbb{R}, x = -y\}. C'est une droite, engendrée par (1,−1)(1, -1) ;
  • pour λ=3\lambda = 3, on résout AX=3XAX = 3X, donc on résout {x+2y=3x2x+y=3y⇔x=y\begin{cases} x + 2y = 3x \\ 2x + y = 3y \end{cases} \Leftrightarrow x = y. Ainsi, Ker(A−3Id)=Vect(1,1)\text{Ker}(A - 3 I_d) = \text{Vect}(1, 1) : c'est une droite.

Corollaire 4. Soit EE un espace vectoriel de dimension nn et soit u∈L(E)u \in \mathcal{L}(E). Alors uu a au plus nn valeurs propres distinctes.

Démonstration. Puisque les valeurs propres sont les racines de χu\chi_u et que deg(χu)=ndeg(\chi_u) = n, on déduit le corollaire. □\square

En particulier card(sp(u))<+∞\text{card}(sp(u)) < +\infty, pareil pour A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}).

Définition 5. Soit λ∈sp(u)\lambda \in sp(u). On appelle "multiplicité" de λ\lambda, notée parfois m(λ)m(\lambda) ou mλm_\lambda, l'ordre de λ\lambda en tant que racine du polynôme χu\chi_u. Ainsi mλm_\lambda est le plus grand nombre tel que χu(X)=(X−λ)mλQ(X)\chi_u(X) = (X - \lambda)^{m_\lambda} Q(X). Si mλ=1m_\lambda = 1, alors on dit que λ\lambda est racine simple.

Proposition 6. Soit u∈L(E)u \in \mathcal{L}(E) et λ∈sp(u)\lambda \in sp(u). Alors on a

1≤dim⁡(Ker(u−λId))≤mλ1 \leq \dim(\text{Ker}(u - \lambda I_d)) \leq m_\lambda

Démonstration. (À savoir.) Puisque λ∈sp(u)\lambda \in sp(u), on a Ker(u−λId)≠{0}\text{Ker}(u - \lambda I_d) \neq \{0\}. Donc sa dimension est au moins 1. On note d0=dim⁡(Ker(u−λId))d_0 = \dim(\text{Ker}(u - \lambda I_d)). Soit (e1,…,ed0)(e_1, \dots, e_{d_0}) une base de ce sev. On la complète en une base BB de EE. Dans cette base, la matrice de uu est :

(λ0……bb0λ0…bb⋮⋮⋱⋮cc0000cc)=(λId0B0C)\begin{pmatrix} \lambda & 0 & \dots & \dots & b & b \\ 0 & \lambda & 0 & \dots & b & b \\ \vdots & \vdots & \ddots & \vdots & c & c \\ 0 & 0 & 0 & 0 & c & c \end{pmatrix} = \begin{pmatrix} \lambda I_{d_0} & B \\ 0 & C \end{pmatrix}

On a alors χu(X)=∣(λ−X)Id0B0C−XIn−d0∣=(λ−X)d0Q(X)\chi_u(X) = \begin{vmatrix} (\lambda - X) I_{d_0} & B \\ 0 & C - X I_{n - d_0} \end{vmatrix} = (\lambda - X)^{d_0} Q(X). Or χu=(X−λ)mλQ~(X)\chi_u = (X - \lambda)^{m_\lambda} \tilde{Q}(X), où mλm_\lambda est le plus grand nombre possible pour cette écriture. Ainsi on a d0≤mλd_0 \leq m_\lambda. □\square

Corollaire 7. Si λ\lambda est simple, alors dim⁡(Ker(u−λId))=1\dim(\text{Ker}(u - \lambda I_d)) = 1.

Démonstration. Par P6 on a 1≤d0≤11 \leq d_0 \leq 1. Donc d0=1d_0 = 1. □\square

Exemple. Pour la matrice A=(1221)A = \begin{pmatrix} 1 & 2 \\ 2 & 1 \end{pmatrix}, on avait m(−1)=dim⁡(Ker(A+I2))=1m(-1) = \dim(\text{Ker}(A + I_2)) = 1 et m(3)=dim⁡(Ker(A−3I2))=1m(3) = \dim(\text{Ker}(A - 3 I_2)) = 1 ; ainsi les 2 valeurs propres sont simples.

Diagonalisation​

Définition 1. Soit f∈L(E)f \in \mathcal{L}(E). On dit que ff est diagonalisable s'il existe une base BB de EE telle que matB(f)mat_B(f) est diagonale. Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}). On dit que AA est diagonalisable si AA est semblable à une matrice diagonale, c'est-à-dire ∃P∈GLn(K)\exists P \in GL_n(\mathbb{K}) et DD diagonale telles que A=PDP−1A = P D P^{-1}. Diagonaliser AA veut donc dire "trouver DD et PP".

Proposition 2. Les points suivants sont équivalents :

  1. ff est diagonalisable ;
  2. il existe une base de EE formée de vp⃗\vec{vp} de ff ;
  3. la somme des sepsep de ff vaut EE ;
  4. la somme des dimensions des sepsep vaut dim⁡(E)\dim(E).

Démonstration. On montre que 1 est équivalent à 2. Soit BB une base telle que matB(f)mat_B(f) est diagonale. On note B=(e1,…,en)B = (e_1, \dots, e_n) et (λ10⋱0λn)\begin{pmatrix} \lambda_1 & & 0 \\ & \ddots & \\ 0 & & \lambda_n \end{pmatrix} la matrice. Par lecture, on a f(ei)=λieif(e_i) = \lambda_i e_i. Donc eie_i est un vp⃗\vec{vp} de ff. Réciproquement : si (e1,…,en)(e_1, \dots, e_n) est une base de vp⃗\vec{vp}, associés à des vpvp λ1,…,λn\lambda_1, \dots, \lambda_n, dans cette base la matrice de ff est (λ10⋱0λn)\begin{pmatrix} \lambda_1 & & 0 \\ & \ddots & \\ 0 & & \lambda_n \end{pmatrix}. On montre que 2 implique 3. Soit (e1,…,en)(e_1, \dots, e_n) une base de vp⃗\vec{vp}. Alors chaque eie_i appartient à ⨁λ∈sp(f)Ker(f−λId)\displaystyle\bigoplus_{\lambda \in sp(f)} \text{Ker}(f - \lambda I_d). Donc Vect(e1,…,en)⊂⨁Ker(f−λId)\text{Vect}(e_1, \dots, e_n) \subset \bigoplus \text{Ker}(f - \lambda I_d). Donc E=⨁Ker(f−λId)E = \bigoplus \text{Ker}(f - \lambda I_d). On va montrer que 3 implique 4. On a vu que la somme des sepsep est directe. Ainsi dim⁡(E)=dim⁡(⨁λ∈sp(f)Ker(f−λId))=∑λ∈sp(f)dim⁡(Ker(f−λId))\dim(E) = \dim\Big(\displaystyle\bigoplus_{\lambda \in sp(f)} \text{Ker}(f - \lambda I_d)\Big) = \sum_{\lambda \in sp(f)} \dim(\text{Ker}(f - \lambda I_d)). On va montrer que 4 implique 2. Supposons que ∑λ∈sp(f)dim⁡(Ker(f−λId))=dim⁡(E)\displaystyle\sum_{\lambda \in sp(f)} \dim(\text{Ker}(f - \lambda I_d)) = \dim(E). Puisque la somme est directe, si BiB_i est une base de Ker(f−λi)\text{Ker}(f - \lambda_i), pour λi∈sp(f)\lambda_i \in sp(f), on a ⋃Bi\bigcup B_i est une base de EE. Donc 2. □\square

Exercice. Pour notre matrice AA, elle est diagonalisable car ∑λ∈sp(A)dim⁡(Ker(A−λI2))=1+1=dim⁡(R2)\displaystyle\sum_{\lambda \in sp(A)} \dim(\text{Ker}(A - \lambda I_2)) = 1 + 1 = \dim(\mathbb{R}^2).

Théorème III.3. (Condition nécessaire et suffisante de diagonalisation) ff est diagonalisable si et seulement si

{χf est scindeˊ sur K∀λ∈sp(f),dim⁡(Ker(f−λId))=mλ\begin{cases} \chi_f \text{ est scindé sur } \mathbb{K} \\ \forall \lambda \in sp(f), \dim(\text{Ker}(f - \lambda I_d)) = m_\lambda \end{cases}

Rappel. On dit que PP est scindé sur K\mathbb{K} s'il est de la forme ∏(X−λi)mi\prod (X - \lambda_i)^{m_i}.

Exemple. (X−1)X(X+2)(X - 1) X (X + 2) est scindé (à racines simples), X2+1X^2 + 1 n'est pas scindé sur R\mathbb{R}.

Démonstration. Supposons ff diagonalisable. Soit (λ10⋱0λn)\begin{pmatrix} \lambda_1 & & 0 \\ & \ddots & \\ 0 & & \lambda_n \end{pmatrix} sa matrice de diagonalisation. Alors χf(X)=∣λ1−X0⋱0λn−X∣=∏i=1n(λi−X)\chi_f(X) = \begin{vmatrix} \lambda_1 - X & & 0 \\ & \ddots & \\ 0 & & \lambda_n - X \end{vmatrix} = \displaystyle\prod_{i=1}^{n} (\lambda_i - X). De plus, par P2, n=∑λ∈sp(f)dim⁡(Ker(f−λId))≤∑λ∈sp(f)m(λ)=nn = \displaystyle\sum_{\lambda \in sp(f)} \dim(\text{Ker}(f - \lambda I_d)) \leq \sum_{\lambda \in sp(f)} m(\lambda) = n. ∑(m(λ)−dim⁡(Ker(f−λId)))=0\sum (m(\lambda) - \dim(\text{Ker}(f - \lambda I_d))) = 0. Or par II-P6, ces termes sont supérieurs à 0. Donc ils sont nuls, et dim⁡(Ker(f−λId))=m(λ)\dim(\text{Ker}(f - \lambda I_d)) = m(\lambda). □\square

Réciproque. Puisque ∑λ∈sp(f)dim⁡(Ker(f−λId))=∑λ∈sp(f)m(λ)=deg(χf)=n\displaystyle\sum_{\lambda \in sp(f)} \dim(\text{Ker}(f - \lambda I_d)) = \sum_{\lambda \in sp(f)} m(\lambda) = deg(\chi_f) = n.

Corollaire 4. Si ff a nn valeurs propres distinctes (où n=dim⁡(E)n = \dim(E)), alors ces vpvp sont simples et ff est diagonalisable.

Démonstration. ∑λ∈sp(f)dim⁡(sep(λ))=1+⋯+1⏟n fois=dim⁡(E)\displaystyle\sum_{\lambda \in sp(f)} \dim(sep(\lambda)) = \underbrace{1 + \dots + 1}_{n \text{ fois}} = \dim(E). Donc on applique P2. □\square

Exercice. Soit TT une matrice triangulaire, avec des coefficients diagonaux distincts 2 à 2. Alors elle est diagonalisable.

Solution. Supposons que T=(λ1……⋱0…λn)T = \begin{pmatrix} \lambda_1 & \dots & \dots \\ & \ddots & \\ 0 & \dots & \lambda_n \end{pmatrix}. On a χT(X)=det(T−XIn)=∣λ1−X…λ1n⋱0…λn−X∣=∏i=1n(λi−X)\chi_T(X) = det(T - X I_n) = \begin{vmatrix} \lambda_1 - X & \dots & \lambda_{1n} \\ & \ddots & \\ 0 & \dots & \lambda_n - X \end{vmatrix} = \displaystyle\prod_{i=1}^{n} (\lambda_i - X). Donc les éléments diagonaux sont des vpvp.

Exercice. Soit A=(0100)A = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} ; est-elle diagonalisable ? Alors χA(X)=∣−X10−X∣=X2\chi_A(X) = \begin{vmatrix} -X & 1 \\ 0 & -X \end{vmatrix} = X^2. Donc AA n'a qu'une seule vp : 0. On a 2 méthodes de réponse :

  1. On calcule sep(0)sep(0) : on résout AX=0AX = 0, on a donc {y=00=0\begin{cases} y = 0 \\ 0 = 0 \end{cases}. Donc dim⁡(sep(0))=1\dim(sep(0)) = 1. Or 0 est de multiplicité 2. Ainsi AA est non diagonalisable.
  2. Puisque 0 est la seule vp, si AA était diagonalisable elle serait semblable à (0000)\begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}, c'est-à-dire A=POP−1=0A = P O P^{-1} = 0. Absurde car A≠0A \neq 0.

Polynôme d'endomorphisme​

Soit P∈Rn[X]P \in \mathbb{R}_n[X] un polynôme. On le note P(X)=∑k=0NakXkP(X) = \displaystyle\sum_{k=0}^{N} a_k X^k.

  • Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}). On pose P(A)=∑k=0NakAkP(A) = \displaystyle\sum_{k=0}^{N} a_k A^k. Cela définit un élément de Mn(K)\mathcal{M}_n(\mathbb{K}). On parle d'un polynôme de matrice.
  • Soit EE un espace vectoriel de dimension nn, et f∈L(E)f \in \mathcal{L}(E). On pose P(f)=∑k=0NakfkP(f) = \displaystyle\sum_{k=0}^{N} a_k f^k, avec la convention f0=Idf^0 = I_d. On parle d'un polynôme d'endomorphisme. On rappelle que fk=f∘⋯∘f⏟k foisf^k = \underbrace{f \circ \dots \circ f}_{k \text{ fois}}.

Remarque. Si AA est une matrice de ff dans une base BB, alors P(A)P(A) est une matrice de P(f)P(f) dans la base BB.

Exemple. Si P(X)=X−1P(X) = X - 1, alors P(f)=f−IdP(f) = f - I_d.

Proposition 2. On dit que Φ:{P↦P(f)K[X]→L(E)\Phi : \begin{cases} P \mapsto P(f) \\ \mathbb{K}[X] \to \mathcal{L}(E) \end{cases} est un morphisme d'algèbre unitaire. Cela veut dire :

  • (P+αQ)(f)=P(f)+αQ(f)(P + \alpha Q)(f) = P(f) + \alpha Q(f)
  • (PQ)(f)=P(f)∘Q(f)(PQ)(f) = P(f) \circ Q(f)
  • 1(f)=Id1(f) = Id

Démonstration. Évident ou admis. □\square

Rappel. Si f∈L(E)f \in \mathcal{L}(E) et g∈L(E)g \in \mathcal{L}(E), alors en général, f∘g≠g∘ff \circ g \neq g \circ f. Si AA et BB sont dans Mn(K)\mathcal{M}_n(\mathbb{K}), alors en général AB≠BAAB \neq BA.

Vocabulaire. On dit que ff et gg commutent lorsque f∘g=g∘ff \circ g = g \circ f. Pareil pour les 2 matrices.

Proposition 3. Si ff et gg commutent, alors tout polynôme en ff commute avec tout polynôme en gg.

Démonstration. Soit PP et QQ deux polynômes, P(X)=∑0NakXkP(X) = \displaystyle\sum_{0}^{N} a_k X^k et Q(X)=∑0MbkXkQ(X) = \displaystyle\sum_{0}^{M} b_k X^k. On part de f∘g=g∘ff \circ g = g \circ f. On compose par ff à gauche, donc f2∘g=f∘f∘g=f∘g∘f=g∘f∘f=g∘f2f^2 \circ g = f \circ f \circ g = f \circ g \circ f = g \circ f \circ f = g \circ f^2. De même, on prouve par récurrence que fk∘g=g∘fkf^k \circ g = g \circ f^k : on suppose cela vrai au rang kk : fk∘g=g∘fk⇒f∘fk∘g=f∘g∘fk=g∘fk+1f^k \circ g = g \circ f^k \Rightarrow f \circ f^k \circ g = f \circ g \circ f^k = g \circ f^{k+1}. La propriété est donc vraie. On a P(f)∘g=(∑0Nakfk)∘g=∑0Nak(fk∘g)=∑0Nakg∘fk=g∘(∑0Nakfk)=g∘P(f)P(f) \circ g = \Big(\displaystyle\sum_{0}^{N} a_k f^k\Big) \circ g = \sum_{0}^{N} a_k (f^k \circ g) = \sum_{0}^{N} a_k g \circ f^k = g \circ \Big(\sum_{0}^{N} a_k f^k\Big) = g \circ P(f). On a, en inversant les rôles, Q(g)∘P(f)=P(f)∘Q(g)Q(g) \circ P(f) = P(f) \circ Q(g). □\square

Proposition 3 bis. Pareil pour les matrices.

Définition 4. Soit f∈L(E)f \in \mathcal{L}(E). On dit que ff annule PP lorsque P(f)=0P(f) = 0. On dit aussi que PP est annulateur de ff.

Exemple. Soit ss une symétrie. Soit P(X)=X2−1P(X) = X^2 - 1. Alors P(s)=s2−Id=0P(s) = s^2 - I_d = 0. Donc ss annule PP. Soit pp un projecteur. Soit Q(X)=X2−XQ(X) = X^2 - X. Alors Q(p)=p2−p=p∘p−p=0Q(p) = p^2 - p = p \circ p - p = 0. Donc pp annule QQ.

Proposition 5. Soit λ∈sp(f)\lambda \in sp(f) et xx un vp⃗\vec{vp} associé. Alors P(f)x=P(λ)xP(f) x = P(\lambda) x.

Démonstration. (À savoir.) On a f(x)=λxf(x) = \lambda x. Donc f(f(x))=f(λx)⇔f∘f(x)=λf(x)=λλx=λ2xf(f(x)) = f(\lambda x) \Leftrightarrow f \circ f(x) = \lambda f(x) = \lambda \lambda x = \lambda^2 x. On itère : ∀k∈N,fk(x)=λkx\forall k \in \mathbb{N}, f^k(x) = \lambda^k x, puis par combinaison linéaire : (∑0Nakfk)(x)=∑0N(akλkx)=p(λ)x\Big(\displaystyle\sum_{0}^{N} a_k f^k\Big)(x) = \sum_{0}^{N} (a_k \lambda^k x) = p(\lambda) x. □\square

Corollaire 6. Supposons que ff annule PP : alors sp(f)⊂P−1({0})sp(f) \subset P^{-1}(\{0\}) (c'est-à-dire que les valeurs propres de ff sont contenues dans les racines de PP).

Démonstration. Si P(f)=0P(f) = 0, par P5, p(λ)x=0p(\lambda) x = 0 lorsque λ∈sp(f)\lambda \in sp(f). Or x⃗\vec{x} est un vp⃗⇒x≠0\vec{vp} \Rightarrow x \neq 0. Donc p(λ)=0p(\lambda) = 0. □\square

Exemple. Si ss est une symétrie, quelles sont les valeurs propres possibles ? On a que P(X)=X2−1P(X) = X^2 - 1 est annulateur de ss. Donc sp(f)⊂P−1({0})sp(f) \subset P^{-1}(\{0\}). Or P(X)=0⇔X2−1=0⇔X=1P(X) = 0 \Leftrightarrow X^2 - 1 = 0 \Leftrightarrow X = 1 ou X=−1X = -1. Donc les vpvp d'une symétrie peuvent être 1 ou −1-1.

Exemple. Si pp est un projecteur, sp(p)⊂{0,1}sp(p) \subset \{0, 1\} (à faire).

Remarque. Le C6 donne un critère pour chercher les vpvp de ff. Mais C6 n'est qu'une inclusion ; par exemple si f=Idf = I_d et P(X)=X(X−1)P(X) = X(X - 1), alors sp(f)={1}sp(f) = \{1\}, et P(f)=Id(Id−Id)=0P(f) = I_d(I_d - I_d) = 0, et P−1({0})={0,1}P^{-1}(\{0\}) = \{0, 1\}.

Théorème III.7. Soit EE un espace vectoriel de dimension finie. Pour que ff soit diagonalisable, il faut et il suffit qu'il existe un polynôme annulateur scindé à racines simples.

Démonstration. Soit ff diagonalisable ; alors dans une base BB de vp⃗\vec{vp}, MatB(f)=(λ10⋱0λn)Mat_B(f) = \begin{pmatrix} \lambda_1 & & 0 \\ & \ddots & \\ 0 & & \lambda_n \end{pmatrix}. Ici les λi\lambda_i peuvent se répéter (ils ne sont pas différents deux à deux). Soient (μi)i=1,…,n(\mu_i)_{i=1,\dots,n} les vpvp (sans répétitions). On pose P(X)=∏i=1d(X−μi)P(X) = \displaystyle\prod_{i=1}^{d} (X - \mu_i). Alors PP est scindé à racines simples. Soit MM la matrice ci-dessus, alors

P(M)=(0⋱μ2−μ1⋱μd−μ1)(μ1−μ20μ3−μ2⋱μd−μ2)⋯(μ1−μd⋱0⋱0)=0P(M) = \begin{pmatrix} 0 & & & & \\ & \ddots & & & \\ & & \mu_2 - \mu_1 & & \\ & & & \ddots & \\ & & & & \mu_d - \mu_1 \end{pmatrix} \begin{pmatrix} \mu_1 - \mu_2 & & & & \\ & 0 & & & \\ & & \mu_3 - \mu_2 & & \\ & & & \ddots & \\ & & & & \mu_d - \mu_2 \end{pmatrix} \cdots \begin{pmatrix} \mu_1 - \mu_d & & & & \\ & \ddots & & & \\ & & 0 & & \\ & & & \ddots & \\ & & & & 0 \end{pmatrix} = 0

La réciproque est admise et repose sur le lemme des noyaux. □\square

Exercice. Soit A=(01…0⋮⋱⋱⋮0⋮⋱110…0)A = \begin{pmatrix} 0 & 1 & \dots & 0 \\ \vdots & \ddots & \ddots & \vdots \\ 0 & \vdots & \ddots & 1 \\ 1 & 0 & \dots & 0 \end{pmatrix}. Alors An=?A^n = ? Donc... on en déduit les vpvp de AA.

Application. Projecteurs et symétries sont diagonalisables.

Exemple des symétries. Soit ss tel que s2=Ids^2 = I_d. Soit P(X)=X2−1P(X) = X^2 - 1. Alors PP est annulateur de ss. Or P(X)=(X−1)(X+1)P(X) = (X - 1)(X + 1) est scindé à racines simples. Donc {1 et −1 sont les vp possibles de ss est diagonalisable\begin{cases} 1 \text{ et } -1 \text{ sont les } vp \text{ possibles de } s \\ s \text{ est diagonalisable} \end{cases}. On a donc E=⨁λ∈sp(s)sep(λ)=Ker(s−Id)⨁Ker(s+Id)E = \displaystyle\bigoplus_{\lambda \in sp(s)} sep(\lambda) = \text{Ker}(s - I_d) \bigoplus \text{Ker}(s + I_d), si on suppose que 1 et −1-1 sont valeurs propres.

Remarque. Si sp(s)={1}sp(s) = \{1\}, s=Ids = I_d (car ss a pour matrice InI_n). Si sp(s)={−1}sp(s) = \{-1\}, s=−Ids = -I_d (car ss a pour matrice −In-I_n). Si 1 et −1-1 sont vpvp, alors : Ker(s−Id)={x∈E,sx=x}\text{Ker}(s - Id) = \{x \in E, sx = x\} est le plan de symétrie et Ker(s+Id)={x∈E,sx=−x}\text{Ker}(s + I_d) = \{x \in E, sx = -x\} est la direction de la symétrie.

Théorème III.8. (Théorème de Cayley-Hamilton) Soit f∈L(E)f \in \mathcal{L}(E) et χf\chi_f son polynôme caractéristique. Alors χf(f)=0\chi_f(f) = 0.

Démonstration. Ce n'est pas χf(f)=det(f−fId)=det(0)=0\chi_f(f) = det(f - f I_d) = det(0) = 0, car χf(X):=det(f−XId)\chi_f(X) := det(f - X I_d). On ne peut pas remplacer XX par ff ici. □\square

Une application des polynômes annulateurs. Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}) et PP annulateur de AA. Alors on peut calculer l'inverse de AA via PP :

  • si P=χAP = \chi_A et si P(0)=0P(0) = 0, alors AA n'est pas inversible, car P(0)=det(A)P(0) = det(A) ;
  • si PP est un polynôme tel que P(0)≠0P(0) \neq 0, alors AA est inversible. En effet si P(X)=∑k=0dakXkP(X) = \displaystyle\sum_{k=0}^{d} a_k X^k, alors P(A)=0⇔a0Id+a1A+⋯+adAd=0⇔−1a0(a1A+⋯+adAd)=Id⇔A[−1a0(a1Id+⋯+adAd−1)]=IdP(A) = 0 \Leftrightarrow a_0 I_d + a_1 A + \dots + a_d A^d = 0 \Leftrightarrow -\dfrac{1}{a_0}(a_1 A + \dots + a_d A^d) = I_d \Leftrightarrow A\Big[-\dfrac{1}{a_0}(a_1 I_d + \dots + a_d A^{d-1})\Big] = I_d. Donc A∈GLn(K)A \in GL_n(\mathbb{K}), et A−1=−1a0(a1Id+⋯+adAd−1)A^{-1} = -\dfrac{1}{a_0}(a_1 I_d + \dots + a_d A^{d-1}). Ceci permet d'exprimer A−1A^{-1} en fonction de puissances de AA. C'est pratique si on connaît χA\chi_A.

Trigonalisation​

Définition. Un endomorphisme est trigonalisable s'il existe une base dans laquelle sa matrice est triangulaire. Une matrice est trigonalisable si elle est semblable à une matrice triangulaire : AA trigonalisable ⇔A=PTP−1\Leftrightarrow A = P T P^{-1}, où P∈GLn(K)P \in GL_n(\mathbb{K}) et TT triangulaire.

Remarque. Une matrice diagonalisable est trigonalisable (réciproque fausse).

Théorème III.3. Une matrice (ou l'endomorphisme qu'elle représente) est trigonalisable si et seulement si son polynôme caractéristique est scindé.

Démonstration. Soit AA une matrice.

  • ⇒\Rightarrow : ∃P∈GLn\exists P \in GL_n et T=(t11⋱0tnn)T = \begin{pmatrix} t_{11} & & \\ & \ddots & \\ 0 & & t_{nn} \end{pmatrix} tels que A=PTP−1A = P T P^{-1}. Or 2 matrices semblables ont le même polynôme caractéristique. Ainsi χA(X)=det(T−XIn)=∣t11−X⋱0tnn−X∣=∏i=1n(tii−X)\chi_A(X) = det(T - X I_n) = \begin{vmatrix} t_{11} - X & & \\ & \ddots & \\ 0 & & t_{nn} - X \end{vmatrix} = \displaystyle\prod_{i=1}^{n} (t_{ii} - X) est scindé.
  • ⇐\Leftarrow : on le montre par récurrence sur nn.
    • Initialisation n=1n = 1. Les matrices sont des scalaires. Donc elles sont diagonales.
    • Hérédité. Supposons que la propriété est vraie au rang nn. Soit A∈Mn+1(K)A \in \mathcal{M}_{n+1}(\mathbb{K}) une matrice telle que χA\chi_A est scindé. Puisque χA\chi_A est scindé, il admet au moins une racine λ1\lambda_1. C'est donc une valeur propre de AA. Soit x1x_1 un vecteur propre associé. On complète ce vecteur en une base de Rn+1\mathbb{R}^{n+1}. Soit PP la matrice de passage depuis la base canonique vers cette base. Alors P−1AP=(λ1l…l0a1…a1⋮⋮⋮⋮0a1a1a1)P^{-1} A P = \begin{pmatrix} \lambda_1 & l & \dots & l \\ 0 & a_1 & \dots & a_1 \\ \vdots & \vdots & \vdots & \vdots \\ 0 & a_1 & a_1 & a_1 \end{pmatrix}, où L∈M1,nL \in \mathcal{M}_{1,n} et A1∈MnA_1 \in \mathcal{M}_n. Cette dernière matrice étant triangulaire par blocs, on a χA(X)=∣λ1−XL0A1−XIn∣=(λ1−X)χA1(X)\chi_A(X) = \begin{vmatrix} \lambda_1 - X & L \\ 0 & A_1 - X I_n \end{vmatrix} = (\lambda_1 - X) \chi_{A_1}(X). Or χA1\chi_{A_1} est scindé car χA\chi_A l'est. Donc A1A_1 est trigonalisable par hypothèse de récurrence : ∃T1\exists T_1 triangulaire, ∃Q∈GLn\exists Q \in GL_n, A1=QT1Q−1A_1 = Q T_1 Q^{-1}. Soit R=(10…00q…q⋮⋮⋮⋮0q…q)R = \begin{pmatrix} 1 & 0 & \dots & 0 \\ 0 & q & \dots & q \\ \vdots & \vdots & \vdots & \vdots \\ 0 & q & \dots & q \end{pmatrix}. Alors (λ1L0A1)=R(λ1LQ0T1)R−1\begin{pmatrix} \lambda_1 & L \\ 0 & A_1 \end{pmatrix} = R \begin{pmatrix} \lambda_1 & LQ \\ 0 & T_1 \end{pmatrix} R^{-1} (calcul à vérifier). Or la matrice (λ1LQ0T1)\begin{pmatrix} \lambda_1 & LQ \\ 0 & T_1 \end{pmatrix} est triangulaire. AA est semblable à cette matrice, elle est donc trigonalisable.

□\square

Corollaire 3. Toute matrice est trigonalisable sur C\mathbb{C}.

Démonstration. χA\chi_A est un polynôme scindé par le théorème de d'Alembert. On applique ensuite T2. □\square

Exemple. Les matrices nilpotentes. On dit que A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}) est nilpotente s'il existe k∈Nk \in \mathbb{N} tel que Ak=0A^k = 0. Par exemple : A=(0100)A = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}. Alors A2=(0000)A^2 = \begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}. Soit P(X)=X2P(X) = X^2. Alors P(A)=0P(A) = 0, donc PP est annulateur de AA. Donc sp(A)⊂P−1({0})={0}sp(A) \subset P^{-1}(\{0\}) = \{0\}. Donc la seule valeur propre de AA est 0. Si AA était diagonalisable, elle serait semblable à 0 (la matrice diagonale avec les valeurs propres sur la diagonale). Or AA semblable à 0 donne A=P−10P=0A = P^{-1} 0 P = 0. Or A≠0A \neq 0. Donc AA est non diagonalisable. On peut retrouver cela en analysant sep(0)sep(0) (déjà fait). Cette analyse est valable pour toute matrice nilpotente.

Applications aux suites récurrentes linéaires​

But. Que dire des suites suivantes : un+2=un+1+unu_{n+2} = u_{n+1} + u_n ? Et {un+1=un+3vnvn+1=2un−vn\begin{cases} u_{n+1} = u_n + 3 v_n \\ v_{n+1} = 2 u_n - v_n \end{cases} ?

Puissances des matrices​

Soit A∈Mn(K)A \in \mathcal{M}_n(\mathbb{K}). On veut calculer AkA^k. Voici 2 cas où on a une méthode :

  • Si AA est diagonalisable, A=PDP−1A = P D P^{-1}, où D=(λ10⋱0λn)D = \begin{pmatrix} \lambda_1 & & 0 \\ & \ddots & \\ 0 & & \lambda_n \end{pmatrix}. Alors on a : A2=PDP−1PDP−1=PD2P−1A^2 = P D P^{-1} P D P^{-1} = P D^2 P^{-1}. De même, par récurrence évidente, Ak=PDkP−1A^k = P D^k P^{-1}. De plus Dk=(λ1k0⋱0λnk)D^k = \begin{pmatrix} \lambda_1^k & & 0 \\ & \ddots & \\ 0 & & \lambda_n^k \end{pmatrix}. Ceci permet de calculer AkA^k.
  • Si AA est de la forme A=αId+NA = \alpha I_d + N, où NN est nilpotente. Soit pp tel que Np≠0N^p \neq 0 et Np+1=0N^{p+1} = 0. Alors Ak=(αIn+N)k=∑j=0k(kj)(αIn)k−jNj=∑j=0p(kj)αk−jNjA^k = (\alpha I_n + N)^k = \displaystyle\sum_{j=0}^{k} \binom{k}{j} (\alpha I_n)^{k-j} N^j = \sum_{j=0}^{p} \binom{k}{j} \alpha^{k-j} N^j.

Remarque. Si AA et BB sont deux matrices, (A+B)2=A2+B2+AB+BA(A + B)^2 = A^2 + B^2 + AB + BA : pas de binôme en général, sauf si AA et BB commutent.

Exemple. Si A=(210021002)A = \begin{pmatrix} 2 & 1 & 0 \\ 0 & 2 & 1 \\ 0 & 0 & 2 \end{pmatrix}. Notons que A=2I3+NA = 2 I_3 + N, où N=(010001000)N = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} ; alors N3=0N^3 = 0. Ainsi

Ak=(2I3+N)k=∑j=02(kj)αk−jNj=2kI3+(k1)2k−1N+(k2)2k−2N2=(2kk2k−1k(k−1)22k−202kk2k−1002k)A^k = (2 I_3 + N)^k = \sum_{j=0}^{2} \binom{k}{j} \alpha^{k-j} N^j = 2^k I_3 + \binom{k}{1} 2^{k-1} N + \binom{k}{2} 2^{k-2} N^2 = \begin{pmatrix} 2^k & k 2^{k-1} & \dfrac{k(k-1)}{2} 2^{k-2} \\ 0 & 2^k & k 2^{k-1} \\ 0 & 0 & 2^k \end{pmatrix}

Systèmes de suites récurrentes​

On considère pp suites v1,…,vpv^1, \dots, v^p où elles sont "couplées" par le système : {un+11=a11un1+⋯+a1punp⋮un+1p=ap1un1+⋯+appunp\begin{cases} u_{n+1}^1 = a_{11} u_n^1 + \dots + a_{1p} u_n^p \\ \vdots \\ u_{n+1}^p = a_{p1} u_n^1 + \dots + a_{pp} u_n^p \end{cases}. On pose Xn=(un1⋮unp)∈RpX_n = \begin{pmatrix} u_n^1 \\ \vdots \\ u_n^p \end{pmatrix} \in \mathbb{R}^p. Alors le système se réécrit Xn+1=AXnX_{n+1} = A X_n, où A=(a11…a1p⋮⋮ap1…app)A = \begin{pmatrix} a_{11} & \dots & a_{1p} \\ \vdots & & \vdots \\ a_{p1} & \dots & a_{pp} \end{pmatrix}. C'est une suite "géométrique" : on a X2=AX1=AAX0=A2X0X_2 = A X_1 = A A X_0 = A^2 X_0, puis par récurrence directe, Xn=AnX0X_n = A^n X_0. On peut avoir des formules pour les (uni)(u_n^i) si on sait calculer AkA^k.

Suites récurrentes d'ordre pp​

Problématique. On a une seule suite (un)n≥0(u_n)_{n \geq 0} qui vérifie :

un+p=∑i=0p−1aiun+iu_{n+p} = \sum_{i=0}^{p-1} a_i u_{n+i}

Exemple. un+2=un+1+unu_{n+2} = u_{n+1} + u_n ; ici p=2p = 2, a0=a1=1a_0 = a_1 = 1.

Résolution. On va transformer ce problème "linéaire scalaire et d'ordre pp" en un problème "vectoriel d'ordre 1". Pour cela, on pose Xn=(un⋮un+p−1)∈KpX_n = \begin{pmatrix} u_n \\ \vdots \\ u_{n+p-1} \end{pmatrix} \in \mathbb{K}^p. On a alors :

Xn+1=(un+1⋮un+p)=(010…0001…0⋱0……01a0………ap−1)(un⋮un+p−1)=AXnX_{n+1} = \begin{pmatrix} u_{n+1} \\ \vdots \\ u_{n+p} \end{pmatrix} = \begin{pmatrix} 0 & 1 & 0 & \dots & 0 \\ 0 & 0 & 1 & \dots & 0 \\ & & & \ddots & \\ 0 & \dots & \dots & 0 & 1 \\ a_0 & \dots & \dots & \dots & a_{p-1} \end{pmatrix} \begin{pmatrix} u_n \\ \vdots \\ u_{n+p-1} \end{pmatrix} = A X_n

On a à nouveau Xn=AnX0X_n = A^n X_0. La matrice AA a une forme particulière.

Définition - Proposition. Soit P(X)=Xp−ap−1Xp−1−⋯−a0P(X) = X^p - a_{p-1} X^{p-1} - \dots - a_0. Alors AA est la matrice "compagnon" de PP. On a ainsi χA(X)=(−1)pP(X)\chi_A(X) = (-1)^p P(X).

Démonstration. (Exercice classique.) On a χA(X)=det(A−XIp)=∣−X1…00⋱⋱00…⋱1a0…ap−2ap−1−X∣\chi_A(X) = det(A - X I_p) = \begin{vmatrix} -X & 1 & \dots & 0 \\ 0 & \ddots & \ddots & 0 \\ 0 & \dots & \ddots & 1 \\ a_0 & \dots & a_{p-2} & a_{p-1} - X \end{vmatrix}. On fait C1←C1+XC2C_1 \leftarrow C_1 + X C_2 : la première colonne devient

(0−X20⋮0a0+Xa1)\begin{pmatrix} 0 \\ -X^2 \\ 0 \\ \vdots \\ 0 \\ a_0 + X a_1 \end{pmatrix}

On ajoute ensuite X2C3X^2 C_3 à C1C_1 :

(00−X3⋮a0+a1X+a2X2)\begin{pmatrix} 0 \\ 0 \\ -X^3 \\ \vdots \\ a_0 + a_1 X + a_2 X^2 \end{pmatrix}

Ainsi de suite, on trouve en ajoutant XC2+⋯+Xp−1CpX C_2 + \dots + X^{p-1} C_p à C1C_1 :

χA(X)=∣010⋮−X⋱00⋱1a0+a1X+⋯+(ap−1−X)Xp−1……ap−1−X∣=(−1)p+2P(X)×∣10−X⋱⋱⋱0−X1∣\chi_A(X) = \begin{vmatrix} 0 & 1 & & 0 \\ \vdots & -X & \ddots & \\ 0 & 0 & \ddots & 1 \\ a_0 + a_1 X + \dots + (a_{p-1} - X) X^{p-1} & \dots & \dots & a_{p-1} - X \end{vmatrix} = (-1)^{p+2} P(X) \times \begin{vmatrix} 1 & & & 0 \\ -X & \ddots & & \\ & \ddots & \ddots & \\ 0 & & -X & 1 \end{vmatrix}

Et donc χA(X)=(−1)pP(X)\chi_A(X) = (-1)^p P(X). □\square

Ainsi, pour diagonaliser AA afin de calculer AnA^n, on doit résoudre P(X)=0P(X) = 0 pour chercher les valeurs propres.

Lien avec l'an dernier.

un+2=aun+1+bunu_{n+2} = a u_{n+1} + b u_n

a pour équation r2=ar+br^2 = ar + b. Les racines sont en fait les valeurs propres de (01ba)\begin{pmatrix} 0 & 1 \\ b & a \end{pmatrix}. De sorte que (01ba)n=P(λ1n00λ2n)P−1\begin{pmatrix} 0 & 1 \\ b & a \end{pmatrix}^n = P \begin{pmatrix} \lambda_1^n & 0 \\ 0 & \lambda_2^n \end{pmatrix} P^{-1}, où les λi\lambda_i sont ses racines.