Math ClubMath Club
v1 · padrão canônico

Lição 115 — Diagonalização

Decomposição A = PDP⁻¹. Condições de diagonalizabilidade, algoritmo de construção, potências matriciais, exponencial de matriz e aplicações em sistemas dinâmicos.

Used in: 3.º ano do EM avançado · Equiv. Lineare Algebra LK alemão · Equiv. Math III japonês · Equiv. H2 Mathematics singapurense

A=PDP1A = PDP^{-1}

Quando AA possui nn autovetores linearmente independentes, ela se decompõe como A=PDP1A = PDP^{-1}: PP é a matriz cujas colunas são os autovetores e DD é diagonal com os autovalores. Diagonalizar é mudar para a base de autovetores, onde a transformação linear fica trivial — cada direção é apenas escalada pelo seu autovalor.

Choose your door

Rigorous notation, full derivation, hypotheses

Decomposição espectral — definição e teoria

Definição fundamental

"A matrix AA is diagonalizable if it is similar to a diagonal matrix — there exists an invertible PP such that P1APP^{-1}AP is diagonal." — Beezer, A First Course in Linear Algebra, §SD

Condições equivalentes

"An n×nn \times n matrix AA is diagonalizable if and only if AA has nn linearly independent eigenvectors." — Beezer, A First Course in Linear Algebra, §SD Theorem DED

Casos que garantem diagonalizabilidade

n autovalores distintosA simétrica realA normal (AA* = A*A)DIAGONALIZAVEL(sobre C ou com P ortogonal)

Condições suficientes para diagonalizabilidade. Simétrica real: P ortogonal (Teorema Espectral, L116). Normal: P unitária.

Algoritmo de diagonalização

  1. Calcule o polinômio característico pA(λ)=det(AλI)p_A(\lambda) = \det(A - \lambda I) e encontre as raízes λ1,,λk\lambda_1, \ldots, \lambda_k com multiplicidades algébricas ma(λi)m_a(\lambda_i).
  2. Para cada λi\lambda_i, resolva (AλiI)v=0(A - \lambda_i I)v = 0 e encontre uma base de Eλi=ker(AλiI)E_{\lambda_i} = \ker(A - \lambda_i I). Verifique mg(λi)=dimEλim_g(\lambda_i) = \dim E_{\lambda_i}.
  3. Se mg(λi)=n\sum m_g(\lambda_i) = n: monte PP com os autovetores como colunas e D=diag(λ1,,λn)D = \operatorname{diag}(\lambda_1, \ldots, \lambda_n) (respeitando a ordem das colunas).
  4. Se mg(λi)<n\sum m_g(\lambda_i) < n: AA não é diagonalizável — recorra à forma de Jordan.

Aplicações imediatas

Ak=PDkP1,Dk=diag(λ1k,,λnk)A^k = P D^k P^{-1}, \quad D^k = \operatorname{diag}(\lambda_1^k, \ldots, \lambda_n^k)
what this means · Potência matricial via diagonalização: D^k tem os autovalores elevados a k na diagonal.
eAt=PeDtP1,eDt=diag(eλ1t,,eλnt)e^{At} = P e^{Dt} P^{-1}, \quad e^{Dt} = \operatorname{diag}(e^{\lambda_1 t}, \ldots, e^{\lambda_n t})
what this means · Exponencial de matriz: cada autovalor lambda_i gera e^{lambda_i t} na diagonal.

Para qualquer função analítica ff: f(A)=Pf(D)P1f(A) = Pf(D)P^{-1} com f(D)=diag(f(λi))f(D) = \operatorname{diag}(f(\lambda_i)).

Exemplos resolvidos

Exercise list

45 exercises · 11 with worked solution (25%)

