Fatore cada uma das matrizes como um produto P D P^{-1}, onde D é diagonal. E=\left[\begin{array}{ccc}2 & 2 & 1 \\\\ 0 & 1 & 2 \\\\ 0 & 0 & -1\end{array}\right]
Passo 1
O que o enunciado pede: fatorar a matriz E como um produto P D P^{-1}, onde D é diagonal e P é invertível.
Informações fornecidas (peguei estes números diretamente do enunciado):
Portanto, o objetivo final é fornecer uma matriz diagonal D, uma matriz invertível P e a inversa P^{-1} tais que E = P\,D\,P^{-1}.]
Passo 2
Noções básicas que vamos usar (mini-aula): eigenvalores e eigenvetores.
Um vetor não-nulo v é um autovetor de uma matriz A se existe um escalar \lambda tal que
A v = \lambda v.
O escalar \lambda é chamado autovalor (ou eigenvalor) associado a v.
Para encontrar autovalores, resolve-se a equação característica \det(A-\lambda I)=0.
Observação prática muito útil: se A é triangular (superior ou inferior), então os autovalores de A são exatamente os elementos da diagonal principal. Aqui E é triangular superior; portanto os autovalores são as entradas diagonais de E.
Aplicando isso ao nosso caso: os autovalores de E são
Agora vamos calcular os autovetores correspondentes, um por um, resolvendo [(E-\lambda I)v=0). Vou explicar cada equação do sistema e escolher valores convenientes para obter vetores fáceis de manipular.
Se v=(x,y,z)^T, o sistema (E-2I)v=0 dá as equações:
\begin{cases} 2y + z = 0,\\[4pt] - y + 2z = 0,\\[4pt] -3z = 0. \end{cases}
Da terceira equação vem z=0\]. Substituindo nas anteriores obtemos \[y=0. A primeira equação não envolve x, logo x é livre; escolhemos x=1 para um vetor simples. Assim um autovetor associado a \lambda_1=2 é
\begin{cases} x + 2y + z = 0,\\[4pt] 2z = 0,\\[4pt] -2z = 0. \end{cases}
Da segunda/terceira equações vem z=0\]. Então a primeira reduz a \[x + 2y = 0, isto é x = -2y. Escolhendo y=1 (escolha conveniente), obtemos o autovetor
Formamos a matriz P com os autovetores como colunas, na mesma ordem dos autovalores escolhidos (2,1,-1). A matriz diagonal D terá esses autovalores na diagonal, na mesma ordem:
Observação: P é triangular superior com diagonal 1,1,3, portanto \det P = 3 \neq 0 e P é invertível.]
Passo 5
Vamos calcular P^{-1}. Por ser P triangular superior, P^{-1} também será triangular superior; uma maneira prática é resolver P x = e_j para cada coluna j (back-substitution). Faço isso coluna a coluna:
Para e_1=(1,0,0)^T, resolver P x = e_1 dá x^{(1)}=(1,0,0)^T.
Para e_2=(0,1,0)^T, resolver P x = e_2 dá x^{(2)}=(2,1,0)^T.
Para e_3=(0,0,1)^T, resolver P x = e_3 dá x^{(3)}=\bigl(\tfrac{5}{3},1,\tfrac{1}{3}\bigr)^T.
Verificação rápida (mini-checagem): como verificação parcial, multiplicando P por P^{-1} obtém-se a identidade (confere-se a construção da inversa), e calcular P D P^{-1} recupera E. Por exemplo, P D =