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 uL(E)u \in \mathcal{L}(E), n=dimEn = \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 x0x \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 uFu|_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 xiEλix_i \in E_{\lambda_i}, sendo o enunciado conhecido para r1r - 1. Aplique uu e subtraia λr\lambda_r vezes a relação:

i=1r1(λ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(Xidu)\chi_u(X) = \det(X\,\mathrm{id} - u) — calculado em qualquer base como det(XInA)\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(tru)Xn1++(1)ndetu.\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 é dimEλ\dim E_\lambda, e 1dimEλmλ1 \leq \dim E_\lambda \leq m_\lambda.

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

Geométrica \leq algébrica: seja d=dimEλ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 (X2)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 (X2X - 2 contra (X2)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=λSpuEλE = \bigoplus_{\lambda \in \operatorname{Sp} u} E_\lambda;
  3. χu\chi_u se decompõe sobre KK e dimEλ=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λ)dimEλ\chi_u = \prod (X - \lambda)^{\dim E_\lambda} se decompõe com multiplicidades correspondentes. Reciprocamente, suponha que χu\chi_u se decomponha com dimEλ=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

λdimEλ=λ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 dimker(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Π)=4k13J+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: 413J+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 PDP1PDP^{-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 é (Xtii)\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 (n1)×(n1)(n-1) \times (n-1) BB, existe uma QQ invertível com Q1BQQ^{-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=(3111)B = \begin{pmatrix}3 & -1\\ 1 & 1\end{pmatrix}: χB=X24X+4=(X2)2\chi_B = X^2 - 4X + 4 = (X - 2)^2, e ker(B2I)=ker(1111)\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)=1e1+2e2,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 e2e_2', e reescalar e2e_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=akXkK[X]P = \sum a_k X^k \in K[X], ponha P(u)=akukL(E)P(u) = \sum a_k u^k \in \mathcal{L}(E). A aplicação PP(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    μuPP(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, μuFμu\mu_{u|_F} \mid \mu_u.

Demonstração. (1) A divisibilidade é a definição de gerador. Se u(x)=λxu(x) = \lambda x, x0x \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(uF)=μu(u)F=0\mu_u(u|_F) = \mu_u(u)|_F = 0, e aplique (1) a uFu|_F.

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

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

μJ=X23X=X(X3):\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=X21\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=(trJ)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=P1P2PrP = P_1 P_2 \cdots P_r com os PiP_i dois a dois coprimos, então

kerP(u)=kerP1(u)kerPr(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 xkerP(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 x2kerP2(u)x_2 \in \ker P_2(u) e, simetricamente, x1kerP1(u)x_1 \in \ker P_1(u): a soma preenche kerP(u)\ker P(u); as duas parcelas ficam dentro de kerP(u)\ker P(u) (PiPP_i \mid P). Diretividade: xkerP1(u)kerP2(u)x \in \ker P_1(u) \cap \ker P_2(u)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 X21=(X1)(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(X1)=1.\frac{1}{2}(X + 1) - \frac12(X - 1) = 1 .

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

π+=A+I2=12(110110002),π=IA2=12(110110000).\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π++etπ=(coshtsinht0sinhtcosht000et).\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=kerP(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 λSpu(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(X1)X(X-1), que se decompõe com raízes simples — diagonalizáveis com espectro {0,1}\subseteq \{0, 1\}, e E=kerpker(pid)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 X21X^2 - 1): diagonalizáveis quando charK2\operatorname{char} K \neq 2, espectro {±1}\subseteq \{\pm 1\}. Um endomorfismo com u3=u2u^3 = u^2 e u2uu^2 \neq u: anulado por X2(X1)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 keru2=keru\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=(010100001),χR=(X1)(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 e1ie2e_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=(X2)2(X3)\chi_D = (X - 2)^2(X - 3) mas μD=(X2)(X3)\mu_D = (X - 2)(X - 3), pois (D2I)(D3I)=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(X3)\chi_{N'} = X^2(X - 3) e μN=X2(X3)\mu_{N'} = X^2(X - 3) — a raiz dupla é genuinamente necessária, pois NN' não é diagonalizável do lado ker\ker (Ne2=e10N'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μun\deg \mu_u \leq n.

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

ud(x)=a0xa1u(x)ad1ud1(x),u^d(x) = -a_0 x - a_1 u(x) - \dots - a_{d-1}u^{d-1}(x),

e ponha Px=Xd+ad1Xd1++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 (C0D)\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(XIC)\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=X25X2\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(A5I)=2IA(A - 5I) = 2I

A1=12(A5I)=(213/21/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,,An1)\operatorname{Vect}(I, A, \dots, A^{n-1})dimK[A]=degμAn\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 é SpA+SpB\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 χ=X22cosθX+1\chi = X^2 - 2\cos\theta\,X + 1, que se decompõe sobre C\C, não sobre R\Rdiagonalizá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 dimEλ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ννdd'\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 jij \leq i. Mostramos por indução que

(Ak)ij=0sempre que ji+k1,(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)iaj,(A^{k+1})_{ij} = \sum_{\ell} (A^k)_{i\ell}\,a_{\ell j} ,

e cada termo se anula: ou i+k1\ell \leq i + k - 1 (o primeiro fator é 00 por indução), ou i+k\ell \geq i + k, e nesse caso ji+kj \leq i + k \leq \ell mata o segundo fator. Em k=nk = n a condição ji+n1j \leq i + n - 1 vale para todos i,jni, 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=N1Nr,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 ν=ud\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í: νmaxmi=0\nu^{\max m_i} = 0 em cada parcela, logo ν\nu é nilpotente.

Unicidade. Seja u=d+νu = d' + \nu' outro par desse tipo. Como dd' 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 ddd - 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+k1=j=0k+k1(k+k1j)νj(ν)k+k1j,(\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 jkj \geq k' (primeiro fator nulo), ou k+k1jkk + 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=dd = d', ν=ν\nu = \nu'.

Exemplo 3.27 (Potências e exponenciais)

A=(3111)A = \begin{pmatrix} 3 & 1\\ -1 & 1\end{pmatrix}: χA=X24X+4=(X2)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=A2I=(1111)N = A - 2I = \begin{pmatrix} 1 & 1\\ -1 & -1\end{pmatrix}, N2=0N^2 = 0. Então

Ak=(2I+N)k=2kI+k2k1N,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=X22X3=(X3)(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=(1111)P = \begin{pmatrix} 1 & 1\\ 1 & -1\end{pmatrix}P1AP=diag(3,1)P^{-1}AP = \operatorname{diag}(3, -1).

B=JIB = 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á P1BP=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=(X1)2\chi_C = (X-1)^2, único autovalor 11; ker(CI)=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 (X1)2(X-1)^2 e CIC \neq I, logo μC=(X1)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 u25u+6id=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.

X25X+6=(X2)(X3)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=2idu = 2\,\mathrm{id}), {3}\{3\} (u=3idu = 3\,\mathrm{id}) ou {2,3}\{2, 3\}.

Potências: procure uk=akid+bkuu^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=3k2kb_k = 3^k - 2^k, ak=32k23ka_k = 3\cdot2^k - 2\cdot 3^k:

uk=(32k23k)id+(3k2k)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 uFu|_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(uF)=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+Fn1F_{n+1} = F_n + F_{n-1}):

Fn=φnψn5,φ=1+52, ψ=152.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=X2X1\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φψ5=1\frac{\varphi - \psi}{\sqrt5} = 1.)

Exercício 3.6 ★★

Sejam uL(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 (detu2=(detu)20\det u^2 = (\det u)^2 \neq 0), logo todos os μi0\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μiid)=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+k2k1N+(k2)2k2N2=(2kk2k1(k2)2k202kk2k1002k),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 AMn(C)A \in \mathcal{M}_n(\C) com Ak=IA^k = I para algum k1k \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.

Xk1X^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 Xk1X^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, é PIP1=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 xEλ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 uL(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 uGu|_G produziria um autovetor de uu dentro de GG — contradizendo GF={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 uFu|_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=λ(FEλ)F = \bigoplus_\lambda (F \cap E_\lambda), pois uFu|_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\}: χuG\chi_{u|_G} se decompõe sobre C\C, logo uGu|_G tem um autovetor xGx \in G (Teorema 3.9 ou diretamente a existência de uma raiz); mas todo autovetor de uu está em FF, logo xFG={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 AM2(C)A \in \mathcal{M}_2(\C) com os dois autovalores de módulo <1< 1. Prove que Ak0A^k \to 0 entrada a entrada quando kk \to \infty. (Triangularize: A=P(T)P1A = 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=PTP1A = PTP^{-1}, T=(λc0μ)T = \begin{pmatrix} \lambda & c\\ 0 & \mu\end{pmatrix}, λ,μ<1\abs\lambda, \abs\mu < 1. Então Ak=PTkP1A^k = PT^kP^{-1}, e basta que Tk0T^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,μk0\abs{\lambda}^k, \abs\mu^k \to 0).

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

Exercício 3.12 ★★

Seja uL(Cn)u \in \mathcal{L}(\C^n) com rku=1\operatorname{rk} u = 1 (n2n \geq 2). Mostre que χu=Xn1(Xtru)\chi_u = X^{n-1}(X - \operatorname{tr} u) e que uu é diagonalizável se e somente se tru0\operatorname{tr} u \neq 0. (Lembre do Exercício 2.5 que u2=(tru)uu^2 = (\operatorname{tr} u)\,u.)

Solução

Solução de Exercício 3.12.

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

Se tru0\operatorname{tr} u \neq 0: o autovalor tru\operatorname{tr} u é raiz de χu\chi_u, logo carrega um autovetor; os autoespaços de 00 e de tru\operatorname{tr} u têm dimensões n1n - 1 e 1\geq 1, com soma n\geq n: eles preenchem EE, e uu é diagonalizável (Teorema 3.6). Se tru=0\operatorname{tr} u = 0: pelo Exercício 2.5, u2=(tru)u=0u^2 = (\operatorname{tr} u)u = 0 com u0u \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=ak1un+k1++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 k1k \geq 1, escalares a0,,ak1Ca_0, \dots, a_{k-1} \in \C com a00a_0 \neq 0, o polinômio mônico P=Xkak1Xk1a1Xa0P = X^k - a_{k-1}X^{k-1} - \dots - a_1 X - a_0 e a recorrência

(R) ⁣:un+k=ak1un+k1++a1un+1+a0un(n0).(\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=(0101a0a1ak1)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+k1)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(XIC)\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,λ,,λk1)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 SL(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) é kerP(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 kerP(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 a00a_0 \neq 0).
  3. Para λ0\lambda \neq 0 e m1m \geq 1, mostre que

    ker(Sλid)m={(Q(n)λn)n:QCm1[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,QiCmi1[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+14unu_{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 i2i \geq 2, e un=iciλinu_n = \sum_i c_i \lambda_i^n com c10c_1 \neq 0. Mostre que unc1λ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)=a22b2q(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 an22bn2=(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

    anbn2=1bn(an+2bn)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 3223 - 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 unCnm1ρn\abs{u_n} \leq C\,n^{m-1}\rho^n com m=maximim = \max_i m_i; (b) se existe uma única raiz λ1\lambda_1 de módulo máximo e Q10Q_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=JIA = J - I: usando o espectro de JJ (Exemplo 2.19), mostre que

    (An)ii=2n+2(1)n3,(An)ij=2n(1)n3(ij),(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 limwn+1/wn\lim w_{n+1}/w_n.
  4. (O caminho) Para o grafo caminho 1231 - 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 2n/212^{\,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 k3k \geq 3, seja ω=e2iπ/k\omega = \eu^{2\iu\pi/k} e seja WMk(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,,k10, \dots, k-1).

  1. Mostre que WTW^{\mathsf T} é a matriz companheira de Xk1X^k - 1, deduza χW=μW=Xk1\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,,ω(k1)j)Tf_j = (1, \omega^{-j}, \omega^{-2j}, \dots, \omega^{-(k-1)j})^{\mathsf T}.
  2. Uma matriz circulante é C=c0I+c1W++ck1Wk1C = c_0 I + c_1 W + \dots + c_{k-1}W^{k-1}. Mostre que todas as circulantes comutam, que a base (f0,,fk1)(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,,k1j = 0, \dots, k-1.
  3. Deduza detC=j=0k1c^(ω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+W1)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 1kmxm(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 k1k - 1 primeiras coordenadas de CvnCv_n são un+1,,un+k1u_{n+1}, \dots, u_{n+k-1} (a superdiagonal desloca), e a última é a0un++ak1un+k1a_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(XIkC)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 Dk1D_{k-1} para os coeficientes a1,,ak1a_1, \dots, a_{k-1}; o segundo menor é triangular superior com diagonal 1-1: determinante (1)k1(-1)^{k-1}, com sinal (1)k+1(-1)^{k+1} vindo da posição. A indução sobre kk (base k=1k = 1: Xa0X - a_0) dá

Dk(X)=X(Xk1ak1Xk2a1)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, …, CTek1=ekC^{\mathsf T}e_{k-1} = e_k, de modo que (e1,CTe1,,(CT)k1e1)(e_1, C^{\mathsf T}e_1, \dots, (C^{\mathsf T})^{k-1}e_1) é a base canônica: livre. Um polinômio Q0Q \neq 0 de grau <k< k tem então Q(CT)e10Q(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 degP=k\deg P = k: μC=P\mu_C = P.

3. Para v=(1,λ,,λk1)Tv = (1, \lambda, \dots, \lambda^{k-1})^{\mathsf T}: as linhas 11 a k1k-1 de CvCv dão λ,λ2,,λk1\lambda, \lambda^2, \dots, \lambda^{k-1}, isto é, λ\lambda vezes as primeiras k1k - 1 entradas de vv; a última linha dá mamλm=λkP(λ)=λk=λλk1\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=λinmamλ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,,k1n = 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=X2X6=(X3)(X+2)P = X^2 - X - 6 = (X - 3)(X + 2): solução geral un=A3n+B(2)nu_n = A\,3^n + B(-2)^n. Condições iniciais: A+B=1A + B = 1, 3A2B=83A - 2B = 8: A=2A = 2, B=1B = -1:

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

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

6. P(S)((un))P(S)\bigl((u_n)\bigr) é a sequência nun+kak1un+k1a0unn \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 é kerP(S)\ker P(S), um subespaço. A aplicação kerP(S)Ck\ker P(S) \to \C^k, u(u0,,uk1)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 SL(S)S \in \mathcal{L}(\mathcal{S}). Logo

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

8. Para QC[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 degQ1\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 degQm1\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, 0j<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 λn0\lambda^n \neq 0) o polinômio jcjXj\sum_j c_jX^j a se anular em todo nNn \in \N, logo a ser nulo. Reciprocamente, dimker(Sλ)mm\dim\ker(S - \lambda)^m \leq m: desenvolvendo (Sλ)m=j(mj)(λ)mjSj(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,,um1u_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 degQimi1\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=X24X+4=(X2)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=(1n)2n.u_n = (1 - n)\,2^n .

Confira: u2=4u14u0=4u_2 = 4u_1 - 4u_0 = -4 e (12)4=4(1 - 2)\cdot4 = -4.

11. Escreva un=λ1n(c1+i2ci(λ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 c10c_1 \neq 0: unc1λ1nu_n \sim c_1\lambda_1^n. Em particular un0u_n \neq 0 para nn grande, e

un+1un=λ1c1+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)22(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)=12=1q(a_0, b_0) = 1 - 2 = -1: an22bn2=(1)n+1a_n^2 - 2b_n^2 = (-1)^{n+1}. Estruturalmente: q(a,b)=(a2b)(a+2b)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 a2ba - \sqrt2 b por 121 - \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)(12)=1=detM(1 + \sqrt2)(1 - \sqrt2) = -1 = \det M a cada passo.

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

anbn2=an22bn2bn(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 anbn1a_n \geq b_n \geq 1 (indução: ambos crescem), logo an+2bn(1+2)bn2bna_n + \sqrt2 b_n \geq (1 + \sqrt2)b_n \geq 2b_n. Autovalores de MM: χM=X22X1\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), bnc(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)=3220.1721/(3 + 2\sqrt2) = 3 - 2\sqrt2 \approx 0.172.

14. (a) Da questão 9: uniQi(n)λin(iQi(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)Cinmi1Cinm1\abs{Q_i(n)} \leq C_i n^{m_i - 1} \leq C_i n^{m-1} para n1n \geq 1: some as constantes. (b) Sejam ρ=maxi2λi<λ1\rho' = \max_{i \geq 2}\abs{\lambda_i} < \abs{\lambda_1} e d=degQ1d = \deg Q_1, de coeficiente dominante c0c \neq 0. Então un=Q1(n)λ1n+Rnu_n = Q_1(n)\lambda_1^n + R_n com RnCnm1ρn\abs{R_n} \leq Cn^{m-1}\rho'^n, e

RnQ1(n)λ1n=O(nm1d(ρ/λ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

unQ1(n)λ1ncndλ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=(1n)2nu_n = (1-n)2^n a razão é

(n)2n+1(1n)2n=2n1n2=λ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)iAj=(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=JI=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(δij13),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 iii \to \ell \to i para os dois vizinhos \ell); fora da diagonal (41)/3=1(4 - 1)/3 = 1 (o único passeio iji \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 X2X1X^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 ψn0\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,2cos3π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 n1n \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 22n/2/4=2n/212\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 121231\,2\,1\,2\,3 e 123231\,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)=32n+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=ei1W^{\mathsf T}e_i = e_{i-1} para i1i \geq 1 e WTe0=ek1W^{\mathsf T}e_0 = e_{k-1}; renomeando na ordem e0,e1,e_0, e_1, \dots isso é exatamente a matriz companheira de Xk1X^k - 1 (a0=1a_0 = 1, os demais am=0a_m = 0). Questão 2: χW=χWT=Xk1=μW\chi_{W} = \chi_{W^{\mathsf T}} = X^k - 1 = \mu_{W}. As raízes ωj\omega^j (j=0,,k1j = 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(m1)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,,fk1)(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): detC=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}:

detC=(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+W1)M = \frac12(W + W^{-1}) é uma circulante (W1=Wk1W^{-1} = W^{k-1}), com autovalores 12(ωj+ωj)=cos2π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 ωj1\omega^{-j} \neq 1). Somando as coordenadas de x(0)x^{(0)}: mxm(0)=α0k\sum_m x^{(0)}_m = \alpha_0\,k, logo α0=1kmxm(0)\alpha_0 = \frac1k\sum_m x^{(0)}_m, a média.

24. x(n)=Mnx(0)=jαjcosn(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 j0j \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 kerP(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.