Mathematics · Livro 4 · Bachelor Year 2

Matemática universitária — Graduação 2

Matemática universitária — Graduação 2 · Bachelor Year 2

3Redução de endomorfismos

Para entender um endomorfismo, procure as direções que ele apenas estica. Este capítulo constrói a maquinaria — autovalores, polinômios característico e minimal, o lema da decomposição em núcleos — e suas recompensas: critérios de diagonalização e de triangularização, Cayley–Hamilton, a decomposição de Dunford e o cálculo de potências e exponenciais de que o Capítulo 16 se alimentará. Em todo o capítulo, EE é um KK-espaço vetorial de dimensão finita (K=RK = \R ou C\C) e u∈L(E)u \in \mathcal{L}(E), n=dim⁡En = \dim E.

3.1 Autovalores e autovetores

Definição 3.1

λ∈K\lambda \in K é um autovalor de uu quando u(x)=λxu(x) = \lambda x para algum x≠0x \neq 0 (um autovetor); o autoespaço é Eλ(u)=ker⁡(u−λ id)E_\lambda(u) = \ker(u - \lambda\,\mathrm{id}). O conjunto dos autovalores é o espectro Sp⁡(u)\operatorname{Sp}(u). Um subespaço FF é estável quando u(F)⊆Fu(F) \subseteq F; os autoespaços são estáveis, e os subespaços estáveis permitem endomorfismos induzidos u∣Fu|_F.

Teorema 3.2 (Independência dos autoespaços)

Autovetores associados a autovalores dois a dois distintos formam uma família livre; equivalentemente, a soma dos autoespaços Eλ1+⋯+EλrE_{\lambda_1} + \dots + E_{\lambda_r} (λi\lambda_i distintos) é direta. Em particular, uu tem no máximo nn autovalores.

Demonstração. Por indução sobre rr. Suponha x1+⋯+xr=0x_1 + \dots + x_r = 0 com xi∈Eλix_i \in E_{\lambda_i}, sendo o enunciado conhecido para r−1r - 1. Aplique uu e subtraia λr\lambda_r vezes a relação:

∑i=1r−1(λi−λr) xi=0,\sum_{i=1}^{r-1} (\lambda_i - \lambda_r)\, x_i = 0 ,

logo, por indução, cada (λi−λr)xi=0(\lambda_i - \lambda_r)x_i = 0, isto é, xi=0x_i = 0 para i<ri < r, e então xr=0x_r = 0. Somas diretas de espaços não nulos num espaço de dimensão nn têm no máximo nn parcelas. ∎

A matriz A = psmallmatrix2 & 1\\ 1 & 2 psmallmatrix agindo no plano: o vetor genérico e_1 é tirado de sua reta, mas as direções próprias v_1 = (1,1) e v_2 = (1,-1) são apenas esticadas — por 3 e por 1 (de modo que Av_2 = v_2: a imagem tracejada coincide com v_2). A diagonalização é a mudança para a base (v_1, v_2), na qual A se torna diag(3, 1).
A matriz A=(2112)A = \left(\begin{smallmatrix}2 & 1\\ 1 & 2\end{smallmatrix}\right) agindo no plano: o vetor genérico e1e_1 é tirado de sua reta, mas as direções próprias v1=(1,1)v_1 = (1,1) e v2=(1,−1)v_2 = (1,-1) são apenas esticadas — por 33 e por 11 (de modo que Av2=v2Av_2 = v_2: a imagem tracejada coincide com v2v_2). A diagonalização é a mudança para a base (v1,v2)(v_1, v_2), na qual AA se torna diag⁡(3,1)\operatorname{diag}(3, 1).

Definição 3.3 (Polinômio característico)

χu(X)=det⁡(X id−u)\chi_u(X) = \det(X\,\mathrm{id} - u) — calculado em qualquer base como det⁡(XIn−A)\det(XI_n - A), um polinômio mônico de grau nn, invariante por semelhança (Teorema 2.17). Suas raízes em KK são exatamente os autovalores (λ\lambda é autovalor   ⟺  u−λ id\iff u - \lambda\,\mathrm{id} não injetivo   ⟺  χu(λ)=0\iff \chi_u(\lambda) = 0), e

χu(X)=Xn−(tr⁡u) Xn−1+⋯+(−1)ndet⁡u.\chi_u(X) = X^n - (\operatorname{tr} u)\, X^{n-1} + \dots + (-1)^n \det u .

A multiplicidade algébrica mλm_\lambda de um autovalor é sua multiplicidade como raiz de χu\chi_u; a multiplicidade geométrica é dim⁡Eλ\dim E_\lambda, e 1≤dim⁡Eλ≤mλ1 \leq \dim E_\lambda \leq m_\lambda.

Demonstração dos fatos enunciados. As afirmações sobre os coeficientes: desenvolva det⁡(XI−A)\det(XI - A) pela fórmula das permutações; a permutação identidade contribui com ∏i(X−aii)=Xn−(∑aii)Xn−1+…\prod_i (X - a_{ii}) = X^n - (\sum a_{ii})X^{n-1} + \dots, e toda outra permutação fixa no máximo n−2n - 2 posições diagonais, contribuindo com grau ≤n−2\leq n - 2: os dois coeficientes de topo são os anunciados; X=0X = 0 dá o termo constante det⁡(−A)=(−1)ndet⁡A\det(-A) = (-1)^n\det A.

Geométrica ≤\leq algébrica: seja d=dim⁡Eλd = \dim E_\lambda e complete uma base de EλE_\lambda numa base de EE; a matriz de uu é triangular superior por blocos com bloco superior esquerdo λId\lambda I_d, logo χu(X)=(X−λ)d χ(bloco inferior)(X)\chi_u(X) = (X - \lambda)^d\, \chi_{\text{(bloco inferior)}}(X): a multiplicidade de λ\lambda vale ao menos dd. ∎

Exemplo 3.4 (Mesmo χ\chi, geometria diferente)

As matrizes

(2002)e(2102)\begin{pmatrix}2 & 0\\ 0 & 2\end{pmatrix} \qquad\text{e}\qquad \begin{pmatrix}2 & 1\\ 0 & 2\end{pmatrix}

têm o mesmo polinômio característico (X−2)2(X - 2)^2, o mesmo traço, o mesmo determinante, o mesmo espectro — e no entanto não são semelhantes: a primeira tem E2E_2 de dimensão 22 (multiplicidade geométrica 22), a segunda de dimensão 11. O polinômio característico enxerga apenas multiplicidades algébricas; as dimensões dos autoespaços são o invariante mais fino, e o polinômio minimal arbitra (X−2X - 2 contra (X−2)2(X - 2)^2). Moral para toda discussão de diagonalizabilidade: χ\chi seleciona os candidatos, mas são os núcleos que depositam os votos.

Definição 3.5 (Diagonalizável, triangularizável)

uu é diagonalizável quando EE tem uma base de autovetores (matricialmente: semelhante a uma matriz diagonal); triangularizável quando sua matriz em alguma base é triangular superior.

Teorema 3.6 (Critérios de diagonalizabilidade)

As afirmações seguintes são equivalentes:

  1. uu é diagonalizável;
  2. E=⨁λ∈Sp⁡uEλE = \bigoplus_{\lambda \in \operatorname{Sp} u} E_\lambda;
  3. χu\chi_u se decompõe sobre KK e dim⁡Eλ=mλ\dim E_\lambda = m_\lambda para todo autovalor;
  4. (suficiente, não necessária) χu\chi_u tem nn raízes distintas em KK.

Demonstração. (1   ⟺  \iff 2): uma base de autovetores se separa em bases dos EλE_\lambda e, reciprocamente, concatenar bases das parcelas diretas dá uma base de EE (o Teorema 3.2 torna a soma direta; a igualdade de dimensões faz dela o espaço todo).

(2   ⟺  \iff 3): na base diagonal, χu=∏(X−λ)dim⁡Eλ\chi_u = \prod (X - \lambda)^{\dim E_\lambda} se decompõe com multiplicidades correspondentes. Reciprocamente, suponha que χu\chi_u se decomponha com dim⁡Eλ=mλ\dim E_\lambda = m_\lambda em toda parte; então a soma direta dos autoespaços (direta pelo Teorema 3.2) tem dimensão

∑λdim⁡Eλ=∑λmλ=deg⁡χu=n,\sum_{\lambda}\dim E_\lambda = \sum_{\lambda} m_\lambda = \deg\chi_u = n ,

a igualdade do meio porque o grau de um polinômio que se decompõe é a soma das multiplicidades de suas raízes: a soma é todo o EE. Note onde cada hipótese trabalhou: a decomposição preencheu o grau, a igualdade das multiplicidades preencheu as dimensões.

(4 ⇒\Rightarrow 1): nn autovalores distintos dão nn autovetores independentes (Teorema 3.2): uma base. ∎

Método 3.7 (Decidir a diagonalizabilidade)

Na prática, teste nesta ordem — cada etapa pode encerrar o serviço. (1) Um polinômio anulador com raízes simples e que se decompõe se apresenta sozinho (u2=idu^2 = \mathrm{id}, u2=uu^2 = u, uk=idu^k = \mathrm{id})? Se sim: diagonalizável, sem cálculo algum (Corolário 3.17 adiante). (2) Calcule χu\chi_u; se ele tiver nn raízes distintas em KK: diagonalizável (Teorema 3.6 (4)). (3) Caso contrário, apenas para cada raiz múltipla λ\lambda, compare dim⁡ker⁡(u−λ id)\dim\ker(u - \lambda\,\mathrm{id}) com a multiplicidade mλm_\lambda: qualquer déficit mata a diagonalizabilidade; a igualdade em toda parte a demonstra. Nunca calcule autoespaços de raízes simples (sua dimensão está forçada a ser 11), e nunca triangularize só para decidir.

Exemplo 3.8 (A diagonalização posta a trabalhar)

A=I+J=(211121112)A = I + J = \left(\begin{smallmatrix}2 & 1 & 1\\ 1 & 2 & 1\\ 1 & 1 & 2\end{smallmatrix}\right), com JJ a matriz de uns: de Sp⁡(J)={3,0}\operatorname{Sp}(J) = \{3, 0\} (Exemplo 2.19), Sp⁡(A)={4,1}\operatorname{Sp}(A) = \{4, 1\}, com autoespaços R(1,1,1)\R(1,1,1) e o plano {x+y+z=0}\{x + y + z = 0\}: dimensões 1+2=31 + 2 = 3, diagonalizável (Teorema 3.6 (2)). Potências sem nenhuma matriz de mudança de base: com Π=J/3\Pi = J/3 o projetor sobre R(1,1,1)\R(1,1,1),

A=4 Π+1⋅(I−Π)⟹Ak=4k Π+(I−Π)=4k−13 J+I.A = 4\,\Pi + 1\cdot(I - \Pi) \quad\Longrightarrow\quad A^k = 4^k\,\Pi + (I - \Pi) = \frac{4^k - 1}{3}\,J + I .

(Confira k=1k = 1: 4−13J+I=A\frac{4-1}3 J + I = A.) A lição final: quando os autoespaços estão visíveis, os projetores espectrais calculam potências mais depressa do que PDP−1PDP^{-1} jamais fará — e a fórmula exibe a dinâmica: AkA^k cresce como 4k4^k ao longo de (1,1,1)(1,1,1) e fica parado no plano ortogonal.

Teorema 3.9 (Triangularização)

uu é triangularizável sobre KK se e somente se χu\chi_u se decompõe sobre KK. Em particular, todo endomorfismo de um C\C-espaço vetorial é triangularizável.

Demonstração. (⇒\Rightarrow) O polinômio característico de uma matriz triangular é ∏(X−tii)\prod(X - t_{ii}): decompõe-se.

(⇐\Leftarrow) Indução sobre nn. Como χu\chi_u se decompõe, ele tem uma raiz λ\lambda: escolha um autovetor e1e_1. Numa base que começa por e1e_1, a matriz é (λ∗0B)\begin{pmatrix} \lambda & \ast\\ 0 & B\end{pmatrix}, e χu=(X−λ)χB\chi_u = (X - \lambda)\chi_B: χB\chi_B também se decompõe. Pela hipótese de indução aplicada à matriz (n−1)×(n−1)(n-1) \times (n-1) BB, existe uma QQ invertível com Q−1BQQ^{-1}BQ triangular superior; conjugar a matriz inteira por (100Q)\begin{pmatrix}1 & 0\\ 0 & Q\end{pmatrix} a triangulariza. ∎