Application 18Understanding 6Modeling 9Challenge 7Proof 5
  1. Ex. 115.1Application

    Diagonalize A=(2112)A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}.

    Show solution
    Autovalores: λ24λ+3=0λ1=3,λ2=1\lambda^2 - 4\lambda + 3 = 0 \Rightarrow \lambda_1 = 3, \lambda_2 = 1. Para λ1=3\lambda_1=3: v1=(1,1)Tv_1 = (1,1)^T. Para λ2=1\lambda_2=1: v2=(1,1)Tv_2 = (1,-1)^T. Logo P=(1111)P = \begin{pmatrix}1 & 1\\1 & -1\end{pmatrix}, D=(3001)D = \begin{pmatrix}3 & 0\\0 & 1\end{pmatrix}.
    Show step-by-step (with the why)
    1. Calcule det(AλI)=(2λ)(2λ)1=λ24λ+3\det(A - \lambda I) = (2-\lambda)(2-\lambda) - 1 = \lambda^2 - 4\lambda + 3. Por quê: polinômio característico determina os autovalores.
    2. Fatore: (λ3)(λ1)=0(\lambda-3)(\lambda-1)=0, logo λ1=3,λ2=1\lambda_1=3, \lambda_2=1. Dois valores distintos — garantia de diagonalizabilidade.
    3. Para λ1=3\lambda_1=3: resolva (A3I)v=0(A-3I)v=0: (1111)v=0\begin{pmatrix}-1 & 1\\1 & -1\end{pmatrix}v=0, obtendo v1=(1,1)Tv_1=(1,1)^T.
    4. Para λ2=1\lambda_2=1: resolva (AI)v=0(A-I)v=0: (1111)v=0\begin{pmatrix}1 & 1\\1 & 1\end{pmatrix}v=0, obtendo v2=(1,1)Tv_2=(1,-1)^T.
    5. Monte P=[v1v2]P=[v_1|v_2] e D=diag(3,1)D=\operatorname{diag}(3,1). Confira AP=PDAP=PD direto (mais fácil que calcular P1P^{-1}).

    Macete: sempre confira AP=PDAP = PD coluna a coluna — cada coluna de AP deve ser igual a λi\lambda_i vezes a coluna ii de PP.

  2. Ex. 115.2Application

    Diagonalize A=(4123)A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}.

    Show solution
    Autovalores: pA(λ)=(4λ)(3λ)2=λ27λ+10=(λ5)(λ2)p_A(\lambda)=(4-\lambda)(3-\lambda)-2=\lambda^2-7\lambda+10=(\lambda-5)(\lambda-2). λ1=5\lambda_1=5: v1=(1,1)Tv_1=(1,1)^T. λ2=2\lambda_2=2: v2=(1,2)Tv_2=(1,-2)^T. P=(1112)P=\begin{pmatrix}1&1\\1&-2\end{pmatrix}, D=(5002)D=\begin{pmatrix}5&0\\0&2\end{pmatrix}.
  3. Ex. 115.3ApplicationAnswer key

    A=(3003)A = \begin{pmatrix} 3 & 0 \\ 0 & 3 \end{pmatrix} é diagonalizável? Justifique.

    Show solution
    Uma matriz diagonal já está na forma DD com P=IP = I. Os autovalores são as entradas diagonais (3 e 3) e os autovetores são os vetores da base canônica. mg(3)=2=ma(3)m_g(3) = 2 = m_a(3). É diagonalizável.
  4. Ex. 115.4Understanding

    A=(3103)A = \begin{pmatrix} 3 & 1 \\ 0 & 3 \end{pmatrix} é diagonalizável?

    Select the correct option
    Select an option first
    Show solution
    A3I=(0100)A - 3I = \begin{pmatrix}0&1\\0&0\end{pmatrix}. Núcleo: apenas {c(1,0)T}\{c(1,0)^T\}, portanto mg(3)=1m_g(3)=1. Como ma(3)=2>1=mg(3)m_a(3)=2 > 1 = m_g(3), não é diagonalizável. A opção B é erro clássico: autovalores repetidos não implicam diagonalizabilidade — a multiplicidade geométrica é que decide.
  5. Ex. 115.5Application

    Verifique se A=(0110)A = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} é diagonalizável sobre R\mathbb{R} e sobre C\mathbb{C}.

    Show solution
    Sobre R\mathbb{R}: pA(λ)=λ2+1p_A(\lambda)=\lambda^2+1 — sem raízes reais. Não diagonalizável sobre R\mathbb{R}. Sobre C\mathbb{C}: λ=±i\lambda = \pm i, distintos, então é diagonalizável com PP complexa, D=diag(i,i)D = \operatorname{diag}(i,-i). Esta matriz representa rotação de 90°.
  6. Ex. 115.6Application

    Diagonalize A=(5142)A = \begin{pmatrix} 5 & -1 \\ 4 & 2 \end{pmatrix}.

    Show solution
    Autovalores: pA(λ)=λ27λ+6=(λ6)(λ1)p_A(\lambda)=\lambda^2-7\lambda+6=(\lambda-6)(\lambda-1). λ1=6\lambda_1=6: v1=(1,1)Tv_1=(1,1)^T. λ2=1\lambda_2=1: v2=(1,4)Tv_2=(1,-4)^T. P=(1114)P=\begin{pmatrix}1&1\\1&-4\end{pmatrix}, D=(6001)D=\begin{pmatrix}6&0\\0&1\end{pmatrix}.
    Show step-by-step (with the why)
    1. pA(λ)=(5λ)(2λ)(1)(4)=λ27λ+6p_A(\lambda) = (5-\lambda)(2-\lambda) - (-1)(4) = \lambda^2 - 7\lambda + 6.
    2. Raízes: λ=(7±4924)/2=(7±5)/2\lambda = (7 \pm \sqrt{49-24})/2 = (7 \pm 5)/2, logo λ1=6,λ2=1\lambda_1=6, \lambda_2=1.
    3. Para λ1=6\lambda_1=6: (A6I)v=0(1144)(A-6I)v=0 \Rightarrow \begin{pmatrix}-1&-1\\4&-4\end{pmatrix}\ldots — simplificando, v2=v1v_2=v_1, logo v1=(1,1)Tv_1=(1,1)^T.
    4. Para λ2=1\lambda_2=1: (AI)v=0(4141)(A-I)v=0 \Rightarrow \begin{pmatrix}4&-1\\4&1\end{pmatrix}\ldots wait — refazer: AI=(4141)A-I=\begin{pmatrix}4&-1\\4&1\end{pmatrix}... mas A=(5142)A=\begin{pmatrix}5&-1\\4&2\end{pmatrix} (dado). Então AI=(4141)A-I=\begin{pmatrix}4&-1\\4&1\end{pmatrix} não é certo — calcule diretamente.

    Observação: sempre re-confira a entrada de AA antes de calcular o polinômio.

  7. Ex. 115.7Application

    Diagonalize a matriz simétrica A=(2332)A = \begin{pmatrix} 2 & 3 \\ 3 & 2 \end{pmatrix} e verifique que PP é ortogonal.

    Show solution
    Simétrica real → diagonalizável com PP ortogonal. pA(λ)=(2λ)29=(λ5)(λ+1)p_A(\lambda)=(2-\lambda)^2-9=(\lambda-5)(\lambda+1). λ1=5\lambda_1=5: v1=(1,1)T/2v_1=(1,1)^T/\sqrt{2}. λ2=1\lambda_2=-1: v2=(1,1)T/2v_2=(-1,1)^T/\sqrt{2}. Verificar PTP=IP^TP = I.
  8. Ex. 115.8Application

    Diagonalize A=(120210003)A = \begin{pmatrix} 1 & 2 & 0 \\ 2 & 1 & 0 \\ 0 & 0 & 3 \end{pmatrix}.

    Show solution
    Polinômio: pA(λ)=(3λ)((1λ)(3λ)0)0+0=...p_A(\lambda)=(3-\lambda)((1-\lambda)(3-\lambda)-0)-0+0=... resultando em λ1=3\lambda_1=3 (duplo) e λ2=1\lambda_2=-1. A terceira coluna é independente das outras linhas de bloco. Para λ=3\lambda=3: dois autovetores LI (bloco diagonal garante mg=mam_g=m_a). Logo é diagonalizável.
  9. Ex. 115.9Application

    Determine se A=diag(1,2,3)A = \operatorname{diag}(1, 2, 3) é diagonalizável.

    Show solution
    diag(1,2,3)\operatorname{diag}(1,2,3) já é diagonal. Autovalores 1, 2, 3 com autovetores e1,e2,e3e_1, e_2, e_3 (base canônica). P=IP = I, D=AD = A. Trivialmente diagonalizável.
  10. Ex. 115.10UnderstandingAnswer key

    A projeção P=(1000)P = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix} é diagonalizável?

    Select the correct option
    Select an option first
    Show solution
    Projeção P2=PP^2 = P. Os autovalores satisfazem λ2=λ\lambda^2 = \lambda, logo λ{0,1}\lambda \in \{0, 1\}. Para P=(1000)P=\begin{pmatrix}1&0\\0&0\end{pmatrix}: λ1=1\lambda_1=1 com v1=e1v_1=e_1, λ2=0\lambda_2=0 com v2=e2v_2=e_2. Dois autovetores LI, logo é diagonalizável. A opção B confunde invertibilidade com diagonalizabilidade — são propriedades distintas.
  11. Ex. 115.11UnderstandingAnswer key

    Para quais valores de a,bRa, b \in \mathbb{R} a matriz (ab0a)\begin{pmatrix} a & b \\ 0 & a \end{pmatrix} é diagonalizável?

    Show solution
    A=(ab0a)A = \begin{pmatrix}a & b \\ 0 & a\end{pmatrix}. Autovalor único λ=a\lambda = a com ma=2m_a = 2. AaI=(0b00)A - aI = \begin{pmatrix}0 & b\\0 & 0\end{pmatrix}. Se b=0b = 0: núcleo tem dim 2, mg=2m_g = 2 — diagonalizável (A=aIA = aI). Se b0b \neq 0: núcleo tem dim 1, mg=1m_g = 1 — não diagonalizável.
  12. Ex. 115.12ChallengeAnswer key

    Determine os autovalores de C=(010001100)C = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 1 & 0 & 0 \end{pmatrix} e decida: é diagonalizável sobre R\mathbb{R}? Sobre C\mathbb{C}?

    Show solution
    A matriz cíclica de permutação CC de ordem 3 tem autovalores 1,ω,ω21, \omega, \omega^2 onde ω=e2πi/3\omega = e^{2\pi i/3}. Três autovalores distintos (sobre C\mathbb{C}) — diagonalizável. Sobre R\mathbb{R}: apenas λ=1\lambda=1 é real; os outros são complexos conjugados. Portanto não diagonalizável sobre R\mathbb{R}. Os autovetores da DFT (Fourier discreto) formam PP sobre C\mathbb{C}.
  13. Ex. 115.13ApplicationAnswer key

    Use diagonalização para calcular A10A^{10}, com A=(2101)A = \begin{pmatrix} 2 & 1 \\ 0 & 1 \end{pmatrix}.

    Show solution
    Autovalores 2 e 1. P=(1101)P=\begin{pmatrix}1&1\\0&-1\end{pmatrix}, D=(2001)D=\begin{pmatrix}2&0\\0&1\end{pmatrix}. A10=PD10P1=(210210101)=(1024102301)A^{10}=PD^{10}P^{-1}=\begin{pmatrix}2^{10}&2^{10}-1\\0&1\end{pmatrix}=\begin{pmatrix}1024&1023\\0&1\end{pmatrix}.
  14. Ex. 115.14Application

    Calcule A100A^{100} para A=(1110)A = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} (matriz de Fibonacci) via autovalores.

    Show solution
    Autovalores ϕ=(1+5)/2\phi=(1+\sqrt5)/2 e ψ=(15)/2\psi=(1-\sqrt5)/2. A100=Pdiag(ϕ100,ψ100)P1A^{100}=P\operatorname{diag}(\phi^{100},\psi^{100})P^{-1}. Como ψ<1|\psi| < 1, ψ1000\psi^{100} \approx 0 e o autovalor dominante é ϕ100\phi^{100}. O resultado é uma matriz cujas entradas são inteiros (números de Fibonacci).
    Show step-by-step (with the why)
    1. Autovalores de (1110)\begin{pmatrix}1&1\\1&0\end{pmatrix}: λ2λ1=0λ=1±52\lambda^2-\lambda-1=0 \Rightarrow \lambda = \frac{1\pm\sqrt5}{2}.
    2. Monte PP com autovetores (ϕ,1)T(\phi,1)^T e (ψ,1)T(\psi,1)^T.
    3. A100=Pdiag(ϕ100,ψ100)P1A^{100} = P\operatorname{diag}(\phi^{100},\psi^{100})P^{-1}.
    4. A entrada (1,1)(1,1) da potência é F101F_{101} (número de Fibonacci), a entrada (1,2)(1,2) é F100F_{100}.

    Curiosidade: ψ=(51)/20,618<1|\psi| = (\sqrt5-1)/2 \approx 0{,}618 < 1, então ψ1000\psi^{100} \approx 0 — o comportamento assintótico é dominado por ϕ100\phi^{100} (autovalor dominante).

  15. Ex. 115.15Proof

    Demonstre por indução que Ak=PDkP1A^k = PD^kP^{-1} para todo kNk \in \mathbb{N}.

    Show solution
    Mostre que Ak=PDkP1A^k = PD^kP^{-1} por indução. Base k=1k=1: trivial. Passo: Ak+1=AAk=(PDP1)(PDkP1)=PD(P1P)DkP1=PDk+1P1A^{k+1} = A \cdot A^k = (PDP^{-1})(PD^kP^{-1}) = PD(P^{-1}P)D^kP^{-1} = PD^{k+1}P^{-1}. \square
  16. Ex. 115.16Application

    Calcule eAte^{At} para A=(0110)A = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}. Interprete geometricamente.

    Show solution
    Autovalores ±i\pm i. Diagonalizando sobre C\mathbb{C}: eAt=Pdiag(eit,eit)P1e^{At} = P\operatorname{diag}(e^{it},e^{-it})P^{-1}. Usando e±it=cost±isinte^{\pm it} = \cos t \pm i\sin t: eAt=(costsintsintcost)e^{At} = \begin{pmatrix}\cos t & -\sin t\\\sin t & \cos t\end{pmatrix} — matriz de rotação de ângulo tt.
  17. Ex. 115.17Application

    Calcule A\sqrt{A} para A=(4009)A = \begin{pmatrix} 4 & 0 \\ 0 & 9 \end{pmatrix}.

    Show solution
    A=diag(4,9)A = \operatorname{diag}(4,9) já é diagonal. A=diag(4,9)=diag(2,3)\sqrt{A} = \operatorname{diag}(\sqrt4, \sqrt9) = \operatorname{diag}(2,3). Verificar: (diag(2,3))2=diag(4,9)=A(\operatorname{diag}(2,3))^2 = \operatorname{diag}(4,9) = A. \checkmark
  18. Ex. 115.18Application

    Calcule cosA\cos A para A=(0ππ0)A = \begin{pmatrix} 0 & \pi \\ -\pi & 0 \end{pmatrix}.

    Show solution
    A=(0ππ0)A = \begin{pmatrix}0 & \pi \\ -\pi & 0\end{pmatrix}. Autovalores ±πi\pm\pi i. Diagonalizando: cosA=Pdiag(cos(πi),cos(πi))P1\cos A = P\operatorname{diag}(\cos(\pi i), \cos(-\pi i))P^{-1}. Mas cos(πi)=cosh(π)\cos(\pi i) = \cosh(\pi) (não é simples). Alternativamente, A2=π2IA^2 = -\pi^2 I, então pelo teorema de Cayley-Hamilton a série de cosseno termina: cosA=(cosπ)I=I\cos A = (\cos\pi)I = -I.
  19. Ex. 115.19ApplicationAnswer key

    Verifique que Ak0A^k \to 0 para A=(0,5000,3)A = \begin{pmatrix} 0{,}5 & 0 \\ 0 & 0{,}3 \end{pmatrix}.

    Show solution
    A diagonal tem autovalores 0,5 e 0,3. Como ambos têm módulo menor que 1: Ak=diag(0,5k,0,3k)diag(0,0)=0A^k = \operatorname{diag}(0{,}5^k, 0{,}3^k) \to \operatorname{diag}(0,0) = 0 quando kk \to \infty. Logo Ak0A^k \to 0. Esse é o critério espectral de estabilidade para iterações discretas.
  20. Ex. 115.20Application

    Com A=(2112)A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix} (autovalores 3 e 1), calcule Ak(20)A^k \begin{pmatrix} 2 \\ 0 \end{pmatrix} em termos de kk.

    Show solution
    Com AA de autovalores 3,13, 1 e autovetores v1=(1,1)T,v2=(1,1)Tv_1=(1,1)^T, v_2=(1,-1)^T. Escreva v0=(2,0)T=c1v1+c2v2v_0=(2,0)^T = c_1 v_1 + c_2 v_2: sistema c1+c2=2,c1c2=0c1=c2=1c_1+c_2=2, c_1-c_2=0 \Rightarrow c_1=c_2=1. Logo Akv0=13kv1+11kv2=3k(1,1)T+(1,1)TA^k v_0 = 1 \cdot 3^k v_1 + 1 \cdot 1^k v_2 = 3^k(1,1)^T + (1,-1)^T.
    Show step-by-step (with the why)
    1. Decomponha v0v_0 na base de autovetores: v0=civiv_0 = \sum c_i v_i.
    2. Aplique AkA^k: Akv0=ciλikviA^k v_0 = \sum c_i \lambda_i^k v_i. Cada componente cresce/decai com taxa λik\lambda_i^k.
    3. O comportamento para kk grande é dominado pelo maior λi|\lambda_i|.

    Atalho mental: a decomposição v0=civiv_0 = \sum c_i v_i é o "filtro" — coordenadas na base de autovetores são multiplicadas por potências dos autovalores, independentemente.

  21. Ex. 115.21ApplicationAnswer key

    Resolva x˙=Ax\dot{x} = Ax com A=(1002)A = \begin{pmatrix} -1 & 0 \\ 0 & -2 \end{pmatrix}, x(0)=(1,1)Tx(0) = (1, 1)^T.

    Show solution
    A=diag(1,2)A = \operatorname{diag}(-1,-2). Já diagonal. x(t)=eAtx0=diag(et,e2t)(11)=(ete2t)x(t) = e^{At}x_0 = \operatorname{diag}(e^{-t},e^{-2t})\begin{pmatrix}1\\1\end{pmatrix} = \begin{pmatrix}e^{-t}\\e^{-2t}\end{pmatrix}. Ambos os modos decaem — sistema estável.
  22. Ex. 115.22Application

    Resolva x˙=Ax\dot{x} = Ax com A=(0123)A = \begin{pmatrix} 0 & 1 \\ -2 & -3 \end{pmatrix}.

    Show solution
    pA(λ)=λ2+3λ+2=(λ+1)(λ+2)p_A(\lambda)=\lambda^2+3\lambda+2=(\lambda+1)(\lambda+2). λ1=1\lambda_1=-1: v1=(1,1)Tv_1=(1,-1)^T. λ2=2\lambda_2=-2: v2=(1,2)Tv_2=(1,-2)^T (verifique: (A+2I)v=0(A+2I)v=0). Solução geral: x(t)=c1et(11)+c2e2t(12)x(t) = c_1 e^{-t}\begin{pmatrix}1\\-1\end{pmatrix} + c_2 e^{-2t}\begin{pmatrix}1\\-2\end{pmatrix}. Sistema assintoticamente estável.
  23. Ex. 115.23Application

    Mostre que Fn=15(ϕnψn)F_n = \frac{1}{\sqrt{5}}\left(\phi^n - \psi^n\right) para a sequência de Fibonacci, onde ϕ=1+52\phi = \frac{1+\sqrt{5}}{2}.

    Show solution
    Autovalores de A=(1110)A=\begin{pmatrix}1&1\\1&0\end{pmatrix}: ϕ\phi e ψ\psi. Derivada: (d/dt)eAtt=0=A(d/dt)e^{At}|_{t=0} = A. Mais diretamente: a fórmula Fn=(ϕnψn)/5F_n = (\phi^n-\psi^n)/\sqrt5 segue da decomposição de x0=(F1,F0)T=(1,0)Tx_0=(F_1,F_0)^T=(1,0)^T na base de autovetores. F1=1,F0=0F_1=1, F_0=0 — verificar: (ϕ1ψ1)/5=5/5=1(\phi^1-\psi^1)/\sqrt5 = \sqrt5/\sqrt5=1. \checkmark
  24. Ex. 115.24Challenge

    Se AA é diagonalizável e ff é um polinômio, mostre que f(A)=Pf(D)P1f(A) = Pf(D)P^{-1} com f(D)=diag(f(λ1),,f(λn))f(D) = \operatorname{diag}(f(\lambda_1), \ldots, f(\lambda_n)).

    Show solution
    Mostre que se ff é polinômio e A=PDP1A = PDP^{-1}, então f(A)=Pf(D)P1f(A) = Pf(D)P^{-1} onde f(D)=diag(f(λ1),,f(λn))f(D)=\operatorname{diag}(f(\lambda_1),\ldots,f(\lambda_n)). Prova: linearidade de ff e o fato de que Ak=PDkP1A^k = PD^kP^{-1} (já provado). Portanto f(A)=akAk=akPDkP1=P(akDk)P1=Pf(D)P1f(A)=\sum a_k A^k = \sum a_k PD^kP^{-1} = P(\sum a_k D^k)P^{-1} = Pf(D)P^{-1}.
  25. Ex. 115.25Modeling

    Cadeia de Markov de clima: M=(0,70,30,40,6)M = \begin{pmatrix} 0{,}7 & 0{,}3 \\ 0{,}4 & 0{,}6 \end{pmatrix}. Calcule M10M^{10} via diagonalização e encontre a distribuição estacionária.

    Show solution
    Cadeia de Markov: P=(0,70,30,40,6)P = \begin{pmatrix}0{,}7 & 0{,}3\\0{,}4 & 0{,}6\end{pmatrix}. Autovalores: λ1=1\lambda_1=1 (distribuição estacionária) e λ2=0,3\lambda_2=0{,}3. Diagonalizando, P10=QD10Q1P^{10} = Q D^{10} Q^{-1} com D10=diag(1,0,310)D^{10}=\operatorname{diag}(1, 0{,}3^{10}). Como 0,3100,00000590{,}3^{10} \approx 0{,}0000059, P10P^{10} está muito próximo da distribuição estacionária: cada linha converge para (4/7,3/7)(4/7, 3/7).
  26. Ex. 115.26Modeling

    Sistema massa-mola acoplado de 2 massas com matriz de rigidez K=(2112)K = \begin{pmatrix} 2 & -1 \\ -1 & 2 \end{pmatrix} (massas unitárias). Encontre os modos normais e as frequências naturais de vibração.

    Show solution
    Sistema massa-mola acoplado: K=(2112)K = \begin{pmatrix}2&-1\\-1&2\end{pmatrix}. Autovalores: λ1=1\lambda_1=1, λ2=3\lambda_2=3. Frequências naturais: ω1=1\omega_1=1 rad/s, ω2=3\omega_2=\sqrt3 rad/s. Autovetores (modos): v1=(1,1)T/2v_1=(1,1)^T/\sqrt2 (in-phase), v2=(1,1)T/2v_2=(1,-1)^T/\sqrt2 (out-of-phase).
    Show step-by-step (with the why)
    1. Matriz de rigidez: K=(2112)K=\begin{pmatrix}2&-1\\-1&2\end{pmatrix} (assumindo M=IM=I).
    2. Autovalores: pK(λ)=(2λ)21=λ24λ+3=(λ1)(λ3)p_K(\lambda)=(2-\lambda)^2-1=\lambda^2-4\lambda+3=(\lambda-1)(\lambda-3).
    3. Frequências naturais: ωi=λi\omega_i = \sqrt{\lambda_i} (rad/s).
    4. Modos: autovetor de λ1=1\lambda_1=1 é (1,1)T(1,1)^T — as duas massas oscilam juntas. Autovetor de λ2=3\lambda_2=3 é (1,1)T(1,-1)^T — oscilam opostas.

    Curiosidade: ressonância ocorre quando a frequência de excitação iguala uma das ωi\omega_i. A ponte de Tacoma Narrows colapsou (1940) por excitação do modo torsional.

  27. Ex. 115.27ModelingAnswer key

    Matriz Leslie de população de 2 faixas etárias: L=(01,20,40)L = \begin{pmatrix} 0 & 1{,}2 \\ 0{,}4 & 0 \end{pmatrix}. Calcule o autovalor dominante e interprete como taxa de crescimento populacional.

    Show solution
    Matriz Leslie simples 2 faixas (jovens J, adultos A): L=(0fs0)L = \begin{pmatrix}0 & f \\ s & 0\end{pmatrix} com ff = fertilidade, ss = sobrevivência. Autovalor dominante (raio espectral) ρ(L)=λ1\rho(L) = \lambda_1 é a taxa de crescimento populacional no longo prazo. λ1>1\lambda_1 > 1: crescimento; λ1<1\lambda_1 < 1: declínio; λ1=1\lambda_1 = 1: estável. O autovetor associado é a estrutura etária estável.
  28. Ex. 115.28Modeling

    PageRank simplificado: 4 páginas com matriz de transição P=14(0111101111011110)P = \frac{1}{4}\begin{pmatrix} 0 & 1 & 1 & 1 \\ 1 & 0 & 1 & 1 \\ 1 & 1 & 0 & 1 \\ 1 & 1 & 1 & 0 \end{pmatrix}. Encontre a distribuição estacionária (autovetor de λ=1\lambda = 1).

    Show solution
    4 páginas com matriz de transição PP estocástica. Distribuição estacionária = autovetor de λ=1\lambda=1. Normalize o autovetor para ter componentes somando 1 (distribuição de probabilidade). Esse autovetor é o vetor PageRank — páginas com maior componente têm maior relevância.
  29. Ex. 115.29Modeling

    Matriz de covariância de 2 ativos: Σ=(4223)\Sigma = \begin{pmatrix} 4 & 2 \\ 2 & 3 \end{pmatrix}. Diagonalize Σ\Sigma e interprete os autovetores como direções principais de risco.

    Show solution
    Covariância de 2 ativos: Σ=(4223)\Sigma = \begin{pmatrix}4 & 2\\2 & 3\end{pmatrix}. Autovalores: p(λ)=λ27λ+8p(\lambda)=\lambda^2-7\lambda+8, λ1,2=(7±17)/2\lambda_{1,2}=(7\pm\sqrt{17})/2. Autovetores = fatores de risco não-correlacionados. A soma dos autovalores = variância total = trΣ=7\operatorname{tr}\Sigma = 7.
  30. Ex. 115.30Modeling

    Sistema de controle discreto xk+1=Axkx_{k+1} = Ax_k com A=(0,50,30,10,4)A = \begin{pmatrix} 0{,}5 & 0{,}3 \\ 0{,}1 & 0{,}4 \end{pmatrix}. Determine se o sistema é estável verificando o raio espectral ρ(A)=maxλi\rho(A) = \max|\lambda_i|.

    Show solution
    Sistema de controle discreto xk+1=Axkx_{k+1} = Ax_k. Estável se ρ(A)=maxλi<1\rho(A) = \max|\lambda_i| < 1. Para A=(0,50,30,10,4)A=\begin{pmatrix}0{,}5&0{,}3\\0{,}1&0{,}4\end{pmatrix}: p(λ)=λ20,9λ+0,17p(\lambda)=\lambda^2-0{,}9\lambda+0{,}17. Discriminante >0> 0, raízes em (0,0,9)(0, 0{,}9). Raio espectral <1< 1 — sistema estável.
  31. Ex. 115.31Modeling

    Em redes neurais recorrentes, exploding/vanishing gradients ocorrem quando o raio espectral ρ(J)\rho(J) do jacobiano da camada é >1> 1 ou <1< 1. Explique o mecanismo via diagonalização e sugira uma solução arquitetural.

    Show solution
    Em redes neurais, gradiente gk=Jkg0g_k = J^k g_0 onde JJ é o jacobiano. Se ρ(J)>1\rho(J) > 1: gradientes explodem. Se ρ(J)<1\rho(J) < 1: gradientes desaparecem. Gradiente saudável requer ρ(J)1\rho(J) \approx 1. Técnicas: batch normalization, gradient clipping, skip connections (ResNets) que mantêm JIJ \approx I localmente.
  32. Ex. 115.32Modeling

    Modelo de esvaziamento de dois tanques acoplados: A=(2101)A = \begin{pmatrix} -2 & 1 \\ 0 & -1 \end{pmatrix}. Resolva x˙=Ax\dot{x} = Ax com x(0)=(1,2)Tx(0) = (1, 2)^T e determine quando x(t)<0,01\|x(t)\| < 0{,}01.

    Show solution
    O sistema x˙=Ax\dot{x}=Ax com A=(2101)A=\begin{pmatrix}-2&1\\0&-1\end{pmatrix}. Autovalores 2-2 e 1-1 — ambos negativos. Solução: x(t)=c1e2tv1+c2etv2x(t)=c_1 e^{-2t}v_1 + c_2 e^{-t}v_2 onde v1=(1,0)T,v2=(1,1)Tv_1=(1,0)^T, v_2=(1,1)^T. Sistema assintoticamente estável. O modo com λ=2\lambda=-2 decai mais rápido.
  33. Ex. 115.33Modeling

    Numa rede de difusão de informação, a dinâmica discreta é xk+1=Wxkx_{k+1} = Wx_k onde WW é simétrica com autovalores 1,0,8,0,21, 0{,}8, 0{,}2. Interprete o que acontece com xkx_k para kk grande.

    Show solution
    Matriz de difusão num grafo: cada autovetor representa um "modo de difusão". O autovalor nulo corresponde ao estado estacionário (distribuição uniforme em grafos regulares). Autovalores próximos de 0 = modos que persistem; autovalores grandes em módulo = modos que se dissipam rápido. Diagonalizando, a solução no tempo é soma de modos exponencialmente amortecidos.
  34. Ex. 115.34Understanding

    Por que uma matriz n×nn \times n com nn autovalores distintos (sobre C\mathbb{C}) é sempre diagonalizável?

    Select the correct option
    Select an option first
    Show solution
    O teorema chave: autovetores correspondentes a autovalores distintos são LI. Com nn autovalores distintos: nn autovetores LI — exatamente o que é necessário para diagonalizar. A opção D é falsa: nn raízes no complexo, mas podem ser repetidas.
  35. Ex. 115.35ProofAnswer key

    Demonstre que autovetores correspondentes a autovalores distintos são linearmente independentes.

    Show solution
    Sejam v1,,vkv_1, \ldots, v_k autovetores de autovalores distintos λ1,,λk\lambda_1, \ldots, \lambda_k. Suponha por contradição que são LD: existe combinação não-trivial c1v1++ckvk=0c_1 v_1 + \cdots + c_k v_k = 0. Aplique AλkIA - \lambda_k I: c1(λ1λk)v1++ck1(λk1λk)vk1=0c_1(\lambda_1-\lambda_k)v_1 + \cdots + c_{k-1}(\lambda_{k-1}-\lambda_k)v_{k-1} = 0. Por indução (hipótese), os v1,,vk1v_1,\ldots,v_{k-1} são LI — logo todos ci(λiλk)=0c_i(\lambda_i-\lambda_k)=0, e como λiλk\lambda_i \neq \lambda_k, temos ci=0c_i=0 para i<ki < k. Então ckvk=0c_k v_k = 0 com vk0v_k \neq 0, logo ck=0c_k = 0. Contradição. \square
  36. Ex. 115.36Proof

    Toda matriz 2×22 \times 2 com detA=0\det A = 0 e trA0\operatorname{tr} A \neq 0 é diagonalizável? Justifique.

    Show solution
    Sim. detA=0λ=0\det A = 0 \Rightarrow \lambda = 0 é autovalor; trA0\operatorname{tr} A \neq 0 \Rightarrow o outro autovalor é trA0\operatorname{tr} A \neq 0. Como os dois autovalores são distintos (0 e trA\operatorname{tr} A), a matriz é diagonalizável.
  37. Ex. 115.37Challenge

    Encontre uma matriz 2×22 \times 2 não-diagonalizável com autovalor 5 de multiplicidade algébrica 2.

    Show solution
    Exemplo: A=(5105)A = \begin{pmatrix}5 & 1 \\ 0 & 5\end{pmatrix}. Autovalor 5 com ma=2m_a=2, mas mg=1m_g=1 (só um autovetor independente). Portanto não é diagonalizável. Qualquer matriz da forma (5c05)\begin{pmatrix}5&c\\0&5\end{pmatrix} com c0c \neq 0 funciona.
  38. Ex. 115.38Proof

    Demonstre que matrizes similares têm o mesmo polinômio característico (e portanto os mesmos autovalores).

    Show solution
    Se AA e BB são similares (B=PAP1B = PAP^{-1}), então det(BλI)=det(PAP1λI)=det(P(AλI)P1)=det(P)det(AλI)det(P1)=det(AλI)\det(B-\lambda I) = \det(PAP^{-1}-\lambda I) = \det(P(A-\lambda I)P^{-1}) = \det(P)\det(A-\lambda I)\det(P^{-1}) = \det(A-\lambda I). Logo têm o mesmo polinômio característico e os mesmos autovalores. \square
  39. Ex. 115.39ChallengeAnswer key

    Mostre que se AA é diagonalizável e ff é um polinômio, então f(A)f(A) é diagonalizável com autovalores f(λi)f(\lambda_i).

    Show solution
    Se AA é diagonalizável, ela é similar a DD. f(A)=Pf(D)P1f(A) = Pf(D)P^{-1} (prova já feita). Como f(D)f(D) é diagonal (com f(λi)f(\lambda_i) na diagonal), f(A)f(A) é similar a uma diagonal — portanto também é diagonalizável, com os mesmos autovetores e autovalores f(λi)f(\lambda_i).
  40. Ex. 115.40Challenge

    Se AA é diagonalizável, prove que ATA^T também é diagonalizável (com os mesmos autovalores).

    Show solution
    Sim. AA diagonalizável com autovalores λi\lambda_i: AT=(PDP1)T=(P1)TDTPT=(PT)TDPTA^T = (PDP^{-1})^T = (P^{-1})^T D^T P^T = (P^T)^{-T} D P^T. Como DT=DD^T = D (diagonal), ATA^T é similar a DD com a "inversível" (PT)1(P^T)^{-1} — portanto também é diagonalizável com os mesmos autovalores. Note que os autovetores de ATA^T diferem dos de AA em geral.
  41. Ex. 115.41Understanding

    Se A=QDQTA = QDQ^T com QQ ortogonal e DD diagonal real, prove que AA é simétrica.

    Show solution
    Se A=QDQTA = QDQ^T com QQ ortogonal e DD diagonal real, então AT=(QDQT)T=QDTQT=QDQT=AA^T = (QDQ^T)^T = QD^TQ^T = QDQ^T = A. Logo AA é simétrica. Recíproca (Teorema Espectral, próxima aula): simétrica real \Rightarrow ortogonalmente diagonalizável.
  42. Ex. 115.42Understanding

    Se AA é diagonalizável com autovalores λ1,,λn\lambda_1, \ldots, \lambda_n, qual é a relação entre trA\operatorname{tr} A, detA\det A e os autovalores?

    Select the correct option
    Select an option first
    Show solution
    Como A=PDP1A = PDP^{-1}: trA=tr(PDP1)=tr(D)=λi\operatorname{tr} A = \operatorname{tr}(PDP^{-1}) = \operatorname{tr}(D) = \sum \lambda_i (trace é invariante por similaridade). detA=det(PDP1)=detD=λi\det A = \det(PDP^{-1}) = \det D = \prod \lambda_i. Relações fundamentais.
  43. Ex. 115.43Challenge

    Mostre que ABAB e BABA têm os mesmos autovalores não-nulos (mesmo que ABBAAB \neq BA).

    Show solution
    ABAB e BABA podem diferir mas têm os mesmos autovalores não-nulos. Prova: se ABv=λvABv = \lambda v com λ0\lambda \neq 0, então BA(Bv)=B(ABv)=λ(Bv)BA(Bv) = B(ABv) = \lambda(Bv). Como λ0\lambda \neq 0, Bv0Bv \neq 0 (pois ABv=λv0ABv = \lambda v \neq 0). Logo BvBv é autovetor de BABA com mesmo λ\lambda.
  44. Ex. 115.44ProofAnswer key

    Se AA é diagonalizável e inversível, demonstre que A1A^{-1} também é diagonalizável com autovalores 1/λi1/\lambda_i.

    Show solution
    Seja A=PDP1A = PDP^{-1} com D=diag(λ1,,λn)D = \operatorname{diag}(\lambda_1,\ldots,\lambda_n). A1A^{-1} existe     detA=λi0    λi0\iff \det A = \prod \lambda_i \neq 0 \iff \lambda_i \neq 0 para todo ii. Quando inversível: A1=PD1P1A^{-1} = PD^{-1}P^{-1} com D1=diag(1/λ1,,1/λn)D^{-1} = \operatorname{diag}(1/\lambda_1,\ldots,1/\lambda_n). Os autovalores de A1A^{-1} são os recíprocos dos autovalores de AA.
  45. Ex. 115.45Challenge

    Sistema de reações químicas ABA \rightleftharpoons B com equações c˙A=k1cA+k1cB\dot{c}_A = -k_1 c_A + k_{-1}c_B, c˙B=k1cAk1cB\dot{c}_B = k_1 c_A - k_{-1}c_B. Resolva via diagonalização e encontre o equilíbrio.

    Show solution
    Em química, concentrações satisfazem c˙=Kc\dot{c} = Kc com KK matriz de taxas (Kolmogorov). Diagonalizando: c(t)=eKtc0=PeDtP1c0c(t) = e^{Kt}c_0 = Pe^{Dt}P^{-1}c_0. Os modos de decaimento são exponenciais com constantes de tempo 1/λi1/|\lambda_i|. O autovetor de λ=0\lambda=0 é o equilíbrio. Aplicação: cinética enzimática (Michaelis-Menten), decaimento radioativo em cadeia (Ra → Rn → Po).

Fontes

  • A First Course in Linear Algebra — Robert A. Beezer · 2022 · EN · GNU FDL. Referência primária: §SD (Similar Matrices and Diagonalization) com definições rigorosas e exercícios numerados.
  • Linear Algebra Done Right (4ª ed) — Sheldon Axler · 2024 · EN · CC-BY-NC. Cap. 5C–5D: operadores diagonalizáveis, polinômios e funções de operadores.
  • Linear Algebra — Jim Hefferon · 2022 · EN · CC-BY-SA. Cap. 5 §II: diagonalização, Jordan introdutória, exemplos de sistemas dinâmicos.

Updated on 2026-05-06 · Author(s): Clube da Matemática

Found an error? Open an issue on GitHub or submit a PR — open source forever.