Cálculo de las potencias matriciales de una matriz no diagonalizable.

Question

Definir \begin{equation} A = \begin{pmatrix} \frac{1}{2} &\frac{1}{2} & 0\\ 0& \frac{3}{4} & \frac{1}{4}\\ 0& \frac{1}{4} & \frac{3}{4} \end {pmatrix}. \ end {equation}

Tenga en cuenta que la suma de las dimensiones de los espacios propios de $A$ es solo dos. $A$ no es diagonalizable. ¿Cómo podemos calcular $A^n$ ?

Answer 1

Aquí es una manera diferente de usar un clásico truco, convirtiendo el problema en un binomio de expansión. En efecto, podemos escribir :

$$A=\frac12(I+B) \ \text{where} \ B:=\begin{pmatrix}0&1&0 \\0&1/2&1/2\\0&1/2&1/2\end{pmatrix}$$

donde la matriz de $B$ tiene la siguiente particularidad

$$B^n=C \ \text{for all} \ n>1 \ \text{where} \ C:=\begin{pmatrix}0&1/2&1/2\\0&1/2&1/2\\0&1/2&1/2\end{pmatrix}$$

Por lo tanto

$$A^n = \dfrac{1}{2^n}\left(I+\binom{n}{1}B+\binom{n}{2}B^2+\binom{n}{3}B^3+\cdots+\binom{n}{n}B^n\right)$$

$$A^n = \dfrac{1}{2^n}\left(I+nB+\binom{n}{2}C+\binom{n}{3}C+\cdots+\binom{n}{n}C\right)\tag{1}$$

Como es bien sabido, $\sum_{k=0}^n \binom{n}{k}=2^n$, la reducción de la (1) a :

$$A^n = \dfrac{1}{2^n}\left(I+nB+(2^n-n-1)C\right)$$

Basta ahora para reemplazar a $B$ e $C$ por su expresión

$$A^n = \dfrac{1}{2^n}\left(\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}+n\begin{pmatrix}0&1&0 \\0&1/2&1/2\\0&1/2&1/2\end{pmatrix}+(2^n-n-1)\begin{pmatrix}0&1/2&1/2\\0&1/2&1/2\\0&1/2&1/2\end{pmatrix}\right)$$

para obtener el resultado (coincidiendo con la suya !).

Answer 2

Tenga en cuenta que su matriz $A$ tiene vectores propios generalizados

\begin{equation}v_1=\begin{pmatrix}1 \\ 0 \\ 0\end{pmatrix}, v_2=\begin{pmatrix}0 \\ 2 \\ -2\end{pmatrix}, v_3=\begin{pmatrix}0\\ 0 \\ 1\end{pmatrix}.\end{equation}

Así, junto al Jordán de descomposición, $A=\big(v_1,v_2,v_3\big)J\big(v_1,v_2,v_3\big)^{-1}$, donde

\begin{equation}J=\begin{pmatrix}\frac12 & 1 & 0\\0 & \frac12 & 0\\0 & 0 & 1\end{pmatrix}.\end{equation}

El problema del cálculo de $A^n$ se reduce así para el cálculo de $J^n$. Deje $a_{ij}^{(n)}$ denotar la entrada de $J^n$ en la $i$-ésima fila y $j$-ésima columna.

El producto de un arbitrario $3\times3$-matriz con $J$ está dada por: \begin{equation} \begin{pmatrix} a&b&c\\d&e&f\\g&h&i \end{pmatrix} J = \begin{pmatrix} \frac a2&a+\frac b2&c\\\frac d2&d+\frac e2&f\\\frac g2&g+\frac h2&i \end{pmatrix}. \end{equation}

Podemos deducir que, para todos los $n\in\Bbb N$: \begin{align} a_{11}^{(n)}&=a_{22}^{(n)}=\frac1{2^n}, \\a_{21}^{(n)}&=a_{31}^{(n)}=0,\\ a_{13}^{(n)}&=a_{23}^{(n)}=a_{32}^{(n)}=0, \\ a_{33}^{(n)}&=1,\\ a_{12}^{(n+1)}&=a_{11}^{(n)}+\frac{a_{12}^{(n)}}2=\frac1{2^n}+\frac{a_{12}^{(n)}}2. \end{align}

Por lo tanto, todos los $a_{ij}^{(n)}$ son explícitamente conocido excepto por $a_{12}^{(n)}$. Tenga en cuenta que, por la última ecuación, \begin{equation}a_{12}^{(n+1)}=2^{-n}+\frac{a_{12}^{(n)}}2 = 2^{-n}+2^{-n}+\frac{a_{12}^{(n-1)}}4 = \dots = (n+1)\cdot2^{-n}.\end{equation}

Por lo tanto, \begin{equation}J^n=\begin{pmatrix}2^{-n}&n\cdot 2^{1-n} & 0\\0 & 2^{-n} & 0\\0 & 0 & 1\end{pmatrix}.\end{equation}

Y, según algunos cálculos, nos encontramos con que \begin{equation} A^n=\big(v_1,v_2,v_3\big)J^n\big(v_1,v_2,v_3\big)^{-1}= \begin{pmatrix} 2^{-n} & n\cdot 2^{-n-1} - 2^{-n-1} + \frac12 & {1-\frac{n+1}{2^n}\over2}\\ 0 & {2^{-n}+1\over2} & {1-2^{-n}\over2} \\ 0 & {1-2^{-n}\over2} & {2^{-n}+1\over2} \end{pmatrix}. \end{equation}