So the formula to generate the fibonacci sequence in matrix form is: $$ \begin{pmatrix} 1 & 1 \\ 1 & 0 \\ \end{pmatrix}^n = \begin{pmatrix} F_{n+1} & F_n \\ F_n & F_{n-1} \\ \end{pmatrix} $$
So, is there a way to generalize it to a larger matrix, a 4x4 for example?
An attempt at a direct generalization gives the cubes of the Fibonacci numbers:
$\pmatrix{ 3&6&-3&-1\\ 1&0&0&0\\ 0&1&0&0\\ 0&0&1&0 }^n\pmatrix{ 2197&512&125&27\\ 512&125&27&8\\ 125&27&8&1\\ 27&8&1&1 }=\pmatrix{ F_{n+7}^3&F_{n+6}^3&F_{n+5}^3&F_{n+4}^3\\ F_{n+6}^3&F_{n+5}^3&F_{n+4}^3&F_{n+3}^3\\ F_{n+5}^3&F_{n+4}^3&F_{n+3}^3&F_{n+2}^3\\ F_{n+4}^3&F_{n+3}^3&F_{n+2}^3&F_{n+1}^3\\ } $
It's not as clean as the original because there's an "accident" in the original where the coefficients of the recurrence relation are also the first two terms in the series. No such luck here.