Exemplo 3.10 (Triangularizando à mão)

B=(3−111)B = \begin{pmatrix}3 & -1\\ 1 & 1\end{pmatrix}: χB=X2−4X+4=(X−2)2\chi_B = X^2 - 4X + 4 = (X - 2)^2, e ker⁡(B−2I)=ker⁡(1−11−1)\ker(B - 2I) = \ker\left(\begin{smallmatrix}1 & -1\\ 1 & -1\end{smallmatrix}\right) é a reta gerada por e1′=(1,1)e_1' = (1, 1): um único autovalor, um autoespaço de dimensão um — não diagonalizável, mas triangularizável (Teorema 3.9). Complete a base com e2′=(1,0)e_2' = (1, 0) e calcule:

u(e1′)=(2,2)=2e1′,u(e2′)=(3,1)=1⋅e1′+2 e2′,u(e_1') = (2, 2) = 2e_1', \qquad u(e_2') = (3, 1) = 1\cdot e_1' + 2\, e_2' ,

de modo que, na base (e1′,e2′)(e_1', e_2'), a matriz é T=(2102)T = \left(\begin{smallmatrix}2 & 1\\ 0 & 2\end{smallmatrix}\right). A lição final: a diagonal de TT estava forçada (as duas entradas têm de ser o autovalor duplo 22); só a entrada do canto dependia da escolha de e2′e_2', e reescalar e2′e_2' pode torná-la qualquer valor não nulo — o resistente “11” é a sombra da parte nilpotente que Dunford vai isolar.

3.2 Polinômios de um endomorfismo

Definição 3.11

Para P=∑akXk∈K[X]P = \sum a_k X^k \in K[X], ponha P(u)=∑akuk∈L(E)P(u) = \sum a_k u^k \in \mathcal{L}(E). A aplicação P↦P(u)P \mapsto P(u) é um morfismo de álgebras K[X]→L(E)K[X] \to \mathcal{L}(E) (Definição 1.33); seu núcleo {P:P(u)=0}\{P : P(u) = 0\} é um ideal de K[X]K[X], não nulo (a família (id,u,…,un2)(\mathrm{id}, u, \dots, u^{n^2}) é ligada no n2n^2 de dimensão L(E)\mathcal{L}(E)), logo gerado por um único polinômio mônico μu\mu_u: o polinômio minimal (Teorema 1.26).

Proposição 3.12

  1. P(u)=0  ⟺  μu∣PP(u) = 0 \iff \mu_u \mid P; os autovalores de uu são raízes de todo polinômio anulador, e as raízes de μu\mu_u são exatamente os autovalores.
  2. Se FF é estável, μu∣F∣μu\mu_{u|_F} \mid \mu_u.

Demonstração. (1) A divisibilidade é a definição de gerador. Se u(x)=λxu(x) = \lambda x, x≠0x \neq 0, então 0=P(u)(x)=P(λ)x0 = P(u)(x) = P(\lambda)x, logo P(λ)=0P(\lambda) = 0: os autovalores são raízes dos anuladores, em particular de μu\mu_u. Reciprocamente, se λ\lambda é raiz, μu=(X−λ)Q\mu_u = (X - \lambda)Q com Q(u)≠0Q(u) \neq 0 (o grau de μu\mu_u é mínimo): escolha yy com Q(u)(y)≠0Q(u)(y) \neq 0; então (u−λ)(Q(u)(y))=μu(u)(y)=0(u - \lambda)(Q(u)(y)) = \mu_u(u)(y) = 0 exibe o autovetor Q(u)(y)Q(u)(y).

(2) μu(u∣F)=μu(u)∣F=0\mu_u(u|_F) = \mu_u(u)|_F = 0, e aplique (1) a u∣Fu|_F. ∎

Exemplo 3.13 (Polinômios minimais achados à mão)

O polinômio minimal se calcula testando graus sucessivos. Para a matriz de uns J∈M3(R)J \in \mathcal{M}_3(\R): J≠λIJ \neq \lambda I (o grau 11 está fora), e J2=3JJ^2 = 3J, logo

μJ=X2−3X=X(X−3):\mu_J = X^2 - 3X = X(X - 3) :

grau 22, decompõe-se, raízes simples — JJ é diagonalizável com espectro {0,3}\{0, 3\} (Corolário 3.17 adiante), confirmando Exemplo 2.19 sem um único determinante. Para a matriz de troca AA de Exemplo 3.15: A≠±IA \neq \pm I e A2=IA^2 = I dão μA=X2−1\mu_A = X^2 - 1. Nos dois casos o padrão é o mesmo: adivinhe uma identidade de grau baixo a partir da estrutura (o posto um força J2=(tr⁡J) JJ^2 = (\operatorname{tr}J)\,J; uma involução força A2=IA^2 = I) e depois verifique que nenhum divisor próprio anula. Os polinômios minimais são em geral achados, e não calculados a partir de χ\chi.

Teorema 3.14 (Lema da decomposição em núcleos)

Se P=P1P2⋯PrP = P_1 P_2 \cdots P_r com os PiP_i dois a dois coprimos, então

ker⁡P(u)=ker⁡P1(u)⊕⋯⊕ker⁡Pr(u),\ker P(u) = \ker P_1(u) \oplus \dots \oplus \ker P_r(u),

e as projeções sobre as parcelas são polinômios em uu.

Demonstração. Basta tratar o caso r=2r = 2 e fazer indução. Bézout em K[X]K[X] (Teorema 1.26): UP1+VP2=1U P_1 + V P_2 = 1, de modo que, para todo xx,

x=U(u)P1(u)(x)⏟=: x2+V(u)P2(u)(x)⏟=: x1.x = \underbrace{U(u)P_1(u)(x)}_{=:\,x_2} + \underbrace{V(u)P_2(u)(x)}_{=:\,x_1}.

Se x∈ker⁡P(u)x \in \ker P(u): P2(u)(x2)=U(u) P(u)(x)=0P_2(u)(x_2) = U(u)\,P(u)(x) = 0 (polinômios em uu comutam), logo x2∈ker⁡P2(u)x_2 \in \ker P_2(u) e, simetricamente, x1∈ker⁡P1(u)x_1 \in \ker P_1(u): a soma preenche ker⁡P(u)\ker P(u); as duas parcelas ficam dentro de ker⁡P(u)\ker P(u) (Pi∣PP_i \mid P). Diretividade: x∈ker⁡P1(u)∩ker⁡P2(u)x \in \ker P_1(u) \cap \ker P_2(u) dá x=U(u)P1(u)x+V(u)P2(u)x=0x = U(u)P_1(u)x + V(u)P_2(u)x = 0. As fórmulas para x1,x2x_1, x_2 exibem as projeções como V(u)P2(u)V(u)P_2(u) e U(u)P1(u)U(u)P_1(u). ∎

Exemplo 3.15 (O lema dos núcleos com projetores explícitos)

Seja A=(010100001)A = \left(\begin{smallmatrix}0 & 1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 1\end{smallmatrix}\right) (troca as duas primeiras coordenadas). Então A2=IA^2 = I: o polinômio X2−1=(X−1)(X+1)X^2 - 1 = (X - 1)(X + 1) anula AA, seus fatores são coprimos e Bézout é explícito:

12(X+1)−12(X−1)=1.\frac{1}{2}(X + 1) - \frac12(X - 1) = 1 .

Seguindo a demonstração do Teorema 3.14, as projeções sobre ker⁡(A−I)\ker(A - I) e ker⁡(A+I)\ker(A + I) são os polinômios em AA

π+=A+I2=12(110110002),π−=I−A2=12(1−10−110000).\pi_+ = \frac{A + I}{2} = \frac12\begin{pmatrix} 1 & 1 & 0\\ 1 & 1 & 0\\ 0 & 0 & 2\end{pmatrix}, \qquad \pi_- = \frac{I - A}{2} = \frac12\begin{pmatrix} 1 & -1 & 0\\ -1 & 1 & 0\\ 0 & 0 & 0\end{pmatrix}.

Verificação: π++π−=I\pi_+ + \pi_- = I, π+π−=0\pi_+\pi_- = 0, π±2=π±\pi_\pm^2 = \pi_\pm, e as imagens são o plano {x=y}\{x = y\} (vetores simétricos, autovalor 11) e a reta R(1,−1,0)\R(1, -1, 0) (antissimétricos, autovalor −1-1). O lema dos núcleos não é um enunciado de existência: os coeficientes de Bézout são as fórmulas dos projetores.

Exemplo 3.16 (Os projetores também calculam a exponencial)

A mesma matriz de troca, um dividendo adiante. Como A=π+−π−A = \pi_+ - \pi_- com projetores ortogonais no sentido algébrico (π+π−=0\pi_+\pi_- = 0), toda potência obedece a Ak=π++(−1)kπ−A^k = \pi_+ + (-1)^k\pi_-, e a série exponencial se reagrupa por projetor:

etA=∑ktkk!(π++(−1)kπ−)=et π++e−t π−=(cosh⁡tsinh⁡t0sinh⁡tcosh⁡t000et).\eu^{tA} = \sum_k \frac{t^k}{k!}\bigl(\pi_+ + (-1)^k\pi_-\bigr) = \eu^{t}\,\pi_+ + \eu^{-t}\,\pi_- = \begin{pmatrix} \cosh t & \sinh t & 0\\ \sinh t & \cosh t & 0\\ 0 & 0 & \eu^{t} \end{pmatrix}.

(Confira t=0t = 0: a identidade; derivada em 00: AA.) A decomposição espectral converte uma série de matrizes em duas séries escalares — exatamente o mecanismo que o Capítulo 16 vai acionar em todo sistema diagonalizável, e a razão pela qual as funções hiperbólicas governam acoplamentos simétricos.

Corolário 3.17 (Diagonalizabilidade pelo polinômio minimal)

uu é diagonalizável   ⟺  \iff μu\mu_u se decompõe sobre KK com raízes simples   ⟺  \iff algum polinômio anulador de uu se decompõe com raízes simples.

Demonstração. Se P(u)=0P(u) = 0 com P=∏i(X−λi)P = \prod_{i}(X - \lambda_i) (λi\lambda_i distintos), o lema dá E=ker⁡P(u)=⨁iker⁡(u−λi)E = \ker P(u) = \bigoplus_i \ker(u - \lambda_i): uma soma direta de autoespaços, logo uu é diagonalizável (Teorema 3.6). Reciprocamente, um uu diagonalizável é anulado por ∏λ∈Sp⁡u(X−λ)\prod_{\lambda \in \operatorname{Sp}u}(X - \lambda) (que anula cada autoespaço), o qual se decompõe com raízes simples; e μu\mu_u o divide tendo as mesmas raízes (Proposição 3.12): μu\mu_u é exatamente esse produto. ∎

Exemplo 3.18

As projeções satisfazem p2=pp^2 = p: anuladas por X(X−1)X(X-1), que se decompõe com raízes simples — diagonalizáveis com espectro ⊆{0,1}\subseteq \{0, 1\}, e E=ker⁡p⊕ker⁡(p−id)E = \ker p \oplus \ker(p - \mathrm{id}): a análise geométrica do primeiro ano, redemonstrada em uma linha. As simetrias (s2=ids^2 = \mathrm{id}, anulador X2−1X^2 - 1): diagonalizáveis quando char⁡K≠2\operatorname{char} K \neq 2, espectro ⊆{±1}\subseteq \{\pm 1\}. Um endomorfismo com u3=u2u^3 = u^2 e u2≠uu^2 \neq u: anulado por X2(X−1)X^2(X - 1), não necessariamente diagonalizável — o critério o detecta (a raiz dupla 00 precisa ser testada: diagonalizável se e somente se além disso ker⁡u2=ker⁡u\ker u^2 = \ker u).

Exemplo 3.19 (O corpo decide: uma rotação em R3\R^3)

Seja RR o quarto de volta em torno do eixo zz:

R=(0−10100001),χR=(X−1)(X2+1).R = \begin{pmatrix} 0 & -1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 1 \end{pmatrix}, \qquad \chi_R = (X - 1)(X^2 + 1).

Sobre R\R: o único autovalor é 11, com autoespaço igual ao eixo Re3\R e_3 — uma reta de vetores fixos e nenhuma outra redução: RR não é diagonalizável nem triangularizável em M3(R)\mathcal{M}_3(\R) (χR\chi_R não se decompõe). Sobre C\C: três autovalores distintos 1,i,−i1, \iu, -\iu, de modo que RR é diagonalizável, com autovetores e3e_3 e e1∓ie2e_1 \mp \iu e_2. A geometria era audível na álgebra: as rotações do plano não têm direções invariantes reais, e os autovalores complexos ±i\pm\iu de módulo 11 guardam o ângulo (±π2\pm\frac\pi2) que a matriz real só consegue exprimir misturando coordenadas.

Exemplo 3.20 (Minimal contra característico)

Para D=diag⁡(2,2,3)D = \operatorname{diag}(2, 2, 3): χD=(X−2)2(X−3)\chi_D = (X - 2)^2(X - 3) mas μD=(X−2)(X−3)\mu_D = (X - 2)(X - 3), pois (D−2I)(D−3I)=0(D - 2I)(D - 3I) = 0 (verifique na base canônica), ao passo que nenhum dos fatores sozinho anula DD. Para o bloco de deslocamento N=(0100)⊕(3)N = \left(\begin{smallmatrix}0 & 1\\ 0 & 0\end{smallmatrix}\right) \oplus (3), isto é, N′=(010000003)N' = \left(\begin{smallmatrix}0 & 1 & 0\\ 0 & 0 & 0\\ 0 & 0 & 3\end{smallmatrix}\right): χN′=X2(X−3)\chi_{N'} = X^2(X - 3) e μN′=X2(X−3)\mu_{N'} = X^2(X - 3) — a raiz dupla é genuinamente necessária, pois N′N' não é diagonalizável do lado ker⁡\ker (N′e2=e1≠0N'e_2 = e_1 \neq 0). Regra prática: μ\mu e χ\chi têm as mesmas raízes (Proposição 3.12); a multiplicidade em μ\mu mede o tamanho do maior bloco nilpotente, e a de χ\chi a dimensão total do subespaço característico.

Teorema 3.21 (Cayley–Hamilton)

χu(u)=0\chi_u(u) = 0; por consequência, μu∣χu\mu_u \mid \chi_u e deg⁡μu≤n\deg \mu_u \leq n.

Demonstração. Fixe x≠0x \neq 0 e seja dd maximal com (x,u(x),…,ud−1(x))(x, u(x), \dots, u^{d-1}(x)) livre; escreva

ud(x)=−a0x−a1u(x)−⋯−ad−1ud−1(x),u^d(x) = -a_0 x - a_1 u(x) - \dots - a_{d-1}u^{d-1}(x),

e ponha Px=Xd+ad−1Xd−1+⋯+a0P_x = X^d + a_{d-1}X^{d-1} + \dots + a_0, de modo que Px(u)(x)=0P_x(u)(x) = 0. Complete a família livre numa base de EE: nela, uu tem a forma por blocos (C∗0D)\begin{pmatrix} C & \ast\\ 0 & D\end{pmatrix}, em que CC é a matriz companheira de PxP_x, cujo polinômio característico é PxP_x (desenvolva det⁡(XI−C)\det(XI - C) ao longo da primeira coluna, por indução sobre dd). Logo χu=Px⋅χD\chi_u = P_x \cdot \chi_D, e

χu(u)(x)=χD(u)(Px(u)(x))=0.\chi_u(u)(x) = \chi_D(u)\bigl(P_x(u)(x)\bigr) = 0 .

O argumento vale para todo xx: χu(u)=0\chi_u(u) = 0. ∎

Exemplo 3.22 (Cayley–Hamilton em ação)

A=(1234)A = \begin{pmatrix} 1 & 2\\ 3 & 4\end{pmatrix}: χA=X2−5X−2\chi_A = X^2 - 5X - 2, logo A2=5A+2IA^2 = 5A + 2I. Toda potência de AA desaba para uma combinação de II e AA:

A4=(5A+2I)2=25A2+20A+4I=145A+54I=(199290435634),A^4 = (5A + 2I)^2 = 25A^2 + 20A + 4I = 145A + 54I = \begin{pmatrix} 199 & 290\\ 435 & 634\end{pmatrix},

e a inversa vem de brinde: A(A−5I)=2IA(A - 5I) = 2I dá

A−1=12(A−5I)=(−213/2−1/2).A^{-1} = \tfrac12(A - 5I) = \begin{pmatrix} -2 & 1\\ 3/2 & -1/2\end{pmatrix}.

A lição final: Cayley–Hamilton comprime toda a álgebra K[A]K[A] em Vect⁡(I,A,…,An−1)\operatorname{Vect}(I, A, \dots, A^{n-1}) — dim⁡K[A]=deg⁡μA≤n\dim K[A] = \deg\mu_A \leq n, por maiores que sejam as potências de que você precise.

Observação 3.23 (Armadilhas comuns)

(i) Autovalores não se somam: Sp⁡(A+B)\operatorname{Sp}(A + B) não é Sp⁡A+Sp⁡B\operatorname{Sp}A + \operatorname{Sp}B, e uma soma de matrizes diagonalizáveis não precisa ser diagonalizável — (1100)+(0001)=(1101)\left(\begin{smallmatrix}1 & 1\\ 0 & 0\end{smallmatrix}\right) + \left(\begin{smallmatrix}0 & 0\\ 0 & 1\end{smallmatrix}\right) = \left(\begin{smallmatrix}1 & 1\\ 0 & 1\end{smallmatrix}\right) é soma de duas matrizes diagonalizáveis (cada uma com autovalores distintos) e não é diagonalizável; só as famílias que comutam se comportam bem (Exercício 3.9). (ii) “χu\chi_u se decompõe” é uma hipótese sobre o corpo: uma rotação plana tem χ=X2−2cos⁡θ X+1\chi = X^2 - 2\cos\theta\,X + 1, que se decompõe sobre C\C, não sobre R\R — diagonalizável em M2(C)\mathcal{M}_2(\C), não triangularizável em M2(R)\mathcal{M}_2(\R). (iii) A desigualdade vai da geométrica ≤\leq à algébrica, nunca ao contrário; testar apenas dim⁡Eλ≥1\dim E_\lambda \geq 1 nada prova sobre a diagonalizabilidade. (iv) μu\mu_u não é χu\chi_u: a igualdade vale exatamente quando cada autovalor tem uma única cadeia de blocos (por exemplo, as matrizes companheiras, no problema de fim de semana deste capítulo); usar χ\chi onde μ\mu é necessário infla todo cálculo de potências. (v) As partes dd e ν\nu de Dunford são polinômios em uu — uma decomposição u=d′+ν′u = d' + \nu' com as propriedades certas mas com d′ν′≠ν′d′d'\nu' \neq \nu'd' não é a de Dunford e nunca é única.

Observação 3.24 (Onde este capítulo é usado)

A redução é o cavalo de batalha do resto do livro: potências e exponenciais de matrizes movem os sistemas diferenciais lineares do Capítulo 16; o teorema espectral do Capítulo 12 é a diagonalização tornada ortogonal; as funções geradoras (Capítulo 23) redemonstram analiticamente as assintóticas de recorrência do problema de fim de semana deste capítulo. No volume do terceiro ano de graduação, o mesmo programa roda em dimensão infinita: a teoria espectral dos operadores compactos autoadjuntos, na qual sequências de autovalores substituem os espectros finitos, e a teoria de Perron–Frobenius das matrizes positivas, que explica por que os autovalores dominantes de problemas de contagem são positivos e simples.

3.3 Nilpotentes e a decomposição de Dunford

Proposição 3.25 (Endomorfismos nilpotentes)

Para uu com χu\chi_u que se decompõe, as afirmações seguintes são equivalentes: un=0u^n = 0; uk=0u^k = 0 para algum kk; Sp⁡(u)={0}\operatorname{Sp}(u) = \{0\}; χu=Xn\chi_u = X^n; uu é triangularizável com diagonal nula. Um endomorfismo nilpotente tem μu=X(ıˊndice de nilpoteˆncia)\mu_u = X^{\text{(índice de nilpotência)}}, e o índice satisfaz ≤n\leq n.

Demonstração. uk=0u^k = 0 faz de todo autovalor uma raiz de XkX^k: espectro {0}\{0\} (não vazio quando χ\chi se decompõe — sobre C\C sempre). Então χu=Xn\chi_u = X^n (todas as raízes nulas) e Cayley–Hamilton dá un=0u^n = 0; a triangularização (Teorema 3.9) põe zeros na diagonal (a diagonal carrega os autovalores). Reciprocamente, seja AA estritamente triangular superior: aij=0a_{ij} = 0 para j≤ij \leq i. Mostramos por indução que

(Ak)ij=0sempre que j≤i+k−1,(A^k)_{ij} = 0 \qquad \text{sempre que } j \leq i + k - 1,

isto é, cada potência empurra a região nula uma diagonal acima. Para k=1k = 1 essa é a hipótese. Para o passo,

(Ak+1)ij=∑ℓ(Ak)iℓ aℓj,(A^{k+1})_{ij} = \sum_{\ell} (A^k)_{i\ell}\,a_{\ell j} ,

e cada termo se anula: ou ℓ≤i+k−1\ell \leq i + k - 1 (o primeiro fator é 00 por indução), ou ℓ≥i+k\ell \geq i + k, e nesse caso j≤i+k≤ℓj \leq i + k \leq \ell mata o segundo fator. Em k=nk = n a condição j≤i+n−1j \leq i + n - 1 vale para todos i,j≤ni, j \leq n: An=0A^n = 0. O polinômio minimal divide XnX^n e a anulação define o índice. ∎

Teorema 3.26 (Decomposição de Dunford)

Suponha que χu\chi_u se decomponha sobre KK (automático para K=CK = \C). Então existe um único par (d,ν)(d, \nu) com

u=d+ν,d diagonalizaˊvel,ν nilpotente,dν=νd,u = d + \nu, \qquad d \text{ diagonalizável}, \quad \nu \text{ nilpotente}, \quad d\nu = \nu d ,

e além disso dd e ν\nu são polinômios em uu.

Demonstração. Existência. Escreva χu=∏i=1r(X−λi)mi\chi_u = \prod_{i=1}^{r} (X - \lambda_i)^{m_i} (λi\lambda_i distintos) e ponha Ni=ker⁡(u−λi)miN_i = \ker(u - \lambda_i)^{m_i}, os subespaços característicos. Por Cayley–Hamilton e pelo lema dos núcleos (Teorema 3.14),

E=N1⊕⋯⊕Nr,E = N_1 \oplus \dots \oplus N_r ,

com projeções πi\pi_i polinomiais em uu; cada NiN_i é estável (polinômios em uu comutam com uu). Defina d=∑iλiπid = \sum_i \lambda_i \pi_i: um polinômio em uu, diagonalizável (ele age como λi\lambda_i em NiN_i, logo EE se decompõe em seus autoespaços). Então ν=u−d\nu = u - d é um polinômio em uu (logo comuta com dd) e, em cada NiN_i, age como u−λiu - \lambda_i, com (u−λi)mi=0(u - \lambda_i)^{m_i} = 0 aí: νmax⁡mi=0\nu^{\max m_i} = 0 em cada parcela, logo ν\nu é nilpotente.

Unicidade. Seja u=d′+ν′u = d' + \nu' outro par desse tipo. Como d′d' e ν′\nu' comutam entre si, eles comutam com u=d′+ν′u = d' + \nu', logo com todo polinômio em uu — em particular com dd e ν\nu. Então d−d′d - d' é diagonalizável (duas aplicações diagonalizáveis que comutam são simultaneamente diagonalizáveis: Exercício 3.9) e é igual a ν′−ν\nu' - \nu, que é nilpotente: se νk=0\nu^k = 0 e ν′k′=0\nu'^{k'} = 0, a comutatividade autoriza o desenvolvimento binomial

(ν′−ν)k+k′−1=∑j=0k+k′−1(k+k′−1j) ν′ j (−ν)k+k′−1−j,(\nu' - \nu)^{k + k' - 1} = \sum_{j=0}^{k+k'-1}\binom{k + k' - 1}{j} \,\nu'^{\,j}\,(-\nu)^{k + k' - 1 - j} ,

no qual todo termo morre: ou j≥k′j \geq k' (primeiro fator nulo), ou k+k′−1−j≥kk + k' - 1 - j \geq k (segundo fator nulo), e uma das duas condições sempre vale. Um nilpotente diagonalizável é nulo (seu espectro é {0}\{0\} e ele é diagonal em alguma base): d=d′d = d', ν=ν′\nu = \nu'. ∎

Exemplo 3.27 (Potências e exponenciais)

A=(31−11)A = \begin{pmatrix} 3 & 1\\ -1 & 1\end{pmatrix}: χA=X2−4X+4=(X−2)2\chi_A = X^2 - 4X + 4 = (X-2)^2, um único autovalor 22, autoespaço de dimensão 11: não diagonalizável. Dunford: D=2ID = 2I, N=A−2I=(11−1−1)N = A - 2I = \begin{pmatrix} 1 & 1\\ -1 & -1\end{pmatrix}, N2=0N^2 = 0. Então

Ak=(2I+N)k=2kI+k 2k−1N,etA=e2t(I+tN),A^k = (2I + N)^k = 2^k I + k\,2^{k-1} N , \qquad \eu^{tA} = \eu^{2t}(I + tN),

pelo binômio de Newton para elementos que comutam, resp. pela série exponencial (Capítulo 16) separada em parcelas que comutam. A redução transforma dinâmica matricial em dinâmica escalar.

Observação 3.28 (Perspectivas dentro deste volume)

A redução é um entroncamento; eis os quatro ramais a observar. No Capítulo 5, normas adaptadas transformam “todos os autovalores de módulo <1< 1” em “alguma norma de operador <1< 1”, fazendo os espectros governarem a convergência de potências e séries. No Capítulo 16, a receita do Exemplo 3.27 torna-se a solução geral de X′=AXX' = AX: Dunford separa etA\eu^{tA} em blocos do tipo polinômio vezes exponencial, e a estabilidade se lê nas partes reais dos autovalores. No Capítulo 12, um produto escalar impõe o que a mera álgebra linear não consegue: as matrizes simétricas tornam-se ortogonalmente diagonalizáveis, sem nenhuma parte nilpotente. E no Capítulo 23, as assintóticas de autovalor dominante do problema de fim de semana deste capítulo reaparecem analiticamente, como a menor singularidade de uma função geradora — duas línguas para uma mesma taxa de crescimento.

3.4 Exercícios

Exercício 3.1 ★

Diagonalize (autovalores, bases dos autoespaços, PP invertível):

A=(1221),B=(011101110).A = \begin{pmatrix} 1 & 2\\ 2 & 1 \end{pmatrix}, \qquad B = \begin{pmatrix} 0 & 1 & 1\\ 1 & 0 & 1\\ 1 & 1 & 0 \end{pmatrix}.
Solução

Solução de Exercício 3.1.

AA: χA=X2−2X−3=(X−3)(X+1)\chi_A = X^2 - 2X - 3 = (X - 3)(X + 1). Autovetores: para 33: (1,1)(1,1); para −1-1: (1,−1)(1,-1). Logo P=(111−1)P = \begin{pmatrix} 1 & 1\\ 1 & -1\end{pmatrix} dá P−1AP=diag⁡(3,−1)P^{-1}AP = \operatorname{diag}(3, -1).

B=J−IB = J - I, em que JJ é a matriz de uns. JJ tem posto 11 com Jv=3vJv = 3v para v=(1,1,1)v = (1,1,1) e Jw=0Jw = 0 no plano x+y+z=0x + y + z = 0: o espectro de BB é {2,−1}\{2, -1\} com autoespaços Vect⁡(1,1,1)\operatorname{Vect}(1,1,1) (dimensão 11) e {x+y+z=0}\{x + y + z = 0\} (dimensão 22, base (1,−1,0),(1,0,−1)(1,-1,0), (1,0,-1)). PP com essas três colunas dá P−1BP=diag⁡(2,−1,−1)P^{-1}BP = \operatorname{diag}(2, -1, -1).

Exercício 3.2 ★

Mostre que C=(1101)C = \begin{pmatrix} 1 & 1\\ 0 & 1\end{pmatrix} não é diagonalizável, de duas maneiras: pelos autoespaços e pelo polinômio minimal.

Solução

Solução de Exercício 3.2.

Autoespaços: χC=(X−1)2\chi_C = (X-1)^2, único autovalor 11; ker⁡(C−I)=ker⁡(0100)\ker(C - I) = \ker\begin{pmatrix} 0&1\\ 0&0\end{pmatrix} é a reta Vect⁡(e1)\operatorname{Vect}(e_1): dimensão 1<2=m11 < 2 = m_1, logo não é diagonalizável (Teorema 3.6).

Polinômio minimal: μC\mu_C divide (X−1)2(X-1)^2 e C≠IC \neq I, logo μC=(X−1)2\mu_C = (X-1)^2: uma raiz dupla, logo não é diagonalizável (Corolário 3.17).

Exercício 3.3 ★

Seja uu tal que u2−5u+6 id=0u^2 - 5u + 6\,\mathrm{id} = 0. Prove que uu é diagonalizável, determine os espectros possíveis e calcule uku^k como combinação de id\mathrm{id} e uu.

Solução

Solução de Exercício 3.3.

X2−5X+6=(X−2)(X−3)X^2 - 5X + 6 = (X-2)(X-3): decompõe-se com raízes simples, logo uu é diagonalizável (Corolário 3.17), com Sp⁡(u)⊆{2,3}\operatorname{Sp}(u) \subseteq \{2, 3\}. Espectros possíveis: {2}\{2\} (u=2 idu = 2\,\mathrm{id}), {3}\{3\} (u=3 idu = 3\,\mathrm{id}) ou {2,3}\{2, 3\}.

Potências: procure uk=ak id+bk uu^k = a_k\,\mathrm{id} + b_k\,u. Nos autoespaços, isso se lê 2k=ak+2bk2^k = a_k + 2b_k e 3k=ak+3bk3^k = a_k + 3b_k: resolvendo, bk=3k−2kb_k = 3^k - 2^k, ak=3⋅2k−2⋅3ka_k = 3\cdot2^k - 2\cdot 3^k:

uk=(3⋅2k−2⋅3k) id+(3k−2k) u.u^k = (3\cdot 2^k - 2\cdot 3^k)\,\mathrm{id} + (3^k - 2^k)\, u .

(Válido para os três espectros: as identidades valem autovalor a autovalor.)

Exercício 3.4 ★★

Seja uu diagonalizável e FF um subespaço estável. Prove que u∣Fu|_F é diagonalizável (restrinja um polinômio anulador com raízes simples e que se decompõe).

Solução

Solução de Exercício 3.4.

uu diagonalizável: P=∏λ(X−λ)P = \prod_{\lambda}(X - \lambda) sobre o espectro anula uu, decompõe-se e tem raízes simples. Então P(u∣F)=P(u)∣F=0P(u|_F) = P(u)|_F = 0: a restrição é anulada por um polinômio que se decompõe com raízes simples, logo é diagonalizável (Corolário 3.17).

Exercício 3.5 ★★

(Fibonacci) Seja A=(1110)A = \begin{pmatrix} 1 & 1\\ 1 & 0\end{pmatrix}. Diagonalize AA sobre R\R e deduza a fórmula de Binet para a sequência de Fibonacci (F0=0F_0 = 0, F1=1F_1 = 1, Fn+1=Fn+Fn−1F_{n+1} = F_n + F_{n-1}):

Fn=φn−ψn5,φ=1+52, ψ=1−52.F_n = \frac{\varphi^n - \psi^n}{\sqrt 5}, \qquad \varphi = \frac{1 + \sqrt5}{2},\ \psi = \frac{1 - \sqrt5}{2}.
Solução

Solução de Exercício 3.5.

χA=X2−X−1\chi_A = X^2 - X - 1, de raízes φ\varphi e ψ\psi (distintas): diagonalizável, com autovetores (φ,1)(\varphi, 1) e (ψ,1)(\psi, 1). A recorrência dá (Fn+1Fn)=An(10)\begin{pmatrix} F_{n+1}\\ F_n \end{pmatrix} = A^n \begin{pmatrix}1\\ 0\end{pmatrix}. Decomponha (1,0)(1, 0) nos autovetores: (1,0)=1φ−ψ((φ,1)−(ψ,1))(1,0) = \frac{1}{\varphi - \psi}\bigl((\varphi, 1) - (\psi, 1)\bigr) com φ−ψ=5\varphi - \psi = \sqrt5. Aplicar AnA^n multiplica cada componente própria pela nn-ésima potência de seu autovalor; lendo a segunda coordenada:

Fn=φn−ψn5.F_n = \frac{\varphi^n - \psi^n}{\sqrt 5} .

(Confira: n=1n = 1 dá φ−ψ5=1\frac{\varphi - \psi}{\sqrt5} = 1.)

Exercício 3.6 ★★

Sejam u∈L(E)u \in \mathcal{L}(E) com u2u^2 diagonalizável e uu invertível (K=CK = \C). Prove que uu é diagonalizável. Dê um contraexemplo quando uu não é invertível.

Solução

Solução de Exercício 3.6.

Seja P=∏i(X−μi)P = \prod_i (X - \mu_i) um anulador de u2u^2 que se decompõe com raízes simples μi\mu_i (o espectro de u2u^2). Como uu é invertível, 00 não é autovalor de u2u^2 (det⁡u2=(det⁡u)2≠0\det u^2 = (\det u)^2 \neq 0), logo todos os μi≠0\mu_i \neq 0. Então

Q(X)=∏i(X2−μi)=∏i(X−μi)(X+μi)Q(X) = \prod_i (X^2 - \mu_i) = \prod_i (X - \sqrt{\mu_i})(X + \sqrt{\mu_i})

anula uu:   Q(u)=∏i(u2−μi id)=P(u2)=0\;Q(u) = \prod_i (u^2 - \mu_i\,\mathrm{id}) = P(u^2) = 0. Suas raízes ±μi\pm \sqrt{\mu_i} (raízes quadradas complexas) são duas a duas distintas porque os μi\mu_i são distintos e não nulos (μi=−μj\sqrt{\mu_i} = -\sqrt{\mu_j} daria μi=μj\mu_i = \mu_j). Decompõe-se com raízes simples: uu é diagonalizável.

Contraexemplo sem invertibilidade: u=(0100)u = \begin{pmatrix} 0 & 1\\ 0 & 0\end{pmatrix}: u2=0u^2 = 0 é diagonalizável, uu não é.

Exercício 3.7 ★★

Calcule a decomposição de Dunford, AkA^k e etA\eu^{tA} para

A=(210021002).A = \begin{pmatrix} 2 & 1 & 0\\ 0 & 2 & 1\\ 0 & 0 & 2 \end{pmatrix}.
Solução

Solução de Exercício 3.7.

A=2I+NA = 2I + N com NN o deslocamento (Ne2=e1N e_2 = e_1, Ne3=e2Ne_3 = e_2), N3=0N^3 = 0, N2=E13N^2 = E_{13}: isso é a decomposição de Dunford (2I2I diagonal, NN nilpotente, e comutam; a unicidade faz dela a decomposição). Binômio com termos que comutam:

Ak=2kI+k2k−1N+(k2)2k−2N2=(2kk2k−1(k2)2k−202kk2k−1002k),A^k = 2^k I + k 2^{k-1} N + \binom k2 2^{k-2} N^2 = \begin{pmatrix} 2^k & k2^{k-1} & \binom k2 2^{k-2}\\ 0 & 2^k & k2^{k-1}\\ 0 & 0 & 2^k \end{pmatrix},
etA=e2t(I+tN+t22N2)=e2t(1tt2/201t001).\eu^{tA} = \eu^{2t}\Bigl(I + tN + \frac{t^2}{2}N^2\Bigr) = \eu^{2t}\begin{pmatrix} 1 & t & t^2/2\\ 0 & 1 & t\\ 0 & 0 & 1 \end{pmatrix}.

Exercício 3.8 ★★

Seja A∈Mn(C)A \in \mathcal{M}_n(\C) com Ak=IA^k = I para algum k≥1k \geq 1. Prove que AA é diagonalizável e que seus autovalores são raízes kk-ésimas da unidade. Deduza que uma matriz complexa invertível de ordem finita semelhante a uma matriz triangular com diagonal de uns é a identidade.

Solução

Solução de Exercício 3.8.

Xk−1X^k - 1 anula AA e se decompõe sobre C\C com as kk raízes distintas e2iπj/k\eu^{2\iu\pi j/k}: AA é diagonalizável (Corolário 3.17) e seus autovalores, raízes de Xk−1X^k - 1, são raízes kk-ésimas da unidade.

Se além disso AA é semelhante a uma matriz triangular com diagonal de uns: todos os autovalores valem 11, e AA, diagonalizável com único autovalor 11, é P I P−1=IP\,I\,P^{-1} = I.

Exercício 3.9 ★★★

(Diagonalização simultânea) Sejam u,vu, v diagonalizáveis e que comutam. Prove que eles são simultaneamente diagonalizáveis: alguma base diagonaliza os dois. (Cada autoespaço de uu é estável por vv; diagonalize aí as restrições de vv, usando Exercício 3.4.)

Solução

Solução de Exercício 3.9.

Escreva E=⨁λEλ(u)E = \bigoplus_\lambda E_\lambda(u) (Teorema 3.6). Cada Eλ(u)E_\lambda(u) é estável por vv: para x∈Eλx \in E_\lambda, u(v(x))=v(u(x))=λv(x)u(v(x)) = v(u(x)) = \lambda v(x). A restrição de vv a Eλ(u)E_\lambda(u) é diagonalizável (Exercício 3.4): escolha uma base de Eλ(u)E_\lambda(u) formada de autovetores de vv. Concatenando essas bases sobre todos os λ\lambda obtém-se uma base de EE cujos vetores são autovetores de ambos uu (por pertencerem a Eλ(u)E_\lambda(u)) e vv (por construção).

Exercício 3.10 ★★★

Seja u∈L(Cn)u \in \mathcal{L}(\C^n). Prove que uu é diagonalizável se e somente se todo subespaço uu-estável tem um subespaço suplementar uu-estável. (Para ⇐\Leftarrow: aplique a propriedade a F=∑λEλ(u)F = \sum_\lambda E_\lambda(u), a soma de todos os autoespaços; se um suplementar estável GG fosse não nulo, triangularizar u∣Gu|_G produziria um autovetor de uu dentro de GG — contradizendo G∩F={0}G \cap F = \{0\}.)

Solução

Solução de Exercício 3.10.

(⇒\Rightarrow) Seja uu diagonalizável e FF estável. Então u∣Fu|_F é diagonalizável (Exercício 3.4): FF tem uma base de autovetores, que se estende, dentro de cada autoespaço global EλE_\lambda, a uma base de EλE_\lambda (teorema da base incompleta dentro de EλE_\lambda, partindo da parte da base de FF que aí se encontra — note que F=⨁λ(F∩Eλ)F = \bigoplus_\lambda (F \cap E_\lambda), pois u∣Fu|_F é diagonalizável). Os vetores acrescentados geram um suplementar estável (cada um está em algum EλE_\lambda, logo o espaço que geram é estável por uu).

(⇐\Leftarrow) Seja F=∑λEλ(u)F = \sum_\lambda E_\lambda(u) (um subespaço estável) e GG um suplementar estável. Se G≠{0}G \neq \{0\}: χu∣G\chi_{u|_G} se decompõe sobre C\C, logo u∣Gu|_G tem um autovetor x∈Gx \in G (Teorema 3.9 ou diretamente a existência de uma raiz); mas todo autovetor de uu está em FF, logo x∈F∩G={0}x \in F \cap G = \{0\}: contradição. Portanto G={0}G = \{0\} e E=FE = F: os autoespaços preenchem EE, isto é, uu é diagonalizável.

Exercício 3.11 ★★★

(Raio espectral à la Gelfand, 2×22\times2 aperitivo da análise que vem) Seja A∈M2(C)A \in \mathcal{M}_2(\C) com os dois autovalores de módulo <1< 1. Prove que Ak→0A^k \to 0 entrada a entrada quando k→∞k \to \infty. (Triangularize: A=P(T)P−1A = P(T)P^{-1} com TT triangular superior; calcule TkT^k explicitamente — distinga autovalores iguais e distintos — e majore.)

Solução

Solução de Exercício 3.11.

Triangularize: A=PTP−1A = PTP^{-1}, T=(λc0μ)T = \begin{pmatrix} \lambda & c\\ 0 & \mu\end{pmatrix}, ∣λ∣,∣μ∣<1\abs\lambda, \abs\mu < 1. Então Ak=PTkP−1A^k = PT^kP^{-1}, e basta que Tk→0T^k \to 0.

Autovalores distintos: a indução dá

Tk=(λkc λk−μkλ−μ0μk),T^k = \begin{pmatrix} \lambda^k & c\,\dfrac{\lambda^k - \mu^k}{\lambda - \mu}\\[4pt] 0 & \mu^k \end{pmatrix},

e cada entrada tende a 00 (∣λ∣k,∣μ∣k→0\abs{\lambda}^k, \abs\mu^k \to 0).

Autovalores iguais (μ=λ\mu = \lambda): T=λI+cE12T = \lambda I + cE_{12} e Tk=λkI+kλk−1cE12T^k = \lambda^k I + k\lambda^{k-1}cE_{12}; a entrada kλk−1→0k\lambda^{k-1} \to 0, pois ∣λ∣<1\abs\lambda < 1 (a geométrica vence a polinomial). Nos dois casos Tk→0T^k \to 0 entrada a entrada, logo Ak=PTkP−1→0A^k = PT^kP^{-1} \to 0 (a multiplicação matricial por P,P−1P, P^{-1} fixas é contínua nas entradas — cada entrada do produto é uma combinação linear fixa).

Exercício 3.12 ★★

Seja u∈L(Cn)u \in \mathcal{L}(\C^n) com rk⁡u=1\operatorname{rk} u = 1 (n≥2n \geq 2). Mostre que χu=Xn−1(X−tr⁡u)\chi_u = X^{n-1}(X - \operatorname{tr} u) e que uu é diagonalizável se e somente se tr⁡u≠0\operatorname{tr} u \neq 0. (Lembre do Exercício 2.5 que u2=(tr⁡u) uu^2 = (\operatorname{tr} u)\,u.)

Solução

Solução de Exercício 3.12.

ker⁡u\ker u tem dimensão n−1n - 1 (teorema do núcleo e da imagem), logo 00 é autovalor de multiplicidade geométrica n−1n - 1, e χu\chi_u é divisível por Xn−1X^{n-1} (Definição 3.3: geométrica ≤\leq algébrica). Escreva χu=Xn−1(X−α)\chi_u = X^{n-1}(X - \alpha); sendo o coeficiente de Xn−1X^{n-1} igual a −tr⁡u-\operatorname{tr} u, obtemos α=tr⁡u\alpha = \operatorname{tr} u: χu=Xn−1(X−tr⁡u)\chi_u = X^{n-1}(X - \operatorname{tr} u).

Se tr⁡u≠0\operatorname{tr} u \neq 0: o autovalor tr⁡u\operatorname{tr} u é raiz de χu\chi_u, logo carrega um autovetor; os autoespaços de 00 e de tr⁡u\operatorname{tr} u têm dimensões n−1n - 1 e ≥1\geq 1, com soma ≥n\geq n: eles preenchem EE, e uu é diagonalizável (Teorema 3.6). Se tr⁡u=0\operatorname{tr} u = 0: pelo Exercício 2.5, u2=(tr⁡u)u=0u^2 = (\operatorname{tr} u)u = 0 com u≠0u \neq 0: uu é um nilpotente não nulo, e um nilpotente diagonalizável é nulo (Proposição 3.25): não é diagonalizável.

3.5 Problema: recorrências lineares e matrizes companheiras

Uma recorrência linear un+k=ak−1un+k−1+⋯+a0unu_{n+k} = a_{k-1}u_{n+k-1} + \dots + a_0 u_n é uma potência de matriz disfarçada, e a redução a converte em fórmulas fechadas, taxas de crescimento e estimativas de erro. Este problema de fim de semana desenvolve o dicionário — matrizes companheiras de um lado, o operador de deslocamento no espaço das sequências do outro — demonstra o teorema fundamental das recorrências lineares (a solução geral é ∑iQi(n)λin\sum_i Q_i(n)\lambda_i^n sobre as raízes do polinômio característico) e gasta os dividendos em aproximação diofantina de 2\sqrt2, na contagem de passeios e palavras, e num anel de sequências acopladas que só a diagonalização simultânea consegue desembaraçar.

Problema 3.1

Problema de fim de semana — o teorema fundamental das recorrências lineares

Fixe k≥1k \geq 1, escalares a0,…,ak−1∈Ca_0, \dots, a_{k-1} \in \C com a0≠0a_0 \neq 0, o polinômio mônico P=Xk−ak−1Xk−1−⋯−a1X−a0P = X^k - a_{k-1}X^{k-1} - \dots - a_1 X - a_0 e a recorrência

(R) ⁣:un+k=ak−1un+k−1+⋯+a1un+1+a0un(n≥0).(\mathcal R)\colon\quad u_{n+k} = a_{k-1}u_{n+k-1} + \dots + a_1 u_{n+1} + a_0 u_n \qquad (n \geq 0).

A matriz companheira de PP é

C=(01⋱⋱01a0a1⋯ak−1)∈Mk(C).C = \begin{pmatrix} 0 & 1 & & \\ & \ddots & \ddots & \\ & & 0 & 1\\ a_0 & a_1 & \cdots & a_{k-1} \end{pmatrix} \in \mathcal{M}_k(\C).

Parte I — O dicionário companheiro.

  1. Mostre que uma sequência (un)(u_n) satisfaz (R)(\mathcal R) se e somente se os vetores vn=(un,un+1,…,un+k−1)Tv_n = (u_n, u_{n+1}, \dots, u_{n+k-1})^{\mathsf T} satisfazem vn+1=Cvnv_{n+1} = Cv_n, e portanto vn=Cnv0v_n = C^n v_0.
  2. Prove que χC=P\chi_C = P (desenvolva det⁡(XI−C)\det(XI - C) ao longo da primeira coluna e faça indução sobre kk), e depois que também μC=P\mu_C = P (passe a CTC^{\mathsf T}, para o qual e1e_1 é cíclico, e note que uma matriz e sua transposta têm o mesmo polinômio minimal).
  3. Mostre que, para cada raiz λ\lambda de PP, o vetor (1,λ,…,λk−1)T(1, \lambda, \dots, \lambda^{k-1})^{\mathsf T} gera o autoespaço de CC associado a λ\lambda; deduza que todo autoespaço de CC tem dimensão 11 e que CC é diagonalizável se e somente se PP tem kk raízes distintas.
  4. Suponha que PP tenha raízes distintas λ1,…,λk\lambda_1, \dots, \lambda_k. Mostre que as sequências geométricas (λin)n(\lambda_i^n)_n formam uma base do espaço das soluções de (R)(\mathcal R), de modo que toda solução é un=∑iciλinu_n = \sum_i c_i\lambda_i^n para constantes cic_i únicas.
  5. Resolva completamente: un+2=un+1+6unu_{n+2} = u_{n+1} + 6u_n, u0=1u_0 = 1, u1=8u_1 = 8.

Parte II — O operador de deslocamento e o teorema fundamental. Seja S\mathcal{S} o C\C-espaço vetorial de todas as sequências complexas e S∈L(S)S \in \mathcal{L}(\mathcal{S}) o deslocamento, S((un)n)=(un+1)nS\bigl((u_n)_n\bigr) = (u_{n+1})_n.

  1. Mostre que o conjunto solução de (R)(\mathcal R) é ker⁡P(S)\ker P(S) e que ele tem dimensão exatamente kk (leve uma solução a seus valores iniciais).
  2. Explique por que o lema da decomposição em núcleos (Teorema 3.14) se aplica a SS no espaço S\mathcal{S} de dimensão infinita sem nenhuma mudança, e escreva a decomposição resultante de ker⁡P(S)\ker P(S) para P=∏i=1r(X−λi)miP = \prod_{i=1}^{r}(X - \lambda_i)^{m_i} (λi\lambda_i distintos, todos não nulos pois a0≠0a_0 \neq 0).
  3. Para λ≠0\lambda \neq 0 e m≥1m \geq 1, mostre que

    ker⁡ (S−λ id)m={ (Q(n) λn)n:Q∈Cm−1[X] },\ker\,(S - \lambda\,\mathrm{id})^m = \bigl\{\,\bigl(Q(n)\,\lambda^n\bigr)_n : Q \in \C_{m-1}[X]\,\bigr\},

    de dimensão mm. (Calcule (S−λ)(Q(n)λn)=λn+1(ΔQ)(n)(S - \lambda)\bigl(Q(n)\lambda^n\bigr) = \lambda^{n+1}(\Delta Q)(n) com ΔQ=Q(X+1)−Q(X)\Delta Q = Q(X + 1) - Q(X) e use que Δ\Delta baixa o grau; para a dimensão, majore-a por mm pelos valores iniciais.)

  4. (O teorema fundamental das recorrências lineares) Conclua: se P=∏i=1r(X−λi)miP = \prod_{i=1}^{r}(X - \lambda_i)^{m_i} com os λi\lambda_i distintos e não nulos, as soluções de (R)(\mathcal R) são exatamente as sequências

    un=∑i=1rQi(n) λin,Qi∈Cmi−1[X],u_n = \sum_{i=1}^{r} Q_i(n)\,\lambda_i^n, \qquad Q_i \in \C_{m_i - 1}[X],

    com polinômios QiQ_i unicamente determinados.

  5. Resolva completamente: un+2=4un+1−4unu_{n+2} = 4u_{n+1} - 4u_n, u0=1u_0 = 1, u1=0u_1 = 0, e confira a resposta em u2u_2.

Parte III — Raízes dominantes e dividendos diofantinos.

  1. Suponha as raízes simples com ∣λ1∣>∣λi∣\abs{\lambda_1} > \abs{\lambda_i} para i≥2i \geq 2, e un=∑iciλinu_n = \sum_i c_i \lambda_i^n com c1≠0c_1 \neq 0. Mostre que un∼c1λ1nu_n \sim c_1\lambda_1^n e un+1/un→λ1u_{n+1}/u_n \to \lambda_1.
  2. (Pell) Defina an+1=an+2bna_{n+1} = a_n + 2b_n, bn+1=an+bnb_{n+1} = a_n + b_n, a0=b0=1a_0 = b_0 = 1. Mostre que q(a,b)=a2−2b2q(a, b) = a^2 - 2b^2 satisfaz q(an+1,bn+1)=−q(an,bn)q(a_{n+1}, b_{n+1}) = -q(a_n, b_n), logo an2−2bn2=(−1)n+1a_n^2 - 2b_n^2 = (-1)^{n+1}; relacione isso com o determinante de M=(1211)M = \left(\begin{smallmatrix}1 & 2\\ 1 & 1\end{smallmatrix}\right).
  3. Deduza a estimativa de erro

    ∣anbn−2∣=1bn (an+2 bn)≤12bn2,\Bigl|\frac{a_n}{b_n} - \sqrt2\Bigr| = \frac{1}{b_n\,(a_n + \sqrt2\,b_n)} \leq \frac{1}{2b_n^2},

    e mostre que ela decai geometricamente com razão 3−223 - 2\sqrt2 (ache os autovalores de MM e o crescimento de bnb_n).

  4. (Crescimento geral) A partir da questão 9, prove: (a) se toda raiz satisfaz ∣λi∣≤ρ\abs{\lambda_i} \leq \rho, então ∣un∣≤C nm−1ρn\abs{u_n} \leq C\,n^{m-1}\rho^n com m=max⁡imim = \max_i m_i; (b) se existe uma única raiz λ1\lambda_1 de módulo máximo e Q1≠0Q_1 \neq 0, então un+1/un→λ1u_{n+1}/u_n \to \lambda_1 — confira isso na solução da questão 10.

Parte IV — Contando passeios e palavras. Para um grafo finito com conjunto de vértices {1,…,N}\{1, \dots, N\}, a matriz de adjacência AA tem Aij=1A_{ij} = 1 se ijij é uma aresta, e 00 caso contrário.

  1. Prove que (An)ij(A^n)_{ij} é o número de passeios de comprimento nn de ii até jj (sequências de nn arestas, cada passo ao longo de uma aresta).
  2. (O triângulo) Para o grafo completo com 33 vértices, A=J−IA = J - I: usando o espectro de JJ (Exemplo 2.19), mostre que

    (An)ii=2n+2(−1)n3,(An)ij=2n−(−1)n3(i≠j),(A^n)_{ii} = \frac{2^n + 2(-1)^n}{3}, \qquad (A^n)_{ij} = \frac{2^n - (-1)^n}{3} \quad (i \neq j),

    e confira os dois em n=2n = 2 listando os passeios.

  3. (Palavras sem 1111) Seja wnw_n o número de palavras binárias de comprimento nn sem dois 11 consecutivos. Codifique as palavras por sua última letra para obter uma matriz de transferência, mostre que wn+2=wn+1+wnw_{n+2} = w_{n+1} + w_n, deduza wn=Fn+2w_n = F_{n+2} (Fibonacci, Exercício 3.5) e dê a taxa de crescimento lim⁡wn+1/wn\lim w_{n+1}/w_n.
  4. (O caminho) Para o grafo caminho 1−2−31 - 2 - 3, mostre que os autovalores de AA são 2,0,−2\sqrt2, 0, -\sqrt2 com autovetores (1,±2,1)(1, \pm\sqrt2, 1) e (1,0,−1)(1, 0, -1), e deduza que o número de passeios de comprimento nn de uma ponta à outra é ((2)n+(−2)n)/4\bigl((\sqrt2)^n + (-\sqrt2)^n\bigr)/4: zero para nn ímpar, e 2 n/2−12^{\,n/2 - 1} para nn par. Confira em n=4n = 4.
  5. (Fórmula do traço) Mostre que o número total de passeios fechados de comprimento nn (todos os pontos de partida) é tr⁡(An)=∑iλin\operatorname{tr}(A^n) = \sum_i \lambda_i^n, e verifique-o no triângulo.

Parte V — Um anel de sequências: diagonalização simultânea. Fixe k≥3k \geq 3, seja ω=e2iπ/k\omega = \eu^{2\iu\pi/k} e seja W∈Mk(C)W \in \mathcal{M}_k(\C) o deslocamento cíclico: Wei=ei+1W e_i = e_{i+1} (índices módulo kk, colunas indexadas por 0,…,k−10, \dots, k-1).

  1. Mostre que WTW^{\mathsf T} é a matriz companheira de Xk−1X^k - 1, deduza χW=μW=Xk−1\chi_W = \mu_W = X^k - 1 e que WW é diagonalizável com os kk autovalores simples ωj\omega^j e autovetores fj=(1,ω−j,ω−2j,…,ω−(k−1)j)Tf_j = (1, \omega^{-j}, \omega^{-2j}, \dots, \omega^{-(k-1)j})^{\mathsf T}.
  2. Uma matriz circulante é C=c0I+c1W+⋯+ck−1Wk−1C = c_0 I + c_1 W + \dots + c_{k-1}W^{k-1}. Mostre que todas as circulantes comutam, que a base (f0,…,fk−1)(f_0, \dots, f_{k-1}) diagonaliza todas elas simultaneamente, e que os autovalores de CC são c^(ωj)=∑mcmωjm\widehat c(\omega^j) = \sum_m c_m \omega^{jm}, j=0,…,k−1j = 0, \dots, k-1.
  3. Deduza det⁡C=∏j=0k−1c^(ωj)\det C = \prod_{j=0}^{k-1} \widehat c(\omega^j) e verifique que k=3k = 3 recupera a fatoração de Exercício 2.8.
  4. (A média do colar) Seja x(n+1)=Mx(n)x^{(n+1)} = Mx^{(n)} com M=12(W+W−1)M = \frac12(W + W^{-1}): cada um de kk números dispostos em anel é substituído pela média de seus dois vizinhos. Mostre que os autovalores de MM são cos⁡(2πj/k)\cos(2\pi j/k) e que o coeficiente de x(0)x^{(0)} em f0f_0 é a média 1k∑mxm(0)\frac1k\sum_m x^{(0)}_m (some as coordenadas dos fjf_j).
  5. Conclua: para kk ímpar, x(n)x^{(n)} converge para o vetor constante cujo valor é a média dos valores iniciais; para k=4k = 4, exiba o autovalor responsável pela não convergência e a obstrução exata (um coeficiente de média alternada que deve se anular).
  6. (Síntese) Em uma frase cada: como a matriz companheira converte a análise de (R)(\mathcal R) em redução; onde o lema da decomposição em núcleos não precisou de dimensão finita; por que os autovalores dominantes governam as taxas de crescimento e o erro diofantino; por que as potências da matriz de adjacência contam passeios; e o que se ganha com matrizes que comutam. Nomeie os dois cumes: o teorema fundamental das recorrências lineares e — para as matrizes positivas da Parte IV, no volume do terceiro ano de graduação — o teorema de Perron–Frobenius.
Solução

Solução de Problema 3.1.

1. As k−1k - 1 primeiras coordenadas de CvnCv_n são un+1,…,un+k−1u_{n+1}, \dots, u_{n+k-1} (a superdiagonal desloca), e a última é a0un+⋯+ak−1un+k−1a_0 u_n + \dots + a_{k-1}u_{n+k-1}. Assim vn+1=Cvnv_{n+1} = Cv_n vale para todo nn se e somente se as últimas coordenadas coincidem para todo nn, isto é, se e somente se (R)(\mathcal R) vale. Iterando, vn=Cnv0v_n = C^nv_0.

2. Desenvolva Dk(X)=det⁡(XIk−C)D_k(X) = \det(XI_k - C) ao longo da primeira coluna: as duas entradas não nulas são XX (posição (1,1)(1,1)) e −a0-a_0 (posição (k,1)(k,1)). O primeiro menor tem a forma Dk−1D_{k-1} para os coeficientes a1,…,ak−1a_1, \dots, a_{k-1}; o segundo menor é triangular superior com diagonal −1-1: determinante (−1)k−1(-1)^{k-1}, com sinal (−1)k+1(-1)^{k+1} vindo da posição. A indução sobre kk (base k=1k = 1: X−a0X - a_0) dá

Dk(X)=X(Xk−1−ak−1Xk−2−⋯−a1)−a0=P(X).D_k(X) = X\bigl(X^{k-1} - a_{k-1}X^{k-2} - \dots - a_1\bigr) - a_0 = P(X).

Para μC\mu_C: como Q(CT)=Q(C)TQ(C^{\mathsf T}) = Q(C)^{\mathsf T} para qualquer polinômio, CC e CTC^{\mathsf T} têm os mesmos anuladores, logo o mesmo polinômio minimal. Para CTC^{\mathsf T}: as colunas se leem CTe1=e2C^{\mathsf T}e_1 = e_2, …, CTek−1=ekC^{\mathsf T}e_{k-1} = e_k, de modo que (e1,CTe1,…,(CT)k−1e1)(e_1, C^{\mathsf T}e_1, \dots, (C^{\mathsf T})^{k-1}e_1) é a base canônica: livre. Um polinômio Q≠0Q \neq 0 de grau <k< k tem então Q(CT)e1≠0Q(C^{\mathsf T})e_1 \neq 0 (é uma combinação não trivial de vetores da base): deg⁡μ≥k\deg\mu \geq k. Como μ∣χ=P\mu \mid \chi = P com deg⁡P=k\deg P = k: μC=P\mu_C = P.

3. Para v=(1,λ,…,λk−1)Tv = (1, \lambda, \dots, \lambda^{k-1})^{\mathsf T}: as linhas 11 a k−1k-1 de CvCv dão λ,λ2,…,λk−1\lambda, \lambda^2, \dots, \lambda^{k-1}, isto é, λ\lambda vezes as primeiras k−1k - 1 entradas de vv; a última linha dá ∑mamλm=λk−P(λ)=λk=λ⋅λk−1\sum_m a_m\lambda^m = \lambda^k - P(\lambda) = \lambda^k = \lambda\cdot\lambda^{k-1}. Logo Cv=λvCv = \lambda v. Reciprocamente, as equações (Cx)i=λxi(Cx)_i = \lambda x_i para i<ki < k leem-se xi+1=λxix_{i+1} = \lambda x_i: todo autovetor é proporcional a vv — todo autoespaço tem dimensão exatamente 11. Diagonalizável se e somente se as dimensões dos autoespaços somam kk (Teorema 3.6), se e somente se há kk autovalores distintos, se e somente se PP tem kk raízes distintas (os autovalores são as raízes de χC=P\chi_C = P).

4. Cada (λin)n(\lambda_i^n)_n resolve (R)(\mathcal R): λin+k=λin λik=λin∑mamλim\lambda_i^{n+k} = \lambda_i^n\,\lambda_i^k = \lambda_i^n\sum_m a_m\lambda_i^m. Liberdade: uma combinação nula ∑iciλin=0\sum_i c_i\lambda_i^n = 0 para n=0,…,k−1n = 0, \dots, k-1 é um sistema de Vandermonde (Exercício 2.11) nos cic_i: todos os ci=0c_i = 0. O espaço das soluções tem dimensão kk (questão 6, cuja demonstração é elementar e independente): kk soluções livres formam uma base, e as coordenadas são únicas.

5. P=X2−X−6=(X−3)(X+2)P = X^2 - X - 6 = (X - 3)(X + 2): solução geral un=A 3n+B(−2)nu_n = A\,3^n + B(-2)^n. Condições iniciais: A+B=1A + B = 1, 3A−2B=83A - 2B = 8: A=2A = 2, B=−1B = -1:

un=2⋅3n−(−2)n.u_n = 2\cdot 3^n - (-2)^n .

(Confira: u2=u1+6u0=14u_2 = u_1 + 6u_0 = 14 e 2⋅9−4=142\cdot9 - 4 = 14.)

6. P(S)((un))P(S)\bigl((u_n)\bigr) é a sequência n↦un+k−ak−1un+k−1−⋯−a0unn \mapsto u_{n+k} - a_{k-1}u_{n+k-1} - \dots - a_0u_n: ela se anula se e somente se (R)(\mathcal R) vale, de modo que o conjunto solução é ker⁡P(S)\ker P(S), um subespaço. A aplicação ker⁡P(S)→Ck\ker P(S) \to \C^k, u↦(u0,…,uk−1)u \mapsto (u_0, \dots, u_{k-1}), é linear, injetiva (a recorrência determina uk,uk+1,…u_{k}, u_{k+1}, \dots a partir dos kk primeiros valores, por indução) e sobrejetiva (defina unu_n recursivamente a partir de dados iniciais quaisquer): dimensão kk.

7. A demonstração do Teorema 3.14 usa apenas: a identidade de Bézout em C[X]\C[X] e o fato de que polinômios num endomorfismo fixo comutam. Nenhum dos dois menciona a dimensão do espaço ambiente: o lema vale palavra por palavra para S∈L(S)S \in \mathcal{L}(\mathcal{S}). Logo

ker⁡P(S)=⨁i=1rker⁡ (S−λi id)mi.\ker P(S) = \bigoplus_{i=1}^{r} \ker\,(S - \lambda_i\,\mathrm{id})^{m_i}.

8. Para Q∈C[X]Q \in \C[X]: (S−λ)(Q(n)λn)n(S - \lambda)\bigl(Q(n)\lambda^n\bigr)_n tem nn-ésimo termo Q(n+1)λn+1−λQ(n)λn=λn+1(ΔQ)(n)Q(n{+}1)\lambda^{n+1} - \lambda Q(n)\lambda^n = \lambda^{n+1}(\Delta Q)(n), com ΔQ=Q(X+1)−Q(X)\Delta Q = Q(X{+}1) - Q(X) de grau deg⁡Q−1\deg Q - 1 (os termos dominantes se cancelam). Iterando, (S−λ)m(Q(n)λn)=(λn+m(ΔmQ)(n))n(S - \lambda)^m\bigl(Q(n)\lambda^n\bigr) = \bigl(\lambda^{n+m}(\Delta^m Q)(n)\bigr)_n, e ΔmQ=0\Delta^m Q = 0 quando deg⁡Q≤m−1\deg Q \leq m - 1: o conjunto da direita está contido no núcleo. Ele é um subespaço de dimensão mm: as sequências (njλn)n(n^j\lambda^n)_n, 0≤j<m0 \leq j < m, são livres, pois ∑jcjnjλn=0\sum_j c_j n^j\lambda^n = 0 para todo nn força (dividindo por λn≠0\lambda^n \neq 0) o polinômio ∑jcjXj\sum_j c_jX^j a se anular em todo n∈Nn \in \N, logo a ser nulo. Reciprocamente, dim⁡ker⁡(S−λ)m≤m\dim\ker(S - \lambda)^m \leq m: desenvolvendo (S−λ)m=∑j(mj)(−λ)m−jSj(S - \lambda)^m = \sum_j \binom mj(-\lambda)^{m-j}S^j, a equação (S−λ)mu=0(S - \lambda)^m u = 0 é uma recorrência linear de ordem mm (coeficiente dominante 11), de modo que uu fica determinada por u0,…,um−1u_0, \dots, u_{m-1} como na questão 6. A igualdade das dimensões conclui.

9. Combine as questões 7 e 8: toda solução se decompõe de modo único como soma de elementos dos ker⁡(S−λi)mi\ker(S - \lambda_i)^{m_i}, isto é, un=∑iQi(n)λinu_n = \sum_i Q_i(n)\lambda_i^n com deg⁡Qi≤mi−1\deg Q_i \leq m_i - 1; os QiQ_i são únicos porque a decomposição é direta e, dentro de cada parcela, os coeficientes de QiQ_i são coordenadas na base (njλin)j(n^j\lambda_i^n)_j (questão 8). Verificação de bom senso nas dimensões: ∑imi=k\sum_i m_i = k.

10. P=X2−4X+4=(X−2)2P = X^2 - 4X + 4 = (X - 2)^2: soluções (a+bn)2n(a + bn)2^n. Dados iniciais: a=1a = 1, 2(a+b)=02(a + b) = 0, logo b=−1b = -1:

un=(1−n) 2n.u_n = (1 - n)\,2^n .

Confira: u2=4u1−4u0=−4u_2 = 4u_1 - 4u_0 = -4 e (1−2)⋅4=−4(1 - 2)\cdot4 = -4.

11. Escreva un=λ1n(c1+∑i≥2ci(λi/λ1)n)u_n = \lambda_1^n\bigl(c_1 + \sum_{i\geq2} c_i(\lambda_i/\lambda_1)^n\bigr); cada razão tem módulo <1< 1, logo o colchete tende a c1≠0c_1 \neq 0: un∼c1λ1nu_n \sim c_1\lambda_1^n. Em particular un≠0u_n \neq 0 para nn grande, e

un+1un=λ1 c1+o(1)c1+o(1)⟶λ1.\frac{u_{n+1}}{u_n} = \lambda_1\,\frac{c_1 + o(1)}{c_1 + o(1)} \longrightarrow \lambda_1 .

12. Calcule:

q(an+1,bn+1)=(an+2bn)2−2(an+bn)2=−an2+2bn2=−q(an,bn).q(a_{n+1}, b_{n+1}) = (a_n + 2b_n)^2 - 2(a_n + b_n)^2 = -a_n^2 + 2b_n^2 = -q(a_n, b_n).

Com q(a0,b0)=1−2=−1q(a_0, b_0) = 1 - 2 = -1: an2−2bn2=(−1)n+1a_n^2 - 2b_n^2 = (-1)^{n+1}. Estruturalmente: q(a,b)=(a−2 b)(a+2 b)q(a, b) = (a - \sqrt2\,b)(a + \sqrt2\,b) e a aplicação linear MM multiplica o fator a+2ba + \sqrt2 b por 1+21 + \sqrt2 e o fator a−2ba - \sqrt2 b por 1−21 - \sqrt2 (calcule: an+1+2bn+1=(1+2)(an+2bn)a_{n+1} + \sqrt2 b_{n+1} = (1 + \sqrt2)(a_n + \sqrt2 b_n)); o produto é multiplicado por (1+2)(1−2)=−1=det⁡M(1 + \sqrt2)(1 - \sqrt2) = -1 = \det M a cada passo.

13. Como an2−2bn2=(an−2bn)(an+2bn)=(−1)n+1a_n^2 - 2b_n^2 = (a_n - \sqrt2 b_n)(a_n + \sqrt2 b_n) = (-1)^{n+1},

∣anbn−2∣=∣an2−2bn2∣bn(an+2bn)=1bn(an+2bn)≤12bn2,\Bigl|\frac{a_n}{b_n} - \sqrt2\Bigr| = \frac{\abs{a_n^2 - 2b_n^2}}{b_n(a_n + \sqrt2 b_n)} = \frac{1}{b_n(a_n + \sqrt2 b_n)} \leq \frac1{2b_n^2},

usando an≥bn≥1a_n \geq b_n \geq 1 (indução: ambos crescem), logo an+2bn≥(1+2)bn≥2bna_n + \sqrt2 b_n \geq (1 + \sqrt2)b_n \geq 2b_n. Autovalores de MM: χM=X2−2X−1\chi_M = X^2 - 2X - 1, de raízes 1±21 \pm \sqrt2; como (a0,b0)(a_0, b_0) tem componente não nula no autovetor dominante (todas as entradas positivas), bn∼c(1+2)nb_n \sim c(1 + \sqrt2)^n com c>0c > 0 (questão 11). Logo o erro é ≍(1+2)−2n=(3+22)−n\asymp (1 + \sqrt2)^{-2n} = (3 + 2\sqrt2)^{-n}: decaimento geométrico de razão 1/(3+22)=3−22≈0.1721/(3 + 2\sqrt2) = 3 - 2\sqrt2 \approx 0.172.

14. (a) Da questão 9: ∣un∣≤∑i∣Qi(n)∣∣λi∣n≤(∑i∣Qi(n)∣)ρn\abs{u_n} \leq \sum_i \abs{Q_i(n)}\abs{\lambda_i}^n \leq \bigl(\sum_i \abs{Q_i(n)}\bigr)\rho^n, e cada ∣Qi(n)∣≤Cinmi−1≤Cinm−1\abs{Q_i(n)} \leq C_i n^{m_i - 1} \leq C_i n^{m-1} para n≥1n \geq 1: some as constantes. (b) Sejam ρ′=max⁡i≥2∣λi∣<∣λ1∣\rho' = \max_{i \geq 2}\abs{\lambda_i} < \abs{\lambda_1} e d=deg⁡Q1d = \deg Q_1, de coeficiente dominante c≠0c \neq 0. Então un=Q1(n)λ1n+Rnu_n = Q_1(n)\lambda_1^n + R_n com ∣Rn∣≤Cnm−1ρ′n\abs{R_n} \leq Cn^{m-1}\rho'^n, e

RnQ1(n)λ1n=O(nm−1−d(ρ′/∣λ1∣)n)⟶0\frac{R_n}{Q_1(n)\lambda_1^n} = O\Bigl(n^{m-1-d} \bigl(\rho'/\abs{\lambda_1}\bigr)^n\Bigr) \longrightarrow 0

(a geométrica vence a polinomial). Logo

un∼Q1(n) λ1n∼c ndλ1n,un+1un⟶λ1(pois Q1(n+1)/Q1(n)→1).u_n \sim Q_1(n)\,\lambda_1^n \sim c\,n^d\lambda_1^n, \qquad \frac{u_{n+1}}{u_n} \longrightarrow \lambda_1 \quad\text{(pois } Q_1(n{+}1)/Q_1(n) \to 1\text{)}.

Confira na questão 10: para un=(1−n)2nu_n = (1-n)2^n a razão é

(−n)2n+1(1−n)2n=2 −n1−n⟶2=λ1.\frac{(-n)2^{n+1}}{(1-n)2^n} = 2\,\frac{-n}{1-n} \longrightarrow 2 = \lambda_1 .

15. Indução sobre nn. Para n=1n = 1, AijA_{ij} conta os passeios de comprimento 11. Passo: um passeio de comprimento n+1n + 1 de ii a jj é um passeio de comprimento nn de ii até algum vértice ℓ\ell seguido de uma aresta ℓj\ell j:

#{passeios}=∑ℓ(An)iℓAℓj=(An+1)ij.\#\{\text{passeios}\} = \sum_{\ell} (A^n)_{i\ell}A_{\ell j} = (A^{n+1})_{ij}.

16. J=3ΠJ = 3\Pi, em que Π=J/3\Pi = J/3 é a projeção sobre Vect⁡(1,1,1)\operatorname{Vect}(1,1,1) ao longo do plano x+y+z=0x + y + z = 0 (Π2=Π\Pi^2 = \Pi, pois J2=3JJ^2 = 3J). Então A=J−I=2Π−(I−Π)A = J - I = 2\Pi - (I - \Pi) e, como Π\Pi e I−ΠI - \Pi são projeções complementares,

An=2n Π+(−1)n(I−Π),i.e.(An)ij=2n3+(−1)n(δij−13),A^n = 2^n\,\Pi + (-1)^n (I - \Pi), \qquad\text{i.e.}\qquad (A^n)_{ij} = \frac{2^n}3 + (-1)^n\Bigl(\delta_{ij} - \frac13\Bigr),

o que dá as duas fórmulas exibidas. Em n=2n = 2: diagonal (4+2)/3=2(4 + 2)/3 = 2 (passeios i→ℓ→ii \to \ell \to i para os dois vizinhos ℓ\ell); fora da diagonal (4−1)/3=1(4 - 1)/3 = 1 (o único passeio i→ℓ→ji \to \ell \to j pelo terceiro vértice).

17. Sejam wn(0),wn(1)w_n^{(0)}, w_n^{(1)} o número de palavras admissíveis de comprimento nn terminadas em 00, resp. em 11. Acrescentando uma letra: um 00 pode seguir qualquer coisa, um 11 só pode seguir um 00:

(wn+1(0)wn+1(1))=(1110)(wn(0)wn(1)).\begin{pmatrix} w_{n+1}^{(0)}\\ w_{n+1}^{(1)}\end{pmatrix} = \begin{pmatrix} 1 & 1\\ 1 & 0\end{pmatrix} \begin{pmatrix} w_n^{(0)}\\ w_n^{(1)}\end{pmatrix}.

Somando, wn+2=wn+1+wnw_{n+2} = w_{n+1} + w_n (ou: condicione na primeira letra). Com w1=2w_1 = 2, w2=3w_2 = 3: wn=Fn+2w_n = F_{n+2} por indução (F3=2F_3 = 2, F4=3F_4 = 3, mesma recorrência). Crescimento: as raízes de X2−X−1X^2 - X - 1 são φ>∣ψ∣\varphi > \abs\psi (Exercício 3.5), e a componente em φ\varphi é não nula (os wnw_n são positivos e ψn→0\psi^n \to 0), logo a questão 11 dá wn+1/wn→φ=1+52w_{n+1}/w_n \to \varphi = \frac{1 + \sqrt5}2.

18. A=(010101010)A = \left(\begin{smallmatrix} 0&1&0\\ 1&0&1\\ 0&1&0\end{smallmatrix}\right). Verificação:

A(1,±2,1)T=(±2,2,±2)T=±2 (1,±2,1)T,A(1,0,−1)T=0:A(1, \pm\sqrt2, 1)^{\mathsf T} = (\pm\sqrt2, 2, \pm\sqrt2)^{\mathsf T} = \pm\sqrt2\,(1, \pm\sqrt2, 1)^{\mathsf T}, \qquad A(1, 0, -1)^{\mathsf T} = 0 :

autovalores 2,−2,0\sqrt2, -\sqrt2, 0 (=2cos⁡π4,2cos⁡3π4,2cos⁡π2= 2\cos\frac\pi4, 2\cos\frac{3\pi}4, 2\cos\frac\pi2). Decomponha e1e_1 na base de autovetores e leia a terceira coordenada, ou use a simetria: com v±=(1,±2,1)v_\pm = (1, \pm\sqrt2, 1), v0=(1,0,−1)v_0 = (1, 0, -1), verifica-se e1=14v++14v−+12v0e_1 = \frac14 v_+ + \frac14 v_- + \frac12 v_0, de modo que, para n≥1n \geq 1,

(An)13=(14(2)nv++14(−2)nv−+0) ⁣3=(2)n+(−2)n4,(A^n)_{13} = \Bigl(\tfrac14(\sqrt2)^n v_+ + \tfrac14(-\sqrt2)^n v_- + 0\Bigr)_{\!3} = \frac{(\sqrt2)^n + (-\sqrt2)^n}{4},

zero para nn ímpar (grafo bipartido: as pontas estão a distância par), e 2⋅2n/2/4=2n/2−12\cdot 2^{n/2}/4 = 2^{n/2 - 1} para nn par. Em n=4n = 4: 21=22^{1} = 2, o que corresponde aos dois passeios 1 2 1 2 31\,2\,1\,2\,3 e 1 2 3 2 31\,2\,3\,2\,3.

19. Os passeios fechados de comprimento nn a partir de ii são (An)ii(A^n)_{ii}; somando sobre ii obtém-se tr⁡(An)\operatorname{tr}(A^n). Triangularizando AA (sobre C\C), AnA^n é triangular com diagonal λin\lambda_i^n: tr⁡(An)=∑iλin\operatorname{tr}(A^n) = \sum_i\lambda_i^n. Triângulo: tr⁡(An)=3 2n+2(−1)n3=2n+2(−1)n=2n+(−1)n+(−1)n\operatorname{tr}(A^n) = 3\,\frac{2^n + 2(-1)^n}3 = 2^n + 2(-1)^n = 2^n + (-1)^n + (-1)^n: o espectro {2,−1,−1}\{2, -1, -1\}, coerente com a questão 16.

20. As colunas de WTW^{\mathsf T}: WTei=ei−1W^{\mathsf T}e_i = e_{i-1} para i≥1i \geq 1 e WTe0=ek−1W^{\mathsf T}e_0 = e_{k-1}; renomeando na ordem e0,e1,…e_0, e_1, \dots isso é exatamente a matriz companheira de Xk−1X^k - 1 (a0=1a_0 = 1, os demais am=0a_m = 0). Questão 2: χW=χWT=Xk−1=μW\chi_{W} = \chi_{W^{\mathsf T}} = X^k - 1 = \mu_{W}. As raízes ωj\omega^j (j=0,…,k−1j = 0, \dots, k-1) são as kk raízes kk-ésimas distintas da unidade: WW é diagonalizável (questão 3, ou Exercício 3.8: Wk=IW^k = I). Autovetores: Wfj=∑mω−jmem+1=∑m′ω−j(m′−1)em′=ωjfjWf_j = \sum_m \omega^{-jm}e_{m+1} = \sum_{m'}\omega^{-j(m'-1)}e_{m'} = \omega^j f_j.

21. As circulantes são polinômios em WW, e polinômios numa matriz fixa comutam entre si. Cada fjf_j é autovetor de toda potência: Wmfj=ωjmfjW^m f_j = \omega^{jm}f_j, logo

Cfj=∑mcmωjmfj=c^(ωj) fj:Cf_j = \sum_m c_m\omega^{jm} f_j = \widehat c(\omega^j)\,f_j :

a base (f0,…,fk−1)(f_0, \dots, f_{k-1}) (livre: Vandermonde nos ω−j\omega^{-j} distintos, Exercício 2.11) diagonaliza todas as circulantes de uma só vez, com os autovalores enunciados.

22. O determinante é o produto dos autovalores (diagonalize): det⁡C=∏jc^(ωj)\det C = \prod_{j}\widehat c(\omega^j). Para k=3k = 3, c0=ac_0 = a, c1=bc_1 = b, c2=cc_2 = c e ω=j=e2iπ/3\omega = j = \eu^{2\iu\pi/3}:

det⁡C=(a+b+c)(a+bj+cj2)(a+bj2+cj4),\det C = (a + b + c)(a + bj + cj^2)(a + bj^2 + cj^4),

e j4=jj^4 = j: exatamente a fatoração de Exercício 2.8.

23. M=12(W+W−1)M = \frac12(W + W^{-1}) é uma circulante (W−1=Wk−1W^{-1} = W^{k-1}), com autovalores 12(ωj+ω−j)=cos⁡2πjk\frac12(\omega^j + \omega^{-j}) = \cos\frac{2\pi j}k na mesma base fjf_j. Coordenadas: escreva x(0)=∑jαjfjx^{(0)} = \sum_j \alpha_j f_j. As coordenadas de fjf_j somam ∑mω−jm\sum_m \omega^{-jm}, o que vale kk para j=0j = 0 e 00 caso contrário (soma geométrica de razão ω−j≠1\omega^{-j} \neq 1). Somando as coordenadas de x(0)x^{(0)}: ∑mxm(0)=α0 k\sum_m x^{(0)}_m = \alpha_0\,k, logo α0=1k∑mxm(0)\alpha_0 = \frac1k\sum_m x^{(0)}_m, a média.

24. x(n)=Mnx(0)=∑jαjcos⁡n(2πjk)fjx^{(n)} = M^nx^{(0)} = \sum_j \alpha_j\cos^n\bigl(\tfrac{2\pi j}k\bigr)f_j. Para kk ímpar, ∣cos⁡(2πj/k)∣<1\abs{\cos(2\pi j/k)} < 1 para todo j≠0j \neq 0 (o ângulo nunca vale 00 nem π\pi), de modo que todos os termos exceto j=0j = 0 tendem a 00: x(n)→α0f0x^{(n)} \to \alpha_0 f_0, o vetor constante igual à média — promediar num anel ímpar iguala tudo. Para k=4k = 4 os autovalores são 1,0,−1,01, 0, -1, 0: o termo j=2j = 2 α2(−1)nf2\alpha_2(-1)^nf_2 com f2=(1,−1,1,−1)Tf_2 = (1, -1, 1, -1)^{\mathsf T} oscila para sempre. A obstrução é a média alternada: multiplicando as coordenadas de x(0)x^{(0)} por (−1)m(-1)^m e somando, o mesmo cálculo de soma geométrica dá ∑m(−1)mxm(0)=4α2\sum_m (-1)^mx^{(0)}_m = 4\alpha_2: o processo converge se e somente se x0(0)−x1(0)+x2(0)−x3(0)=0x^{(0)}_0 - x^{(0)}_1 + x^{(0)}_2 - x^{(0)}_3 = 0, e então converge para a média.

25. A matriz companheira converte uma recorrência escalar de ordem kk numa recorrência vetorial de primeira ordem, de modo que as fórmulas fechadas viram enunciados sobre CnC^n — o terreno de casa da redução (questões 1–5). O lema da decomposição em núcleos é álgebra polinomial pura (Bézout mais comutação), logo ele reparte ker⁡P(S)\ker P(S) mesmo com S\mathcal{S} de dimensão infinita (questões 7–9). Os autovalores dominantes governam o crescimento porque toda outra contribuição é geometricamente desprezível após a normalização — e é também por isso que o erro de Pell decai como o quadrado da raiz dominante (questões 11–14). As potências da matriz de adjacência contam passeios porque a multiplicação de matrizes soma sobre os vértices intermediários, de modo que os espectros contam passeios fechados (questões 15–19). Matrizes que comutam compartilham uma base de autovetores, e uma única base de Fourier diagonaliza toda a álgebra das circulantes de um golpe só (questões 20–24). Cumes: o teorema fundamental das recorrências lineares (questão 9); e, para matrizes não negativas, a razão pela qual raízes dominantes como φ\varphi ou 1+21 + \sqrt2 são automaticamente reais, positivas e simples é o teorema de Perron–Frobenius, demonstrado no volume do terceiro ano de graduação.

Termos definidos neste capítulo

Ver todos os 395 termos do glossário